Page 1 Source Listing LWR88 2025-03-12 18:22 /tmp/ifortyAX8ZD.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 2025-03-12 18:22 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" 99 INCLUDE "mpp.h" 100 # 1 "./sp.h" 1 101 # 4 102 103 # 77 "LWR88.F" 2 104 C----------------------------------------------------------------------- 105 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 106 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 107 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 108 C BANDTA FOR DEFINITION 109 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 110 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE 111 C BDCOMB FOR DEFINITION 112 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 113 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 114 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 115 C FROM THE ABOVE PARAMETERS. 174 PARAMETER (L=LM) 175 PARAMETER (NBLW=163,NBLX=47,NBLY=15) 176 PARAMETER (NBLM=NBLY-1) 177 PARAMETER (LP1=L+1,LP2=L+2,LP3=L+3) 178 PARAMETER (LM1=L-1,LM2=L-2,LM3=L-3) 179 PARAMETER (LL=2*L,LLP1=LL+1,LLP2=LL+2,LLP3=LL+3) 180 PARAMETER (LLM1=LL-1,LLM2=LL-2,LLM3=LL-3) 181 PARAMETER (LP1M=LP1*LP1,LP1M1=LP1M-1) 182 PARAMETER (LP1V=LP1*(1+2*L/2)) 183 PARAMETER (LP121=LP1*NBLY) 184 PARAMETER (LL3P=3*L+2) 185 PARAMETER (NB=12) 186 PARAMETER (INLTE=3,INLTEP=INLTE+1,NNLTE=56) 187 PARAMETER (NB1=NB-1) 188 PARAMETER (KO2=12) 189 PARAMETER (KO21=KO2+1,KO2M=KO2-1) 190 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 191 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 192 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 193 C BANDTA FOR DEFINITION 194 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 195 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE Page 3 Source Listing LWR88 2025-03-12 18:22 LWR88.F 196 C BDCOMB FOR DEFINITION 197 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 198 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 199 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 200 C FROM THE ABOVE PARAMETERS. 201 C THE FOLLOWING COMMON BLOCKS CONTAIN PRETABULATED CO2 TRANSMISSION 202 C FUNCTIONS, EVALUATED USING THE METHODS OF FELS AND 203 C SCHWARZKOPF (1981) AND SCHWARZKOPF AND FELS (1985), 204 C***COMMON CO2BD3 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE 205 C AND PRESSURE DERIVATIVES FOR THE 560-800 CM-1 BAND. ALSO INCLUDED 206 C ARE THE STANDARD TEMPERATURES AND THE WEIGHTING FUNCTION. THESE 207 C DATA ARE IN BLOCK DATA BD3: 208 C CO251 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 209 C WITH P(SFC)=1013.25 MB 210 C CO258 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 211 C WITH P(SFC)= ^810 MB 212 C CDT51 = FIRST TEMPERATURE DERIVATIVE OF CO251 213 C CDT58 = FIRST TEMPERATURE DERIVATIVE OF CO258 214 C C2D51 = SECOND TEMPERATURE DERIVATIVE OF CO251 215 C C2D58 = SECOND TEMPERATURE DERIVATIVE OF CO251 216 C CO2M51 = TRANSMISSION FCTNS FOR T0 FOR ADJACENT PRESSURE 217 C LEVELS, WITH NO PRESSURE QUADRATURE. USED FOR 218 C NEARBY LAYER COMPUTATIONS. P(SFC)=1013.25 MB 219 C CO2M58 = SAME AS CO2M51,WITH P(SFC)= ^810 MB 220 C CDTM51 = FIRST TEMPERATURE DERIVATIVE OF CO2M51 221 C CDTM58 = FIRST TEMPERATURE DERIVATIVE OF CO2M58 222 C C2DM51 = SECOND TEMPERATURE DERIVATIVE OF CO2M51 223 C C2DM58 = SECOND TEMPERATURE DERIVATIVE OF CO2M58 224 C STEMP = STANDARD TEMPERATURES FOR MODEL PRESSURE LEVEL 225 C STRUCTURE WITH P(SFC)=1013.25 MB 226 C GTEMP = WEIGHTING FUNCTION FOR MODEL PRESSURE LEVEL 227 C STRUCTURE WITH P(SFC)=1013.25 MB. 228 C B0 = TEMP. COEFFICIENT USED FOR CO2 TRANS. FCTN. 229 C CORRECTION FOR T(K). (SEE REF. 4 AND BD3) 230 C B1 = TEMP. COEFFICIENT, USED ALONG WITH B0 231 C B2 = TEMP. COEFFICIENT, USED ALONG WITH B0 232 C B3 = TEMP. COEFFICIENT, USED ALONG WITH B0 233 C 234 COMMON/CO2BD3/CO251(LP1,LP1),CO258(LP1,LP1),CDT51(LP1,LP1), 235 1 CDT58(LP1,LP1),C2D51(LP1,LP1),C2D58(LP1,LP1),CO2M51(L), 236 2 CO2M58(L),CDTM51(L),CDTM58(L),C2DM51(L),C2DM58(L), 237 3 STEMP(LP1),GTEMP(LP1),B0,B1,B2,B3 238 C 239 C***COMMON CO2BD2 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE 240 C AND PRESSURE DERIVATIVES FOR THE 560-670 CM-1 PART OF THE 15 UM 241 C CO2 BAND. THESE DATA ARE IN BLOCK DATA BD2. 242 C CO231 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 243 C WITH P(SFC)=1013.25 MB 244 C CO238 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 245 C WITH P(SFC)= ^810 MB 246 C CDT31 = FIRST TEMPERATURE DERIVATIVE OF CO231 247 C CDT38 = FIRST TEMPERATURE DERIVATIVE OF CO238 248 C C2D31 = SECOND TEMPERATURE DERIVATIVE OF CO231 249 C C2D38 = SECOND TEMPERATURE DERIVATIVE OF CO231 250 C 251 COMMON / CO2BD2 / CO231(LP1),CO238(LP1),CDT31(LP1), 252 1 CDT38(LP1),C2D31(LP1),C2D38(LP1) Page 4 Source Listing LWR88 2025-03-12 18:22 LWR88.F 253 C 254 C***COMMON CO2BD4 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE 255 C AND PRESSURE DERIVATIVES FOR THE 670-800 CM-1 PART OF THE 15 UM 256 C CO2 BAND. THESE DATA ARE IN BLOCK DATA BD4. 257 C CO271 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 258 C WITH P(SFC)=1013.25 MB 259 C CO278 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 260 C WITH P(SFC)= ^810 MB 261 C CDT71 = FIRST TEMPERATURE DERIVATIVE OF CO271 262 C CDT78 = FIRST TEMPERATURE DERIVATIVE OF CO278 263 C C2D71 = SECOND TEMPERATURE DERIVATIVE OF CO271 264 C C2D78 = SECOND TEMPERATURE DERIVATIVE OF CO271 265 C 266 COMMON / CO2BD4 / CO271(LP1),CO278(LP1),CDT71(LP1), 267 1 CDT78(LP1),C2D71(LP1),C2D78(LP1) 268 C 269 C***COMMON CO2BD5 CONTAINS CO2 TRANSMISSION FUNCTIONS FOR THE 2270- 270 C 2380 PART OF THE 4.3 UM CO2 BAND. THESE DATA ARE IN BLOCK DATA BD5. 271 C CO211 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE) 272 C WITH P(SFC)=1013.25 MB 273 C CO218 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE) 274 C WITH P(SFC)= ^810 MB 275 C 276 COMMON / CO2BD5 / CO211(LP1),CO218(LP1) 277 C 278 C COMMON BLOCK BANDTA CONTAINS RANDOM BAND PARAMETERS FOR THE LW 279 C CALCULATIONS USING 10 CM-1 WIDE BANDS.THE 15 UM CO2 COMPLEX 280 C IS 2 BANDS,560-670 AND 670-800 CM-1. OZONE COEFFICIENTS ARE 281 C IN 3 BANDS,670-800 (14.1 UM),990-1070 AND 1070-1200 (9.6 UM). 282 C THE (NBLW) BANDS NOW INCLUDE: 283 C 56 BANDS, 10 CM-1 WIDE 0 - 560 CM-1 284 C 2 BANDS, 15 UM COMPLEX 560 - 670 CM-1 285 C 670 - 800 CM-1 286 C 3 "CONTINUUM" BANDS 800 - 900 CM-1 287 C 900 - 990 CM-1 288 C 1070 - 1200 CM-1 289 C 1 BAND FOR 9.6 UM BAND 990 - 1070 CM-1 290 C 100 BANDS, 10 CM-1 WIDE 1200 - 2200 CM-1 291 C 1 BAND FOR 4.3 UM SRC 2270 - 2380 CM-1 292 C THUS NBLW PRESENTLY EQUALS 163 293 C ALL BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 294 C 295 C ARNDM = RANDOM "A" PARAMETER FOR (NBLW) BANDS 296 C BRNDM = RANDOM "B" PARAMETER FOR (NBLW) BANDS 297 C BETAD = CONTINUUM COEFFICIENTS FOR (NBLW) BANDS 298 C AP,BP = CAPPHI COEFFICIENTS FOR (NBLW) BANDS 299 C ATP,BTP = CAPPSI COEFFICIENTS FOR (NBLW) BANDS 300 C BANDLO = LOWEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 301 C BANDHI = HIGHEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 302 C AO3RND = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 303 C BANDS 304 C BO3RND = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 305 C BANDS 306 C AB15 = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS 307 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 308 C DATA FOR ARNDM,BRNDM,AP,BP,ATP,BTP,AO3RND,BO3RND ARE OBTAINED BY 309 C USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS ARE FROM Page 5 Source Listing LWR88 2025-03-12 18:22 LWR88.F 310 C ROBERTS (1976). 311 COMMON / BANDTA / ARNDM(NBLW),BRNDM(NBLW),BETAD(NBLW),AP(NBLW), 312 1 BP(NBLW),ATP(NBLW),BTP(NBLW),BANDLO(NBLW), 313 2 BANDHI(NBLW),AO3RND(3),BO3RND(3),AB15(2) 314 C 315 C COMMON BLOCK BDWIDE CONTAINS RANDOM BAND PARAMETERS FOR SPECIFIC 316 C WIDE BANDS. AT PRESENT,THE INFORMATION CONSISTS OF 1) RANDOM 317 C MODEL PARAMETERS FOR THE 15 UM BAND,560-800 CM-1; 2) THE 318 C CONTINUUM COEFFICIENT FOR THE 800-990,1070-1200 CM-1 BAND 319 C SPECIFICALLY: 320 C AWIDE = RANDOM "A" PARAMETER FOR BAND 321 C BWIDE = RANDOM "B" PARAMETER FOR BAND 322 C BETAWD = CONTINUUM COEFFICIENTS FOR BAND 323 C APWD,BPWD = CAPPHI COEFFICIENTS FOR BAND 324 C ATPWD,BTPWD = CAPPSI COEFFICIENTS FOR BAND 325 C BDLOWD = LOWEST FREQUENCY IN EACH FREQ BAND 326 C BDHIWD = HIGHEST FREQUENCY IN EACH FREQ BAND 327 C AB15WD = THE PRODUCT ARNDM*BRNDM FOR THE ONE BAND 328 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 329 C BETINW = CONT.COEFFICIENT FOR A SPECIFIED WIDE 330 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 331 C SKO2D = 1./BETINW, USED IN SPA88 FOR CONT. COEFFS 332 C SKC1R = BETAWD/BETINW, USED FOR CONT. COEFF. FOR 333 C 15 UM BAND IN FST88 334 C SKO3R = RATIO OF CONT. COEFF. FOR 9.9 UM BAND TO 335 C BETINW, USED FOR 9.6 UM CONT COEFF IN FST88 336 C DATA FOR AWIDE,BWIDE,APWD,BPWD,ATPWD,BTPWD,AO3WD,BO3WD ARE 337 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 338 C ARE FROM ROBERTS (1976). 339 COMMON / BDWIDE / AWIDE,BWIDE,BETAWD, 340 1 APWD,BPWD,ATPWD,BTPWD, 341 2 BDLOWD,BDHIWD,BETINW, 342 3 AB15WD,SKO2D,SKC1R,SKO3R 343 C 344 C COMMON BLOCK BDCOMB CONTAINS RANDOM BAND PARAMETERS FOR THE LW 345 C CALCULATIONS USING COMBINED WIDE FREQUENCY BANDS BETWEEN 160 AND 346 C 1200 CM-1,AS WELL AS THE 2270-2380 BAND FOR SOURCE CALC. 347 C BANDS 1-8: COMBINED WIDE FREQUENCY BANDS FOR 160-560 CM-1 348 C BANDS 9-14: FREQUENCY BANDS,AS IN BANDTA (NARROW BANDS) 349 C FOR 560-1200 CM-1 350 C BAND 15: FREQUENCY BAND 2270-2380 CM-1,USED FOR SOURCE 351 C CALCULATION ONLY 352 C THUS NBLY PRESENTLY EQUALS 15 353 C 354 C BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 355 C ACOMB = RANDOM "A" PARAMETER FOR (NBLY) BANDS 356 C BCOMB = RANDOM "B" PARAMETER FOR (NBLY) BANDS 357 C BETACM = CONTINUUM COEFFICIENTS FOR (NBLY) BANDS 358 C APCM,BPCM = CAPPHI COEFFICIENTS FOR (NBLY) BANDS 359 C ATPCM,BTPCM = CAPPSI COEFFICIENTS FOR (NBLY) BANDS 360 C BDLOCM = LOWEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 361 C BDHICM = HIGHEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 362 C AO3CM = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 363 C BANDS 364 C BO3CM = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 365 C BANDS 366 C AB15CM = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS Page 6 Source Listing LWR88 2025-03-12 18:22 LWR88.F 367 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 368 C BETINC = CONT.COEFFICIENT FOR A SPECIFIED WIDE 369 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 370 C IBAND = INDEX NO OF THE 40 WIDE BANDS USED IN 371 C COMBINED WIDE BAND CALCULATIONS. IN OTHER 372 C WORDS,INDEX TELLING WHICH OF THE 40 WIDE 373 C BANDS BETWEEN 160-560 CM-1 ARE INCLUDED IN 374 C EACH OF THE FIRST 8 COMBINED WIDE BANDS 375 C DATA FOR ACOMB,BCOMB,APCM,BPCM,ATPCM,BTPCM,AO3CM,BO3CM ARE 376 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 377 C ARE FROM ROBERTS (1976). IBAND INDEX VALUES ARE OBTAINED BY 378 C EXPERIMENTATION. 379 COMMON / BDCOMB / IBAND(40),ACOMB(NBLY),BCOMB(NBLY), 380 1 BETACM(NBLY),APCM(NBLY),BPCM(NBLY),ATPCM(NBLY), 381 2 BTPCM(NBLY),BDLOCM(NBLY),BDHICM(NBLY),BETINC, 382 3 AO3CM(3),BO3CM(3),AB15CM(2) 383 C 384 C 385 DIMENSION PRESS(IDIM1:IDIM2,LP1),TEMP(IDIM1:IDIM2,LP1), 386 1 RH2O(IDIM1:IDIM2,L),QO3(IDIM1:IDIM2,L) 387 DIMENSION CLDFAC(IDIM1:IDIM2,LP1,LP1),CAMT(IDIM1:IDIM2,LP1) 388 DIMENSION NCLDS(IDIM1:IDIM2),KTOP(IDIM1:IDIM2,LP1), 389 1 KBTM(IDIM1:IDIM2,LP1) 390 DIMENSION HEATRA(IDIM1:IDIM2,L),GRNFLX(IDIM1:IDIM2), 391 1 TOPFLX(IDIM1:IDIM2) 392 DIMENSION DELP2(IDIM1:IDIM2,L) 393 C 394 DIMENSION QH2O(IDIM1:IDIM2,L),T(IDIM1:IDIM2,LP1) 395 DIMENSION P(IDIM1:IDIM2,LP1),DELP(IDIM1:IDIM2,L) 396 DIMENSION CO21(IDIM1:IDIM2,LP1,LP1),CO2NBL(IDIM1:IDIM2,L) 397 DIMENSION CO2SP1(IDIM1:IDIM2,LP1),CO2SP2(IDIM1:IDIM2,LP1) 398 DIMENSION VAR1(IDIM1:IDIM2,L),VAR2(IDIM1:IDIM2,L), 399 1 VAR3(IDIM1:IDIM2,L),VAR4(IDIM1:IDIM2,L) 400 DIMENSION CNTVAL(IDIM1:IDIM2,LP1) 401 DIMENSION TOTO3(IDIM1:IDIM2,LP1),TPHIO3(IDIM1:IDIM2,LP1), 402 1 TOTPHI(IDIM1:IDIM2,LP1) 403 DIMENSION TOTVO2(IDIM1:IDIM2,LP1),EMX1(IDIM1:IDIM2), 404 1 EMX2(IDIM1:IDIM2),EMPL(IDIM1:IDIM2,LLP1) 405 C 406 DIMENSION CO2R(IDIM1:IDIM2,LP1),DIFT(IDIM1:IDIM2,LP1) 407 DIMENSION CO2R1(IDIM1:IDIM2,LP1),DCO2D1(IDIM1:IDIM2,LP1) 408 DIMENSION D2CD21(IDIM1:IDIM2,LP1),D2CD22(IDIM1:IDIM2,LP1) 409 DIMENSION CO2R2(IDIM1:IDIM2,LP1),DCO2D2(IDIM1:IDIM2,LP1) 410 DIMENSION CO2MR(IDIM1:IDIM2,L),CO2MD(IDIM1:IDIM2,L), 411 1 CO2M2D(IDIM1:IDIM2,L) 412 DIMENSION TDAV(IDIM1:IDIM2,LP1),TSTDAV(IDIM1:IDIM2,LP1), 413 1 VV(IDIM1:IDIM2,L),VSUM3(IDIM1:IDIM2,LP1),VSUM1(IDIM1:IDIM2), 414 2 VSUM2(IDIM1:IDIM2) 415 DIMENSION A1(IDIM1:IDIM2),A2(IDIM1:IDIM2) 416 DIMENSION DCO2DT(IDIM1:IDIM2,LP1),D2CDT2(IDIM1:IDIM2,LP1) 417 C 418 DIMENSION TEXPSL(IDIM1:IDIM2,LP1),TLSQU(IDIM1:IDIM2,LP1) 419 DIMENSION VSUM4(IDIM1:IDIM2,L) 420 EQUIVALENCE (VSUM3,TLSQU,TEXPSL) 421 EQUIVALENCE (VV,VSUM4) 422 c 423 C Page 7 Source Listing LWR88 2025-03-12 18:22 LWR88.F 424 C****COMPUTE FLUX PRESSURES (P) AND DIFFERENCES (DELP2,DELP) 425 C****COMPUTE FLUX LEVEL TEMPERATURES (T) AND CONTINUUM TEMPERATURE 426 C CORRECTIONS (TEXPSL) 427 DO 103 K=2,L 428 DO 103 I=MYIS,MYIE 429 P(I,K)=HAF*(PRESS(I,K-1)+PRESS(I,K)) 430 T(I,K)=HAF*(TEMP(I,K-1)+TEMP(I,K)) 431 103 CONTINUE 432 DO 105 I=MYIS,MYIE 433 P(I,1)=ZERO 434 P(I,LP1)=PRESS(I,LP1) 435 T(I,1)=TEMP(I,1) 436 T(I,LP1)=TEMP(I,LP1) 437 105 CONTINUE 438 DO 107 K=1,L 439 DO 107 I=MYIS,MYIE 440 DELP2(I,K)=P(I,K+1)-P(I,K) 441 DELP(I,K)=ONE/DELP2(I,K) 442 107 CONTINUE 443 C****COMPUTE ARGUMENT FOR CONT.TEMP.COEFF. 444 C (THIS IS 1800.(1./TEMP-1./296.)) 445 DO 125 K=1,LP1 446 DO 125 I=MYIS,MYIE 447 TEXPSL(I,K)=H18E3/TEMP(I,K)-H6P08108 448 C...THEN TAKE EXPONENTIAL 449 TEXPSL(I,K)=EXP(TEXPSL(I,K)) 450 125 CONTINUE 451 C***COMPUTE OPTICAL PATHS FOR H2O AND O3, USING THE DIFFUSIVITY 452 C APPROXIMATION FOR THE ANGULAR INTEGRATION (1.66). OBTAIN THE 453 C UNWEIGHTED VALUES(VAR1,VAR3) AND THE WEIGHTED VALUES(VAR2,VAR4). 454 C THE QUANTITIES H3M4(.0003) AND H3M3(.003) APPEARING IN THE VAR2 AND 455 C VAR4 EXPRESSIONS ARE THE APPROXIMATE VOIGT CORRECTIONS FOR H2O AND 456 C O3,RESPECTIVELY. 457 C 458 DO 131 K=1,L 459 DO 131 I=MYIS,MYIE 460 QH2O(I,K)=RH2O(I,K)*DIFFCTR 461 C---VV IS THE LAYER-MEAN PRESSURE (IN ATM),WHICH IS NOT THE SAME AS 462 C THE LEVEL PRESSURE (PRESS) 463 VV(I,K)=HAF*(P(I,K+1)+P(I,K))*P0INV 464 VAR1(I,K)=DELP2(I,K)*QH2O(I,K)*GINV 465 VAR3(I,K)=DELP2(I,K)*QO3(I,K)*DIFFCTR*GINV 466 VAR2(I,K)=VAR1(I,K)*(VV(I,K)+H3M4) 467 VAR4(I,K)=VAR3(I,K)*(VV(I,K)+H3M3) 468 C COMPUTE OPTICAL PATH FOR THE H2O CONTINUUM, USING ROBERTS COEFFS. 469 C (BETINW),AND TEMP. CORRECTION (TEXPSL). THE DIFFUSIVITY FACTOR 470 C (WHICH CANCELS OUT IN THIS EXPRESSION) IS ASSUMED TO BE 1.66. THE 471 C USE OF THE DIFFUSIVITY FACTOR HAS BEEN SHOWN TO BE A SIGNIFICANT 472 C SOURCE OF ERROR IN THE CONTINUUM CALCS.,BUT THE TIME PENALTY OF 473 C AN ANGULAR INTEGRATION IS SEVERE. 474 C 475 CNTVAL(I,K)=TEXPSL(I,K)*RH2O(I,K)*VAR2(I,K)*BETINW/ 476 1 (RH2O(I,K)+RATH2OMW) 477 131 CONTINUE 478 C COMPUTE SUMMED OPTICAL PATHS FOR H2O,O3 AND CONTINUUM 479 DO 201 I=MYIS,MYIE 480 TOTPHI(I,1)=ZERO Page 8 Source Listing LWR88 2025-03-12 18:22 LWR88.F 481 TOTO3(I,1)=ZERO 482 TPHIO3(I,1)=ZERO 483 TOTVO2(I,1)=ZERO 484 201 CONTINUE 485 DO 203 K=2,LP1 486 DO 203 I=MYIS,MYIE 487 TOTPHI(I,K)=TOTPHI(I,K-1)+VAR2(I,K-1) 488 TOTO3(I,K)=TOTO3(I,K-1)+VAR3(I,K-1) 489 TPHIO3(I,K)=TPHIO3(I,K-1)+VAR4(I,K-1) 490 TOTVO2(I,K)=TOTVO2(I,K-1)+CNTVAL(I,K-1) 491 203 CONTINUE 492 C---EMX1 IS THE ADDITIONAL PRESSURE-SCALED MASS FROM PRESS(L) TO 493 C P(L). IT IS USED IN NEARBY LAYER AND EMISS CALCULATIONS. 494 C---EMX2 IS THE ADDITIONAL PRESSURE-SCALED MASS FROM PRESS(L) TO 495 C P(LP1). IT IS USED IN CALCULATIONS BETWEEN FLUX LEVELS L AND LP1. 496 C 497 DO 801 I=MYIS,MYIE 498 EMX1(I)=QH2O(I,L)*PRESS(I,L)*(PRESS(I,L)-P(I,L))*GP0INV 499 EMX2(I)=QH2O(I,L)*PRESS(I,L)*(P(I,LP1)-PRESS(I,L))*GP0INV 500 801 CONTINUE 501 C---EMPL IS THE PRESSURE SCALED MASS FROM P(K) TO PRESS(K) (INDEX 2-LP1) 502 C OR TO PRESS(K+1) (INDEX LP2-LL) 503 DO 811 K=1,L 504 DO 811 I=MYIS,MYIE 505 EMPL(I,K+1)=QH2O(I,K)*P(I,K+1)*(P(I,K+1)-PRESS(I,K))*GP0INV 506 811 CONTINUE 507 DO 812 K=1,LM1 508 DO 812 I=MYIS,MYIE 509 EMPL(I,LP2+K-1)=QH2O(I,K+1)*P(I,K+1)*(PRESS(I,K+1)-P(I,K+1)) 510 1 *GP0INV 511 812 CONTINUE 512 DO 821 I=MYIS,MYIE 513 EMPL(I,1)=VAR2(I,L) 514 EMPL(I,LLP1)=EMPL(I,LL) 515 821 CONTINUE 516 C***COMPUTE WEIGHTED TEMPERATURE (TDAV) AND PRESSURE (TSTDAV) INTEGRALS 517 C FOR USE IN OBTAINING TEMP. DIFFERENCE BET. SOUNDING AND STD. 518 C TEMP. SOUNDING (DIFT) 519 DO 161 I=MYIS,MYIE 520 TSTDAV(I,1)=ZERO 521 TDAV(I,1)=ZERO 522 161 CONTINUE 523 DO 162 K=1,LP1 524 DO 162 I=MYIS,MYIE 525 VSUM3(I,K)=TEMP(I,K)-STEMP(K) 526 162 CONTINUE 527 DO 163 K=1,L 528 DO 165 I=MYIS,MYIE 529 VSUM2(I)=GTEMP(K)*DELP2(I,K) 530 VSUM1(I)=VSUM2(I)*VSUM3(I,K) 531 TSTDAV(I,K+1)=TSTDAV(I,K)+VSUM2(I) 532 TDAV(I,K+1)=TDAV(I,K)+VSUM1(I) 533 165 CONTINUE 534 163 CONTINUE 535 C 536 C****EVALUATE COEFFICIENTS FOR CO2 PRESSURE INTERPOLATION (A1,A2) 537 DO 171 I=MYIS,MYIE Page 9 Source Listing LWR88 2025-03-12 18:22 LWR88.F 538 A1(I)=(PRESS(I,LP1)-P0XZP8)/P0XZP2 539 A2(I)=(P0-PRESS(I,LP1))/P0XZP2 540 171 CONTINUE 541 C***PERFORM CO2 PRESSURE INTERPOLATION ON ALL INPUTTED TRANSMISSION 542 C FUNCTIONS AND TEMP. DERIVATIVES 543 C---SUCCESSIVELY COMPUTING CO2R,DCO2DT AND D2CDT2 IS DONE TO SAVE 544 C STORAGE (AT A SLIGHT LOSS IN COMPUTATION TIME) 545 DO 184 K=1,LP1 546 DO 184 I=MYIS,MYIE 547 CO2R1(I,K)=A1(I)*CO231(K)+A2(I)*CO238(K) 548 D2CD21(I,K)=H1M3*(A1(I)*C2D31(K)+A2(I)*C2D38(K)) 549 DCO2D1(I,K)=H1M2*(A1(I)*CDT31(K)+A2(I)*CDT38(K)) 550 CO2R2(I,K)=A1(I)*CO271(K)+A2(I)*CO278(K) 551 D2CD22(I,K)=H1M3*(A1(I)*C2D71(K)+A2(I)*C2D78(K)) 552 DCO2D2(I,K)=H1M2*(A1(I)*CDT71(K)+A2(I)*CDT78(K)) 553 184 CONTINUE 554 DO 190 K=1,L 555 DO 190 I=MYIS,MYIE 556 CO2MR(I,K)=A1(I)*CO2M51(K)+A2(I)*CO2M58(K) 557 CO2MD(I,K)=H1M2*(A1(I)*CDTM51(K)+A2(I)*CDTM58(K)) 558 CO2M2D(I,K)=H1M3*(A1(I)*C2DM51(K)+A2(I)*C2DM58(K)) 559 190 CONTINUE 560 C***COMPUTE CO2 TEMPERATURE INTERPOLATIONS FOR ALL BANDS,USING DIFT 561 C 562 C THE CASE WHERE K=1 IS HANDLED FIRST. WE ARE NOW REPLACING 563 C 3-DIMENSIONAL ARRAYS BY 2-D ARRAYS, TO SAVE SPACE. THUS THIS 564 C CALCULATION IS FOR (I,KP,1) 565 DO 211 KP=2,LP1 566 DO 211 I=MYIS,MYIE 567 DIFT(I,KP)=TDAV(I,KP)/TSTDAV(I,KP) 568 211 CONTINUE 569 DO 212 I=MYIS,MYIE 570 CO21(I,1,1)=1.0 571 CO2SP1(I,1)=1.0 572 CO2SP2(I,1)=1.0 573 212 CONTINUE 574 DO 215 KP=2,LP1 575 DO 215 I=MYIS,MYIE 576 C---CALCULATIONS FOR KP>1 FOR K=1 577 CO2R(I,KP)=A1(I)*CO251(KP,1)+A2(I)*CO258(KP,1) 578 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(KP,1)+A2(I)*CDT58(KP,1)) 579 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(KP,1)+A2(I)*C2D58(KP,1)) 580 CO21(I,KP,1)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 581 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 582 C---CALCULATIONS FOR (EFFECTIVELY) KP=1,K>KP. THESE USE THE 583 C SAME VALUE OF DIFT DUE TO SYMMETRY 584 CO2R(I,KP)=A1(I)*CO251(1,KP)+A2(I)*CO258(1,KP) 585 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(1,KP)+A2(I)*CDT58(1,KP)) 586 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(1,KP)+A2(I)*C2D58(1,KP)) 587 CO21(I,1,KP)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 588 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 589 215 CONTINUE 590 C THE TRANSMISSION FUNCTIONS USED IN SPA88 MAY BE COMPUTED NOW. 591 C---(IN THE 250 LOOP,DIFT REALLY SHOULD BE (I,1,K), BUT DIFT IS 592 C INVARIANT WITH RESPECT TO K,KP,AND SO (I,1,K)=(I,K,1)) 593 DO 250 K=2,LP1 594 DO 250 I=MYIS,MYIE Page 10 Source Listing LWR88 2025-03-12 18:22 LWR88.F 595 CO2SP1(I,K)=CO2R1(I,K)+DIFT(I,K)*(DCO2D1(I,K)+HAF*DIFT(I,K)* 596 1 D2CD21(I,K)) 597 CO2SP2(I,K)=CO2R2(I,K)+DIFT(I,K)*(DCO2D2(I,K)+HAF*DIFT(I,K)* 598 1 D2CD22(I,K)) 599 250 CONTINUE 600 C 601 C NEXT THE CASE WHEN K=2...L 602 DO 220 K=2,L 603 DO 222 KP=K+1,LP1 604 DO 222 I=MYIS,MYIE 605 DIFT(I,KP)=(TDAV(I,KP)-TDAV(I,K))/ 606 1 (TSTDAV(I,KP)-TSTDAV(I,K)) 607 CO2R(I,KP)=A1(I)*CO251(KP,K)+A2(I)*CO258(KP,K) 608 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(KP,K)+A2(I)*CDT58(KP,K)) 609 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(KP,K)+A2(I)*C2D58(KP,K)) 610 CO21(I,KP,K)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 611 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 612 CO2R(I,KP)=A1(I)*CO251(K,KP)+A2(I)*CO258(K,KP) 613 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(K,KP)+A2(I)*CDT58(K,KP)) 614 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(K,KP)+A2(I)*C2D58(K,KP)) 615 CO21(I,K,KP)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+ 616 1 HAF*DIFT(I,KP)*D2CDT2(I,KP)) 617 222 CONTINUE 618 220 CONTINUE 619 C FINALLY THE CASE WHEN K=KP,K=2..LP1 620 DO 206 K=2,LP1 621 DO 206 I=MYIS,MYIE 622 DIFT(I,K)=HAF*(VSUM3(I,K)+VSUM3(I,K-1)) 623 CO2R(I,K)=A1(I)*CO251(K,K)+A2(I)*CO258(K,K) 624 DCO2DT(I,K)=H1M2*(A1(I)*CDT51(K,K)+A2(I)*CDT58(K,K)) 625 D2CDT2(I,K)=H1M3*(A1(I)*C2D51(K,K)+A2(I)*C2D58(K,K)) 626 CO21(I,K,K)=CO2R(I,K)+DIFT(I,K)*(DCO2DT(I,K)+ 627 1 HAF*DIFT(I,K)*D2CDT2(I,K)) 628 206 CONTINUE 629 C--- WE AREN'T DOING NBL TFS ON THE 100 CM-1 BANDS . 630 DO 260 K=1,L 631 DO 260 I=MYIS,MYIE 632 CO2NBL(I,K)=CO2MR(I,K)+VSUM3(I,K)*(CO2MD(I,K)+HAF* 633 1 VSUM3(I,K)*CO2M2D(I,K)) 634 260 CONTINUE 635 C***COMPUTE TEMP. COEFFICIENT BASED ON T(K) (SEE REF.2) 636 DO 264 K=1,LP1 637 DO 264 I=MYIS,MYIE 638 IF (T(I,K).LE.H25E2) THEN 639 TLSQU(I,K)=B0+(T(I,K)-H25E2)* 640 1 (B1+(T(I,K)-H25E2)* 641 2 (B2+B3*(T(I,K)-H25E2))) 642 ELSE 643 TLSQU(I,K)=B0 644 ENDIF 645 264 CONTINUE 646 C***APPLY TO ALL CO2 TFS 647 DO 280 K=1,LP1 648 DO 282 KP=1,LP1 649 DO 282 I=MYIS,MYIE 650 CO21(I,KP,K)=CO21(I,KP,K)*(ONE-TLSQU(I,KP))+TLSQU(I,KP) 651 282 CONTINUE Page 11 Source Listing LWR88 2025-03-12 18:22 LWR88.F 652 280 CONTINUE 653 DO 284 K=1,LP1 654 DO 286 I=MYIS,MYIE 655 CO2SP1(I,K)=CO2SP1(I,K)*(ONE-TLSQU(I,1))+TLSQU(I,1) 656 CO2SP2(I,K)=CO2SP2(I,K)*(ONE-TLSQU(I,1))+TLSQU(I,1) 657 286 CONTINUE 658 284 CONTINUE 659 DO 288 K=1,L 660 DO 290 I=MYIS,MYIE 661 CO2NBL(I,K)=CO2NBL(I,K)*(ONE-TLSQU(I,K))+TLSQU(I,K) 662 290 CONTINUE 663 288 CONTINUE 664 CALL FST88(HEATRA,GRNFLX,TOPFLX, 665 1 QH2O,PRESS,P,DELP,DELP2,TEMP,T, 666 2 CLDFAC,NCLDS,KTOP,KBTM,CAMT, 667 3 CO21,CO2NBL,CO2SP1,CO2SP2, 668 4 VAR1,VAR2,VAR3,VAR4,CNTVAL, 669 5 TOTO3,TPHIO3,TOTPHI,TOTVO2, 670 6 EMX1,EMX2,EMPL) 671 RETURN 672 END ENTRY POINTS Name lwr88_ Page 12 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 103 Label 431 427,428 105 Label 437 432 107 Label 442 438,439 125 Label 450 445,446 131 Label 477 458,459 161 Label 522 519 162 Label 526 523,524 163 Label 534 527 165 Label 533 528 171 Label 540 537 184 Label 553 545,546 190 Label 559 554,555 201 Label 484 479 203 Label 491 485,486 206 Label 628 620,621 211 Label 568 565,566 212 Label 573 569 215 Label 589 574,575 220 Label 618 602 222 Label 617 603,604 250 Label 599 593,594 260 Label 634 630,631 264 Label 645 636,637 280 Label 652 647 282 Label 651 648,649 284 Label 658 653 286 Label 657 654 288 Label 663 659 290 Label 662 660 801 Label 500 497 811 Label 506 503,504 812 Label 511 507,508 821 Label 515 512 A1 Local 415 R(4) 4 1 71 538,547,548,549,550,551,552,556,55 7,558,577,578,579,584,585,586,607, 608,609,612,613,614,623,624,625 A2 Local 415 R(4) 4 1 71 539,547,548,549,550,551,552,556,55 7,558,577,578,579,584,585,586,607, 608,609,612,613,614,623,624,625 BANDTA Common 311 5900 BDCOMB Common 379 736 BDWIDE Common 339 56 CAMT Dummy 4 R(4) 4 2 3266 ARG,INOUT 666 CLDFAC Dummy 3 R(4) 4 3 150236 ARG,INOUT 666 CNTVAL Local 400 R(4) 4 2 3266 475,490,668 CO21 Local 396 R(4) 4 3 150236 570,580,587,610,615,626,650,667 CO2BD2 Common 251 1104 CO2BD3 Common 234 52248 CO2BD4 Common 266 1104 CO2BD5 Common 276 368 CO2M2D Local 411 R(4) 4 2 3195 558,633 Page 13 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References CO2MD Local 410 R(4) 4 2 3195 557,632 CO2MR Local 410 R(4) 4 2 3195 556,632 CO2NBL Local 396 R(4) 4 2 3195 632,661,667 CO2R Local 406 R(4) 4 2 3266 577,580,584,587,607,610,612,615,62 3,626 CO2R1 Local 407 R(4) 4 2 3266 547,595 CO2R2 Local 409 R(4) 4 2 3266 550,597 CO2SP1 Local 397 R(4) 4 2 3266 571,595,655,667 CO2SP2 Local 397 R(4) 4 2 3266 572,597,656,667 D2CD21 Local 408 R(4) 4 2 3266 548,596 D2CD22 Local 408 R(4) 4 2 3266 551,598 D2CDT2 Local 416 R(4) 4 2 3266 579,581,586,588,609,611,614,616,62 5,627 DCO2D1 Local 407 R(4) 4 2 3266 549,595 DCO2D2 Local 409 R(4) 4 2 3266 552,597 DCO2DT Local 416 R(4) 4 2 3266 578,580,585,587,608,610,613,615,62 4,626 DELP Local 395 R(4) 4 2 3195 441,665 DELP2 Local 392 R(4) 4 2 3195 440,441,464,465,529,665 DIFT Local 406 R(4) 4 2 3266 567,580,581,587,588,595,597,605,61 0,611,615,616,622,626,627 EMPL Local 404 R(4) 4 2 6461 505,509,513,514,670 EMX1 Local 403 R(4) 4 1 71 498,670 EMX2 Local 404 R(4) 4 1 71 499,670 EXP Func 449 scalar 449 FST88 Subr 664 664 GLB_TABLE Common 152 336 SAVE GRNFLX Dummy 2 R(4) 4 1 71 ARG,INOUT 664 HCON Common 34 872 HEATRA Dummy 2 R(4) 4 2 3195 ARG,INOUT 664 I Local 428 I(4) 4 scalar 428,429,430,432,433,434,435,436,43 9,440,441,446,447,449,459,460,463, 464,465,466,467,475,476,479,480,48 1,482,483,486,487,488,489,490,497, 498,499,504,505,508,509,512,513,51 4,519,520,521,524,525,528,529,530, 531,532,537,538,539,546,547,548,54 9,550,551,552,555,556,557,558,566, 567,569,570,571,572,575,577,578,57 9,580,581,584,585,586,587,588,594, 595,596,597,598,604,605,606,607,60 8,609,610,611,612,613,614,615,616, 621,622,623,624,625,626,627,631,63 2,633,637,638,639,640,641,643,649, 650,654,655,656,660,661 IDIM1 Param 95 I(4) 4 scalar 385,386,387,388,389,390,391,392,39 4,395,396,397,398,399,400,401,402, 403,404,406,407,408,409,410,411,41 2,413,414,415,416,418,419 IDIM2 Param 95 I(4) 4 scalar 385,386,387,388,389,390,391,392,39 4,395,396,397,398,399,400,401,402, 403,404,406,407,408,409,410,411,41 2,413,414,415,416,418,419 IGSTL Param 90 I(4) 4 scalar 95 IGSTR Param 90 I(4) 4 scalar 95 Page 14 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References IM Param 78 I(4) 4 scalar 92,95,158,159,160,161,162,166,167, 172 INLTE Param 186 I(4) 4 scalar 186 INLTEP Param 186 I(4) 4 scalar INPES Param 87 I(4) 4 scalar 92,95,142,143,144,153,154 ITAIL Param 92 I(4) 4 scalar JDIM1 Param 96 I(4) 4 scalar JDIM2 Param 96 I(4) 4 scalar JGSTL Param 91 I(4) 4 scalar 96 JGSTR Param 91 I(4) 4 scalar 96 JM Param 78 I(4) 4 scalar 93,96,158,159,160,161,162,166,167, 172 JNPES Param 87 I(4) 4 scalar 93,96,142,143,144,153,154 JTAIL Param 93 I(4) 4 scalar K Local 427 I(4) 4 scalar 427,429,430,438,440,441,445,447,44 9,458,460,463,464,465,466,467,475, 476,485,487,488,489,490,503,505,50 7,509,523,525,527,529,530,531,532, 545,547,548,549,550,551,552,554,55 6,557,558,593,595,596,597,598,602, 603,605,606,607,608,609,610,612,61 3,614,615,620,622,623,624,625,626, 627,630,632,633,636,638,639,640,64 1,643,647,650,653,655,656,659,661 KBTM Dummy 4 I(4) 4 2 3266 ARG,INOUT 666 KO2 Param 188 I(4) 4 scalar 189 KO21 Param 189 I(4) 4 scalar KO2M Param 189 I(4) 4 scalar KP Local 565 I(4) 4 scalar 565,567,574,577,578,579,580,581,58 4,585,586,587,588,603,605,606,607, 608,609,610,611,612,613,614,615,61 6,648,650 KTOP Dummy 4 I(4) 4 2 3266 ARG,INOUT 666 L Param 174 I(4) 4 scalar 177,178,179,182,184,235,236,386,39 0,392,394,395,396,398,399,410,411, 413,419,427,438,458,498,499,503,51 3,527,554,602,630,659 LL Param 179 I(4) 4 scalar 179,180,514 LL3P Param 184 I(4) 4 scalar LLM1 Param 180 I(4) 4 scalar LLM2 Param 180 I(4) 4 scalar LLM3 Param 180 I(4) 4 scalar LLP1 Param 179 I(4) 4 scalar 404,514 LLP2 Param 179 I(4) 4 scalar LLP3 Param 179 I(4) 4 scalar LM Param 78 I(4) 4 scalar 167,174 LM1 Param 178 I(4) 4 scalar 507 LM2 Param 178 I(4) 4 scalar LM3 Param 178 I(4) 4 scalar LP1 Param 177 I(4) 4 scalar 181,182,183,234,235,237,251,252,26 6,267,276,385,387,388,389,394,395, 396,397,400,401,402,403,406,407,40 8,409,412,413,416,418,434,436,445, 485,499,523,538,539,545,565,574,59 3,603,620,636,647,648,653 Page 15 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Object Declared Type Bytes Dimen Elements Attributes References LP121 Param 183 I(4) 4 scalar LP1M Param 181 I(4) 4 scalar 181 LP1M1 Param 181 I(4) 4 scalar LP1V Param 182 I(4) 4 scalar LP2 Param 177 I(4) 4 scalar 509 LP3 Param 177 I(4) 4 scalar LSM Param 78 I(4) 4 scalar LWR88 Subr 2 MAPPINGS Common 171 5024 SAVE MPPCOM Common 116 1724 SAVE NB Param 185 I(4) 4 scalar 187 NB1 Param 187 I(4) 4 scalar NBLM Param 176 I(4) 4 scalar NBLW Param 175 I(4) 4 scalar 311,312,313 NBLX Param 175 I(4) 4 scalar NBLY Param 175 I(4) 4 scalar 176,183,379,380,381 NCLDS Dummy 4 I(4) 4 1 71 ARG,INOUT 666 NNLTE Param 186 I(4) 4 scalar P Local 395 R(4) 4 2 3266 429,433,434,440,463,498,499,505,50 9,665 PHYCON Common 29 80 PRESS Dummy 3 R(4) 4 2 3266 ARG,INOUT 429,434,498,499,505,509,538,539,66 5 QH2O Local 394 R(4) 4 2 3195 460,464,498,499,505,509,665 QO3 Dummy 3 R(4) 4 2 3195 ARG,INOUT 465 RH2O Dummy 3 R(4) 4 2 3195 ARG,INOUT 460,475,476 T Local 394 R(4) 4 2 3266 430,435,436,638,639,640,641,665 TDAV Local 412 R(4) 4 2 3266 521,532,567,605 TEMP Dummy 3 R(4) 4 2 3266 ARG,INOUT 430,435,436,447,525,665 TEMPCOM Common 157 6603768 SAVE TEXPSL Local 418 R(4) 4 2 3266 447,449,475 TLSQU Local 418 R(4) 4 2 3266 639,643,650,655,656,661 TOPFLX Dummy 2 R(4) 4 1 71 ARG,INOUT 664 TOPO Common 165 17478548 SAVE TOTO3 Local 401 R(4) 4 2 3266 481,488,669 TOTPHI Local 402 R(4) 4 2 3266 480,487,669 TOTVO2 Local 403 R(4) 4 2 3266 483,490,669 TPHIO3 Local 401 R(4) 4 2 3266 482,489,669 TSTDAV Local 412 R(4) 4 2 3266 520,531,567,606 VAR1 Local 398 R(4) 4 2 3195 464,466,668 VAR2 Local 398 R(4) 4 2 3195 466,475,487,513,668 VAR3 Local 399 R(4) 4 2 3195 465,467,488,668 VAR4 Local 399 R(4) 4 2 3195 467,489,668 VSUM1 Local 413 R(4) 4 1 71 530,532 VSUM2 Local 414 R(4) 4 1 71 529,530,531 VSUM3 Local 413 R(4) 4 2 3266 525,530,622,632,633 VSUM4 Local 419 R(4) 4 2 3195 VV Local 413 R(4) 4 2 3195 463,466,467 Page 16 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References AB15 R(4) 4 5892 1 2 COM AB15CM R(4) 4 728 1 2 COM AB15WD R(4) 4 40 scalar COM ACOMB R(4) 4 160 1 15 COM AMOLWT R(4) 4 0 scalar COM AO3CM R(4) 4 704 1 3 COM AO3RND R(4) 4 5868 1 3 COM AP R(4) 4 1956 1 163 COM APCM R(4) 4 340 1 15 COM APWD R(4) 4 12 scalar COM ARNDM R(4) 4 0 1 163 COM ATP R(4) 4 3260 1 163 COM ATPCM R(4) 4 460 1 15 COM ATPWD R(4) 4 20 scalar COM AWIDE R(4) 4 0 scalar COM B0 R(4) 4 52232 scalar COM 639,643 B1 R(4) 4 52236 scalar COM 640 B2 R(4) 4 52240 scalar COM 641 B3 R(4) 4 52244 scalar COM 641 BANDHI R(4) 4 5216 1 163 COM BANDLO R(4) 4 4564 1 163 COM BCOMB R(4) 4 220 1 15 COM BDHICM R(4) 4 640 1 15 COM BDHIWD R(4) 4 32 scalar COM BDLOCM R(4) 4 580 1 15 COM BDLOWD R(4) 4 28 scalar COM BETACM R(4) 4 280 1 15 COM BETAD R(4) 4 1304 1 163 COM BETAWD R(4) 4 8 scalar COM BETINC R(4) 4 700 scalar COM BETINW R(4) 4 36 scalar COM 475 BO3CM R(4) 4 716 1 3 COM BO3RND R(4) 4 5880 1 3 COM BP R(4) 4 2608 1 163 COM BPCM R(4) 4 400 1 15 COM BPWD R(4) 4 16 scalar COM BRNDM R(4) 4 652 1 163 COM BTP R(4) 4 3912 1 163 COM BTPCM R(4) 4 520 1 15 COM BTPWD R(4) 4 24 scalar COM BWIDE R(4) 4 4 scalar COM C2D31 R(4) 4 736 1 46 COM 548 C2D38 R(4) 4 920 1 46 COM 548 C2D51 R(4) 4 33856 2 2116 COM 579,586,609,614,625 C2D58 R(4) 4 42320 2 2116 COM 579,586,609,614,625 C2D71 R(4) 4 736 1 46 COM 551 C2D78 R(4) 4 920 1 46 COM 551 C2DM51 R(4) 4 51504 1 45 COM 558 C2DM58 R(4) 4 51684 1 45 COM 558 CDT31 R(4) 4 368 1 46 COM 549 CDT38 R(4) 4 552 1 46 COM 549 CDT51 R(4) 4 16928 2 2116 COM 578,585,608,613,624 CDT58 R(4) 4 25392 2 2116 COM 578,585,608,613,624 Page 17 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References CDT71 R(4) 4 368 1 46 COM 552 CDT78 R(4) 4 552 1 46 COM 552 CDTM51 R(4) 4 51144 1 45 COM 557 CDTM58 R(4) 4 51324 1 45 COM 557 CO211 R(4) 4 0 1 46 COM CO218 R(4) 4 184 1 46 COM CO231 R(4) 4 0 1 46 COM 547 CO238 R(4) 4 184 1 46 COM 547 CO251 R(4) 4 0 2 2116 COM 577,584,607,612,623 CO258 R(4) 4 8464 2 2116 COM 577,584,607,612,623 CO271 R(4) 4 0 1 46 COM 550 CO278 R(4) 4 184 1 46 COM 550 CO2M51 R(4) 4 50784 1 45 COM 556 CO2M58 R(4) 4 50964 1 45 COM 556 CSUBP R(4) 4 4 scalar COM DIFFCTR R(4) 4 8 scalar COM 460,465 EIGHT R(4) 4 20 scalar COM FIFTY R(4) 4 12 scalar COM FIVE R(4) 4 24 scalar COM FOUR R(4) 4 28 scalar COM G R(4) 4 12 scalar COM G2LI I(4) 4 0 1 239 COM G2LJ I(4) 4 1912 1 389 COM GINV R(4) 4 68 scalar COM 464,465 GP0INV R(4) 4 76 scalar COM 498,499,505,510 GRAVDR R(4) 4 16 scalar COM GTEMP R(4) 4 52048 1 46 COM 529 H101M16 R(4) 4 808 scalar COM H102M5 R(4) 4 784 scalar COM H1036E2 R(4) 4 756 scalar COM H114M11 R(4) 4 384 scalar COM H1174M7 R(4) 4 268 scalar COM H11M10 R(4) 4 296 scalar COM H11M11 R(4) 4 388 scalar COM H1224E3 R(4) 4 112 scalar COM H1226E1 R(4) 4 652 scalar COM H128M5 R(4) 4 240 scalar COM H129M2 R(4) 4 732 scalar COM H12M12 R(4) 4 472 scalar COM H12M13 R(4) 4 492 scalar COM H12M31 R(4) 4 560 scalar COM H135M13 R(4) 4 488 scalar COM H1386E2 R(4) 4 748 scalar COM H1439M5 R(4) 4 236 scalar COM H14M10 R(4) 4 292 scalar COM H14M11 R(4) 4 380 scalar COM H14M12 R(4) 4 468 scalar COM H14M14 R(4) 4 508 scalar COM H14M30 R(4) 4 548 scalar COM H15E2 R(4) 4 136 scalar COM H15M11 R(4) 4 376 scalar COM H15M14 R(4) 4 504 scalar COM H15M5 R(4) 4 824 scalar COM H161E1 R(4) 4 800 scalar COM H165E5 R(4) 4 84 scalar COM Page 18 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References H16E1 R(4) 4 792 scalar COM H16M12 R(4) 4 464 scalar COM H181E1 R(4) 4 148 scalar COM H18E1 R(4) 4 152 scalar COM H18E3 R(4) 4 832 scalar COM 447 H18M11 R(4) 4 372 scalar COM H1E11 R(4) 4 72 scalar COM H1E13 R(4) 4 68 scalar COM H1E15 R(4) 4 64 scalar COM H1E4 R(4) 4 96 scalar COM H1E6 R(4) 4 612 scalar COM H1E8 R(4) 4 76 scalar COM H1M10 R(4) 4 300 scalar COM H1M11 R(4) 4 392 scalar COM H1M13 R(4) 4 496 scalar COM H1M16 R(4) 4 856 scalar COM H1M17 R(4) 4 512 scalar COM H1M18 R(4) 4 516 scalar COM H1M19 R(4) 4 520 scalar COM H1M2 R(4) 4 620 scalar COM 549,552,557,578,585,608,613,624 H1M20 R(4) 4 524 scalar COM H1M21 R(4) 4 528 scalar COM H1M22 R(4) 4 532 scalar COM H1M23 R(4) 4 536 scalar COM H1M24 R(4) 4 540 scalar COM H1M3 R(4) 4 208 scalar COM 548,551,558,579,586,609,614,625 H1M4 R(4) 4 220 scalar COM H1M5 R(4) 4 244 scalar COM H1M6 R(4) 4 260 scalar COM H1M60 R(4) 4 584 scalar COM H1M8 R(4) 4 284 scalar COM H1P082 R(4) 4 740 scalar COM H1P25892 R(4) 4 180 scalar COM H1P4 R(4) 4 176 scalar COM H1P41819 R(4) 4 704 scalar COM H1P4387 R(4) 4 172 scalar COM H1P8 R(4) 4 168 scalar COM H2075E3 R(4) 4 108 scalar COM H20788E3 R(4) 4 104 scalar COM H2118M2 R(4) 4 760 scalar COM H21M12 R(4) 4 460 scalar COM H21M31 R(4) 4 556 scalar COM H235M3 R(4) 4 724 scalar COM H23E2 R(4) 4 132 scalar COM H23M10 R(4) 4 288 scalar COM H23M11 R(4) 4 364 scalar COM H24E3 R(4) 4 100 scalar COM H24M11 R(4) 4 360 scalar COM H24M12 R(4) 4 456 scalar COM H25452M6 R(4) 4 256 scalar COM H257M8 R(4) 4 280 scalar COM H25E2 R(4) 4 816 scalar COM 638,639,640,641 H25M31 R(4) 4 552 scalar COM H26E2 R(4) 4 696 scalar COM H26M30 R(4) 4 544 scalar COM Page 19 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References H28571M2 R(4) 4 852 scalar COM H285M4 R(4) 4 216 scalar COM H28E1 R(4) 4 868 scalar COM H28M11 R(4) 4 356 scalar COM H28M12 R(4) 4 452 scalar COM H29316E2 R(4) 4 648 scalar COM H2945E2 R(4) 4 128 scalar COM H29M12 R(4) 4 448 scalar COM H2E2 R(4) 4 688 scalar COM H2E6 R(4) 4 616 scalar COM H2M11 R(4) 4 368 scalar COM H2P5 R(4) 4 164 scalar COM H2P8 R(4) 4 160 scalar COM H2P9 R(4) 4 156 scalar COM H3082E2 R(4) 4 120 scalar COM H3116E1 R(4) 4 656 scalar COM H323M4 R(4) 4 768 scalar COM H327M8 R(4) 4 276 scalar COM H32M11 R(4) 4 348 scalar COM H35E1 R(4) 4 140 scalar COM H35M11 R(4) 4 344 scalar COM H36M13 R(4) 4 484 scalar COM H37412M5 R(4) 4 232 scalar COM H37M11 R(4) 4 340 scalar COM H37M12 R(4) 4 440 scalar COM H38M12 R(4) 4 436 scalar COM H391M7 R(4) 4 264 scalar COM H394M5 R(4) 4 228 scalar COM H3E2 R(4) 4 124 scalar COM H3M11 R(4) 4 352 scalar COM H3M12 R(4) 4 444 scalar COM H3M14 R(4) 4 500 scalar COM H3M3 R(4) 4 804 scalar COM 467 H3M4 R(4) 4 860 scalar COM 466 H3P5 R(4) 4 828 scalar COM H3P6 R(4) 4 144 scalar COM H41666M2 R(4) 4 636 scalar COM H42M11 R(4) 4 336 scalar COM H42M2 R(4) 4 764 scalar COM H44194M2 R(4) 4 700 scalar COM H44871M2 R(4) 4 200 scalar COM H44M11 R(4) 4 332 scalar COM H44M12 R(4) 4 428 scalar COM H451M6 R(4) 4 788 scalar COM H45M12 R(4) 4 424 scalar COM H45M32 R(4) 4 572 scalar COM H46M13 R(4) 4 480 scalar COM H488E4 R(4) 4 92 scalar COM H48M11 R(4) 4 328 scalar COM H4999M6 R(4) 4 252 scalar COM H4E5 R(4) 4 80 scalar COM H4M12 R(4) 4 432 scalar COM H4M33 R(4) 4 576 scalar COM H53M11 R(4) 4 324 scalar COM H559M3 R(4) 4 204 scalar COM Page 20 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References H55M32 R(4) 4 568 scalar COM H5725E4 R(4) 4 88 scalar COM H5E2 R(4) 4 116 scalar COM H625M2 R(4) 4 668 scalar COM H62M12 R(4) 4 416 scalar COM H62M34 R(4) 4 580 scalar COM H658M2 R(4) 4 752 scalar COM H65M12 R(4) 4 412 scalar COM H67390E2 R(4) 4 772 scalar COM H6938M5 R(4) 4 224 scalar COM H69766E5 R(4) 4 720 scalar COM H6M12 R(4) 4 420 scalar COM H6P08108 R(4) 4 836 scalar COM 447 H71E26 R(4) 4 60 scalar COM H72M11 R(4) 4 320 scalar COM H74M12 R(4) 4 408 scalar COM H75826M4 R(4) 4 736 scalar COM H77M11 R(4) 4 316 scalar COM H77M12 R(4) 4 404 scalar COM H7M6 R(4) 4 248 scalar COM H8121E1 R(4) 4 684 scalar COM H82M11 R(4) 4 308 scalar COM H83E26 R(4) 4 56 scalar COM H83M11 R(4) 4 304 scalar COM H8725M8 R(4) 4 272 scalar COM H8M11 R(4) 4 312 scalar COM H8M13 R(4) 4 476 scalar COM H93M12 R(4) 4 400 scalar COM H96M12 R(4) 4 396 scalar COM H987M4 R(4) 4 212 scalar COM H9M32 R(4) 4 564 scalar COM H9P94 R(4) 4 660 scalar COM HAF R(4) 4 44 scalar COM 429,430,463,581,588,595,597,611,61 6,622,627,632 HM13EZ R(4) 4 592 scalar COM HM1597E1 R(4) 4 812 scalar COM HM161E1 R(4) 4 796 scalar COM HM1797E1 R(4) 4 680 scalar COM HM181E1 R(4) 4 604 scalar COM HM19EZ R(4) 4 596 scalar COM HM1E1 R(4) 4 600 scalar COM HM1E2 R(4) 4 608 scalar COM HM1EZ R(4) 4 692 scalar COM HM2M2 R(4) 4 644 scalar COM HM6666M2 R(4) 4 628 scalar COM HM8E1 R(4) 4 864 scalar COM HMP5 R(4) 4 640 scalar COM HMP575 R(4) 4 588 scalar COM HMP66667 R(4) 4 624 scalar COM HMP805 R(4) 4 840 scalar COM HNINETY R(4) 4 4 scalar COM HP1 R(4) 4 196 scalar COM HP118666 R(4) 4 820 scalar COM HP144 R(4) 4 712 scalar COM HP166666 R(4) 4 632 scalar COM Page 21 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References HP219 R(4) 4 708 scalar COM HP228 R(4) 4 672 scalar COM HP26 R(4) 4 728 scalar COM HP369 R(4) 4 192 scalar COM HP3795 R(4) 4 776 scalar COM HP5048 R(4) 4 780 scalar COM HP518 R(4) 4 188 scalar COM HP526315 R(4) 4 848 scalar COM HP6 R(4) 4 664 scalar COM HP602409 R(4) 4 844 scalar COM HP60241 R(4) 4 676 scalar COM HP8 R(4) 4 184 scalar COM HP805 R(4) 4 744 scalar COM HP816 R(4) 4 716 scalar COM HTMG R(4) 4 743768 3 4183695 COM HUNDRED R(4) 4 0 scalar COM IBAND I(4) 4 0 1 40 COM IBROW I(4) 4 404 scalar COM ICHUNKTAB I(4) 4 828 1 21 COM IE_GLB_TABLE I(4) 4 84 1 21 COM IE_LOC_TABLE I(4) 4 660 1 21 COM ILCOL I(4) 4 396 scalar COM ILPAD1 I(4) 4 412 scalar COM ILPAD2 I(4) 4 416 scalar COM ILPAD3 I(4) 4 420 scalar COM ILPAD4 I(4) 4 424 scalar COM ILPAD5 I(4) 4 428 scalar COM INUMQ I(4) 4 1324 1 100 COM IQUILT_GROUP I(4) 4 1320 scalar COM IRCOL I(4) 4 400 scalar COM IRPAD1 I(4) 4 432 scalar COM IRPAD2 I(4) 4 436 scalar COM IRPAD3 I(4) 4 440 scalar COM IRPAD4 I(4) 4 444 scalar COM IRPAD5 I(4) 4 448 scalar COM IS_GLB_TABLE I(4) 4 0 1 21 COM IS_LOC_TABLE I(4) 4 492 1 21 COM ITEMP I(4) 4 5870016 2 91719 COM ITEMP2 I(4) 4 6236892 2 91719 COM ITROW I(4) 4 408 scalar COM JBPAD1 I(4) 4 452 scalar COM JBPAD2 I(4) 4 456 scalar COM JBPAD3 I(4) 4 460 scalar COM JBPAD4 I(4) 4 464 scalar COM JBPAD5 I(4) 4 468 scalar COM JE_GLB_TABLE I(4) 4 252 1 21 COM JE_LOC_TABLE I(4) 4 744 1 21 COM JS_GLB_TABLE I(4) 4 168 1 21 COM JS_LOC_TABLE I(4) 4 576 1 21 COM JTPAD1 I(4) 4 472 scalar COM JTPAD2 I(4) 4 476 scalar COM JTPAD3 I(4) 4 480 scalar COM JTPAD4 I(4) 4 484 scalar COM JTPAD5 I(4) 4 488 scalar COM L2GI I(4) 4 956 1 239 COM Page 22 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References L2GJ I(4) 4 3468 1 389 COM MPI_COMM_COMP I(4) 4 912 scalar COM MPI_COMM_INTER I(4) 4 916 scalar COM MPI_COMM_INTER_ARRAY I(4) 4 920 1 100 COM MYIE I(4) 4 52 scalar COM 428,432,439,446,459,479,486,497,50 4,508,512,519,524,528,537,546,555, 566,569,575,594,604,621,631,637,64 9,654,660 MYIE1 I(4) 4 56 scalar COM MYIE1_P1 I(4) 4 120 scalar COM MYIE1_P2 I(4) 4 124 scalar COM MYIE1_P3 I(4) 4 128 scalar COM MYIE1_P4 I(4) 4 132 scalar COM MYIE2 I(4) 4 60 scalar COM MYIE2_P1 I(4) 4 136 scalar COM MYIE_P1 I(4) 4 100 scalar COM MYIE_P2 I(4) 4 104 scalar COM MYIE_P3 I(4) 4 108 scalar COM MYIE_P4 I(4) 4 112 scalar COM MYIE_P5 I(4) 4 116 scalar COM MYIS I(4) 4 40 scalar COM 428,432,439,446,459,479,486,497,50 4,508,512,519,524,528,537,546,555, 566,569,575,594,604,621,631,637,64 9,654,660 MYIS1 I(4) 4 44 scalar COM MYIS1_P1 I(4) 4 84 scalar COM MYIS1_P2 I(4) 4 88 scalar COM MYIS1_P3 I(4) 4 92 scalar COM MYIS1_P4 I(4) 4 96 scalar COM MYIS2 I(4) 4 48 scalar COM MYIS_P1 I(4) 4 64 scalar COM MYIS_P2 I(4) 4 68 scalar COM MYIS_P3 I(4) 4 72 scalar COM MYIS_P4 I(4) 4 76 scalar COM MYIS_P5 I(4) 4 80 scalar COM MYJE I(4) 4 236 scalar COM MYJE1 I(4) 4 240 scalar COM MYJE1_P1 I(4) 4 280 scalar COM MYJE1_P2 I(4) 4 284 scalar COM MYJE1_P3 I(4) 4 288 scalar COM MYJE1_P4 I(4) 4 292 scalar COM MYJE2 I(4) 4 244 scalar COM MYJE2_P1 I(4) 4 296 scalar COM MYJE2_P2 I(4) 4 300 scalar COM MYJE2_P3 I(4) 4 304 scalar COM MYJE2_P4 I(4) 4 308 scalar COM MYJE3 I(4) 4 248 scalar COM MYJE3_P4 I(4) 4 312 scalar COM MYJE4 I(4) 4 252 scalar COM MYJE4_P1 I(4) 4 316 scalar COM MYJE4_P4 I(4) 4 320 scalar COM MYJE5 I(4) 4 256 scalar COM MYJE5_P1 I(4) 4 324 scalar COM MYJE5_P2 I(4) 4 328 scalar COM MYJE_P1 I(4) 4 260 scalar COM Page 23 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References MYJE_P2 I(4) 4 264 scalar COM MYJE_P3 I(4) 4 268 scalar COM MYJE_P4 I(4) 4 272 scalar COM MYJE_P5 I(4) 4 276 scalar COM MYJS I(4) 4 140 scalar COM MYJS1 I(4) 4 144 scalar COM MYJS1_P1 I(4) 4 184 scalar COM MYJS1_P2 I(4) 4 188 scalar COM MYJS1_P3 I(4) 4 192 scalar COM MYJS1_P4 I(4) 4 196 scalar COM MYJS2 I(4) 4 148 scalar COM MYJS2_P1 I(4) 4 200 scalar COM MYJS2_P2 I(4) 4 204 scalar COM MYJS2_P3 I(4) 4 208 scalar COM MYJS2_P4 I(4) 4 212 scalar COM MYJS3 I(4) 4 152 scalar COM MYJS3_P4 I(4) 4 216 scalar COM MYJS4 I(4) 4 156 scalar COM MYJS4_P1 I(4) 4 220 scalar COM MYJS4_P4 I(4) 4 224 scalar COM MYJS5 I(4) 4 160 scalar COM MYJS5_P1 I(4) 4 228 scalar COM MYJS5_P2 I(4) 4 232 scalar COM MYJS_P1 I(4) 4 164 scalar COM MYJS_P2 I(4) 4 168 scalar COM MYJS_P3 I(4) 4 172 scalar COM MYJS_P4 I(4) 4 176 scalar COM MYJS_P5 I(4) 4 180 scalar COM MYPE I(4) 4 0 scalar COM MY_E I(4) 4 336 scalar COM MY_IE_GLB I(4) 4 12 scalar COM MY_IE_LOC I(4) 4 28 scalar COM MY_IS_GLB I(4) 4 8 scalar COM MY_IS_LOC I(4) 4 24 scalar COM MY_JE_GLB I(4) 4 20 scalar COM MY_JE_LOC I(4) 4 36 scalar COM MY_JS_GLB I(4) 4 16 scalar COM MY_JS_LOC I(4) 4 32 scalar COM MY_N I(4) 4 332 scalar COM MY_NE I(4) 4 348 scalar COM MY_NEB I(4) 4 364 1 8 COM MY_NW I(4) 4 360 scalar COM MY_S I(4) 4 340 scalar COM MY_SE I(4) 4 352 scalar COM MY_SW I(4) 4 356 scalar COM MY_W I(4) 4 344 scalar COM NPES I(4) 4 4 scalar COM O3DIFCTR R(4) 4 20 scalar COM ONE R(4) 4 40 scalar COM 441,650,655,656,661 P0 R(4) 4 24 scalar COM 539 P0INV R(4) 4 72 scalar COM 463 P0X2 R(4) 4 36 scalar COM P0XZP2 R(4) 4 28 scalar COM 538,539 P0XZP8 R(4) 4 32 scalar COM 538 QUARTR R(4) 4 48 scalar COM Page 24 Source Listing LWR88 2025-03-12 18:22 Symbol Table LWR88.F Name Type Bytes Offset Dimen Elements Attributes References RADCON R(4) 4 40 scalar COM RADCON1 R(4) 4 64 scalar COM RATCO2MW R(4) 4 56 scalar COM RATH2OMW R(4) 4 60 scalar COM 476 RGAS R(4) 4 44 scalar COM RGASSP R(4) 4 48 scalar COM SECPDA R(4) 4 52 scalar COM SIXTY R(4) 4 8 scalar COM SKC1R R(4) 4 48 scalar COM SKO2D R(4) 4 44 scalar COM SKO3R R(4) 4 52 scalar COM STEMP R(4) 4 51864 1 46 COM 525 TEMP1 R(4) 4 0 2 91719 COM TEMP10 R(4) 4 3301884 2 91719 COM TEMP11 R(4) 4 3668760 2 91719 COM TEMP12 R(4) 4 4035636 2 91719 COM TEMP13 R(4) 4 4402512 2 91719 COM TEMP14 R(4) 4 4769388 2 91719 COM TEMP15 R(4) 4 5136264 2 91719 COM TEMP16 R(4) 4 5503140 2 91719 COM TEMP2 R(4) 4 366876 2 91719 COM TEMP2X R(4) 4 0 2 92971 COM TEMP3 R(4) 4 733752 2 91719 COM TEMP4 R(4) 4 1100628 2 91719 COM TEMP5 R(4) 4 1467504 2 91719 COM TEMP6 R(4) 4 1834380 2 91719 COM TEMP7 R(4) 4 2201256 2 91719 COM TEMP8 R(4) 4 2568132 2 91719 COM TEMP9 R(4) 4 2935008 2 91719 COM TEN R(4) 4 16 scalar COM THREE R(4) 4 32 scalar COM TTVG R(4) 4 371884 2 92971 COM TWO R(4) 4 36 scalar COM ZERO R(4) 4 52 scalar COM 433,480,481,482,483,520,521 Page 25 Source Listing LWR88 2025-03-12 18:22 Subprograms/Common Blocks LWR88.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BANDTA Common 311 5900 BDCOMB Common 379 736 BDWIDE Common 339 56 CO2BD2 Common 251 1104 CO2BD3 Common 234 52248 CO2BD4 Common 266 1104 CO2BD5 Common 276 368 GLB_TABLE Common 152 336 SAVE HCON Common 34 872 LWR88 Subr 2 MAPPINGS Common 171 5024 SAVE MPPCOM Common 116 1724 SAVE PHYCON Common 29 80 TEMPCOM Common 157 6603768 SAVE TOPO Common 165 17478548 SAVE 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 Page 26 Source Listing LWR88 2025-03-12 18:22 LWR88.F -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 -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 Page 27 Source Listing LWR88 2025-03-12 18:22 LWR88.F -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 -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 : LWR88.lst no -o COMPILER: Intel(R) Fortran 19.1-1655