Page 1 Source Listing PRODQ2 2014-12-17 20:47 /tmpnwprd/ifortjz44UD.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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 Entry Points PRODQ2.F ENTRY POINTS Name prodq2_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 150 Label 138 67 A1 Param 22 R(4) 4 scalar 41,42,43,46,47,48,49,51,52,53,54,5 5 A2 Param 22 R(4) 4 scalar 41,42,44,46,47,49,51,52,53,55 ADEN Local 88 R(4) 4 scalar 88,92,104,117 ADNH Param 47 R(4) 4 scalar 88 ADNM Param 46 R(4) 4 scalar 88 AEQH Param 53 R(4) 4 scalar 57,71 AEQM Param 51 R(4) 4 scalar 57,71 AEQU Local 71 R(4) 4 scalar 71,74 AMAX1 Func 110 scalar 110,125,131,140 ANMH Param 42 R(4) 4 scalar 85 ANMM Param 41 R(4) 4 scalar 85 ANUM Local 85 R(4) 4 scalar 85,92,93,119 ARHS Local 92 R(4) 4 scalar 92,106,120 B1 Param 23 R(4) 4 scalar 39,51,53,54,55,140 B2 Param 23 R(4) 4 scalar 41,46,47,49,51,52,53,55 BDEN Local 89 R(4) 4 scalar 89,92,104,117 BDNH Param 49 R(4) 4 scalar 89 BDNM Param 48 R(4) 4 scalar 89 BEQH Param 55 R(4) 4 scalar 72 BEQM Param 54 R(4) 4 scalar 72 BEQU Local 72 R(4) 4 scalar 72,74 BETA Param 20 R(4) 4 scalar 20 BNMH Param 44 R(4) 4 scalar 86 BNMM Param 43 R(4) 4 scalar 86 BNUM Local 86 R(4) 4 scalar 86,92,94,119 BRHS Local 93 R(4) 4 scalar 93,106,120 BTG Param 20 R(4) 4 scalar 41,42,44,46,47,49,51,52,53,55 C1 Param 24 R(4) 4 scalar 41,43,51,54 CDEN Local 90 R(4) 4 scalar 90,104,117 CRHS Local 94 R(4) 4 scalar 94,106,120 DLOQ1 Local 96 R(4) 4 scalar 96,108,123 DTQ2 Dummy 2 R(4) 4 scalar ARG,INOUT 108,123 EL Dummy 2 R(4) 4 1 44 ARG,INOUT 81,96,130 ELOQ11 Local 99 R(4) 4 scalar 99,100,106,108 ELOQ12 Local 108 R(4) 4 scalar 108,110,112,113,120,123 ELOQ13 Local 123 R(4) 4 scalar 123,125,127 ELOQ21 Local 98 R(4) 4 scalar 98,99,100,101,102,104 ELOQ22 Local 112 R(4) 4 scalar 112,113,114,115,117,119 ELOQ31 Local 100 R(4) 4 scalar 100,102,106 ELOQ32 Local 113 R(4) 4 scalar 113,115,120 ELOQ41 Local 101 R(4) 4 scalar 101,104 ELOQ42 Local 114 R(4) 4 scalar 114,117,119 Page 5 Source Listing PRODQ2 2014-12-17 20:47 Symbol Table PRODQ2.F Name Object Declared Type Bytes Dimen Elements Attributes References ELOQ51 Local 102 R(4) 4 scalar 102,106 ELOQ52 Local 115 R(4) 4 scalar 115,120 ELOQN Local 127 R(4) 4 scalar 127,129,130 EPS1 Param 16 R(4) 4 scalar 110,125,129 EPS2 Param 16 R(4) 4 scalar 78 EPSGH Param 57 R(4) 4 scalar 77 EPSL Param 16 R(4) 4 scalar 81 EPSQ2 Param 16 R(4) 4 scalar 80,131,133,140 EPSTRB Param 16 R(4) 4 scalar 76 EQOL2 Local 74 R(4) 4 scalar 74,78,98 EXP Func 108 scalar 108,123 G Param 20 R(4) 4 scalar 20 GAM1 Param 21 R(4) 4 scalar GH Dummy 2 R(4) 4 1 44 ARG,INOUT 69 GHL Local 69 R(4) 4 scalar 69,71,72,76,77,85,86,88,89 GM Dummy 2 R(4) 4 1 44 ARG,INOUT 68 GML Local 68 R(4) 4 scalar 68,71,72,76,77,85,86,88,89 IDIM1 Param 19 I(4) 4 scalar IDIM2 Param 19 I(4) 4 scalar IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 IM Param 3 I(4) 4 scalar 16,19 INPES Param 11 I(4) 4 scalar 16,19 ITAIL Param 16 I(4) 4 scalar JDIM1 Param 20 I(4) 4 scalar JDIM2 Param 20 I(4) 4 scalar JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20 JNPES Param 11 I(4) 4 scalar 17,20 JTAIL Param 17 I(4) 4 scalar L Local 67 I(4) 4 scalar 67,68,69,80,81,96,130,131,133 LM Param 3 I(4) 4 scalar 13,60 LM1 Param 13 I(4) 4 scalar 62 LMHK Dummy 2 I(4) 4 scalar ARG,INOUT 65,140 LMHM Local 65 I(4) 4 scalar 65,67 LSM Param 3 I(4) 4 scalar PRODQ2 Subr 2 PRT Param 21 R(4) 4 scalar Q2 Dummy 2 R(4) 4 1 45 ARG,INOUT 80,96,130,131,133,140 RB1 Param 39 R(4) 4 scalar 119 RDEN1 Local 104 R(4) 4 scalar 104,106 RDEN2 Local 117 R(4) 4 scalar 117,119,120 REQU Param 57 R(4) 4 scalar 77 RHS2 Local 119 R(4) 4 scalar 119,121 RHSP1 Local 106 R(4) 4 scalar 106,108 RHSP2 Local 120 R(4) 4 scalar 120,121,123 RHST2 Local 121 R(4) 4 scalar 121,123 SQRT Func 74 scalar 74,96,99 USTAR Dummy 2 R(4) 4 scalar ARG,INOUT 140 Page 6 Source Listing PRODQ2 2014-12-17 20:47 Subprograms/Common Blocks PRODQ2.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References PRODQ2 Subr 2 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -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 Page 7 Source Listing PRODQ2 2014-12-17 20:47 PRODQ2.F 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 no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : PRODQ2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100