Page 1 Source Listing PRODQ2
2025-03-12 18:21 /tmp/ifortfekGfL.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-----------------------------------------------------------------------
38 P A R A M E T E R
39 &(LM1=LM-1)
40 C-----------------------------------------------------------------------
41 P A R A M E T E R
42 &(EPSQ2=0.2,EPSL=0.32,EPSTRB=1.E-24,EPS1=1.E-12,EPS2=0.)
43 C-----------------------------------------------------------------------
44 P A R A M E T E R
45 C-----------------------------------------------------------------------
46 &(G=9.8,BETA=1./270.,BTG=BETA*G
47 &,PRT=1.0,GAM1=0.2222222222222222222
48 &,A1=0.659888514560862645,A2=0.6574209922667784586
49 &,B1=11.87799326209552761,B2=7.226971804046074028
50 &,C1=0.000830955950095854396)
51 C-----------------------------------------------------------------------
52 CUN &(G=9.8,BETA=1./270.,BTG=BETA*G
53 CUN &,PRT=1.0,GAM1=0.2222222222222222222
54 CUN &,A1=0.3310949523016403346,A2=0.8273378704055731278
55 CUN &,B1=5.959709141429526024,B2=3.626088092074591135
56 CUN &,C1=-0.3330651924968952113)
57 C-----------------------------------------------------------------------
58 CMY &(G=9.8,BETA=1./270.,BTG=BETA*G
59 CMY &,PRT=0.8,GAM1=0.2222222222222222222
60 CMY &,A1=0.9222222350809054114,A2=0.7350190142719400952
61 CMY &,B1=16.60000023145629741,B2=10.10000014082581951
62 CMY &,C1=0.0805318118080613468)
63 C-----------------------------------------------------------------------
64 P A R A M E T E R
65 &(RB1=1./B1
66 C--------------COEFFICIENTS OF THE TERMS IN THE NUMERATOR---------------
67 &,ANMM=-3.*A1*A2*(3.*A2+3.*B2*C1+18.*A1*C1-B2)*BTG
68 &,ANMH=-9.*A1*A2*A2*BTG*BTG
69 &,BNMM= A1*(1.-3.*C1)
70 &,BNMH= -A2*BTG
71 C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR-------------
72 &,ADNM=18.*A1*A1*A2*(B2-3.*A2)*BTG
73 &,ADNH= 9.*A1*A2*A2*(12.*A1+3.*B2)*BTG*BTG
74 &,BDNM= 6.*A1*A1
75 &,BDNH= 3.*A2*(7.*A1+B2)*BTG
76 C--------------COEFFICIENTS OF THE EQUILIBRIUM EQUATION-----------------
77 &,AEQM= 3.*A1*A2*B1*(3.*A2+3.*B2*C1+18.*A1*C1-B2)*BTG
78 & +18.*A1*A1*A2*(B2-3.*A2)*BTG
79 &,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:21 PRODQ2.F
80 &,BEQM=-A1*B1*(1.-3.*C1)+6.*A1*A1
81 &,BEQH= A2*B1*BTG+3.*A2*(7.*A1+B2)*BTG
82 C--------------FORBIDDEN TURBULENCE AREA--------------------------------
83 &,REQU=-AEQH/AEQM*1.02,EPSGH=1.E-9)
84 C-----------------------------------------------------------------------
85 D I M E N S I O N
86 & Q2 (LM)
87 D I M E N S I O N
88 & GM (LM1),GH (LM1),EL (LM1)
89 C-----------------------------------------------------------------------
90 C***********************************************************************
91 LMHM=LMHK-1
92 C
93 DO 150 L=1,LMHM
94 GML=GM(L)
95 GHL=GH(L)
96 C--------------COEFFICIENTS OF THE EQUILIBRIUM EQUATION-----------------
97 AEQU=(AEQM*GML+AEQH*GHL)*GHL
98 BEQU= BEQM*GML+BEQH*GHL
99 C--------------EQUILIBRIUM SOLUTION FOR L/Q-----------------------------
100 EQOL2=-0.5*BEQU+SQRT(BEQU*BEQU*0.25-AEQU)
101 C--------------IS THERE PRODUCTION/DISSIPATION ?------------------------
102 IF((GML+GHL*GHL.LE.EPSTRB )
103 & .OR.(GHL.GE.EPSGH.AND.GML/GHL.LE.REQU)
104 & .OR.(EQOL2.LE.EPS2) ) THEN
105 C--------------NO TURBULENCE--------------------------------------------
106 Q2(L)=EPSQ2
107 EL(L)=EPSL
108 C--------------END OF THE NO TURBULENCE BRANCH--------------------------
109 ELSE
110 C--------------COEFFICIENTS OF THE TERMS IN THE NUMERATOR---------------
111 ANUM=(ANMM*GML+ANMH*GHL)*GHL
112 BNUM= BNMM*GML+BNMH*GHL
113 C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR-------------
114 ADEN=(ADNM*GML+ADNH*GHL)*GHL
115 BDEN= BDNM*GML+BDNH*GHL
116 CDEN= 1.
117 C--------------COEFFICIENTS OF THE NUMERATOR OF THE LINEARIZED EQ.------
118 ARHS=-(ANUM*BDEN-BNUM*ADEN)*2.
119 BRHS=- ANUM*4.
120 CRHS=- BNUM*2.
121 C--------------INITIAL VALUE OF L/Q-------------------------------------
122 DLOQ1=EL(L)/SQRT(Q2(L))
123 C--------------FIRST ITERATION FOR L/Q, RHS=0---------------------------
124 ELOQ21=1./EQOL2
125 ELOQ11=SQRT(ELOQ21)
126 ELOQ31=ELOQ21*ELOQ11
127 ELOQ41=ELOQ21*ELOQ21
128 ELOQ51=ELOQ21*ELOQ31
129 C--------------1./DENOMINATOR-------------------------------------------
130 RDEN1=1./(ADEN*ELOQ41+BDEN*ELOQ21+CDEN)
131 C--------------D(RHS)/D(L/Q)--------------------------------------------
132 RHSP1= (ARHS*ELOQ51+BRHS*ELOQ31+CRHS*ELOQ11)*RDEN1*RDEN1
133 C--------------FIRST-GUESS SOLUTION-------------------------------------
134 ELOQ12=ELOQ11+(DLOQ1-ELOQ11)*EXP(RHSP1*DTQ2)
135 C-----------------------------------------------------------------------
136 ELOQ12=AMAX1(ELOQ12,EPS1)
Page 3 Source Listing PRODQ2
2025-03-12 18:21 PRODQ2.F
137 C--------------SECOND ITERATION FOR L/Q---------------------------------
138 ELOQ22=ELOQ12*ELOQ12
139 ELOQ32=ELOQ22*ELOQ12
140 ELOQ42=ELOQ22*ELOQ22
141 ELOQ52=ELOQ22*ELOQ32
142 C--------------1./DENOMINATOR-------------------------------------------
143 RDEN2=1./(ADEN*ELOQ42+BDEN*ELOQ22+CDEN)
144 C-----------------------------------------------------------------------
145 RHS2 =-(ANUM*ELOQ42+BNUM*ELOQ22)*RDEN2+RB1
146 RHSP2= (ARHS*ELOQ52+BRHS*ELOQ32+CRHS*ELOQ12)*RDEN2*RDEN2
147 RHST2=RHS2/RHSP2
148 C--------------CORRECTED SOLUTION---------------------------------------
149 ELOQ13=ELOQ12-RHST2+(RHST2+DLOQ1-ELOQ12)*EXP(RHSP2*DTQ2)
150 C-----------------------------------------------------------------------
151 ELOQ13=AMAX1(ELOQ13,EPS1)
152 C--------------TWO ITERATIONS IS ENOUGH IN MOST CASES ...---------------
153 ELOQN=ELOQ13
154 C-----------------------------------------------------------------------
155 IF(ELOQN.GT.EPS1)THEN
156 Q2(L)=EL(L)*EL(L)/(ELOQN*ELOQN)
157 Q2(L)=AMAX1(Q2(L),EPSQ2)
158 ELSE
159 Q2(L)=EPSQ2
160 ENDIF
161 C--------------END OF TURBULENT BRANCH----------------------------------
162 ENDIF
163 C--------------END OF PRODUCTION/DISSIPATION LOOP-----------------------
164 150 CONTINUE
165 C--------------LOWER BOUNDARY CONDITION FOR Q2--------------------------
166 Q2(LMHK)=AMAX1(B1**(2./3.)*USTAR*USTAR,EPSQ2)
167 C-----------------------------------------------------------------------
168 RETURN
169 END
Page 4 Source Listing PRODQ2
2025-03-12 18:21 Entry Points PRODQ2.F
ENTRY POINTS
Name
prodq2_
SYMBOL CROSS REFERENCE
Name Object Declared Type Bytes Dimen Elements Attributes References
150 Label 164 93
A1 Param 48 R(4) 4 scalar 67,68,69,72,73,74,75,77,78,79,80,8
1
A2 Param 48 R(4) 4 scalar 67,68,70,72,73,75,77,78,79,81
ADEN Local 114 R(4) 4 scalar 114,118,130,143
ADNH Param 73 R(4) 4 scalar 114
ADNM Param 72 R(4) 4 scalar 114
AEQH Param 79 R(4) 4 scalar 83,97
AEQM Param 77 R(4) 4 scalar 83,97
AEQU Local 97 R(4) 4 scalar 97,100
AMAX1 Func 136 scalar 136,151,157,166
ANMH Param 68 R(4) 4 scalar 111
ANMM Param 67 R(4) 4 scalar 111
ANUM Local 111 R(4) 4 scalar 111,118,119,145
ARHS Local 118 R(4) 4 scalar 118,132,146
B1 Param 49 R(4) 4 scalar 65,77,79,80,81,166
B2 Param 49 R(4) 4 scalar 67,72,73,75,77,78,79,81
BDEN Local 115 R(4) 4 scalar 115,118,130,143
BDNH Param 75 R(4) 4 scalar 115
BDNM Param 74 R(4) 4 scalar 115
BEQH Param 81 R(4) 4 scalar 98
BEQM Param 80 R(4) 4 scalar 98
BEQU Local 98 R(4) 4 scalar 98,100
BETA Param 46 R(4) 4 scalar 46
BNMH Param 70 R(4) 4 scalar 112
BNMM Param 69 R(4) 4 scalar 112
BNUM Local 112 R(4) 4 scalar 112,118,120,145
BRHS Local 119 R(4) 4 scalar 119,132,146
BTG Param 46 R(4) 4 scalar 67,68,70,72,73,75,77,78,79,81
C1 Param 50 R(4) 4 scalar 67,69,77,80
CDEN Local 116 R(4) 4 scalar 116,130,143
CRHS Local 120 R(4) 4 scalar 120,132,146
DLOQ1 Local 122 R(4) 4 scalar 122,134,149
DTQ2 Dummy 3 R(4) 4 scalar ARG,INOUT 134,149
EL Dummy 3 R(4) 4 1 44 ARG,INOUT 107,122,156
ELOQ11 Local 125 R(4) 4 scalar 125,126,132,134
ELOQ12 Local 134 R(4) 4 scalar 134,136,138,139,146,149
ELOQ13 Local 149 R(4) 4 scalar 149,151,153
ELOQ21 Local 124 R(4) 4 scalar 124,125,126,127,128,130
ELOQ22 Local 138 R(4) 4 scalar 138,139,140,141,143,145
ELOQ31 Local 126 R(4) 4 scalar 126,128,132
ELOQ32 Local 139 R(4) 4 scalar 139,141,146
ELOQ41 Local 127 R(4) 4 scalar 127,130
ELOQ42 Local 140 R(4) 4 scalar 140,143,145
Page 5 Source Listing PRODQ2
2025-03-12 18:21 Symbol Table PRODQ2.F
Name Object Declared Type Bytes Dimen Elements Attributes References
ELOQ51 Local 128 R(4) 4 scalar 128,132
ELOQ52 Local 141 R(4) 4 scalar 141,146
ELOQN Local 153 R(4) 4 scalar 153,155,156
EPS1 Param 42 R(4) 4 scalar 136,151,155
EPS2 Param 42 R(4) 4 scalar 104
EPSGH Param 83 R(4) 4 scalar 103
EPSL Param 42 R(4) 4 scalar 107
EPSQ2 Param 42 R(4) 4 scalar 106,157,159,166
EPSTRB Param 42 R(4) 4 scalar 102
EQOL2 Local 100 R(4) 4 scalar 100,104,124
EXP Func 134 scalar 134,149
G Param 46 R(4) 4 scalar 46
GAM1 Param 47 R(4) 4 scalar
GH Dummy 3 R(4) 4 1 44 ARG,INOUT 95
GHL Local 95 R(4) 4 scalar 95,97,98,102,103,111,112,114,115
GM Dummy 3 R(4) 4 1 44 ARG,INOUT 94
GML Local 94 R(4) 4 scalar 94,97,98,102,103,111,112,114,115
IDIM1 Param 34 I(4) 4 scalar
IDIM2 Param 34 I(4) 4 scalar
IGSTL Param 29 I(4) 4 scalar 34
IGSTR Param 29 I(4) 4 scalar 34
IM Param 18 I(4) 4 scalar 31,34
INPES Param 26 I(4) 4 scalar 31,34
ITAIL Param 31 I(4) 4 scalar
JDIM1 Param 35 I(4) 4 scalar
JDIM2 Param 35 I(4) 4 scalar
JGSTL Param 30 I(4) 4 scalar 35
JGSTR Param 30 I(4) 4 scalar 35
JM Param 18 I(4) 4 scalar 32,35
JNPES Param 26 I(4) 4 scalar 32,35
JTAIL Param 32 I(4) 4 scalar
L Local 93 I(4) 4 scalar 93,94,95,106,107,122,156,157,159
LM Param 18 I(4) 4 scalar 39,86
LM1 Param 39 I(4) 4 scalar 88
LMHK Dummy 3 I(4) 4 scalar ARG,INOUT 91,166
LMHM Local 91 I(4) 4 scalar 91,93
LSM Param 18 I(4) 4 scalar
PRODQ2 Subr 3
PRT Param 47 R(4) 4 scalar
Q2 Dummy 3 R(4) 4 1 45 ARG,INOUT 106,122,156,157,159,166
RB1 Param 65 R(4) 4 scalar 145
RDEN1 Local 130 R(4) 4 scalar 130,132
RDEN2 Local 143 R(4) 4 scalar 143,145,146
REQU Param 83 R(4) 4 scalar 103
RHS2 Local 145 R(4) 4 scalar 145,147
RHSP1 Local 132 R(4) 4 scalar 132,134
RHSP2 Local 146 R(4) 4 scalar 146,147,149
RHST2 Local 147 R(4) 4 scalar 147,149
SQRT Func 100 scalar 100,122,125
USTAR Dummy 3 R(4) 4 scalar ARG,INOUT 166
Page 6 Source Listing PRODQ2
2025-03-12 18:21 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:21 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:21 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