Page 1 Source Listing MIXLEN 2014-12-17 20:46 /tmpnwprd/ifortQqwMqK.i 1 # 1 "MIXLEN.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE MIXLEN(LMHK,LPBL,HPBL,U,V,T,Q,Q2,APE,Z,GM,GH,EL) 4 C ****************************************************************** 5 C * * 6 C * LEVEL 2.5 MIXING LENGTH * 7 C * * 8 C ****************************************************************** 9 c 10 c July 1997: Modified to restore averaging of layer values of l, 11 c ELL, a la Mesinger 1993, Res. Act. Atmos. Ocean. Mod., No. 18, 12 c 4.36-4.38; 13 c 14 c A problem removed which may have led to the Òabove PBLÓ scheme 15 c for preliminary EL to be used inadvertently at times and places 16 c where the PBL scheme was supposed to have been used 17 c 18 C----------------------------------------------------------------------- 19 INCLUDE "parmeta" 42 INCLUDE "mpp.h" 43 # 1 "./sp.h" 1 44 # 4 45 46 # 21 "MIXLEN.F" 2 47 C----------------------------------------------------------------------- 106 P A R A M E T E R 107 &(LP1=LM+1,LM1=LM-1) 108 C----------------------------------------------------------------------- 109 P A R A M E T E R 110 &(EPSQ2=0.2,EPSL=0.32,EPSRU=1.E-12,EPSRS=1.E-12 111 &,FH=1.01,ALPH=.20,VKRM=.40,ELFC=0.23 ,EL0MAX=1000.) 112 C----------------------------------------------------------------------- 113 P A R A M E T E R 114 C----------------------------------------------------------------------- 115 &(G=9.8,BETA=1./270.,BTG=BETA*G 116 &,PRT=1.0,GAM1=0.2222222222222222222 117 &,A1=0.659888514560862645,A2=0.6574209922667784586 118 &,B1=11.87799326209552761,B2=7.226971804046074028 119 &,C1=0.000830955950095854396) 120 C----------------------------------------------------------------------- 121 CUN &(G=9.8,BETA=1./270.,BTG=BETA*G 122 CUN &,PRT=1.0,GAM1=0.2222222222222222222 123 CUN &,A1=0.3310949523016403346,A2=0.8273378704055731278 124 CUN &,B1=5.959709141429526024,B2=3.626088092074591135 125 CUN &,C1=-0.3330651924968952113) 126 C----------------------------------------------------------------------- 127 CMY &(G=9.8,BETA=1./270.,BTG=BETA*G 128 CMY &,PRT=0.8,GAM1=0.2222222222222222222 129 CMY &,A1=0.9222222350809054114,A2=0.7350190142719400952 130 CMY &,B1=16.60000023145629741,B2=10.10000014082581951 131 CMY &,C1=0.0805318118080613468) 132 C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR------------- 133 P A R A M E T E R 134 &(ADNM=18.*A1*A1*A2*(B2-3.*A2)*BTG 135 &,ADNH= 9.*A1*A2*A2*(12.*A1+3.*B2)*BTG*BTG 136 &,BDNM= 6.*A1*A1 137 &,BDNH= 3.*A2*(7.*A1+B2)*BTG Page 2 Source Listing MIXLEN 2014-12-17 20:46 MIXLEN.F 138 C--------------FREE TERM IN THE EQUILIBRIUM EQUATION FOR (L/Q)**2------- 139 &,AEQM= 3.*A1*A2*B1*(3.*A2+3.*B2*C1+18.*A1*C1-B2)*BTG 140 & +18.*A1*A1*A2*(B2-3.*A2)*BTG 141 &,AEQH= 9.*A1*A2*A2*B1*BTG*BTG+9.*A1*A2*A2*(12.*A1+3.*B2)*BTG*BTG 142 C--------------FORBIDDEN TURBULENCE AREA-------------------------------- 143 &,REQU=-AEQH/AEQM*1.02,EPSGH=1.E-9,EPSGM=REQU*EPSGH 144 C--------------NEAR ISOTROPY FOR SHEAR TURBULENCE, WW/Q2 LOWER LIMIT---- 145 &,UBRYL=(18.*REQU*A1*A1*A2*B2*C1*BTG+9.*A1*A2*A2*B2*BTG*BTG) 146 & /(REQU*ADNM+ADNH) 147 &,UBRY=(1.+EPSRS)*UBRYL 148 &,UBRY3=3.*UBRY 149 &,AUBM=54.*A1*A1*A2*B2*C1*BTG -ADNM*UBRY3 150 &,AUBH=27.*A1*A2*A2*B2*BTG*BTG-ADNH*UBRY3 151 &,BUBM=18.*A1*A1*C1 -BDNM*UBRY3 152 &,BUBH=(9.*A1*A2+3.*A2*B2)*BTG-BDNH*UBRY3 153 &,CUBR=1. - UBRY3,RCUBR=1./CUBR) 154 C----------------------------------------------------------------------- 155 D I M E N S I O N 156 & U (LM),V (LM) 157 &,T (LM),Q (LM),Q2 (LM) 158 D I M E N S I O N 159 & GM (LM1),GH (LM1),EL (LM1) 160 &,APE (LM ) 161 &,Z (LP1) 162 C----------------------------------------------------------------------- 163 D I M E N S I O N 164 & ELM (LM1),ELL (LM ) 165 &,Q1 (LM ),THV (LM ) 166 C----------------------------------------------------------------------- 167 C*********************************************************************** 168 LMHM=LMHK-1 169 LMHP=LMHK+1 170 C--------------FIND THE HEIGHT OF THE PBL------------------------------- 171 LPBL=LMHK 172 DO 100 IVI=1,LMHK 173 L=LMHK-IVI 174 IF(Q2(L).LE.EPSQ2*FH)THEN 175 CVVVVVVVVVVVVVV NOT NECESSARY IF DRIVEN BY TURBL VVVVVVVVVVVVVVVVVVVVVVV 176 C Q2(L)=EPSQ2 177 CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 178 LPBL=L 179 GO TO 110 180 ENDIF 181 100 CONTINUE 182 LPBL=1 183 C--------------THE HEIGHT OF THE PBL------------------------------------ 184 110 HPBL=Z(LPBL)-Z(LMHP) 185 C----------------------------------------------------------------------- 186 DO 120 L=1,LMHK 187 Q1(L)=0. 188 THV(L)=(0.608*Q(L)+1.)*T(L)*APE(L) 189 120 CONTINUE 190 C----------------------------------------------------------------------- 191 DO 130 L=LPBL,LMHK 192 Q1(L)=SQRT(Q2(L)) 193 130 CONTINUE 194 C----------------------------------------------------------------------- Page 3 Source Listing MIXLEN 2014-12-17 20:46 MIXLEN.F 195 SZQ=0. 196 SQ =0. 197 DO 140 L=1,LMHM 198 QDZL=(Q1(L)+Q1(L+1))*(Z(L+1)-Z(L+2)) 199 SZQ=(Z(L+1)+Z(L+2)-Z(LMHP)-Z(LMHP))*QDZL+SZQ 200 SQ=QDZL+SQ 201 RDZ=2./(Z(L)-Z(L+2)) 202 GML=((U(L)-U(L+1))**2+(V(L)-V(L+1))**2)*RDZ*RDZ 203 GM(L)=AMAX1(GML,EPSGM) 204 GHL=(THV(L)-THV(L+1))*RDZ 205 IF(ABS(GHL).LE.EPSGH) GHL=EPSGH 206 GH(L)=GHL 207 140 CONTINUE 208 C--------------COMPUTATION OF ASYMPTOTIC L IN BLACKADAR FORMULA--------- 209 EL0=AMIN1(ALPH*SZQ*0.5/SQ,EL0MAX) 210 C----------------------------------------------------------------------- 211 DO 150 L=1,LMHM 212 GML=GM(L) 213 GHL=GH(L) 214 IF(GHL.GE.EPSGH)THEN 215 IF(GML/GHL.LE.REQU)THEN 216 ELM(L)=EPSL 217 ELSE 218 AUBR=(AUBM*GML+AUBH*GHL)*GHL 219 BUBR= BUBM*GML+BUBH*GHL 220 QOL2ST=(-0.5*BUBR+SQRT(BUBR*BUBR*0.25-AUBR*CUBR))*RCUBR 221 QOL2ST=AMAX1(QOL2ST,1.E-8) 222 ELOQ2X=1./QOL2ST 223 ELM(L)=AMAX1(SQRT(ELOQ2X*Q2(L)),EPSL) 224 ENDIF 225 ELSE 226 ADEN=(ADNM*GML+ADNH*GHL)*GHL 227 BDEN= BDNM*GML+BDNH*GHL 228 QOL2UN= -0.5*BDEN+SQRT(BDEN*BDEN*0.25-ADEN) 229 ELOQ2X=1./((1.+EPSRU)*QOL2UN) 230 ELM(L)=AMAX1(SQRT(ELOQ2X*Q2(L)),EPSL) 231 ENDIF 232 150 CONTINUE 233 C----------------------------------------------------------------------- 234 c LPBLM=LPBL-1 235 c DO 160 L=1,LPBLM 236 c EL(L)=AMIN1((Z(L)-Z(L+2))*ELFC,ELM(L)) 237 c160 CONTINUE 238 c DO 165 L=LPBL,LMHM 239 c VKRMZ=(Z(L+1)-Z(LMHP))*VKRM 240 c165 EL(L)=AMIN1(VKRMZ/(VKRMZ/EL0+1.),ELM(L)) 241 C----------------------------------------------------------------------- 242 c Note: LPBL is the EL value of L of the interface starting with the 243 c lowest interface and going up for the first time having Q2 .le. 244 c a specified value 245 c 246 DO 160 L=1,LPBL 247 ELL(L)= (Z(L)-Z(L+1))*ELFC 248 160 CONTINUE 249 IF(LPBL.LT.LMHK)THEN 250 LPBLP=LPBL+1 251 DO 165 L=LPBLP,LMHK Page 4 Source Listing MIXLEN 2014-12-17 20:46 MIXLEN.F 252 VKRMZ=(0.5*(Z(L)+Z(L+1))-Z(LMHP))*VKRM 253 ELL(L)= VKRMZ/(VKRMZ/EL0+1.) 254 165 CONTINUE 255 ENDIF 256 C----------------------------------------------------------------------- 257 DO 260 L=1,LMHM 258 EL(L)=AMIN1(0.5*(ELL(L)+ELL(L+1)),ELM(L)) 259 260 CONTINUE 260 C----------------------------------------------------------------------- 261 RETURN 262 END ENTRY POINTS Name mixlen_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 97 88 110 Label 100 95 120 Label 105 102 130 Label 109 107 140 Label 123 113 150 Label 148 127 160 Label 164 162 165 Label 170 167 260 Label 175 173 A1 Param 33 R(4) 4 scalar 50,51,52,53,55,56,57,61,65,66,67,6 8 A2 Param 33 R(4) 4 scalar 50,51,53,55,56,57,61,65,66,68 ABS Func 121 scalar 121 ADEN Local 142 R(4) 4 scalar 142,144 ADNH Param 51 R(4) 4 scalar 62,66,142 ADNM Param 50 R(4) 4 scalar 62,65,142 AEQH Param 57 R(4) 4 scalar 59 AEQM Param 55 R(4) 4 scalar 59 ALPH Param 27 R(4) 4 scalar 125 AMAX1 Func 119 scalar 119,137,139,146 AMIN1 Func 125 scalar 125,174 APE Dummy 2 R(4) 4 1 45 ARG,INOUT 104 AUBH Param 66 R(4) 4 scalar 134 AUBM Param 65 R(4) 4 scalar 134 AUBR Local 134 R(4) 4 scalar 134,136 B1 Param 34 R(4) 4 scalar 55,57 B2 Param 34 R(4) 4 scalar 50,51,53,55,56,57,61,65,66,68 BDEN Local 143 R(4) 4 scalar 143,144 BDNH Param 53 R(4) 4 scalar 68,143 BDNM Param 52 R(4) 4 scalar 67,143 BETA Param 31 R(4) 4 scalar 31 BTG Param 31 R(4) 4 scalar 50,51,53,55,56,57,61,65,66,68 BUBH Param 68 R(4) 4 scalar 135 Page 5 Source Listing MIXLEN 2014-12-17 20:46 Symbol Table MIXLEN.F Name Object Declared Type Bytes Dimen Elements Attributes References BUBM Param 67 R(4) 4 scalar 135 BUBR Local 135 R(4) 4 scalar 135,136 C1 Param 35 R(4) 4 scalar 55,61,65,67 CUBR Param 69 R(4) 4 scalar 69,136 EL Dummy 2 R(4) 4 1 44 ARG,INOUT 174 EL0 Local 125 R(4) 4 scalar 125,169 EL0MAX Param 27 R(4) 4 scalar 125 ELFC Param 27 R(4) 4 scalar 163 ELL Local 80 R(4) 4 1 45 163,169,174 ELM Local 80 R(4) 4 1 44 132,139,146,174 ELOQ2X Local 138 R(4) 4 scalar 138,139,145,146 EPSGH Param 59 R(4) 4 scalar 59,121,130 EPSGM Param 59 R(4) 4 scalar 119 EPSL Param 26 R(4) 4 scalar 132,139,146 EPSQ2 Param 26 R(4) 4 scalar 90 EPSRS Param 26 R(4) 4 scalar 63 EPSRU Param 26 R(4) 4 scalar 145 FH Param 27 R(4) 4 scalar 90 G Param 31 R(4) 4 scalar 31 G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GAM1 Param 32 R(4) 4 scalar GH Dummy 2 R(4) 4 1 44 ARG,INOUT 122,129 GHL Local 120 R(4) 4 scalar 120,121,122,129,130,131,134,135,14 2,143 GLB_TABLE Common 37 128 GM Dummy 2 R(4) 4 1 44 ARG,INOUT 119,128 GML Local 118 R(4) 4 scalar 118,119,128,131,134,135,142,143 HPBL Dummy 2 R(4) 4 scalar ARG,INOUT 100 HTMG Scalar 52 R(4) 4 3 4183695 COM IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDIM1 Param 19 I(4) 4 scalar IDIM2 Param 19 I(4) 4 scalar IE_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 ILCOL Scalar 22 I(4) 4 scalar COM ILPAD1 Scalar 23 I(4) 4 scalar COM ILPAD2 Scalar 23 I(4) 4 scalar COM ILPAD3 Scalar 23 I(4) 4 scalar COM ILPAD4 Scalar 23 I(4) 4 scalar COM ILPAD5 Scalar 23 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 16,19,43,44,45,46,47,51,52,57 INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INUMQ Scalar 34 I(4) 4 1 100 COM IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IRCOL Scalar 22 I(4) 4 scalar COM IRPAD1 Scalar 24 I(4) 4 scalar COM IRPAD2 Scalar 24 I(4) 4 scalar COM IRPAD3 Scalar 24 I(4) 4 scalar COM IRPAD4 Scalar 24 I(4) 4 scalar COM IRPAD5 Scalar 24 I(4) 4 scalar COM IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM Page 6 Source Listing MIXLEN 2014-12-17 20:46 Symbol Table MIXLEN.F Name Object Declared Type Bytes Dimen Elements Attributes References IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM ITAIL Param 16 I(4) 4 scalar ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITROW Scalar 22 I(4) 4 scalar COM IVI Local 88 I(4) 4 scalar 88,89 JBPAD1 Scalar 25 I(4) 4 scalar COM JBPAD2 Scalar 25 I(4) 4 scalar COM JBPAD3 Scalar 25 I(4) 4 scalar COM JBPAD4 Scalar 25 I(4) 4 scalar COM JBPAD5 Scalar 25 I(4) 4 scalar COM JDIM1 Param 20 I(4) 4 scalar JDIM2 Param 20 I(4) 4 scalar JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM 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,43,44,45,46,47,51,52,57 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JS_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM JTAIL Param 17 I(4) 4 scalar JTPAD1 Scalar 26 I(4) 4 scalar COM JTPAD2 Scalar 26 I(4) 4 scalar COM JTPAD3 Scalar 26 I(4) 4 scalar COM JTPAD4 Scalar 26 I(4) 4 scalar COM JTPAD5 Scalar 26 I(4) 4 scalar COM L Local 89 I(4) 4 scalar 89,90,94,102,103,104,107,108,113,1 14,115,117,118,119,120,122,127,128 ,129,132,139,146,162,163,167,168,1 69,173,174 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LM Param 3 I(4) 4 scalar 23,52,72,73,76,80,81 LM1 Param 23 I(4) 4 scalar 75,80 LMHK Dummy 2 I(4) 4 scalar ARG,INOUT 84,85,87,88,89,102,107,165,167 LMHM Local 84 I(4) 4 scalar 84,113,127,173 LMHP Local 85 I(4) 4 scalar 85,100,115,168 LP1 Param 23 I(4) 4 scalar 77 LPBL Dummy 2 I(4) 4 scalar ARG,INOUT 87,94,98,100,107,162,165,166 LPBLP Local 166 I(4) 4 scalar 166,167 LSM Param 3 I(4) 4 scalar MAPPINGS Common 56 5024 MIXLEN Subr 2 MPI_COMM_COMP Scalar 30 I(4) 4 scalar COM MPI_COMM_INTER Scalar 31 I(4) 4 scalar COM MPI_COMM_INTER_ARRAY Scalar 32 I(4) 4 1 100 COM MPPCOM Common 1 1464 MYIE Scalar 4 I(4) 4 scalar COM MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM MYIE1_P3 Scalar 8 I(4) 4 scalar COM MYIE1_P4 Scalar 8 I(4) 4 scalar COM MYIE2 Scalar 4 I(4) 4 scalar COM Page 7 Source Listing MIXLEN 2014-12-17 20:46 Symbol Table MIXLEN.F Name Object Declared Type Bytes Dimen Elements Attributes References MYIE2_P1 Scalar 9 I(4) 4 scalar COM MYIE_P1 Scalar 7 I(4) 4 scalar COM MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM MYIE_P4 Scalar 7 I(4) 4 scalar COM MYIE_P5 Scalar 7 I(4) 4 scalar COM MYIS Scalar 4 I(4) 4 scalar COM MYIS1 Scalar 4 I(4) 4 scalar COM MYIS1_P1 Scalar 6 I(4) 4 scalar COM MYIS1_P2 Scalar 6 I(4) 4 scalar COM MYIS1_P3 Scalar 6 I(4) 4 scalar COM MYIS1_P4 Scalar 6 I(4) 4 scalar COM MYIS2 Scalar 4 I(4) 4 scalar COM MYIS_P1 Scalar 5 I(4) 4 scalar COM MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM MYIS_P4 Scalar 5 I(4) 4 scalar COM MYIS_P5 Scalar 5 I(4) 4 scalar COM MYJE Scalar 15 I(4) 4 scalar COM MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM MYJE2_P1 Scalar 18 I(4) 4 scalar COM MYJE2_P2 Scalar 18 I(4) 4 scalar COM MYJE2_P3 Scalar 18 I(4) 4 scalar COM MYJE2_P4 Scalar 18 I(4) 4 scalar COM MYJE3 Scalar 15 I(4) 4 scalar COM MYJE3_P4 Scalar 18 I(4) 4 scalar COM MYJE4 Scalar 15 I(4) 4 scalar COM MYJE4_P1 Scalar 19 I(4) 4 scalar COM MYJE4_P4 Scalar 19 I(4) 4 scalar COM MYJE5 Scalar 15 I(4) 4 scalar COM MYJE5_P1 Scalar 19 I(4) 4 scalar COM MYJE5_P2 Scalar 19 I(4) 4 scalar COM MYJE_P1 Scalar 16 I(4) 4 scalar COM MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM MYJE_P4 Scalar 16 I(4) 4 scalar COM MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM MYJS1_P3 Scalar 12 I(4) 4 scalar COM MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM MYJS2_P1 Scalar 13 I(4) 4 scalar COM MYJS2_P2 Scalar 13 I(4) 4 scalar COM MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM MYJS3 Scalar 10 I(4) 4 scalar COM MYJS3_P4 Scalar 13 I(4) 4 scalar COM Page 8 Source Listing MIXLEN 2014-12-17 20:46 Symbol Table MIXLEN.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJS4 Scalar 10 I(4) 4 scalar COM MYJS4_P1 Scalar 14 I(4) 4 scalar COM MYJS4_P4 Scalar 14 I(4) 4 scalar COM MYJS5 Scalar 10 I(4) 4 scalar COM MYJS5_P1 Scalar 14 I(4) 4 scalar COM MYJS5_P2 Scalar 14 I(4) 4 scalar COM MYJS_P1 Scalar 11 I(4) 4 scalar COM MYJS_P2 Scalar 11 I(4) 4 scalar COM MYJS_P3 Scalar 11 I(4) 4 scalar COM MYJS_P4 Scalar 11 I(4) 4 scalar COM MYJS_P5 Scalar 11 I(4) 4 scalar COM MYPE Scalar 2 I(4) 4 scalar COM MY_E Scalar 20 I(4) 4 scalar COM MY_IE_GLB Scalar 2 I(4) 4 scalar COM MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM MY_JS_LOC Scalar 3 I(4) 4 scalar COM MY_N Scalar 20 I(4) 4 scalar COM MY_NE Scalar 21 I(4) 4 scalar COM MY_NEB Scalar 21 I(4) 4 1 8 COM MY_NW Scalar 21 I(4) 4 scalar COM MY_S Scalar 20 I(4) 4 scalar COM MY_SE Scalar 21 I(4) 4 scalar COM MY_SW Scalar 21 I(4) 4 scalar COM MY_W Scalar 20 I(4) 4 scalar COM NPES Scalar 2 I(4) 4 scalar COM PRT Param 32 R(4) 4 scalar Q Dummy 2 R(4) 4 1 45 ARG,INOUT 104 Q1 Local 81 R(4) 4 1 45 103,108,114 Q2 Dummy 2 R(4) 4 1 45 ARG,INOUT 90,108,139,146 QDZL Local 114 R(4) 4 scalar 114,115,116 QOL2ST Local 136 R(4) 4 scalar 136,137,138 QOL2UN Local 144 R(4) 4 scalar 144,145 RCUBR Param 69 R(4) 4 scalar 136 RDZ Local 117 R(4) 4 scalar 117,118,120 REQU Param 59 R(4) 4 scalar 59,61,62,131 SQ Local 112 R(4) 4 scalar 112,116,125 SQRT Func 108 scalar 108,136,139,144,146 SZQ Local 111 R(4) 4 scalar 111,115,125 T Dummy 2 R(4) 4 1 45 ARG,INOUT 104 TEMP1 Scalar 43 R(4) 4 2 91719 COM TEMP10 Scalar 45 R(4) 4 2 91719 COM TEMP11 Scalar 45 R(4) 4 2 91719 COM TEMP12 Scalar 45 R(4) 4 2 91719 COM TEMP13 Scalar 46 R(4) 4 2 91719 COM TEMP14 Scalar 46 R(4) 4 2 91719 COM TEMP15 Scalar 46 R(4) 4 2 91719 COM TEMP16 Scalar 46 R(4) 4 2 91719 COM TEMP2 Scalar 43 R(4) 4 2 91719 COM TEMP2X Scalar 51 R(4) 4 2 92971 COM TEMP3 Scalar 43 R(4) 4 2 91719 COM Page 9 Source Listing MIXLEN 2014-12-17 20:46 Symbol Table MIXLEN.F Name Object Declared Type Bytes Dimen Elements Attributes References TEMP4 Scalar 43 R(4) 4 2 91719 COM TEMP5 Scalar 44 R(4) 4 2 91719 COM TEMP6 Scalar 44 R(4) 4 2 91719 COM TEMP7 Scalar 44 R(4) 4 2 91719 COM TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 THV Local 81 R(4) 4 1 45 104,120 TOPO Common 50 17478548 TTVG Scalar 51 R(4) 4 2 92971 COM U Dummy 2 R(4) 4 1 45 ARG,INOUT 118 UBRY Param 63 R(4) 4 scalar 64 UBRY3 Param 64 R(4) 4 scalar 65,66,67,68,69 UBRYL Param 61 R(4) 4 scalar 63 V Dummy 2 R(4) 4 1 45 ARG,INOUT 118 VKRM Param 27 R(4) 4 scalar 168 VKRMZ Local 168 R(4) 4 scalar 168,169 Z Dummy 2 R(4) 4 1 46 ARG,INOUT 100,114,115,117,163,168 Page 10 Source Listing MIXLEN 2014-12-17 20:46 Subprograms/Common Blocks MIXLEN.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 37 128 MAPPINGS Common 56 5024 MIXLEN Subr 2 MPPCOM Common 1 1464 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 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 Page 11 Source Listing MIXLEN 2014-12-17 20:46 MIXLEN.F 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 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 : MIXLEN.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100