Page 1 Source Listing PRODQ2 2025-03-12 18:22 /tmp/ifortOWTPPh.i 1 # 1 "PRODQ2.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE PRODQ2(LMHK,DTQ2,USTAR,GM,GH,EL,Q2) 4 C ****************************************************************** 5 C * * 6 C * LEVEL 2.5 Q2 PRODUCTION/DISSIPATION * 7 C * * 8 C ****************************************************************** 9 C----------------------------------------------------------------------- 10 INCLUDE "parmeta" 11 # 1 "./sp.h" 1 12 # 4 13 14 # 11 "PRODQ2.F" 2 15 C----------------------------------------------------------------------- 39 P A R A M E T E R 40 &(LM1=LM-1) 41 C----------------------------------------------------------------------- 42 P A R A M E T E R 43 &(EPSQ2=0.2,EPSL=0.32,EPSTRB=1.E-24,EPS1=1.E-12,EPS2=0.) 44 C----------------------------------------------------------------------- 45 P A R A M E T E R 46 C----------------------------------------------------------------------- 47 &(G=9.8,BETA=1./270.,BTG=BETA*G 48 &,PRT=1.0,GAM1=0.2222222222222222222 49 &,A1=0.659888514560862645,A2=0.6574209922667784586 50 &,B1=11.87799326209552761,B2=7.226971804046074028 51 &,C1=0.000830955950095854396) 52 C----------------------------------------------------------------------- 53 CUN &(G=9.8,BETA=1./270.,BTG=BETA*G 54 CUN &,PRT=1.0,GAM1=0.2222222222222222222 55 CUN &,A1=0.3310949523016403346,A2=0.8273378704055731278 56 CUN &,B1=5.959709141429526024,B2=3.626088092074591135 57 CUN &,C1=-0.3330651924968952113) 58 C----------------------------------------------------------------------- 59 CMY &(G=9.8,BETA=1./270.,BTG=BETA*G 60 CMY &,PRT=0.8,GAM1=0.2222222222222222222 61 CMY &,A1=0.9222222350809054114,A2=0.7350190142719400952 62 CMY &,B1=16.60000023145629741,B2=10.10000014082581951 63 CMY &,C1=0.0805318118080613468) 64 C----------------------------------------------------------------------- 65 P A R A M E T E R 66 &(RB1=1./B1 67 C--------------COEFFICIENTS OF THE TERMS IN THE NUMERATOR--------------- 68 &,ANMM=-3.*A1*A2*(3.*A2+3.*B2*C1+18.*A1*C1-B2)*BTG 69 &,ANMH=-9.*A1*A2*A2*BTG*BTG 70 &,BNMM= A1*(1.-3.*C1) 71 &,BNMH= -A2*BTG 72 C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR------------- 73 &,ADNM=18.*A1*A1*A2*(B2-3.*A2)*BTG 74 &,ADNH= 9.*A1*A2*A2*(12.*A1+3.*B2)*BTG*BTG 75 &,BDNM= 6.*A1*A1 76 &,BDNH= 3.*A2*(7.*A1+B2)*BTG 77 C--------------COEFFICIENTS OF THE EQUILIBRIUM EQUATION----------------- 78 &,AEQM= 3.*A1*A2*B1*(3.*A2+3.*B2*C1+18.*A1*C1-B2)*BTG 79 & +18.*A1*A1*A2*(B2-3.*A2)*BTG 80 &,AEQH= 9.*A1*A2*A2*B1*BTG*BTG+9.*A1*A2*A2*(12.*A1+3.*B2)*BTG*BTG Page 2 Source Listing PRODQ2 2025-03-12 18:22 PRODQ2.F 81 &,BEQM=-A1*B1*(1.-3.*C1)+6.*A1*A1 82 &,BEQH= A2*B1*BTG+3.*A2*(7.*A1+B2)*BTG 83 C--------------FORBIDDEN TURBULENCE AREA-------------------------------- 84 &,REQU=-AEQH/AEQM*1.02,EPSGH=1.E-9) 85 C----------------------------------------------------------------------- 86 D I M E N S I O N 87 & Q2 (LM) 88 D I M E N S I O N 89 & GM (LM1),GH (LM1),EL (LM1) 90 C----------------------------------------------------------------------- 91 C*********************************************************************** 92 LMHM=LMHK-1 93 C 94 DO 150 L=1,LMHM 95 GML=GM(L) 96 GHL=GH(L) 97 C--------------COEFFICIENTS OF THE EQUILIBRIUM EQUATION----------------- 98 AEQU=(AEQM*GML+AEQH*GHL)*GHL 99 BEQU= BEQM*GML+BEQH*GHL 100 C--------------EQUILIBRIUM SOLUTION FOR L/Q----------------------------- 101 EQOL2=-0.5*BEQU+SQRT(BEQU*BEQU*0.25-AEQU) 102 C--------------IS THERE PRODUCTION/DISSIPATION ?------------------------ 103 IF((GML+GHL*GHL.LE.EPSTRB ) 104 & .OR.(GHL.GE.EPSGH.AND.GML/GHL.LE.REQU) 105 & .OR.(EQOL2.LE.EPS2) ) THEN 106 C--------------NO TURBULENCE-------------------------------------------- 107 Q2(L)=EPSQ2 108 EL(L)=EPSL 109 C--------------END OF THE NO TURBULENCE BRANCH-------------------------- 110 ELSE 111 C--------------COEFFICIENTS OF THE TERMS IN THE NUMERATOR--------------- 112 ANUM=(ANMM*GML+ANMH*GHL)*GHL 113 BNUM= BNMM*GML+BNMH*GHL 114 C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR------------- 115 ADEN=(ADNM*GML+ADNH*GHL)*GHL 116 BDEN= BDNM*GML+BDNH*GHL 117 CDEN= 1. 118 C--------------COEFFICIENTS OF THE NUMERATOR OF THE LINEARIZED EQ.------ 119 ARHS=-(ANUM*BDEN-BNUM*ADEN)*2. 120 BRHS=- ANUM*4. 121 CRHS=- BNUM*2. 122 C--------------INITIAL VALUE OF L/Q------------------------------------- 123 DLOQ1=EL(L)/SQRT(Q2(L)) 124 C--------------FIRST ITERATION FOR L/Q, RHS=0--------------------------- 125 ELOQ21=1./EQOL2 126 ELOQ11=SQRT(ELOQ21) 127 ELOQ31=ELOQ21*ELOQ11 128 ELOQ41=ELOQ21*ELOQ21 129 ELOQ51=ELOQ21*ELOQ31 130 C--------------1./DENOMINATOR------------------------------------------- 131 RDEN1=1./(ADEN*ELOQ41+BDEN*ELOQ21+CDEN) 132 C--------------D(RHS)/D(L/Q)-------------------------------------------- 133 RHSP1= (ARHS*ELOQ51+BRHS*ELOQ31+CRHS*ELOQ11)*RDEN1*RDEN1 134 C--------------FIRST-GUESS SOLUTION------------------------------------- 135 ELOQ12=ELOQ11+(DLOQ1-ELOQ11)*EXP(RHSP1*DTQ2) 136 C----------------------------------------------------------------------- 137 ELOQ12=AMAX1(ELOQ12,EPS1) Page 3 Source Listing PRODQ2 2025-03-12 18:22 PRODQ2.F 138 C--------------SECOND ITERATION FOR L/Q--------------------------------- 139 ELOQ22=ELOQ12*ELOQ12 140 ELOQ32=ELOQ22*ELOQ12 141 ELOQ42=ELOQ22*ELOQ22 142 ELOQ52=ELOQ22*ELOQ32 143 C--------------1./DENOMINATOR------------------------------------------- 144 RDEN2=1./(ADEN*ELOQ42+BDEN*ELOQ22+CDEN) 145 C----------------------------------------------------------------------- 146 RHS2 =-(ANUM*ELOQ42+BNUM*ELOQ22)*RDEN2+RB1 147 RHSP2= (ARHS*ELOQ52+BRHS*ELOQ32+CRHS*ELOQ12)*RDEN2*RDEN2 148 RHST2=RHS2/RHSP2 149 C--------------CORRECTED SOLUTION--------------------------------------- 150 ELOQ13=ELOQ12-RHST2+(RHST2+DLOQ1-ELOQ12)*EXP(RHSP2*DTQ2) 151 C----------------------------------------------------------------------- 152 ELOQ13=AMAX1(ELOQ13,EPS1) 153 C--------------TWO ITERATIONS IS ENOUGH IN MOST CASES ...--------------- 154 ELOQN=ELOQ13 155 C----------------------------------------------------------------------- 156 IF(ELOQN.GT.EPS1)THEN 157 Q2(L)=EL(L)*EL(L)/(ELOQN*ELOQN) 158 Q2(L)=AMAX1(Q2(L),EPSQ2) 159 ELSE 160 Q2(L)=EPSQ2 161 ENDIF 162 C--------------END OF TURBULENT BRANCH---------------------------------- 163 ENDIF 164 C--------------END OF PRODUCTION/DISSIPATION LOOP----------------------- 165 150 CONTINUE 166 C--------------LOWER BOUNDARY CONDITION FOR Q2-------------------------- 167 Q2(LMHK)=AMAX1(B1**(2./3.)*USTAR*USTAR,EPSQ2) 168 C----------------------------------------------------------------------- 169 RETURN 170 END Page 4 Source Listing PRODQ2 2025-03-12 18:22 Entry Points PRODQ2.F ENTRY POINTS Name prodq2_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 150 Label 165 94 A1 Param 49 R(4) 4 scalar 68,69,70,73,74,75,76,78,79,80,81,8 2 A2 Param 49 R(4) 4 scalar 68,69,71,73,74,76,78,79,80,82 ADEN Local 115 R(4) 4 scalar 115,119,131,144 ADNH Param 74 R(4) 4 scalar 115 ADNM Param 73 R(4) 4 scalar 115 AEQH Param 80 R(4) 4 scalar 84,98 AEQM Param 78 R(4) 4 scalar 84,98 AEQU Local 98 R(4) 4 scalar 98,101 AMAX1 Func 137 scalar 137,152,158,167 ANMH Param 69 R(4) 4 scalar 112 ANMM Param 68 R(4) 4 scalar 112 ANUM Local 112 R(4) 4 scalar 112,119,120,146 ARHS Local 119 R(4) 4 scalar 119,133,147 B1 Param 50 R(4) 4 scalar 66,78,80,81,82,167 B2 Param 50 R(4) 4 scalar 68,73,74,76,78,79,80,82 BDEN Local 116 R(4) 4 scalar 116,119,131,144 BDNH Param 76 R(4) 4 scalar 116 BDNM Param 75 R(4) 4 scalar 116 BEQH Param 82 R(4) 4 scalar 99 BEQM Param 81 R(4) 4 scalar 99 BEQU Local 99 R(4) 4 scalar 99,101 BETA Param 47 R(4) 4 scalar 47 BNMH Param 71 R(4) 4 scalar 113 BNMM Param 70 R(4) 4 scalar 113 BNUM Local 113 R(4) 4 scalar 113,119,121,146 BRHS Local 120 R(4) 4 scalar 120,133,147 BTG Param 47 R(4) 4 scalar 68,69,71,73,74,76,78,79,80,82 C1 Param 51 R(4) 4 scalar 68,70,78,81 CDEN Local 117 R(4) 4 scalar 117,131,144 CRHS Local 121 R(4) 4 scalar 121,133,147 DLOQ1 Local 123 R(4) 4 scalar 123,135,150 DTQ2 Dummy 3 R(4) 4 scalar ARG,INOUT 135,150 EL Dummy 3 R(4) 4 1 44 ARG,INOUT 108,123,157 ELOQ11 Local 126 R(4) 4 scalar 126,127,133,135 ELOQ12 Local 135 R(4) 4 scalar 135,137,139,140,147,150 ELOQ13 Local 150 R(4) 4 scalar 150,152,154 ELOQ21 Local 125 R(4) 4 scalar 125,126,127,128,129,131 ELOQ22 Local 139 R(4) 4 scalar 139,140,141,142,144,146 ELOQ31 Local 127 R(4) 4 scalar 127,129,133 ELOQ32 Local 140 R(4) 4 scalar 140,142,147 ELOQ41 Local 128 R(4) 4 scalar 128,131 ELOQ42 Local 141 R(4) 4 scalar 141,144,146 Page 5 Source Listing PRODQ2 2025-03-12 18:22 Symbol Table PRODQ2.F Name Object Declared Type Bytes Dimen Elements Attributes References ELOQ51 Local 129 R(4) 4 scalar 129,133 ELOQ52 Local 142 R(4) 4 scalar 142,147 ELOQN Local 154 R(4) 4 scalar 154,156,157 EPS1 Param 43 R(4) 4 scalar 137,152,156 EPS2 Param 43 R(4) 4 scalar 105 EPSGH Param 84 R(4) 4 scalar 104 EPSL Param 43 R(4) 4 scalar 108 EPSQ2 Param 43 R(4) 4 scalar 107,158,160,167 EPSTRB Param 43 R(4) 4 scalar 103 EQOL2 Local 101 R(4) 4 scalar 101,105,125 EXP Func 135 scalar 135,150 G Param 47 R(4) 4 scalar 47 GAM1 Param 48 R(4) 4 scalar GH Dummy 3 R(4) 4 1 44 ARG,INOUT 96 GHL Local 96 R(4) 4 scalar 96,98,99,103,104,112,113,115,116 GM Dummy 3 R(4) 4 1 44 ARG,INOUT 95 GML Local 95 R(4) 4 scalar 95,98,99,103,104,112,113,115,116 IDIM1 Param 35 I(4) 4 scalar IDIM2 Param 35 I(4) 4 scalar IGSTL Param 30 I(4) 4 scalar 35 IGSTR Param 30 I(4) 4 scalar 35 IM Param 18 I(4) 4 scalar 32,35 INPES Param 27 I(4) 4 scalar 32,35 ITAIL Param 32 I(4) 4 scalar JDIM1 Param 36 I(4) 4 scalar JDIM2 Param 36 I(4) 4 scalar JGSTL Param 31 I(4) 4 scalar 36 JGSTR Param 31 I(4) 4 scalar 36 JM Param 18 I(4) 4 scalar 33,36 JNPES Param 27 I(4) 4 scalar 33,36 JTAIL Param 33 I(4) 4 scalar L Local 94 I(4) 4 scalar 94,95,96,107,108,123,157,158,160 LM Param 18 I(4) 4 scalar 40,87 LM1 Param 40 I(4) 4 scalar 89 LMHK Dummy 3 I(4) 4 scalar ARG,INOUT 92,167 LMHM Local 92 I(4) 4 scalar 92,94 LSM Param 18 I(4) 4 scalar PRODQ2 Subr 3 PRT Param 48 R(4) 4 scalar Q2 Dummy 3 R(4) 4 1 45 ARG,INOUT 107,123,157,158,160,167 RB1 Param 66 R(4) 4 scalar 146 RDEN1 Local 131 R(4) 4 scalar 131,133 RDEN2 Local 144 R(4) 4 scalar 144,146,147 REQU Param 84 R(4) 4 scalar 104 RHS2 Local 146 R(4) 4 scalar 146,148 RHSP1 Local 133 R(4) 4 scalar 133,135 RHSP2 Local 147 R(4) 4 scalar 147,148,150 RHST2 Local 148 R(4) 4 scalar 148,150 SQRT Func 101 scalar 101,123,126 USTAR Dummy 3 R(4) 4 scalar ARG,INOUT 167 Page 6 Source Listing PRODQ2 2025-03-12 18:22 Subprograms/Common Blocks PRODQ2.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References PRODQ2 Subr 3 COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -align norec32byte -align norec64byte -align noarray8byte -align noarray16byte -align noarray32byte -align noarray64byte -align noarray128byte -align noarray256byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume nobuffered_stdout -assume byterecl -assume nocontiguous_assumed_shape -assume nocontiguous_pointer -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_complex_align -assume old_unit_star -assume old_inquire_recl -assume old_ldout_format -assume old_ldout_zero -assume noold_logical_assign -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume noprotect_allocates -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume std_minus0_rounding -assume nostd_mod_proc_name -assume std_value -assume realloc_lhs -assume underscore -assume no2underscores -assume norecursion no -auto -auto_scalar no -bintext -ccdefault default -check noarg_temp_created -check noassume -check nobounds -check nocontiguous -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check noshape -check nostack -check nouninitialized -check noudio_iostat -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1910 -D __INTEL_COMPILER_UPDATE=3 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D __amd64 -D __amd64__ -D __INTEL_COMPILER_BUILD_DATE=20200925 -D __INTEL_OFFLOAD Page 7 Source Listing PRODQ2 2025-03-12 18:22 PRODQ2.F -D __MMX__ -D __SSE__ -D __SSE_MATH__ -D __SSE2__ -D __SSE2_MATH__ -D __SSE3__ -D __SSSE3__ -D __SSE4_1__ -D __SSE4_2__ -D __POPCNT__ -D __PCLMUL__ -D __AES__ -D __AVX__ -D __F16C__ -D __AVX_I__ -D __RDRND__ -D __FMA__ -D __FP_FAST_FMA -D __FP_FAST_FMAF -D __BMI__ -D __LZCNT__ -D __AVX2__ -D __haswell -D __haswell__ -D __tune_haswell__ -D __core_avx2 -D __core_avx2__ -D __tune_core_avx2__ -D __CRAY_X86_ROME -D __CRAYXT_COMPUTE_LINUX_TARGET -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg -init noarrays -init nohuge -init noinfinity -init nominus_huge -init nominus_infinity -init nominus_tiny -init nonan -init nosnan -init notiny -init nozero no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude no -o -offload-build=host -openmp-simd -O2 no -pad_source -real_size 32 no -recursive -reentrancy threaded -vec=simd Page 8 Source Listing PRODQ2 2025-03-12 18:22 PRODQ2.F -show nofullpath -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w nodeclarations -w noexternals -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w noshape -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage no -wrap-margins -includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, .f90,./.f90,/opt/cray/pe/mpich/8.1.12/ofi/intel/19.0/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.f90,/usr/lib64/gcc/x86_64-suse-linux/7/include/.f90, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.f90,/usr/include/.f90,/usr/include/.f90,/usr/include/.f90 -list filename : PRODQ2.lst no -o COMPILER: Intel(R) Fortran 19.1-1655