Page 1 Source Listing LWR88 2014-12-17 20:46 /tmpnwprd/ifort2A2wNW.i 1 # 1 "LWR88.F" 2 SUBROUTINE LWR88(HEATRA,GRNFLX,TOPFLX, 3 1 PRESS,TEMP,RH2O,QO3,CLDFAC, 4 2 CAMT,NCLDS,KTOP,KBTM) 5 C--------------------------------------------------------------------- 6 C SUBROUTINE LWR88 COMPUTES TEMPERATURE-CORRECTED CO2 TRANSMISSION 7 C FUNCTIONS AND ALSO COMPUTES THE PRESSURE GRID AND LAYER OPTICAL 8 C PATHS. 9 C INPUTS: (COMMON BLOCKS) 10 C CLDFAC CLDCOM 11 C PRESS,TEMP,RH2O,QO3 RADISW 12 C CAMT,NCLDS,KTOP,KBTM RADISW 13 C CO251,CO258,CDT51,CDT58 CO2BD3 14 C C2D51,C2D58,CO2M51,CO2M58 CO2BD3 15 C CDTM51,CDTM58,C2DM51,C2DM58 CO2BD3 16 C STEMP,GTEMP CO2BD3 17 C CO231,CO238,CDT31,CDT38 CO2BD2 18 C C2D31,C2D38 CO2BD2 19 C CO271,CO278,CDT71,CDT78 CO2BD4 20 C C2D71,C2D78 CO2BD4 21 C BETINW BDWIDE 22 C OUTPUTS: 23 C HEATRA,GRNFLX,TOPFLX LWOUT 24 C CALLED BY: 25 C RADMN OR INPUT ROUTINE OF MODEL 26 C CALLS: 27 C FST88 28 C 29 COMMON/PHYCON/AMOLWT,CSUBP,DIFFCTR,G,GRAVDR,O3DIFCTR,P0, 30 * P0XZP2,P0XZP8,P0X2,RADCON,RGAS,RGASSP,SECPDA 31 COMMON/PHYCON/RATCO2MW,RATH2OMW 32 COMMON/PHYCON/RADCON1 33 COMMON/PHYCON/GINV,P0INV,GP0INV 34 COMMON/HCON/HUNDRED,HNINETY,SIXTY,FIFTY,TEN,EIGHT,FIVE, 35 * FOUR,THREE,TWO,ONE,HAF,QUARTR,ZERO 36 COMMON/HCON/H83E26,H71E26,H1E15,H1E13,H1E11,H1E8,H4E5, 37 * H165E5,H5725E4,H488E4,H1E4,H24E3,H20788E3, 38 * H2075E3,H1224E3,H5E2,H3082E2,H3E2,H2945E2, 39 * H23E2,H15E2,H35E1,H3P6,H181E1,H18E1,H2P9,H2P8, 40 * H2P5,H1P8,H1P4387,H1P4,H1P25892,HP8,HP518, 41 * HP369,HP1 42 COMMON/HCON/H44871M2,H559M3,H1M3,H987M4,H285M4,H1M4, 43 * H6938M5,H394M5,H37412M5,H1439M5,H128M5,H1M5, 44 * H7M6,H4999M6,H25452M6,H1M6,H391M7,H1174M7, 45 * H8725M8,H327M8,H257M8,H1M8,H23M10,H14M10, 46 * H11M10,H1M10,H83M11,H82M11,H8M11,H77M11, 47 * H72M11,H53M11,H48M11,H44M11,H42M11,H37M11, 48 * H35M11,H32M11,H3M11,H28M11,H24M11,H23M11, 49 * H2M11,H18M11,H15M11,H14M11,H114M11,H11M11, 50 * H1M11,H96M12,H93M12,H77M12,H74M12,H65M12, 51 * H62M12,H6M12,H45M12,H44M12,H4M12,H38M12, 52 * H37M12,H3M12,H29M12,H28M12,H24M12,H21M12, 53 * H16M12,H14M12,H12M12,H8M13,H46M13,H36M13, 54 * H135M13,H12M13,H1M13,H3M14,H15M14,H14M14, 55 * H1M17,H1M18,H1M19,H1M20,H1M21,H1M22,H1M23, 56 * H1M24,H26M30,H14M30,H25M31,H21M31,H12M31, 57 * H9M32,H55M32,H45M32,H4M33,H62M34,H1M60 Page 2 Source Listing LWR88 2014-12-17 20:46 LWR88.F 58 COMMON/HCON/HMP575,HM13EZ,HM19EZ,HM1E1,HM181E1,HM1E2 59 COMMON/HCON/H1E6,H2E6,H1M2,HMP66667,HM6666M2,HP166666, 60 * H41666M2,HMP5,HM2M2,H29316E2,H1226E1,H3116E1, 61 * H9P94,HP6,H625M2,HP228,HP60241,HM1797E1, 62 * H8121E1,H2E2,HM1EZ,H26E2,H44194M2,H1P41819 63 COMMON/HCON/HP219,HP144,HP816,H69766E5,H235M3,HP26, 64 * H129M2,H75826M4,H1P082,HP805,H1386E2, 65 * H658M2,H1036E2,H2118M2,H42M2,H323M4, 66 * H67390E2,HP3795,HP5048,H102M5,H451M6 67 COMMON/HCON/H16E1,HM161E1,H161E1,H3M3,H101M16, 68 * HM1597E1,H25E2,HP118666,H15M5,H3P5,H18E3, 69 * H6P08108,HMP805,HP602409,HP526315, 70 * H28571M2,H1M16 71 COMMON/HCON/H3M4 72 COMMON/HCON/HM8E1 73 COMMON/HCON/H28E1 74 C----------------------------------------------------------------------- 75 INCLUDE "parmeta" 98 INCLUDE "mpp.h" 99 # 1 "./sp.h" 1 100 # 4 101 102 # 77 "LWR88.F" 2 103 C----------------------------------------------------------------------- 104 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 105 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 106 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 107 C BANDTA FOR DEFINITION 108 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 109 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE 110 C BDCOMB FOR DEFINITION 111 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 112 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 113 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 114 C FROM THE ABOVE PARAMETERS. 173 PARAMETER (L=LM) 174 PARAMETER (NBLW=163,NBLX=47,NBLY=15) 175 PARAMETER (NBLM=NBLY-1) 176 PARAMETER (LP1=L+1,LP2=L+2,LP3=L+3) 177 PARAMETER (LM1=L-1,LM2=L-2,LM3=L-3) 178 PARAMETER (LL=2*L,LLP1=LL+1,LLP2=LL+2,LLP3=LL+3) 179 PARAMETER (LLM1=LL-1,LLM2=LL-2,LLM3=LL-3) 180 PARAMETER (LP1M=LP1*LP1,LP1M1=LP1M-1) 181 PARAMETER (LP1V=LP1*(1+2*L/2)) 182 PARAMETER (LP121=LP1*NBLY) 183 PARAMETER (LL3P=3*L+2) 184 PARAMETER (NB=12) 185 PARAMETER (INLTE=3,INLTEP=INLTE+1,NNLTE=56) 186 PARAMETER (NB1=NB-1) 187 PARAMETER (KO2=12) 188 PARAMETER (KO21=KO2+1,KO2M=KO2-1) 189 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 190 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 191 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 192 C BANDTA FOR DEFINITION 193 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 194 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE Page 3 Source Listing LWR88 2014-12-17 20:46 LWR88.F 195 C BDCOMB FOR DEFINITION 196 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 197 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 198 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 199 C FROM THE ABOVE PARAMETERS. 200 C THE FOLLOWING COMMON BLOCKS CONTAIN PRETABULATED CO2 TRANSMISSION 201 C FUNCTIONS, EVALUATED USING THE METHODS OF FELS AND 202 C SCHWARZKOPF (1981) AND SCHWARZKOPF AND FELS (1985), 203 C***COMMON CO2BD3 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE 204 C AND PRESSURE DERIVATIVES FOR THE 560-800 CM-1 BAND. ALSO INCLUDED 205 C ARE THE STANDARD TEMPERATURES AND THE WEIGHTING FUNCTION. THESE 206 C DATA ARE IN BLOCK DATA BD3: 207 C CO251 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 208 C WITH P(SFC)=1013.25 MB 209 C CO258 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 210 C WITH P(SFC)= ^810 MB 211 C CDT51 = FIRST TEMPERATURE DERIVATIVE OF CO251 212 C CDT58 = FIRST TEMPERATURE DERIVATIVE OF CO258 213 C C2D51 = SECOND TEMPERATURE DERIVATIVE OF CO251 214 C C2D58 = SECOND TEMPERATURE DERIVATIVE OF CO251 215 C CO2M51 = TRANSMISSION FCTNS FOR T0 FOR ADJACENT PRESSURE 216 C LEVELS, WITH NO PRESSURE QUADRATURE. USED FOR 217 C NEARBY LAYER COMPUTATIONS. P(SFC)=1013.25 MB 218 C CO2M58 = SAME AS CO2M51,WITH P(SFC)= ^810 MB 219 C CDTM51 = FIRST TEMPERATURE DERIVATIVE OF CO2M51 220 C CDTM58 = FIRST TEMPERATURE DERIVATIVE OF CO2M58 221 C C2DM51 = SECOND TEMPERATURE DERIVATIVE OF CO2M51 222 C C2DM58 = SECOND TEMPERATURE DERIVATIVE OF CO2M58 223 C STEMP = STANDARD TEMPERATURES FOR MODEL PRESSURE LEVEL 224 C STRUCTURE WITH P(SFC)=1013.25 MB 225 C GTEMP = WEIGHTING FUNCTION FOR MODEL PRESSURE LEVEL 226 C STRUCTURE WITH P(SFC)=1013.25 MB. 227 C B0 = TEMP. COEFFICIENT USED FOR CO2 TRANS. FCTN. 228 C CORRECTION FOR T(K). (SEE REF. 4 AND BD3) 229 C B1 = TEMP. COEFFICIENT, USED ALONG WITH B0 230 C B2 = TEMP. COEFFICIENT, USED ALONG WITH B0 231 C B3 = TEMP. COEFFICIENT, USED ALONG WITH B0 232 C 233 COMMON/CO2BD3/CO251(LP1,LP1),CO258(LP1,LP1),CDT51(LP1,LP1), 234 1 CDT58(LP1,LP1),C2D51(LP1,LP1),C2D58(LP1,LP1),CO2M51(L), 235 2 CO2M58(L),CDTM51(L),CDTM58(L),C2DM51(L),C2DM58(L), 236 3 STEMP(LP1),GTEMP(LP1),B0,B1,B2,B3 237 C 238 C***COMMON CO2BD2 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE 239 C AND PRESSURE DERIVATIVES FOR THE 560-670 CM-1 PART OF THE 15 UM 240 C CO2 BAND. THESE DATA ARE IN BLOCK DATA BD2. 241 C CO231 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 242 C WITH P(SFC)=1013.25 MB 243 C CO238 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 244 C WITH P(SFC)= ^810 MB 245 C CDT31 = FIRST TEMPERATURE DERIVATIVE OF CO231 246 C CDT38 = FIRST TEMPERATURE DERIVATIVE OF CO238 247 C C2D31 = SECOND TEMPERATURE DERIVATIVE OF CO231 248 C C2D38 = SECOND TEMPERATURE DERIVATIVE OF CO231 249 C 250 COMMON / CO2BD2 / CO231(LP1),CO238(LP1),CDT31(LP1), 251 1 CDT38(LP1),C2D31(LP1),C2D38(LP1) Page 4 Source Listing LWR88 2014-12-17 20:46 LWR88.F 252 C 253 C***COMMON CO2BD4 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE 254 C AND PRESSURE DERIVATIVES FOR THE 670-800 CM-1 PART OF THE 15 UM 255 C CO2 BAND. THESE DATA ARE IN BLOCK DATA BD4. 256 C CO271 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 257 C WITH P(SFC)=1013.25 MB 258 C CO278 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 259 C WITH P(SFC)= ^810 MB 260 C CDT71 = FIRST TEMPERATURE DERIVATIVE OF CO271 261 C CDT78 = FIRST TEMPERATURE DERIVATIVE OF CO278 262 C C2D71 = SECOND TEMPERATURE DERIVATIVE OF CO271 263 C C2D78 = SECOND TEMPERATURE DERIVATIVE OF CO271 264 C 265 COMMON / CO2BD4 / CO271(LP1),CO278(LP1),CDT71(LP1), 266 1 CDT78(LP1),C2D71(LP1),C2D78(LP1) 267 C 268 C***COMMON CO2BD5 CONTAINS CO2 TRANSMISSION FUNCTIONS FOR THE 2270- 269 C 2380 PART OF THE 4.3 UM CO2 BAND. THESE DATA ARE IN BLOCK DATA BD5. 270 C CO211 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 271 C WITH P(SFC)=1013.25 MB 272 C CO218 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 273 C WITH P(SFC)= ^810 MB 274 C 275 COMMON / CO2BD5 / CO211(LP1),CO218(LP1) 276 C 277 C COMMON BLOCK BANDTA CONTAINS RANDOM BAND PARAMETERS FOR THE LW 278 C CALCULATIONS USING 10 CM-1 WIDE BANDS.THE 15 UM CO2 COMPLEX 279 C IS 2 BANDS,560-670 AND 670-800 CM-1. OZONE COEFFICIENTS ARE 280 C IN 3 BANDS,670-800 (14.1 UM),990-1070 AND 1070-1200 (9.6 UM). 281 C THE (NBLW) BANDS NOW INCLUDE: 282 C 56 BANDS, 10 CM-1 WIDE 0 - 560 CM-1 283 C 2 BANDS, 15 UM COMPLEX 560 - 670 CM-1 284 C 670 - 800 CM-1 285 C 3 "CONTINUUM" BANDS 800 - 900 CM-1 286 C 900 - 990 CM-1 287 C 1070 - 1200 CM-1 288 C 1 BAND FOR 9.6 UM BAND 990 - 1070 CM-1 289 C 100 BANDS, 10 CM-1 WIDE 1200 - 2200 CM-1 290 C 1 BAND FOR 4.3 UM SRC 2270 - 2380 CM-1 291 C THUS NBLW PRESENTLY EQUALS 163 292 C ALL BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 293 C 294 C ARNDM = RANDOM "A" PARAMETER FOR (NBLW) BANDS 295 C BRNDM = RANDOM "B" PARAMETER FOR (NBLW) BANDS 296 C BETAD = CONTINUUM COEFFICIENTS FOR (NBLW) BANDS 297 C AP,BP = CAPPHI COEFFICIENTS FOR (NBLW) BANDS 298 C ATP,BTP = CAPPSI COEFFICIENTS FOR (NBLW) BANDS 299 C BANDLO = LOWEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 300 C BANDHI = HIGHEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 301 C AO3RND = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 302 C BANDS 303 C BO3RND = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 304 C BANDS 305 C AB15 = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS 306 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 307 C DATA FOR ARNDM,BRNDM,AP,BP,ATP,BTP,AO3RND,BO3RND ARE OBTAINED BY 308 C USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS ARE FROM Page 5 Source Listing LWR88 2014-12-17 20:46 LWR88.F 309 C ROBERTS (1976). 310 COMMON / BANDTA / ARNDM(NBLW),BRNDM(NBLW),BETAD(NBLW),AP(NBLW), 311 1 BP(NBLW),ATP(NBLW),BTP(NBLW),BANDLO(NBLW), 312 2 BANDHI(NBLW),AO3RND(3),BO3RND(3),AB15(2) 313 C 314 C COMMON BLOCK BDWIDE CONTAINS RANDOM BAND PARAMETERS FOR SPECIFIC 315 C WIDE BANDS. AT PRESENT,THE INFORMATION CONSISTS OF 1) RANDOM 316 C MODEL PARAMETERS FOR THE 15 UM BAND,560-800 CM-1; 2) THE 317 C CONTINUUM COEFFICIENT FOR THE 800-990,1070-1200 CM-1 BAND 318 C SPECIFICALLY: 319 C AWIDE = RANDOM "A" PARAMETER FOR BAND 320 C BWIDE = RANDOM "B" PARAMETER FOR BAND 321 C BETAWD = CONTINUUM COEFFICIENTS FOR BAND 322 C APWD,BPWD = CAPPHI COEFFICIENTS FOR BAND 323 C ATPWD,BTPWD = CAPPSI COEFFICIENTS FOR BAND 324 C BDLOWD = LOWEST FREQUENCY IN EACH FREQ BAND 325 C BDHIWD = HIGHEST FREQUENCY IN EACH FREQ BAND 326 C AB15WD = THE PRODUCT ARNDM*BRNDM FOR THE ONE BAND 327 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 328 C BETINW = CONT.COEFFICIENT FOR A SPECIFIED WIDE 329 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 330 C SKO2D = 1./BETINW, USED IN SPA88 FOR CONT. COEFFS 331 C SKC1R = BETAWD/BETINW, USED FOR CONT. COEFF. FOR 332 C 15 UM BAND IN FST88 333 C SKO3R = RATIO OF CONT. COEFF. FOR 9.9 UM BAND TO 334 C BETINW, USED FOR 9.6 UM CONT COEFF IN FST88 335 C DATA FOR AWIDE,BWIDE,APWD,BPWD,ATPWD,BTPWD,AO3WD,BO3WD ARE 336 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 337 C ARE FROM ROBERTS (1976). 338 COMMON / BDWIDE / AWIDE,BWIDE,BETAWD, 339 1 APWD,BPWD,ATPWD,BTPWD, 340 2 BDLOWD,BDHIWD,BETINW, 341 3 AB15WD,SKO2D,SKC1R,SKO3R 342 C 343 C COMMON BLOCK BDCOMB CONTAINS RANDOM BAND PARAMETERS FOR THE LW 344 C CALCULATIONS USING COMBINED WIDE FREQUENCY BANDS BETWEEN 160 AND 345 C 1200 CM-1,AS WELL AS THE 2270-2380 BAND FOR SOURCE CALC. 346 C BANDS 1-8: COMBINED WIDE FREQUENCY BANDS FOR 160-560 CM-1 347 C BANDS 9-14: FREQUENCY BANDS,AS IN BANDTA (NARROW BANDS) 348 C FOR 560-1200 CM-1 349 C BAND 15: FREQUENCY BAND 2270-2380 CM-1,USED FOR SOURCE 350 C CALCULATION ONLY 351 C THUS NBLY PRESENTLY EQUALS 15 352 C 353 C BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 354 C ACOMB = RANDOM "A" PARAMETER FOR (NBLY) BANDS 355 C BCOMB = RANDOM "B" PARAMETER FOR (NBLY) BANDS 356 C BETACM = CONTINUUM COEFFICIENTS FOR (NBLY) BANDS 357 C APCM,BPCM = CAPPHI COEFFICIENTS FOR (NBLY) BANDS 358 C ATPCM,BTPCM = CAPPSI COEFFICIENTS FOR (NBLY) BANDS 359 C BDLOCM = LOWEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 360 C BDHICM = HIGHEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 361 C AO3CM = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 362 C BANDS 363 C BO3CM = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 364 C BANDS 365 C AB15CM = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS Page 6 Source Listing LWR88 2014-12-17 20:46 LWR88.F 366 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 367 C BETINC = CONT.COEFFICIENT FOR A SPECIFIED WIDE 368 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 369 C IBAND = INDEX NO OF THE 40 WIDE BANDS USED IN 370 C COMBINED WIDE BAND CALCULATIONS. IN OTHER 371 C WORDS,INDEX TELLING WHICH OF THE 40 WIDE 372 C BANDS BETWEEN 160-560 CM-1 ARE INCLUDED IN 373 C EACH OF THE FIRST 8 COMBINED WIDE BANDS 374 C DATA FOR ACOMB,BCOMB,APCM,BPCM,ATPCM,BTPCM,AO3CM,BO3CM ARE 375 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 376 C ARE FROM ROBERTS (1976). IBAND INDEX VALUES ARE OBTAINED BY 377 C EXPERIMENTATION. 378 COMMON / BDCOMB / IBAND(40),ACOMB(NBLY),BCOMB(NBLY), 379 1 BETACM(NBLY),APCM(NBLY),BPCM(NBLY),ATPCM(NBLY), 380 2 BTPCM(NBLY),BDLOCM(NBLY),BDHICM(NBLY),BETINC, 381 3 AO3CM(3),BO3CM(3),AB15CM(2) 382 C 383 C 384 DIMENSION PRESS(IDIM1:IDIM2,LP1),TEMP(IDIM1:IDIM2,LP1), 385 1 RH2O(IDIM1:IDIM2,L),QO3(IDIM1:IDIM2,L) 386 DIMENSION CLDFAC(IDIM1:IDIM2,LP1,LP1),CAMT(IDIM1:IDIM2,LP1) 387 DIMENSION NCLDS(IDIM1:IDIM2),KTOP(IDIM1:IDIM2,LP1), 388 1 KBTM(IDIM1:IDIM2,LP1) 389 DIMENSION HEATRA(IDIM1:IDIM2,L),GRNFLX(IDIM1:IDIM2), 390 1 TOPFLX(IDIM1:IDIM2) 391 DIMENSION DELP2(IDIM1:IDIM2,L) 392 C 393 DIMENSION QH2O(IDIM1:IDIM2,L),T(IDIM1:IDIM2,LP1) 394 DIMENSION P(IDIM1:IDIM2,LP1),DELP(IDIM1:IDIM2,L) 395 DIMENSION CO21(IDIM1:IDIM2,LP1,LP1),CO2NBL(IDIM1:IDIM2,L) 396 DIMENSION CO2SP1(IDIM1:IDIM2,LP1),CO2SP2(IDIM1:IDIM2,LP1) 397 DIMENSION VAR1(IDIM1:IDIM2,L),VAR2(IDIM1:IDIM2,L), 398 1 VAR3(IDIM1:IDIM2,L),VAR4(IDIM1:IDIM2,L) 399 DIMENSION CNTVAL(IDIM1:IDIM2,LP1) 400 DIMENSION TOTO3(IDIM1:IDIM2,LP1),TPHIO3(IDIM1:IDIM2,LP1), 401 1 TOTPHI(IDIM1:IDIM2,LP1) 402 DIMENSION TOTVO2(IDIM1:IDIM2,LP1),EMX1(IDIM1:IDIM2), 403 1 EMX2(IDIM1:IDIM2),EMPL(IDIM1:IDIM2,LLP1) 404 C 405 DIMENSION CO2R(IDIM1:IDIM2,LP1),DIFT(IDIM1:IDIM2,LP1) 406 DIMENSION CO2R1(IDIM1:IDIM2,LP1),DCO2D1(IDIM1:IDIM2,LP1) 407 DIMENSION D2CD21(IDIM1:IDIM2,LP1),D2CD22(IDIM1:IDIM2,LP1) 408 DIMENSION CO2R2(IDIM1:IDIM2,LP1),DCO2D2(IDIM1:IDIM2,LP1) 409 DIMENSION CO2MR(IDIM1:IDIM2,L),CO2MD(IDIM1:IDIM2,L), 410 1 CO2M2D(IDIM1:IDIM2,L) 411 DIMENSION TDAV(IDIM1:IDIM2,LP1),TSTDAV(IDIM1:IDIM2,LP1), 412 1 VV(IDIM1:IDIM2,L),VSUM3(IDIM1:IDIM2,LP1),VSUM1(IDIM1:IDIM2), 413 2 VSUM2(IDIM1:IDIM2) 414 DIMENSION A1(IDIM1:IDIM2),A2(IDIM1:IDIM2) 415 DIMENSION DCO2DT(IDIM1:IDIM2,LP1),D2CDT2(IDIM1:IDIM2,LP1) 416 C 417 DIMENSION TEXPSL(IDIM1:IDIM2,LP1),TLSQU(IDIM1:IDIM2,LP1) 418 DIMENSION VSUM4(IDIM1:IDIM2,L) 419 EQUIVALENCE (VSUM3,TLSQU,TEXPSL) 420 EQUIVALENCE (VV,VSUM4) 421 c 422 C Page 7 Source Listing LWR88 2014-12-17 20:46 LWR88.F 423 C****COMPUTE FLUX PRESSURES (P) AND DIFFERENCES (DELP2,DELP) 424 C****COMPUTE FLUX LEVEL TEMPERATURES (T) AND CONTINUUM TEMPERATURE 425 C CORRECTIONS (TEXPSL) 426 DO 103 K=2,L 427 DO 103 I=MYIS,MYIE 428 P(I,K)=HAF*(PRESS(I,K-1)+PRESS(I,K)) 429 T(I,K)=HAF*(TEMP(I,K-1)+TEMP(I,K)) 430 103 CONTINUE 431 DO 105 I=MYIS,MYIE 432 P(I,1)=ZERO 433 P(I,LP1)=PRESS(I,LP1) 434 T(I,1)=TEMP(I,1) 435 T(I,LP1)=TEMP(I,LP1) 436 105 CONTINUE 437 DO 107 K=1,L 438 DO 107 I=MYIS,MYIE 439 DELP2(I,K)=P(I,K+1)-P(I,K) 440 DELP(I,K)=ONE/DELP2(I,K) 441 107 CONTINUE 442 C****COMPUTE ARGUMENT FOR CONT.TEMP.COEFF. 443 C (THIS IS 1800.(1./TEMP-1./296.)) 444 DO 125 K=1,LP1 445 DO 125 I=MYIS,MYIE 446 TEXPSL(I,K)=H18E3/TEMP(I,K)-H6P08108 447 C...THEN TAKE EXPONENTIAL 448 TEXPSL(I,K)=EXP(TEXPSL(I,K)) 449 125 CONTINUE 450 C***COMPUTE OPTICAL PATHS FOR H2O AND O3, USING THE DIFFUSIVITY 451 C APPROXIMATION FOR THE ANGULAR INTEGRATION (1.66). OBTAIN THE 452 C UNWEIGHTED VALUES(VAR1,VAR3) AND THE WEIGHTED VALUES(VAR2,VAR4). 453 C THE QUANTITIES H3M4(.0003) AND H3M3(.003) APPEARING IN THE VAR2 AND 454 C VAR4 EXPRESSIONS ARE THE APPROXIMATE VOIGT CORRECTIONS FOR H2O AND 455 C O3,RESPECTIVELY. 456 C 457 DO 131 K=1,L 458 DO 131 I=MYIS,MYIE 459 QH2O(I,K)=RH2O(I,K)*DIFFCTR 460 C---VV IS THE LAYER-MEAN PRESSURE (IN ATM),WHICH IS NOT THE SAME AS 461 C THE LEVEL PRESSURE (PRESS) 462 VV(I,K)=HAF*(P(I,K+1)+P(I,K))*P0INV 463 VAR1(I,K)=DELP2(I,K)*QH2O(I,K)*GINV 464 VAR3(I,K)=DELP2(I,K)*QO3(I,K)*DIFFCTR*GINV 465 VAR2(I,K)=VAR1(I,K)*(VV(I,K)+H3M4) 466 VAR4(I,K)=VAR3(I,K)*(VV(I,K)+H3M3) 467 C COMPUTE OPTICAL PATH FOR THE H2O CONTINUUM, USING ROBERTS COEFFS. 468 C (BETINW),AND TEMP. CORRECTION (TEXPSL). THE DIFFUSIVITY FACTOR 469 C (WHICH CANCELS OUT IN THIS EXPRESSION) IS ASSUMED TO BE 1.66. THE 470 C USE OF THE DIFFUSIVITY FACTOR HAS BEEN SHOWN TO BE A SIGNIFICANT 471 C SOURCE OF ERROR IN THE CONTINUUM CALCS.,BUT THE TIME PENALTY OF 472 C AN ANGULAR INTEGRATION IS SEVERE. 473 C 474 CNTVAL(I,K)=TEXPSL(I,K)*RH2O(I,K)*VAR2(I,K)*BETINW/ 475 1 (RH2O(I,K)+RATH2OMW) 476 131 CONTINUE 477 C COMPUTE SUMMED OPTICAL PATHS FOR H2O,O3 AND CONTINUUM 478 DO 201 I=MYIS,MYIE 479 TOTPHI(I,1)=ZERO Page 8 Source Listing LWR88 2014-12-17 20:46 LWR88.F 480 TOTO3(I,1)=ZERO 481 TPHIO3(I,1)=ZERO 482 TOTVO2(I,1)=ZERO 483 201 CONTINUE 484 DO 203 K=2,LP1 485 DO 203 I=MYIS,MYIE 486 TOTPHI(I,K)=TOTPHI(I,K-1)+VAR2(I,K-1) 487 TOTO3(I,K)=TOTO3(I,K-1)+VAR3(I,K-1) 488 TPHIO3(I,K)=TPHIO3(I,K-1)+VAR4(I,K-1) 489 TOTVO2(I,K)=TOTVO2(I,K-1)+CNTVAL(I,K-1) 490 203 CONTINUE 491 C---EMX1 IS THE ADDITIONAL PRESSURE-SCALED MASS FROM PRESS(L) TO 492 C P(L). IT IS USED IN NEARBY LAYER AND EMISS CALCULATIONS. 493 C---EMX2 IS THE ADDITIONAL PRESSURE-SCALED MASS FROM PRESS(L) TO 494 C P(LP1). IT IS USED IN CALCULATIONS BETWEEN FLUX LEVELS L AND LP1. 495 C 496 DO 801 I=MYIS,MYIE 497 EMX1(I)=QH2O(I,L)*PRESS(I,L)*(PRESS(I,L)-P(I,L))*GP0INV 498 EMX2(I)=QH2O(I,L)*PRESS(I,L)*(P(I,LP1)-PRESS(I,L))*GP0INV 499 801 CONTINUE 500 C---EMPL IS THE PRESSURE SCALED MASS FROM P(K) TO PRESS(K) (INDEX 2-LP1) 501 C OR TO PRESS(K+1) (INDEX LP2-LL) 502 DO 811 K=1,L 503 DO 811 I=MYIS,MYIE 504 EMPL(I,K+1)=QH2O(I,K)*P(I,K+1)*(P(I,K+1)-PRESS(I,K))*GP0INV 505 811 CONTINUE 506 DO 812 K=1,LM1 507 DO 812 I=MYIS,MYIE 508 EMPL(I,LP2+K-1)=QH2O(I,K+1)*P(I,K+1)*(PRESS(I,K+1)-P(I,K+1)) 509 1 *GP0INV 510 812 CONTINUE 511 DO 821 I=MYIS,MYIE 512 EMPL(I,1)=VAR2(I,L) 513 EMPL(I,LLP1)=EMPL(I,LL) 514 821 CONTINUE 515 C***COMPUTE WEIGHTED TEMPERATURE (TDAV) AND PRESSURE (TSTDAV) INTEGRALS 516 C FOR USE IN OBTAINING TEMP. DIFFERENCE BET. SOUNDING AND STD. 517 C TEMP. SOUNDING (DIFT) 518 DO 161 I=MYIS,MYIE 519 TSTDAV(I,1)=ZERO 520 TDAV(I,1)=ZERO 521 161 CONTINUE 522 DO 162 K=1,LP1 523 DO 162 I=MYIS,MYIE 524 VSUM3(I,K)=TEMP(I,K)-STEMP(K) 525 162 CONTINUE 526 DO 163 K=1,L 527 DO 165 I=MYIS,MYIE 528 VSUM2(I)=GTEMP(K)*DELP2(I,K) 529 VSUM1(I)=VSUM2(I)*VSUM3(I,K) 530 TSTDAV(I,K+1)=TSTDAV(I,K)+VSUM2(I) 531 TDAV(I,K+1)=TDAV(I,K)+VSUM1(I) 532 165 CONTINUE 533 163 CONTINUE 534 C 535 C****EVALUATE COEFFICIENTS FOR CO2 PRESSURE INTERPOLATION (A1,A2) 536 DO 171 I=MYIS,MYIE Page 9 Source Listing LWR88 2014-12-17 20:46 LWR88.F 537 A1(I)=(PRESS(I,LP1)-P0XZP8)/P0XZP2 538 A2(I)=(P0-PRESS(I,LP1))/P0XZP2 539 171 CONTINUE 540 C***PERFORM CO2 PRESSURE INTERPOLATION ON ALL INPUTTED TRANSMISSION 541 C FUNCTIONS AND TEMP. DERIVATIVES 542 C---SUCCESSIVELY COMPUTING CO2R,DCO2DT AND D2CDT2 IS DONE TO SAVE 543 C STORAGE (AT A SLIGHT LOSS IN COMPUTATION TIME) 544 DO 184 K=1,LP1 545 DO 184 I=MYIS,MYIE 546 CO2R1(I,K)=A1(I)*CO231(K)+A2(I)*CO238(K) 547 D2CD21(I,K)=H1M3*(A1(I)*C2D31(K)+A2(I)*C2D38(K)) 548 DCO2D1(I,K)=H1M2*(A1(I)*CDT31(K)+A2(I)*CDT38(K)) 549 CO2R2(I,K)=A1(I)*CO271(K)+A2(I)*CO278(K) 550 D2CD22(I,K)=H1M3*(A1(I)*C2D71(K)+A2(I)*C2D78(K)) 551 DCO2D2(I,K)=H1M2*(A1(I)*CDT71(K)+A2(I)*CDT78(K)) 552 184 CONTINUE 553 DO 190 K=1,L 554 DO 190 I=MYIS,MYIE 555 CO2MR(I,K)=A1(I)*CO2M51(K)+A2(I)*CO2M58(K) 556 CO2MD(I,K)=H1M2*(A1(I)*CDTM51(K)+A2(I)*CDTM58(K)) 557 CO2M2D(I,K)=H1M3*(A1(I)*C2DM51(K)+A2(I)*C2DM58(K)) 558 190 CONTINUE 559 C***COMPUTE CO2 TEMPERATURE INTERPOLATIONS FOR ALL BANDS,USING DIFT 560 C 561 C THE CASE WHERE K=1 IS HANDLED FIRST. WE ARE NOW REPLACING 562 C 3-DIMENSIONAL ARRAYS BY 2-D ARRAYS, TO SAVE SPACE. THUS THIS 563 C CALCULATION IS FOR (I,KP,1) 564 DO 211 KP=2,LP1 565 DO 211 I=MYIS,MYIE 566 DIFT(I,KP)=TDAV(I,KP)/TSTDAV(I,KP) 567 211 CONTINUE 568 DO 212 I=MYIS,MYIE 569 CO21(I,1,1)=1.0 570 CO2SP1(I,1)=1.0 571 CO2SP2(I,1)=1.0 572 212 CONTINUE 573 DO 215 KP=2,LP1 574 DO 215 I=MYIS,MYIE 575 C---CALCULATIONS FOR KP>1 FOR K=1 576 CO2R(I,KP)=A1(I)*CO251(KP,1)+A2(I)*CO258(KP,1) 577 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(KP,1)+A2(I)*CDT58(KP,1)) 578 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(KP,1)+A2(I)*C2D58(KP,1)) 579 CO21(I,KP,1)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 580 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 581 C---CALCULATIONS FOR (EFFECTIVELY) KP=1,K>KP. THESE USE THE 582 C SAME VALUE OF DIFT DUE TO SYMMETRY 583 CO2R(I,KP)=A1(I)*CO251(1,KP)+A2(I)*CO258(1,KP) 584 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(1,KP)+A2(I)*CDT58(1,KP)) 585 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(1,KP)+A2(I)*C2D58(1,KP)) 586 CO21(I,1,KP)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 587 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 588 215 CONTINUE 589 C THE TRANSMISSION FUNCTIONS USED IN SPA88 MAY BE COMPUTED NOW. 590 C---(IN THE 250 LOOP,DIFT REALLY SHOULD BE (I,1,K), BUT DIFT IS 591 C INVARIANT WITH RESPECT TO K,KP,AND SO (I,1,K)=(I,K,1)) 592 DO 250 K=2,LP1 593 DO 250 I=MYIS,MYIE Page 10 Source Listing LWR88 2014-12-17 20:46 LWR88.F 594 CO2SP1(I,K)=CO2R1(I,K)+DIFT(I,K)*(DCO2D1(I,K)+HAF*DIFT(I,K)* 595 1 D2CD21(I,K)) 596 CO2SP2(I,K)=CO2R2(I,K)+DIFT(I,K)*(DCO2D2(I,K)+HAF*DIFT(I,K)* 597 1 D2CD22(I,K)) 598 250 CONTINUE 599 C 600 C NEXT THE CASE WHEN K=2...L 601 DO 220 K=2,L 602 DO 222 KP=K+1,LP1 603 DO 222 I=MYIS,MYIE 604 DIFT(I,KP)=(TDAV(I,KP)-TDAV(I,K))/ 605 1 (TSTDAV(I,KP)-TSTDAV(I,K)) 606 CO2R(I,KP)=A1(I)*CO251(KP,K)+A2(I)*CO258(KP,K) 607 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(KP,K)+A2(I)*CDT58(KP,K)) 608 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(KP,K)+A2(I)*C2D58(KP,K)) 609 CO21(I,KP,K)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 610 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 611 CO2R(I,KP)=A1(I)*CO251(K,KP)+A2(I)*CO258(K,KP) 612 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(K,KP)+A2(I)*CDT58(K,KP)) 613 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(K,KP)+A2(I)*C2D58(K,KP)) 614 CO21(I,K,KP)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 615 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 616 222 CONTINUE 617 220 CONTINUE 618 C FINALLY THE CASE WHEN K=KP,K=2..LP1 619 DO 206 K=2,LP1 620 DO 206 I=MYIS,MYIE 621 DIFT(I,K)=HAF*(VSUM3(I,K)+VSUM3(I,K-1)) 622 CO2R(I,K)=A1(I)*CO251(K,K)+A2(I)*CO258(K,K) 623 DCO2DT(I,K)=H1M2*(A1(I)*CDT51(K,K)+A2(I)*CDT58(K,K)) 624 D2CDT2(I,K)=H1M3*(A1(I)*C2D51(K,K)+A2(I)*C2D58(K,K)) 625 CO21(I,K,K)=CO2R(I,K)+DIFT(I,K)*(DCO2DT(I,K)+ 626 1 HAF*DIFT(I,K)*D2CDT2(I,K)) 627 206 CONTINUE 628 C--- WE AREN'T DOING NBL TFS ON THE 100 CM-1 BANDS . 629 DO 260 K=1,L 630 DO 260 I=MYIS,MYIE 631 CO2NBL(I,K)=CO2MR(I,K)+VSUM3(I,K)*(CO2MD(I,K)+HAF* 632 1 VSUM3(I,K)*CO2M2D(I,K)) 633 260 CONTINUE 634 C***COMPUTE TEMP. COEFFICIENT BASED ON T(K) (SEE REF.2) 635 DO 264 K=1,LP1 636 DO 264 I=MYIS,MYIE 637 IF (T(I,K).LE.H25E2) THEN 638 TLSQU(I,K)=B0+(T(I,K)-H25E2)* 639 1 (B1+(T(I,K)-H25E2)* 640 2 (B2+B3*(T(I,K)-H25E2))) 641 ELSE 642 TLSQU(I,K)=B0 643 ENDIF 644 264 CONTINUE 645 C***APPLY TO ALL CO2 TFS 646 DO 280 K=1,LP1 647 DO 282 KP=1,LP1 648 DO 282 I=MYIS,MYIE 649 CO21(I,KP,K)=CO21(I,KP,K)*(ONE-TLSQU(I,KP))+TLSQU(I,KP) 650 282 CONTINUE Page 11 Source Listing LWR88 2014-12-17 20:46 LWR88.F 651 280 CONTINUE 652 DO 284 K=1,LP1 653 DO 286 I=MYIS,MYIE 654 CO2SP1(I,K)=CO2SP1(I,K)*(ONE-TLSQU(I,1))+TLSQU(I,1) 655 CO2SP2(I,K)=CO2SP2(I,K)*(ONE-TLSQU(I,1))+TLSQU(I,1) 656 286 CONTINUE 657 284 CONTINUE 658 DO 288 K=1,L 659 DO 290 I=MYIS,MYIE 660 CO2NBL(I,K)=CO2NBL(I,K)*(ONE-TLSQU(I,K))+TLSQU(I,K) 661 290 CONTINUE 662 288 CONTINUE 663 CALL FST88(HEATRA,GRNFLX,TOPFLX, 664 1 QH2O,PRESS,P,DELP,DELP2,TEMP,T, 665 2 CLDFAC,NCLDS,KTOP,KBTM,CAMT, 666 3 CO21,CO2NBL,CO2SP1,CO2SP2, 667 4 VAR1,VAR2,VAR3,VAR4,CNTVAL, 668 5 TOTO3,TPHIO3,TOTPHI,TOTVO2, 669 6 EMX1,EMX2,EMPL) 670 RETURN 671 END ENTRY POINTS Name lwr88_ Page 12 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 103 Label 346 342,343 105 Label 352 347 107 Label 357 353,354 125 Label 365 360,361 131 Label 392 373,374 161 Label 437 434 162 Label 441 438,439 163 Label 449 442 165 Label 448 443 171 Label 455 452 184 Label 468 460,461 190 Label 474 469,470 201 Label 399 394 203 Label 406 400,401 206 Label 543 535,536 211 Label 483 480,481 212 Label 488 484 215 Label 504 489,490 220 Label 533 517 222 Label 532 518,519 250 Label 514 508,509 260 Label 549 545,546 264 Label 560 551,552 280 Label 567 562 282 Label 566 563,564 284 Label 573 568 286 Label 572 569 288 Label 578 574 290 Label 577 575 801 Label 415 412 811 Label 421 418,419 812 Label 426 422,423 821 Label 430 427 A1 Local 330 R(4) 4 1 45 453,462,463,464,465,466,467,471,47 2,473,492,493,494,499,500,501,522, 523,524,527,528,529,538,539,540 A2 Local 330 R(4) 4 1 45 454,462,463,464,465,466,467,471,47 2,473,492,493,494,499,500,501,522, 523,524,527,528,529,538,539,540 AB15 Scalar 228 R(4) 4 1 2 COM AB15CM Scalar 297 R(4) 4 1 2 COM AB15WD Scalar 257 R(4) 4 scalar COM ACOMB Scalar 294 R(4) 4 1 15 COM AMOLWT Scalar 28 R(4) 4 scalar COM AO3CM Scalar 297 R(4) 4 1 3 COM AO3RND Scalar 228 R(4) 4 1 3 COM AP Scalar 226 R(4) 4 1 163 COM APCM Scalar 295 R(4) 4 1 15 COM APWD Scalar 255 R(4) 4 scalar COM ARNDM Scalar 226 R(4) 4 1 163 COM ATP Scalar 227 R(4) 4 1 163 COM Page 13 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References ATPCM Scalar 295 R(4) 4 1 15 COM ATPWD Scalar 255 R(4) 4 scalar COM AWIDE Scalar 254 R(4) 4 scalar COM B0 Scalar 152 R(4) 4 scalar COM 554,558 B1 Scalar 152 R(4) 4 scalar COM 555 B2 Scalar 152 R(4) 4 scalar COM 556 B3 Scalar 152 R(4) 4 scalar COM 556 BANDHI Scalar 228 R(4) 4 1 163 COM BANDLO Scalar 227 R(4) 4 1 163 COM BANDTA Common 226 5900 BCOMB Scalar 294 R(4) 4 1 15 COM BDCOMB Common 294 736 BDHICM Scalar 296 R(4) 4 1 15 COM BDHIWD Scalar 256 R(4) 4 scalar COM BDLOCM Scalar 296 R(4) 4 1 15 COM BDLOWD Scalar 256 R(4) 4 scalar COM BDWIDE Common 254 56 BETACM Scalar 295 R(4) 4 1 15 COM BETAD Scalar 226 R(4) 4 1 163 COM BETAWD Scalar 254 R(4) 4 scalar COM BETINC Scalar 296 R(4) 4 scalar COM BETINW Scalar 256 R(4) 4 scalar COM 390 BO3CM Scalar 297 R(4) 4 1 3 COM BO3RND Scalar 228 R(4) 4 1 3 COM BP Scalar 227 R(4) 4 1 163 COM BPCM Scalar 295 R(4) 4 1 15 COM BPWD Scalar 255 R(4) 4 scalar COM BRNDM Scalar 226 R(4) 4 1 163 COM BTP Scalar 227 R(4) 4 1 163 COM BTPCM Scalar 296 R(4) 4 1 15 COM BTPWD Scalar 255 R(4) 4 scalar COM BWIDE Scalar 254 R(4) 4 scalar COM C2D31 Scalar 167 R(4) 4 1 46 COM 463 C2D38 Scalar 167 R(4) 4 1 46 COM 463 C2D51 Scalar 150 R(4) 4 2 2116 COM 494,501,524,529,540 C2D58 Scalar 150 R(4) 4 2 2116 COM 494,501,524,529,540 C2D71 Scalar 182 R(4) 4 1 46 COM 466 C2D78 Scalar 182 R(4) 4 1 46 COM 466 C2DM51 Scalar 151 R(4) 4 1 45 COM 473 C2DM58 Scalar 151 R(4) 4 1 45 COM 473 CAMT Dummy 3 R(4) 4 2 2070 ARG,INOUT 581 CDT31 Scalar 166 R(4) 4 1 46 COM 464 CDT38 Scalar 167 R(4) 4 1 46 COM 464 CDT51 Scalar 149 R(4) 4 2 2116 COM 493,500,523,528,539 CDT58 Scalar 150 R(4) 4 2 2116 COM 493,500,523,528,539 CDT71 Scalar 181 R(4) 4 1 46 COM 467 CDT78 Scalar 182 R(4) 4 1 46 COM 467 CDTM51 Scalar 151 R(4) 4 1 45 COM 472 CDTM58 Scalar 151 R(4) 4 1 45 COM 472 CLDFAC Dummy 2 R(4) 4 3 95220 ARG,INOUT 581 CNTVAL Local 315 R(4) 4 2 2070 390,405,583 CO21 Local 311 R(4) 4 3 95220 485,495,502,525,530,541,565,582 CO211 Scalar 191 R(4) 4 1 46 COM CO218 Scalar 191 R(4) 4 1 46 COM CO231 Scalar 166 R(4) 4 1 46 COM 462 Page 14 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References CO238 Scalar 166 R(4) 4 1 46 COM 462 CO251 Scalar 149 R(4) 4 2 2116 COM 492,499,522,527,538 CO258 Scalar 149 R(4) 4 2 2116 COM 492,499,522,527,538 CO271 Scalar 181 R(4) 4 1 46 COM 465 CO278 Scalar 181 R(4) 4 1 46 COM 465 CO2BD2 Common 166 1104 CO2BD3 Common 149 52248 CO2BD4 Common 181 1104 CO2BD5 Common 191 368 CO2M2D Local 326 R(4) 4 2 2025 473,548 CO2M51 Scalar 150 R(4) 4 1 45 COM 471 CO2M58 Scalar 151 R(4) 4 1 45 COM 471 CO2MD Local 325 R(4) 4 2 2025 472,547 CO2MR Local 325 R(4) 4 2 2025 471,547 CO2NBL Local 311 R(4) 4 2 2025 547,576,582 CO2R Local 321 R(4) 4 2 2070 492,495,499,502,522,525,527,530,53 8,541 CO2R1 Local 322 R(4) 4 2 2070 462,510 CO2R2 Local 324 R(4) 4 2 2070 465,512 CO2SP1 Local 312 R(4) 4 2 2070 486,510,570,582 CO2SP2 Local 312 R(4) 4 2 2070 487,512,571,582 CSUBP Scalar 28 R(4) 4 scalar COM D2CD21 Local 323 R(4) 4 2 2070 463,511 D2CD22 Local 323 R(4) 4 2 2070 466,513 D2CDT2 Local 331 R(4) 4 2 2070 494,496,501,503,524,526,529,531,54 0,542 DCO2D1 Local 322 R(4) 4 2 2070 464,510 DCO2D2 Local 324 R(4) 4 2 2070 467,512 DCO2DT Local 331 R(4) 4 2 2070 493,495,500,502,523,525,528,530,53 9,541 DELP Local 310 R(4) 4 2 2025 356,580 DELP2 Local 307 R(4) 4 2 2025 355,356,379,380,444,580 DIFFCTR Scalar 28 R(4) 4 scalar COM 375,380 DIFT Local 321 R(4) 4 2 2070 482,495,496,502,503,510,512,520,52 5,526,530,531,537,541,542 EIGHT Scalar 33 R(4) 4 scalar COM EMPL Local 319 R(4) 4 2 4095 420,424,428,429,585 EMX1 Local 318 R(4) 4 1 45 413,585 EMX2 Local 319 R(4) 4 1 45 414,585 EXP Func 364 scalar 364 FIFTY Scalar 33 R(4) 4 scalar COM FIVE Scalar 33 R(4) 4 scalar COM FOUR Scalar 34 R(4) 4 scalar COM FST88 Subr 579 579 G Scalar 28 R(4) 4 scalar COM G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GINV Scalar 32 R(4) 4 scalar COM 379,380 GLB_TABLE Common 37 128 GP0INV Scalar 32 R(4) 4 scalar COM 413,414,420,425 GRAVDR Scalar 28 R(4) 4 scalar COM GRNFLX Dummy 1 R(4) 4 1 45 ARG,INOUT 579 GTEMP Scalar 152 R(4) 4 1 46 COM 444 H101M16 Scalar 66 R(4) 4 scalar COM H102M5 Scalar 65 R(4) 4 scalar COM Page 15 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References H1036E2 Scalar 64 R(4) 4 scalar COM H114M11 Scalar 48 R(4) 4 scalar COM H1174M7 Scalar 43 R(4) 4 scalar COM H11M10 Scalar 45 R(4) 4 scalar COM H11M11 Scalar 48 R(4) 4 scalar COM H1224E3 Scalar 37 R(4) 4 scalar COM H1226E1 Scalar 59 R(4) 4 scalar COM H128M5 Scalar 42 R(4) 4 scalar COM H129M2 Scalar 63 R(4) 4 scalar COM H12M12 Scalar 52 R(4) 4 scalar COM H12M13 Scalar 53 R(4) 4 scalar COM H12M31 Scalar 55 R(4) 4 scalar COM H135M13 Scalar 53 R(4) 4 scalar COM H1386E2 Scalar 63 R(4) 4 scalar COM H1439M5 Scalar 42 R(4) 4 scalar COM H14M10 Scalar 44 R(4) 4 scalar COM H14M11 Scalar 48 R(4) 4 scalar COM H14M12 Scalar 52 R(4) 4 scalar COM H14M14 Scalar 53 R(4) 4 scalar COM H14M30 Scalar 55 R(4) 4 scalar COM H15E2 Scalar 38 R(4) 4 scalar COM H15M11 Scalar 48 R(4) 4 scalar COM H15M14 Scalar 53 R(4) 4 scalar COM H15M5 Scalar 67 R(4) 4 scalar COM H161E1 Scalar 66 R(4) 4 scalar COM H165E5 Scalar 36 R(4) 4 scalar COM H16E1 Scalar 66 R(4) 4 scalar COM H16M12 Scalar 52 R(4) 4 scalar COM H181E1 Scalar 38 R(4) 4 scalar COM H18E1 Scalar 38 R(4) 4 scalar COM H18E3 Scalar 67 R(4) 4 scalar COM 362 H18M11 Scalar 48 R(4) 4 scalar COM H1E11 Scalar 35 R(4) 4 scalar COM H1E13 Scalar 35 R(4) 4 scalar COM H1E15 Scalar 35 R(4) 4 scalar COM H1E4 Scalar 36 R(4) 4 scalar COM H1E6 Scalar 58 R(4) 4 scalar COM H1E8 Scalar 35 R(4) 4 scalar COM H1M10 Scalar 45 R(4) 4 scalar COM H1M11 Scalar 49 R(4) 4 scalar COM H1M13 Scalar 53 R(4) 4 scalar COM H1M16 Scalar 69 R(4) 4 scalar COM H1M17 Scalar 54 R(4) 4 scalar COM H1M18 Scalar 54 R(4) 4 scalar COM H1M19 Scalar 54 R(4) 4 scalar COM H1M2 Scalar 58 R(4) 4 scalar COM 464,467,472,493,500,523,528,539 H1M20 Scalar 54 R(4) 4 scalar COM H1M21 Scalar 54 R(4) 4 scalar COM H1M22 Scalar 54 R(4) 4 scalar COM H1M23 Scalar 54 R(4) 4 scalar COM H1M24 Scalar 55 R(4) 4 scalar COM H1M3 Scalar 41 R(4) 4 scalar COM 463,466,473,494,501,524,529,540 H1M4 Scalar 41 R(4) 4 scalar COM H1M5 Scalar 42 R(4) 4 scalar COM H1M6 Scalar 43 R(4) 4 scalar COM Page 16 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References H1M60 Scalar 56 R(4) 4 scalar COM H1M8 Scalar 44 R(4) 4 scalar COM H1P082 Scalar 63 R(4) 4 scalar COM H1P25892 Scalar 39 R(4) 4 scalar COM H1P4 Scalar 39 R(4) 4 scalar COM H1P41819 Scalar 61 R(4) 4 scalar COM H1P4387 Scalar 39 R(4) 4 scalar COM H1P8 Scalar 39 R(4) 4 scalar COM H2075E3 Scalar 37 R(4) 4 scalar COM H20788E3 Scalar 36 R(4) 4 scalar COM H2118M2 Scalar 64 R(4) 4 scalar COM H21M12 Scalar 51 R(4) 4 scalar COM H21M31 Scalar 55 R(4) 4 scalar COM H235M3 Scalar 62 R(4) 4 scalar COM H23E2 Scalar 38 R(4) 4 scalar COM H23M10 Scalar 44 R(4) 4 scalar COM H23M11 Scalar 47 R(4) 4 scalar COM H24E3 Scalar 36 R(4) 4 scalar COM H24M11 Scalar 47 R(4) 4 scalar COM H24M12 Scalar 51 R(4) 4 scalar COM H25452M6 Scalar 43 R(4) 4 scalar COM H257M8 Scalar 44 R(4) 4 scalar COM H25E2 Scalar 67 R(4) 4 scalar COM 553,554,555,556 H25M31 Scalar 55 R(4) 4 scalar COM H26E2 Scalar 61 R(4) 4 scalar COM H26M30 Scalar 55 R(4) 4 scalar COM H28571M2 Scalar 69 R(4) 4 scalar COM H285M4 Scalar 41 R(4) 4 scalar COM H28E1 Scalar 72 R(4) 4 scalar COM H28M11 Scalar 47 R(4) 4 scalar COM H28M12 Scalar 51 R(4) 4 scalar COM H29316E2 Scalar 59 R(4) 4 scalar COM H2945E2 Scalar 37 R(4) 4 scalar COM H29M12 Scalar 51 R(4) 4 scalar COM H2E2 Scalar 61 R(4) 4 scalar COM H2E6 Scalar 58 R(4) 4 scalar COM H2M11 Scalar 48 R(4) 4 scalar COM H2P5 Scalar 39 R(4) 4 scalar COM H2P8 Scalar 38 R(4) 4 scalar COM H2P9 Scalar 38 R(4) 4 scalar COM H3082E2 Scalar 37 R(4) 4 scalar COM H3116E1 Scalar 59 R(4) 4 scalar COM H323M4 Scalar 64 R(4) 4 scalar COM H327M8 Scalar 44 R(4) 4 scalar COM H32M11 Scalar 47 R(4) 4 scalar COM H35E1 Scalar 38 R(4) 4 scalar COM H35M11 Scalar 47 R(4) 4 scalar COM H36M13 Scalar 52 R(4) 4 scalar COM H37412M5 Scalar 42 R(4) 4 scalar COM H37M11 Scalar 46 R(4) 4 scalar COM H37M12 Scalar 51 R(4) 4 scalar COM H38M12 Scalar 50 R(4) 4 scalar COM H391M7 Scalar 43 R(4) 4 scalar COM H394M5 Scalar 42 R(4) 4 scalar COM H3E2 Scalar 37 R(4) 4 scalar COM Page 17 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References H3M11 Scalar 47 R(4) 4 scalar COM H3M12 Scalar 51 R(4) 4 scalar COM H3M14 Scalar 53 R(4) 4 scalar COM H3M3 Scalar 66 R(4) 4 scalar COM 382 H3M4 Scalar 70 R(4) 4 scalar COM 381 H3P5 Scalar 67 R(4) 4 scalar COM H3P6 Scalar 38 R(4) 4 scalar COM H41666M2 Scalar 59 R(4) 4 scalar COM H42M11 Scalar 46 R(4) 4 scalar COM H42M2 Scalar 64 R(4) 4 scalar COM H44194M2 Scalar 61 R(4) 4 scalar COM H44871M2 Scalar 41 R(4) 4 scalar COM H44M11 Scalar 46 R(4) 4 scalar COM H44M12 Scalar 50 R(4) 4 scalar COM H451M6 Scalar 65 R(4) 4 scalar COM H45M12 Scalar 50 R(4) 4 scalar COM H45M32 Scalar 56 R(4) 4 scalar COM H46M13 Scalar 52 R(4) 4 scalar COM H488E4 Scalar 36 R(4) 4 scalar COM H48M11 Scalar 46 R(4) 4 scalar COM H4999M6 Scalar 43 R(4) 4 scalar COM H4E5 Scalar 35 R(4) 4 scalar COM H4M12 Scalar 50 R(4) 4 scalar COM H4M33 Scalar 56 R(4) 4 scalar COM H53M11 Scalar 46 R(4) 4 scalar COM H559M3 Scalar 41 R(4) 4 scalar COM H55M32 Scalar 56 R(4) 4 scalar COM H5725E4 Scalar 36 R(4) 4 scalar COM H5E2 Scalar 37 R(4) 4 scalar COM H625M2 Scalar 60 R(4) 4 scalar COM H62M12 Scalar 50 R(4) 4 scalar COM H62M34 Scalar 56 R(4) 4 scalar COM H658M2 Scalar 64 R(4) 4 scalar COM H65M12 Scalar 49 R(4) 4 scalar COM H67390E2 Scalar 65 R(4) 4 scalar COM H6938M5 Scalar 42 R(4) 4 scalar COM H69766E5 Scalar 62 R(4) 4 scalar COM H6M12 Scalar 50 R(4) 4 scalar COM H6P08108 Scalar 68 R(4) 4 scalar COM 362 H71E26 Scalar 35 R(4) 4 scalar COM H72M11 Scalar 46 R(4) 4 scalar COM H74M12 Scalar 49 R(4) 4 scalar COM H75826M4 Scalar 63 R(4) 4 scalar COM H77M11 Scalar 45 R(4) 4 scalar COM H77M12 Scalar 49 R(4) 4 scalar COM H7M6 Scalar 43 R(4) 4 scalar COM H8121E1 Scalar 61 R(4) 4 scalar COM H82M11 Scalar 45 R(4) 4 scalar COM H83E26 Scalar 35 R(4) 4 scalar COM H83M11 Scalar 45 R(4) 4 scalar COM H8725M8 Scalar 44 R(4) 4 scalar COM H8M11 Scalar 45 R(4) 4 scalar COM H8M13 Scalar 52 R(4) 4 scalar COM H93M12 Scalar 49 R(4) 4 scalar COM H96M12 Scalar 49 R(4) 4 scalar COM Page 18 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References H987M4 Scalar 41 R(4) 4 scalar COM H9M32 Scalar 56 R(4) 4 scalar COM H9P94 Scalar 60 R(4) 4 scalar COM HAF Scalar 34 R(4) 4 scalar COM 344,345,378,496,503,510,512,526,53 1,537,542,547 HCON Common 33 872 HEATRA Dummy 1 R(4) 4 2 2025 ARG,INOUT 579 HM13EZ Scalar 57 R(4) 4 scalar COM HM1597E1 Scalar 67 R(4) 4 scalar COM HM161E1 Scalar 66 R(4) 4 scalar COM HM1797E1 Scalar 60 R(4) 4 scalar COM HM181E1 Scalar 57 R(4) 4 scalar COM HM19EZ Scalar 57 R(4) 4 scalar COM HM1E1 Scalar 57 R(4) 4 scalar COM HM1E2 Scalar 57 R(4) 4 scalar COM HM1EZ Scalar 61 R(4) 4 scalar COM HM2M2 Scalar 59 R(4) 4 scalar COM HM6666M2 Scalar 58 R(4) 4 scalar COM HM8E1 Scalar 71 R(4) 4 scalar COM HMP5 Scalar 59 R(4) 4 scalar COM HMP575 Scalar 57 R(4) 4 scalar COM HMP66667 Scalar 58 R(4) 4 scalar COM HMP805 Scalar 68 R(4) 4 scalar COM HNINETY Scalar 33 R(4) 4 scalar COM HP1 Scalar 40 R(4) 4 scalar COM HP118666 Scalar 67 R(4) 4 scalar COM HP144 Scalar 62 R(4) 4 scalar COM HP166666 Scalar 58 R(4) 4 scalar COM HP219 Scalar 62 R(4) 4 scalar COM HP228 Scalar 60 R(4) 4 scalar COM HP26 Scalar 62 R(4) 4 scalar COM HP369 Scalar 40 R(4) 4 scalar COM HP3795 Scalar 65 R(4) 4 scalar COM HP5048 Scalar 65 R(4) 4 scalar COM HP518 Scalar 39 R(4) 4 scalar COM HP526315 Scalar 68 R(4) 4 scalar COM HP6 Scalar 60 R(4) 4 scalar COM HP602409 Scalar 68 R(4) 4 scalar COM HP60241 Scalar 60 R(4) 4 scalar COM HP8 Scalar 39 R(4) 4 scalar COM HP805 Scalar 63 R(4) 4 scalar COM HP816 Scalar 62 R(4) 4 scalar COM HTMG Scalar 52 R(4) 4 3 4183695 COM HUNDRED Scalar 33 R(4) 4 scalar COM I Local 343 I(4) 4 scalar 343,344,345,347,348,349,350,351,35 4,355,356,361,362,364,374,375,378, 379,380,381,382,390,391,394,395,39 6,397,398,401,402,403,404,405,412, 413,414,419,420,423,424,427,428,42 9,434,435,436,439,440,443,444,445, 446,447,452,453,454,461,462,463,46 4,465,466,467,470,471,472,473,481, 482,484,485,486,487,490,492,493,49 4,495,496,499,500,501,502,503,509, 510,511,512,513,519,520,521,522,52 Page 19 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References 3,524,525,526,527,528,529,530,531, 536,537,538,539,540,541,542,546,54 7,548,552,553,554,555,556,558,564, 565,569,570,571,575,576 IBAND Scalar 294 I(4) 4 1 40 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 300,301,302,303,304,305,306,307,30 9,310,311,312,313,314,315,316,317, 318,319,321,322,323,324,325,326,32 7,328,329,330,331,333,334 IDIM2 Param 19 I(4) 4 scalar 300,301,302,303,304,305,306,307,30 9,310,311,312,313,314,315,316,317, 318,319,321,322,323,324,325,326,32 7,328,329,330,331,333,334 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 INLTE Param 101 I(4) 4 scalar 101 INLTEP Param 101 I(4) 4 scalar 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 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 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 Page 20 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 K Local 342 I(4) 4 scalar 342,344,345,353,355,356,360,362,36 4,373,375,378,379,380,381,382,390, 391,400,402,403,404,405,418,420,42 2,424,438,440,442,444,445,446,447, 460,462,463,464,465,466,467,469,47 1,472,473,508,510,511,512,513,517, 518,520,521,522,523,524,525,527,52 8,529,530,535,537,538,539,540,541, 542,545,547,548,551,553,554,555,55 6,558,562,565,568,570,571,574,576 KBTM Dummy 3 I(4) 4 2 2070 ARG,INOUT 581 KO2 Param 103 I(4) 4 scalar 104 KO21 Param 104 I(4) 4 scalar KO2M Param 104 I(4) 4 scalar KP Local 480 I(4) 4 scalar 480,482,489,492,493,494,495,496,49 9,500,501,502,503,518,520,521,522, 523,524,525,526,527,528,529,530,53 1,563,565 KTOP Dummy 3 I(4) 4 2 2070 ARG,INOUT 581 L Param 89 I(4) 4 scalar 92,93,94,97,99,150,151,301,305,307 ,309,310,311,313,314,325,326,328,3 34,342,353,373,413,414,418,428,442 ,469,517,545,574 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LL Param 94 I(4) 4 scalar 94,95,429 LL3P Param 99 I(4) 4 scalar LLM1 Param 95 I(4) 4 scalar LLM2 Param 95 I(4) 4 scalar LLM3 Param 95 I(4) 4 scalar LLP1 Param 94 I(4) 4 scalar 319,429 LLP2 Param 94 I(4) 4 scalar LLP3 Param 94 I(4) 4 scalar LM Param 3 I(4) 4 scalar 52,89 LM1 Param 93 I(4) 4 scalar 422 LM2 Param 93 I(4) 4 scalar LM3 Param 93 I(4) 4 scalar LP1 Param 92 I(4) 4 scalar 96,97,98,149,150,152,166,167,181,1 82,191,300,302,303,304,309,310,311 ,312,315,316,317,318,321,322,323,3 24,327,328,331,333,349,351,360,400 ,414,438,453,454,460,480,489,508,5 18,535,551,562,563,568 LP121 Param 98 I(4) 4 scalar LP1M Param 96 I(4) 4 scalar 96 LP1M1 Param 96 I(4) 4 scalar Page 21 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References LP1V Param 97 I(4) 4 scalar LP2 Param 92 I(4) 4 scalar 424 LP3 Param 92 I(4) 4 scalar LSM Param 3 I(4) 4 scalar LWR88 Subr 1 MAPPINGS Common 56 5024 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 343,347,354,361,374,394,401,412,41 9,423,427,434,439,443,452,461,470, 481,484,490,509,519,536,546,552,56 4,569,575 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 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 343,347,354,361,374,394,401,412,41 9,423,427,434,439,443,452,461,470, 481,484,490,509,519,536,546,552,56 4,569,575 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 Page 22 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 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 NB Param 100 I(4) 4 scalar 102 NB1 Param 102 I(4) 4 scalar NBLM Param 91 I(4) 4 scalar Page 23 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References NBLW Param 90 I(4) 4 scalar 226,227,228 NBLX Param 90 I(4) 4 scalar NBLY Param 90 I(4) 4 scalar 91,98,294,295,296 NCLDS Dummy 3 I(4) 4 1 45 ARG,INOUT 581 NNLTE Param 101 I(4) 4 scalar NPES Scalar 2 I(4) 4 scalar COM O3DIFCTR Scalar 28 R(4) 4 scalar COM ONE Scalar 34 R(4) 4 scalar COM 356,565,570,571,576 P Local 310 R(4) 4 2 2070 344,348,349,355,378,413,414,420,42 4,580 P0 Scalar 28 R(4) 4 scalar COM 454 P0INV Scalar 32 R(4) 4 scalar COM 378 P0X2 Scalar 29 R(4) 4 scalar COM P0XZP2 Scalar 29 R(4) 4 scalar COM 453,454 P0XZP8 Scalar 29 R(4) 4 scalar COM 453 PHYCON Common 28 80 PRESS Dummy 2 R(4) 4 2 2070 ARG,INOUT 344,349,413,414,420,424,453,454,58 0 QH2O Local 309 R(4) 4 2 2025 375,379,413,414,420,424,580 QO3 Dummy 2 R(4) 4 2 2025 ARG,INOUT 380 QUARTR Scalar 34 R(4) 4 scalar COM RADCON Scalar 29 R(4) 4 scalar COM RADCON1 Scalar 31 R(4) 4 scalar COM RATCO2MW Scalar 30 R(4) 4 scalar COM RATH2OMW Scalar 30 R(4) 4 scalar COM 391 RGAS Scalar 29 R(4) 4 scalar COM RGASSP Scalar 29 R(4) 4 scalar COM RH2O Dummy 2 R(4) 4 2 2025 ARG,INOUT 375,390,391 SECPDA Scalar 29 R(4) 4 scalar COM SIXTY Scalar 33 R(4) 4 scalar COM SKC1R Scalar 257 R(4) 4 scalar COM SKO2D Scalar 257 R(4) 4 scalar COM SKO3R Scalar 257 R(4) 4 scalar COM STEMP Scalar 152 R(4) 4 1 46 COM 440 T Local 309 R(4) 4 2 2070 345,350,351,553,554,555,556,580 TDAV Local 327 R(4) 4 2 2070 436,447,482,520 TEMP Dummy 2 R(4) 4 2 2070 ARG,INOUT 345,350,351,362,440,580 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 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 Page 24 Source Listing LWR88 2014-12-17 20:46 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References TEN Scalar 33 R(4) 4 scalar COM TEXPSL Local 333 R(4) 4 2 2070 362,364,390 THREE Scalar 34 R(4) 4 scalar COM TLSQU Local 333 R(4) 4 2 2070 554,558,565,570,571,576 TOPFLX Dummy 1 R(4) 4 1 45 ARG,INOUT 579 TOPO Common 50 17478548 TOTO3 Local 316 R(4) 4 2 2070 396,403,584 TOTPHI Local 317 R(4) 4 2 2070 395,402,584 TOTVO2 Local 318 R(4) 4 2 2070 398,405,584 TPHIO3 Local 316 R(4) 4 2 2070 397,404,584 TSTDAV Local 327 R(4) 4 2 2070 435,446,482,521 TTVG Scalar 51 R(4) 4 2 92971 COM TWO Scalar 34 R(4) 4 scalar COM VAR1 Local 313 R(4) 4 2 2025 379,381,583 VAR2 Local 313 R(4) 4 2 2025 381,390,402,428,583 VAR3 Local 314 R(4) 4 2 2025 380,382,403,583 VAR4 Local 314 R(4) 4 2 2025 382,404,583 VSUM1 Local 328 R(4) 4 1 45 445,447 VSUM2 Local 329 R(4) 4 1 45 444,445,446 VSUM3 Local 328 R(4) 4 2 2070 440,445,537,547,548 VSUM4 Local 334 R(4) 4 2 2025 VV Local 328 R(4) 4 2 2025 378,381,382 ZERO Scalar 34 R(4) 4 scalar COM 348,395,396,397,398,435,436 Page 25 Source Listing LWR88 2014-12-17 20:46 Subprograms/Common Blocks LWR88.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BANDTA Common 226 5900 BDCOMB Common 294 736 BDWIDE Common 254 56 CO2BD2 Common 166 1104 CO2BD3 Common 149 52248 CO2BD4 Common 181 1104 CO2BD5 Common 191 368 GLB_TABLE Common 37 128 HCON Common 33 872 LWR88 Subr 1 MAPPINGS Common 56 5024 MPPCOM Common 1 1464 PHYCON Common 28 80 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__ Page 26 Source Listing LWR88 2014-12-17 20:46 LWR88.F -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 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/, Page 27 Source Listing LWR88 2014-12-17 20:46 LWR88.F /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : LWR88.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100