Page 1 Source Listing FST88 2025-03-12 18:22 /tmp/ifortH7H7AE.i 1 # 1 "FST88.F" 2 C ***************************************************************** 3 C SUBROUTINE FST88 IS THE MAIN COMPUTATION MODULE OF THE 4 C LONG-WAVE RADIATION CODE. IN IT ALL "EMISSIVITY" CALCULATIONS, 5 C INCLUDING CALLS TO TABLE LOOKUP SUBROUTINES. ALSO,AFTER CALLING 6 C SUBROUTINE "SPA88", FINAL COMBINED HEATING RATES AND GROUND 7 C FLUX ARE OBTAINED. 8 C ***************************************************************** 9 C INPUTS: 10 C BETINW,BETAWD,AB15WD BDWIDE 11 C BETAD,BO3RND,AO3RND BANDTA 12 C CLDFAC CLDCOM 13 C QH2O,P,DELP2,DELP,T,VAR1,VAR2, KDACOM 14 C VAR3,VAR4,CNTVAL KDACOM 15 C TOTVO2,TOTO3,TOTPHI,EMPL,EMX1 KDACOM 16 C TPHIO3,EMX2 KDACOM 17 C TEMP,PRESS RADISW 18 C NCLDS,KTOP,KBTM,CAMT RADISW 19 C IND,INDX2,KMAXV,SOURCE,DSRCE TABCOM 20 C SKC1R,SKC3R,KMAXVM,NREP1,NREP2 TABCOM 21 C NST1,NST2,NRP1,NRP2 TABCOM 22 C CO2NBL,CO21 TFCOM 23 C CO2SP1,CO2SP2 TFCOM 24 C OUTPUTS: 25 C HEATRA,GRNFLX,TOPFLX LWOUT 26 C 27 C CALLED BY : RADMN OR MAIN PGM 28 C CALLS : CLO88,E1E288,E3V88,SPA88,NLTE 29 C 30 C PASSED VARIABLES: 31 C IN E3V88: 32 C EMD = E3 FUNCTION FOR H2O LINES (0-560,1200-2200 CM-1) 33 C COMPUTED IN E3V88 34 C TPL = TEMPERATURE INPUT FOR E3 CALCULATION IN E3V88 35 C EMPL = H2O AMOUNT,INPUT FOR E3 CALCULATION IN E3V88 36 C (COMPUTED IN LWR88; STORED IN KDACOM.H) 37 C IN E1E288: 38 C E1CTS1 = E1 FUNCTION FOR THE (I+1)TH LEVEL USING THE 39 C TEMPERATURE OF THE ITH DATA LEVEL,COMPUTED OVER 40 C THE FREQUENCY RANGE 0-560,1200-2200 CM-1. (E1CTS1- 41 C E1CTW1) IS USED IN OBTAINING THE FLUX AT THE TOP 42 C IN THE 0-160,1200-2200 CM-1 RANGE (FLX1E1). 43 C E1CTS2 = E1 FUNCTION FOR THE ITH LEVEL, USING THE TEMP. OF 44 C THE ITH DATA LEVEL,COMPUTED OVER THE FREQUENCY RANGE 45 C 0-560,1200-2200 CM-1. (E1CTS2-E1CTW2) IS ALSO USED 46 C IN OBTAINING THE FLUX AT THE TOP IN THE 0-160,. 47 C 1200-2200 CM-1 RANGE. 48 C E1FLX = E1 FCTN. FOR THE ITH LEVEL,USING THE TEMPERATURE AT 49 C THE TOP OF THE ATMOSPHERE. COMPUTED OVER THE FREQ. 50 C RANGE 0-560,1200-2200 CM-1. USED FOR Q(APPROX) TERM. 51 C (IN COMMON BLOCK TFCOM) 52 C E1CTW1 = LIKE E1CTS1,BUT COMPUTED OVER THE 160-560 CM-1 RANGE 53 C AND USED FOR Q(APPROX,CTS) CALCULATION 54 C E1CTW2 = LIKE E1CTS2,BUT COMPUTED OVER THE 160-560 CM-1 RANGE 55 C AND USED FOR Q(APPROX,CTS) CALCULATION 56 C FXO = TEMPERATURE INDEX USED FOR E1 FUNCTION AND ALSO 57 C USED FOR SOURCE FUNCTION CALC. IN FST88. Page 2 Source Listing FST88 2025-03-12 18:22 FST88.F 58 C DT = TEMP. DIFF.BETWEEN MODEL TEMPS. AND TEMPS. AT 59 C TABULAR VALUES OF E1 AND SOURCE FCTNS. USED IN 60 C FST88 AND IN E1 FUNCTION CALC. 61 C FXOE2 = TEMPERATURE INDEX USED FOR E2 FUNCTION 62 C DTE2 = TEMP. DIFF. BETWEEN MODEL TEMP. AND TEMPS. AT 63 C TABULAR VALUES OF E2 FUNCTION. 64 SUBROUTINE FST88(HEATRA,GRNFLX,TOPFLX, 65 1 QH2O,PRESS,P,DELP,DELP2,TEMP,T, 66 2 CLDFAC,NCLDS,KTOP,KBTM,CAMT, 67 3 CO21,CO2NBL,CO2SP1,CO2SP2, 68 4 VAR1,VAR2,VAR3,VAR4,CNTVAL, 69 5 TOTO3,TPHIO3,TOTPHI,TOTVO2, 70 6 EMX1,EMX2,EMPL) 71 C 72 COMMON/PHYCON/AMOLWT,CSUBP,DIFFCTR,G,GRAVDR,O3DIFCTR,P0, 73 * P0XZP2,P0XZP8,P0X2,RADCON,RGAS,RGASSP,SECPDA 74 COMMON/PHYCON/RATCO2MW,RATH2OMW 75 COMMON/PHYCON/RADCON1 76 COMMON/PHYCON/GINV,P0INV,GP0INV 77 save /PHYCON/ 78 COMMON/HCON/HUNDRED,HNINETY,SIXTY,FIFTY,TEN,EIGHT,FIVE, 79 * FOUR,THREE,TWO,ONE,HAF,QUARTR,ZERO 80 COMMON/HCON/H83E26,H71E26,H1E15,H1E13,H1E11,H1E8,H4E5, 81 * H165E5,H5725E4,H488E4,H1E4,H24E3,H20788E3, 82 * H2075E3,H1224E3,H5E2,H3082E2,H3E2,H2945E2, 83 * H23E2,H15E2,H35E1,H3P6,H181E1,H18E1,H2P9,H2P8, 84 * H2P5,H1P8,H1P4387,H1P4,H1P25892,HP8,HP518, 85 * HP369,HP1 86 COMMON/HCON/H44871M2,H559M3,H1M3,H987M4,H285M4,H1M4, 87 * H6938M5,H394M5,H37412M5,H1439M5,H128M5,H1M5, 88 * H7M6,H4999M6,H25452M6,H1M6,H391M7,H1174M7, 89 * H8725M8,H327M8,H257M8,H1M8,H23M10,H14M10, 90 * H11M10,H1M10,H83M11,H82M11,H8M11,H77M11, 91 * H72M11,H53M11,H48M11,H44M11,H42M11,H37M11, 92 * H35M11,H32M11,H3M11,H28M11,H24M11,H23M11, 93 * H2M11,H18M11,H15M11,H14M11,H114M11,H11M11, 94 * H1M11,H96M12,H93M12,H77M12,H74M12,H65M12, 95 * H62M12,H6M12,H45M12,H44M12,H4M12,H38M12, 96 * H37M12,H3M12,H29M12,H28M12,H24M12,H21M12, 97 * H16M12,H14M12,H12M12,H8M13,H46M13,H36M13, 98 * H135M13,H12M13,H1M13,H3M14,H15M14,H14M14, 99 * H1M17,H1M18,H1M19,H1M20,H1M21,H1M22,H1M23, 100 * H1M24,H26M30,H14M30,H25M31,H21M31,H12M31, 101 * H9M32,H55M32,H45M32,H4M33,H62M34,H1M60 102 COMMON/HCON/HMP575,HM13EZ,HM19EZ,HM1E1,HM181E1,HM1E2 103 COMMON/HCON/H1E6,H2E6,H1M2,HMP66667,HM6666M2,HP166666, 104 * H41666M2,HMP5,HM2M2,H29316E2,H1226E1,H3116E1, 105 * H9P94,HP6,H625M2,HP228,HP60241,HM1797E1, 106 * H8121E1,H2E2,HM1EZ,H26E2,H44194M2,H1P41819 107 COMMON/HCON/HP219,HP144,HP816,H69766E5,H235M3,HP26, 108 * H129M2,H75826M4,H1P082,HP805,H1386E2, 109 * H658M2,H1036E2,H2118M2,H42M2,H323M4, 110 * H67390E2,HP3795,HP5048,H102M5,H451M6 111 COMMON/HCON/H16E1,HM161E1,H161E1,H3M3,H101M16, 112 * HM1597E1,H25E2,HP118666,H15M5,H3P5,H18E3, 113 * H6P08108,HMP805,HP602409,HP526315, 114 * H28571M2,H1M16 Page 3 Source Listing FST88 2025-03-12 18:22 FST88.F 115 COMMON/HCON/H3M4 116 COMMON/HCON/HM8E1 117 COMMON/HCON/H28E1 118 save /HCON/ 119 120 C----------------------------------------------------------------------- 121 INCLUDE "parmeta" 145 INCLUDE "mpp.h" 146 # 1 "./sp.h" 1 147 # 4 148 149 # 123 "FST88.F" 2 150 C----------------------------------------------------------------------- 151 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 152 C IMAX = NO. POINTS ALONG THE LAT. CIRCLE USED IN CALCS. 153 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 154 C***NOTE: THE USER NORMALLY WILL MODIFY ONLY THE IMAX AND L PARAMETERS 155 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 156 C BANDTA FOR DEFINITION 157 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 158 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE 159 C BDCOMB FOR DEFINITION 160 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 161 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 162 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 163 C FROM THE ABOVE PARAMETERS. 222 PARAMETER (L=LM) 223 PARAMETER (IMAX=IM,NCOL=IMAX) 224 PARAMETER (NBLW=163,NBLX=47,NBLY=15) 225 PARAMETER (NBLM=NBLY-1) 226 PARAMETER (LP1=L+1,LP2=L+2,LP3=L+3) 227 PARAMETER (LM1=L-1,LM2=L-2,LM3=L-3) 228 PARAMETER (LL=2*L,LLP1=LL+1,LLP2=LL+2,LLP3=LL+3) 229 PARAMETER (LLM1=LL-1,LLM2=LL-2,LLM3=LL-3) 230 PARAMETER (LP1M=LP1*LP1,LP1M1=LP1M-1) 231 PARAMETER (LP1V=LP1*(1+2*L/2)) 232 PARAMETER (LP121=LP1*NBLY) 233 PARAMETER (LL3P=3*L+2) 234 PARAMETER (NB=12) 235 PARAMETER (INLTE=3,INLTEP=INLTE+1,NNLTE=56) 236 PARAMETER (LP1I=IMAX*LP1,LLP1I=IMAX*LLP1,LL3PI=IMAX*LL3P) 237 PARAMETER (NB1=NB-1) 238 PARAMETER (KO2=12) 239 PARAMETER (KO21=KO2+1,KO2M=KO2-1) 240 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 241 C IMAX = NO. POINTS SENT TO RADFS 242 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 243 C***NOTE: THE USER NORMALLY WILL MODIFY ONLY THE IMAX AND L PARAMETERS 244 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 245 C BANDTA FOR DEFINITION 246 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 247 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE 248 C BDCOMB FOR DEFINITION 249 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 250 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 251 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 252 C FROM THE ABOVE PARAMETERS. Page 4 Source Listing FST88 2025-03-12 18:22 FST88.F 253 C COMMON BLOCK BANDTA CONTAINS RANDOM BAND PARAMETERS FOR THE LW 254 C CALCULATIONS USING 10 CM-1 WIDE BANDS.THE 15 UM CO2 COMPLEX 255 C IS 2 BANDS,560-670 AND 670-800 CM-1. OZONE COEFFICIENTS ARE 256 C IN 3 BANDS,670-800 (14.1 UM),990-1070 AND 1070-1200 (9.6 UM). 257 C THE (NBLW) BANDS NOW INCLUDE: 258 C 56 BANDS, 10 CM-1 WIDE 0 - 560 CM-1 259 C 2 BANDS, 15 UM COMPLEX 560 - 670 CM-1 260 C 670 - 800 CM-1 261 C 3 "CONTINUUM" BANDS 800 - 900 CM-1 262 C 900 - 990 CM-1 263 C 1070 - 1200 CM-1 264 C 1 BAND FOR 9.6 UM BAND 990 - 1070 CM-1 265 C 100 BANDS, 10 CM-1 WIDE 1200 - 2200 CM-1 266 C 1 BAND FOR 4.3 UM SRC 2270 - 2380 CM-1 267 C THUS NBLW PRESENTLY EQUALS 163 268 C ALL BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 269 C 270 C ARNDM = RANDOM "A" PARAMETER FOR (NBLW) BANDS 271 C BRNDM = RANDOM "B" PARAMETER FOR (NBLW) BANDS 272 C BETAD = CONTINUUM COEFFICIENTS FOR (NBLW) BANDS 273 C AP,BP = CAPPHI COEFFICIENTS FOR (NBLW) BANDS 274 C ATP,BTP = CAPPSI COEFFICIENTS FOR (NBLW) BANDS 275 C BANDLO = LOWEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 276 C BANDHI = HIGHEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 277 C AO3RND = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 278 C BANDS 279 C BO3RND = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 280 C BANDS 281 C AB15 = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS 282 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 283 C DATA FOR ARNDM,BRNDM,AP,BP,ATP,BTP,AO3RND,BO3RND ARE OBTAINED BY 284 C USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS ARE FROM 285 C ROBERTS (1976). 286 COMMON / BANDTA / ARNDM(NBLW),BRNDM(NBLW),BETAD(NBLW),AP(NBLW), 287 1 BP(NBLW),ATP(NBLW),BTP(NBLW),BANDLO(NBLW), 288 2 BANDHI(NBLW),AO3RND(3),BO3RND(3),AB15(2) 289 C 290 C COMMON BLOCK BDWIDE CONTAINS RANDOM BAND PARAMETERS FOR SPECIFIC 291 C WIDE BANDS. AT PRESENT,THE INFORMATION CONSISTS OF 1) RANDOM 292 C MODEL PARAMETERS FOR THE 15 UM BAND,560-800 CM-1; 2) THE 293 C CONTINUUM COEFFICIENT FOR THE 800-990,1070-1200 CM-1 BAND 294 C SPECIFICALLY: 295 C AWIDE = RANDOM "A" PARAMETER FOR BAND 296 C BWIDE = RANDOM "B" PARAMETER FOR BAND 297 C BETAWD = CONTINUUM COEFFICIENTS FOR BAND 298 C APWD,BPWD = CAPPHI COEFFICIENTS FOR BAND 299 C ATPWD,BTPWD = CAPPSI COEFFICIENTS FOR BAND 300 C BDLOWD = LOWEST FREQUENCY IN EACH FREQ BAND 301 C BDHIWD = HIGHEST FREQUENCY IN EACH FREQ BAND 302 C AB15WD = THE PRODUCT ARNDM*BRNDM FOR THE ONE BAND 303 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 304 C BETINW = CONT.COEFFICIENT FOR A SPECIFIED WIDE 305 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 306 C SKO2D = 1./BETINW, USED IN SPA88 FOR CONT. COEFFS 307 C SKC1R = BETAWD/BETINW, USED FOR CONT. COEFF. FOR 308 C 15 UM BAND IN FST88 309 C SKO3R = RATIO OF CONT. COEFF. FOR 9.9 UM BAND TO Page 5 Source Listing FST88 2025-03-12 18:22 FST88.F 310 C BETINW, USED FOR 9.6 UM CONT COEFF IN FST88 311 C DATA FOR AWIDE,BWIDE,APWD,BPWD,ATPWD,BTPWD,AO3WD,BO3WD ARE 312 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 313 C ARE FROM ROBERTS (1976). 314 COMMON /BDWIDE/ AWIDE,BWIDE,BETAWD, 315 1 APWD,BPWD,ATPWD,BTPWD, 316 2 BDLOWD,BDHIWD,BETINW, 317 3 AB15WD,SKO2D,SKC1R,SKO3R 318 save /BDWIDE/ 319 C 320 C COMMON BLOCK BDCOMB CONTAINS RANDOM BAND PARAMETERS FOR THE LW 321 C CALCULATIONS USING COMBINED WIDE FREQUENCY BANDS BETWEEN 160 AND 322 C 1200 CM-1,AS WELL AS THE 2270-2380 BAND FOR SOURCE CALC. 323 C BANDS 1-8: COMBINED WIDE FREQUENCY BANDS FOR 160-560 CM-1 324 C BANDS 9-14: FREQUENCY BANDS,AS IN BANDTA (NARROW BANDS) 325 C FOR 560-1200 CM-1 326 C BAND 15: FREQUENCY BAND 2270-2380 CM-1,USED FOR SOURCE 327 C CALCULATION ONLY 328 C THUS NBLY PRESENTLY EQUALS 15 329 C 330 C BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 331 C ACOMB = RANDOM "A" PARAMETER FOR (NBLY) BANDS 332 C BCOMB = RANDOM "B" PARAMETER FOR (NBLY) BANDS 333 C BETACM = CONTINUUM COEFFICIENTS FOR (NBLY) BANDS 334 C APCM,BPCM = CAPPHI COEFFICIENTS FOR (NBLY) BANDS 335 C ATPCM,BTPCM = CAPPSI COEFFICIENTS FOR (NBLY) BANDS 336 C BDLOCM = LOWEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 337 C BDHICM = HIGHEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 338 C AO3CM = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 339 C BANDS 340 C BO3CM = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 341 C BANDS 342 C AB15CM = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS 343 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 344 C BETINC = CONT.COEFFICIENT FOR A SPECIFIED WIDE 345 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 346 C IBAND = INDEX NO OF THE 40 WIDE BANDS USED IN 347 C COMBINED WIDE BAND CALCULATIONS. IN OTHER 348 C WORDS,INDEX TELLING WHICH OF THE 40 WIDE 349 C BANDS BETWEEN 160-560 CM-1 ARE INCLUDED IN 350 C EACH OF THE FIRST 8 COMBINED WIDE BANDS 351 C DATA FOR ACOMB,BCOMB,APCM,BPCM,ATPCM,BTPCM,AO3CM,BO3CM ARE 352 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 353 C ARE FROM ROBERTS (1976). IBAND INDEX VALUES ARE OBTAINED BY 354 C EXPERIMENTATION. 355 COMMON / BDCOMB / IBAND(40),ACOMB(NBLY),BCOMB(NBLY), 356 1 BETACM(NBLY),APCM(NBLY),BPCM(NBLY),ATPCM(NBLY), 357 2 BTPCM(NBLY),BDLOCM(NBLY),BDHICM(NBLY),BETINC, 358 3 AO3CM(3),BO3CM(3),AB15CM(2) 359 save / BDCOMB / 360 C 361 C COMMON BLOCK TABCOM CONTAINS QUANTITIES PRECOMPUTED IN SUBROUTINE 362 C TABLE FOR USE IN THE LONGWAVE RADIATION PROGRAM: 363 C EM1 = E1 FUNCTION, EVALUATED OVER THE 0-560 AND 364 C 1200-2200 CM-1 INTERVALS 365 C EM1WDE = E1 FUNCTION, EVALUATED OVER THE 160-560 CM-1 366 C INTERVAL Page 6 Source Listing FST88 2025-03-12 18:22 FST88.F 367 C TABLE1 = E2 FUNCTION, EVALUATED OVER THE 0-560 AND 368 C 1200-2200 CM-1 INTERVALS 369 C TABLE2 = TEMPERATURE DERIVATIVE OF TABLE1 370 C TABLE3 = MASS DERIVATIVE OF TABLE1 371 C EM3 = E3 FUNCTION, EVALUATED OVER THE 0-560 AND 372 C 1200-2200 CM-1 INTERVALS 373 C SOURCE = PLANCK FUNCTION, EVALUATED AT SPECIFIED TEMPS. FOR 374 C BANDS USED IN CTS CALCULATIONS 375 C DSRCE = TEMPERATURE DERIVATIVE OF SOURCE 376 C IND = INDEX, WITH VALUE IND(I)=I. USED IN FST88 377 C INDX2 = INDEX VALUES USED IN OBTAINING "LOWER TRIANGLE" 378 C ELEMENTS OF AVEPHI,ETC.,IN FST88 379 C KMAXV = INDEX VALUES USED IN OBTAINING "UPPER TRIANGLE" 380 C ELEMENTS OF AVEPHI,ETC.,IN FST88 381 C KMAXVM = KMAXV(L),USED FOR DO LOOP INDICES 382 C 383 COMMON /TABCOM/ IND(IMAX),INDX2(LP1V),KMAXV(LP1), 384 1 KMAXVM 385 COMMON/TABCOM/EM1(28,180),EM1WDE(28,180),TABLE1(28,180), 386 1 TABLE2(28,180),TABLE3(28,180),EM3(28,180),SOURCE(28,NBLY), 387 2 DSRCE(28,NBLY) 388 save /TABCOM/ 389 C 390 DIMENSION QH2O(IDIM1:IDIM2,LP1),PRESS(IDIM1:IDIM2,LP1) 391 DIMENSION P(IDIM1:IDIM2,LP1),DELP(IDIM1:IDIM2,L), 392 & DELP2(IDIM1:IDIM2,L),TEMP(IDIM1:IDIM2,LP1) 393 DIMENSION T(IDIM1:IDIM2,LP1),CLDFAC(IDIM1:IDIM2,LP1,LP1), 394 & CAMT(IDIM1:IDIM2,LP1) 395 DIMENSION NCLDS(IDIM1:IDIM2),KTOP(IDIM1:IDIM2,LP1), 396 & KBTM(IDIM1:IDIM2,LP1) 397 DIMENSION CO21(IDIM1:IDIM2,LP1,LP1),CO2NBL(IDIM1:IDIM2,L) 398 DIMENSION CO2SP1(IDIM1:IDIM2,LP1),CO2SP2(IDIM1:IDIM2,LP1) 399 DIMENSION VAR1(IDIM1:IDIM2,L),VAR2(IDIM1:IDIM2,L), 400 & VAR3(IDIM1:IDIM2,L),VAR4(IDIM1:IDIM2,L) 401 DIMENSION CNTVAL(IDIM1:IDIM2,LP1) 402 DIMENSION HEATRA(IDIM1:IDIM2,L),GRNFLX(IDIM1:IDIM2), 403 & TOPFLX(IDIM1:IDIM2) 404 DIMENSION GXCTS(IDIM1:IDIM2),FLX1E1(IDIM1:IDIM2) 405 DIMENSION AVEPHI(IDIM1:IDIM2,LP1),EMISS(IDIM1:IDIM2,LP1), 406 & EMISSB(IDIM1:IDIM2,LP1) 407 C 408 DIMENSION TOTO3(IDIM1:IDIM2,LP1),TPHIO3(IDIM1:IDIM2,LP1), 409 & TOTPHI(IDIM1:IDIM2,LP1) 410 DIMENSION TOTVO2(IDIM1:IDIM2,LP1),EMX1(IDIM1:IDIM2), 411 & EMX2(IDIM1:IDIM2),EMPL(IDIM1:IDIM2,LLP1) 412 C 413 DIMENSION EXCTS(IDIM1:IDIM2,L),CTSO3(IDIM1:IDIM2,L), 414 & CTS(IDIM1:IDIM2,L),E1FLX(IDIM1:IDIM2,LP1) 415 DIMENSION CO2SP(IDIM1:IDIM2,LP1),TO3SPC(IDIM1:IDIM2,L), 416 & TO3SP(IDIM1:IDIM2,LP1) 417 DIMENSION OSS(IDIM1:IDIM2,LP1),CSS(IDIM1:IDIM2,LP1), 418 & SS1(IDIM1:IDIM2,LP1),SS2(IDIM1:IDIM2,LP1), 419 1 TC(IDIM1:IDIM2,LP1),DTC(IDIM1:IDIM2,LP1) 420 DIMENSION SORC(IDIM1:IDIM2,LP1,NBLY),CSOUR(IDIM1:IDIM2,LP1) 421 CCC 422 DIMENSION AVVO2(IDIM1:IDIM2,LP1),HEATEM(IDIM1:IDIM2,LP1), 423 1 OVER1D(IDIM1:IDIM2,LP1), Page 7 Source Listing FST88 2025-03-12 18:22 FST88.F 424 1 TO31D(IDIM1:IDIM2,LP1),CONT1D(IDIM1:IDIM2,LP1), 425 2 AVMO3(IDIM1:IDIM2,LP1),AVPHO3(IDIM1:IDIM2,LP1), 426 2 C(IDIM1:IDIM2,LLP1),C2(IDIM1:IDIM2,LLP1) 427 DIMENSION ITOP(IDIM1:IDIM2),IBOT(IDIM1:IDIM2), 428 & INDTC(IDIM1:IDIM2) 429 DIMENSION 430 4 DELPTC(IDIM1:IDIM2),PTOP(IDIM1:IDIM2),PBOT(IDIM1:IDIM2), 431 & FTOP(IDIM1:IDIM2), 432 5 FBOT(IDIM1:IDIM2) ,EMSPEC(IDIM1:IDIM2,2) 433 C---DIMENSION OF VARIABLES EQUIVALENCED TO THOSE IN VTEMP--- 434 DIMENSION VTMP3(IDIM1:IDIM2,LP1),DSORC(IDIM1:IDIM2,LP1) 435 DIMENSION ALP(IDIM1:IDIM2,LLP1),CSUB(IDIM1:IDIM2,LLP1), 436 & CSUB2(IDIM1:IDIM2,LLP1) 437 DIMENSION FAC1(IDIM1:IDIM2,LP1) 438 DIMENSION DELPR1(IDIM1:IDIM2,LP1),DELPR2(IDIM1:IDIM2,LP1) 439 DIMENSION EMISDG(IDIM1:IDIM2,LP1),CONTDG(IDIM1:IDIM2,LP1), 440 & TO3DG(IDIM1:IDIM2,LP1) 441 DIMENSION FLXNET(IDIM1:IDIM2,LP1) 442 DIMENSION IXO(IDIM1:IDIM2,LP1) 443 DIMENSION VSUM1(IDIM1:IDIM2,LP1) 444 DIMENSION FLXTHK(IDIM1:IDIM2,LP1) 445 DIMENSION Z1(IDIM1:IDIM2,LP1) 446 C---DIMENSION OF VARIABLES PASSED TO OTHER SUBROUTINES--- 447 C (AND NOT FOUND IN COMMON BLOCKS) 448 DIMENSION E1CTS1(IDIM1:IDIM2,LP1),E1CTS2(IDIM1:IDIM2,L) 449 DIMENSION E1CTW1(IDIM1:IDIM2,LP1),E1CTW2(IDIM1:IDIM2,L) 450 DIMENSION EMD(IDIM1:IDIM2,LLP1),TPL(IDIM1:IDIM2,LLP1) 451 C IT IS POSSIBLE TO EQUIVALENCE EMD,TPL TO THE ABOVE VARIABLES, 452 C AS THEY GET CALLED AT DIFFERENT TIMES 453 DIMENSION FXO(IDIM1:IDIM2,LP1),DT(IDIM1:IDIM2,LP1) 454 DIMENSION FXOE2(IDIM1:IDIM2,LP1),DTE2(IDIM1:IDIM2,LP1) 455 DIMENSION FXOSP(IDIM1:IDIM2,2),DTSP(IDIM1:IDIM2,2) 456 C 457 C DIMENSION OF LOCAL VARIABLES 458 DIMENSION RLOG(IDIM1:IDIM2,L),FLX(IDIM1:IDIM2,LP1) 459 DIMENSION TOTEVV(IDIM1:IDIM2,LP1),CNTTAU(IDIM1:IDIM2,LP1) 460 C 461 EQUIVALENCE (ALP,C,CSUB),(CSUB2,C2) 462 EQUIVALENCE (FAC1,DSORC,OVER1D,DELPR2,FLXNET) 463 EQUIVALENCE (DELPR1,HEATEM) 464 EQUIVALENCE (IXO,AVVO2,FLXTHK,TO3DG) 465 EQUIVALENCE (Z1,AVMO3,CONTDG) 466 EQUIVALENCE (EMISDG,VSUM1,AVPHO3) 467 EQUIVALENCE (EMD(IDIM1,1),E1CTS1(IDIM1,1)), 468 & (EMD(IDIM1,LP2),E1CTS2(IDIM1,1)) 469 EQUIVALENCE (TPL(IDIM1,1),E1CTW1(IDIM1,1)), 470 & (TPL(IDIM1,LP2),E1CTW2(IDIM1,1)) 471 c 472 C 473 C FIRST SECTION IS TABLE LOOKUP FOR SOURCE FUNCTION AND 474 C DERIVATIVE (B AND DB/DT).ALSO,THE NLTE CO2 SOURCE FUNCTION 475 C IS OBTAINED 476 C 477 C---IN CALCS. BELOW, DECREMENTING THE INDEX BY 9 478 C ACCOUNTS FOR THE TABLES BEGINNING AT T=100K. 479 C AT T=100K. 480 DO 101 K=1,LP1 Page 8 Source Listing FST88 2025-03-12 18:22 FST88.F 481 DO 101 I=MYIS,MYIE 482 C---TEMP. INDICES FOR E1,SOURCE 483 VTMP3(I,K)=AINT(TEMP(I,K)*HP1) 484 FXO(I,K)=VTMP3(I,K)-9. 485 DT(I,K)=TEMP(I,K)-TEN*VTMP3(I,K) 486 C---INTEGER INDEX FOR SOURCE (USED IMMEDIATELY) 487 C wne IXO(I,K)=FXO(I,K) 488 IXO(I,K)=max(FXO(I,K), 1.0) 489 101 CONTINUE 490 DO 103 k=1,L 491 DO 103 I=MYIS,MYIE 492 C---TEMP. INDICES FOR E2 (KP=1 LAYER NOT USED IN FLUX CALCULATIONS) 493 VTMP3(I,K)=AINT(T(I,K+1)*HP1) 494 FXOE2(I,K)=VTMP3(I,K)-9. 495 DTE2(I,K)=T(I,K+1)-TEN*VTMP3(I,K) 496 103 CONTINUE 497 C---SPECIAL CASE TO HANDLE KP=LP1 LAYER AND SPECIAL E2 CALCS. 498 DO 105 I=MYIS,MYIE 499 FXOE2(I,LP1)=FXO(I,L) 500 DTE2(I,LP1)=DT(I,L) 501 FXOSP(I,1)=FXOE2(I,LM1) 502 FXOSP(I,2)=FXO(I,LM1) 503 DTSP(I,1)=DTE2(I,LM1) 504 DTSP(I,2)=DT(I,LM1) 505 105 CONTINUE 506 C 507 C---SOURCE FUNCTION FOR COMBINED BAND 1 508 DO 4114 I=MYIS,MYIE 509 DO 4114 K=1,LP1 510 VTMP3(I,K)=SOURCE(IXO(I,K),1) 511 DSORC(I,K)=DSRCE(IXO(I,K),1) 512 4114 CONTINUE 513 DO 4112 K=1,LP1 514 DO 4112 I=MYIS,MYIE 515 SORC(I,K,1)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 516 4112 CONTINUE 517 C---SOURCE FUNCTION FOR COMBINED BAND 2 518 DO 4214 I=MYIS,MYIE 519 DO 4214 K=1,LP1 520 VTMP3(I,K)=SOURCE(IXO(I,K),2) 521 DSORC(I,K)=DSRCE(IXO(I,K),2) 522 4214 CONTINUE 523 DO 4212 K=1,LP1 524 DO 4212 I=MYIS,MYIE 525 SORC(I,K,2)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 526 4212 CONTINUE 527 C---SOURCE FUNCTION FOR COMBINED BAND 3 528 DO 4314 I=MYIS,MYIE 529 DO 4314 K=1,LP1 530 VTMP3(I,K)=SOURCE(IXO(I,K),3) 531 DSORC(I,K)=DSRCE(IXO(I,K),3) 532 4314 CONTINUE 533 DO 4312 K=1,LP1 534 DO 4312 I=MYIS,MYIE 535 SORC(I,K,3)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 536 4312 CONTINUE 537 C---SOURCE FUNCTION FOR COMBINED BAND 4 Page 9 Source Listing FST88 2025-03-12 18:22 FST88.F 538 DO 4414 I=MYIS,MYIE 539 DO 4414 K=1,LP1 540 VTMP3(I,K)=SOURCE(IXO(I,K),4) 541 DSORC(I,K)=DSRCE(IXO(I,K),4) 542 4414 CONTINUE 543 DO 4412 K=1,LP1 544 DO 4412 I=MYIS,MYIE 545 SORC(I,K,4)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 546 4412 CONTINUE 547 C---SOURCE FUNCTION FOR COMBINED BAND 5 548 DO 4514 I=MYIS,MYIE 549 DO 4514 K=1,LP1 550 VTMP3(I,K)=SOURCE(IXO(I,K),5) 551 DSORC(I,K)=DSRCE(IXO(I,K),5) 552 4514 CONTINUE 553 DO 4512 K=1,LP1 554 DO 4512 I=MYIS,MYIE 555 SORC(I,K,5)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 556 4512 CONTINUE 557 C---SOURCE FUNCTION FOR COMBINED BAND 6 558 DO 4614 I=MYIS,MYIE 559 DO 4614 K=1,LP1 560 VTMP3(I,K)=SOURCE(IXO(I,K),6) 561 DSORC(I,K)=DSRCE(IXO(I,K),6) 562 4614 CONTINUE 563 DO 4612 K=1,LP1 564 DO 4612 I=MYIS,MYIE 565 SORC(I,K,6)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 566 4612 CONTINUE 567 C---SOURCE FUNCTION FOR COMBINED BAND 7 568 DO 4714 I=MYIS,MYIE 569 DO 4714 K=1,LP1 570 VTMP3(I,K)=SOURCE(IXO(I,K),7) 571 DSORC(I,K)=DSRCE(IXO(I,K),7) 572 4714 CONTINUE 573 DO 4712 K=1,LP1 574 DO 4712 I=MYIS,MYIE 575 SORC(I,K,7)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 576 4712 CONTINUE 577 C---SOURCE FUNCTION FOR COMBINED BAND 8 578 DO 4814 I=MYIS,MYIE 579 DO 4814 K=1,LP1 580 VTMP3(I,K)=SOURCE(IXO(I,K),8) 581 DSORC(I,K)=DSRCE(IXO(I,K),8) 582 4814 CONTINUE 583 DO 4812 K=1,LP1 584 DO 4812 I=MYIS,MYIE 585 SORC(I,K,8)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 586 4812 CONTINUE 587 C---SOURCE FUNCTION FOR BAND 9 (560-670 CM-1) 588 DO 4914 I=MYIS,MYIE 589 DO 4914 K=1,LP1 590 VTMP3(I,K)=SOURCE(IXO(I,K),9) 591 DSORC(I,K)=DSRCE(IXO(I,K),9) 592 4914 CONTINUE 593 DO 4912 K=1,LP1 594 DO 4912 I=MYIS,MYIE Page 10 Source Listing FST88 2025-03-12 18:22 FST88.F 595 SORC(I,K,9)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 596 4912 CONTINUE 597 C---SOURCE FUNCTION FOR BAND 10 (670-800 CM-1) 598 DO 5014 I=MYIS,MYIE 599 DO 5014 K=1,LP1 600 VTMP3(I,K)=SOURCE(IXO(I,K),10) 601 DSORC(I,K)=DSRCE(IXO(I,K),10) 602 5014 CONTINUE 603 DO 5012 K=1,LP1 604 DO 5012 I=MYIS,MYIE 605 SORC(I,K,10)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 606 5012 CONTINUE 607 C---SOURCE FUNCTION FOR BAND 11 (800-900 CM-1) 608 DO 5114 I=MYIS,MYIE 609 DO 5114 K=1,LP1 610 VTMP3(I,K)=SOURCE(IXO(I,K),11) 611 DSORC(I,K)=DSRCE(IXO(I,K),11) 612 5114 CONTINUE 613 DO 5112 K=1,LP1 614 DO 5112 I=MYIS,MYIE 615 SORC(I,K,11)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 616 5112 CONTINUE 617 C---SOURCE FUNCTION FOR BAND 12 (900-990 CM-1) 618 DO 5214 I=MYIS,MYIE 619 DO 5214 K=1,LP1 620 VTMP3(I,K)=SOURCE(IXO(I,K),12) 621 DSORC(I,K)=DSRCE(IXO(I,K),12) 622 5214 CONTINUE 623 DO 5212 K=1,LP1 624 DO 5212 I=MYIS,MYIE 625 SORC(I,K,12)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 626 5212 CONTINUE 627 C---SOURCE FUNCTION FOR BAND 13 (990-1070 CM-1) 628 DO 5314 I=MYIS,MYIE 629 DO 5314 K=1,LP1 630 VTMP3(I,K)=SOURCE(IXO(I,K),13) 631 DSORC(I,K)=DSRCE(IXO(I,K),13) 632 5314 CONTINUE 633 DO 5312 K=1,LP1 634 DO 5312 I=MYIS,MYIE 635 SORC(I,K,13)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 636 5312 CONTINUE 637 C---SOURCE FUNCTION FOR BAND 14 (1070-1200 CM-1) 638 DO 5414 I=MYIS,MYIE 639 DO 5414 K=1,LP1 640 VTMP3(I,K)=SOURCE(IXO(I,K),14) 641 DSORC(I,K)=DSRCE(IXO(I,K),14) 642 5414 CONTINUE 643 DO 5412 K=1,LP1 644 DO 5412 I=MYIS,MYIE 645 SORC(I,K,14)=VTMP3(I,K)+DT(I,K)*DSORC(I,K) 646 5412 CONTINUE 647 C 648 C THE FOLLOWING SUBROUTINE OBTAINS NLTE SOURCE FUNCTION FOR CO2 649 C 650 C 651 C CALL NLTE Page 11 Source Listing FST88 2025-03-12 18:22 FST88.F 652 C 653 C 654 C---OBTAIN SPECIAL SOURCE FUNCTIONS FOR THE 15 UM BAND (CSOUR) 655 C AND THE WINDOW REGION (SS1) 656 DO 131 K=1,LP1 657 DO 131 I=MYIS,MYIE 658 SS1(I,K)=SORC(I,K,11)+SORC(I,K,12)+SORC(I,K,14) 659 131 CONTINUE 660 DO 143 K=1,LP1 661 DO 143 I=MYIS,MYIE 662 CSOUR(I,K)=SORC(I,K,9)+SORC(I,K,10) 663 143 CONTINUE 664 C 665 C---COMPUTE TEMP**4 (TC) AND VERTICAL TEMPERATURE DIFFERENCES 666 C (OSS,CSS,SS2,DTC). ALL THESE WILL BE USED LATER IN FLUX COMPUTA- 667 C TIONS. 668 C 669 DO 901 K=1,LP1 670 DO 901 I=MYIS,MYIE 671 TC(I,K)=(TEMP(I,K)*TEMP(I,K))**2 672 c if(mype.eq.13.and.i.eq.40) then 673 c print*,'i,k,temp(i,k)=',i,k,temp(i,k) 674 c endif 675 901 CONTINUE 676 DO 903 K=1,L 677 DO 903 I=MYIS,MYIE 678 OSS(I,K+1)=SORC(I,K+1,13)-SORC(I,K,13) 679 CSS(I,K+1)=CSOUR(I,K+1)-CSOUR(I,K) 680 DTC(I,K+1)=TC(I,K+1)-TC(I,K) 681 SS2(I,K+1)=SS1(I,K+1)-SS1(I,K) 682 903 CONTINUE 683 C 684 C 685 C---THE FOLLOWIMG IS A DRASTIC REWRITE OF THE RADIATION CODE TO 686 C (LARGELY) ELIMINATE THREE-DIMENSIONAL ARRAYS. THE CODE WORKS 687 C ON THE FOLLOWING PRINCIPLES: 688 C 689 C LET K = FIXED FLUX LEVEL, KP = VARYING FLUX LEVEL 690 C THEN FLUX(K)=SUM OVER KP : (DELTAB(KP)*TAU(KP,K)) 691 C OVER ALL KP'S, FROM 1 TO LP1. 692 C 693 C WE CAN BREAK DOWN THE CALCULATIONS FOR ALL K'S AS FOLLOWS: 694 C 695 C FOR ALL K'S K=1 TO LP1: 696 C FLUX(K)=SUM OVER KP : (DELTAB(KP)*TAU(KP,K)) (1) 697 C OVER ALL KP'S, FROM K+1 TO LP1 698 C AND 699 C FOR KP FROM K+1 TO LP1: 700 C FLUX(KP) = DELTAB(K)*TAU(K,KP) (2) 701 C 702 C NOW IF TAU(K,KP)=TAU(KP,K) (SYMMETRICAL ARRAYS) 703 C WE CAN COMPUTE A 1-DIMENSIONAL ARRAY TAU1D(KP) FROM 704 C K+1 TO LP1, EACH TIME K IS INCREMENTED. 705 C EQUATIONS (1) AND (2) THEN BECOME: 706 C 707 C TAU1D(KP) = (VALUES FOR TAU(KP,K) AT THE PARTICULAR K) 708 C FLUX(K) = SUM OVER KP : (DELTAB(KP)*TAU1D(KP)) (3) Page 12 Source Listing FST88 2025-03-12 18:22 FST88.F 709 C FLUX(KP) = DELTAB(K)*TAU1D(KP) (4) 710 C 711 C THE TERMS FOR TAU (K,K) AND OTHER SPECIAL TERMS (FOR 712 C NEARBY LAYERS) MUST, OF COURSE, BE HANDLED SEPARATELY, AND 713 C WITH CARE. 714 C 715 C COMPUTE "UPPER TRIANGLE" TRANSMISSION FUNCTIONS FOR 716 C THE 9.6 UM BAND (TO3SP) AND THE 15 UM BAND (OVER1D). ALSO, 717 C THE 718 C STAGE 1...COMPUTE O3 ,OVER TRANSMISSION FCTNS AND AVEPHI 719 C---DO K=1 CALCULATION (FROM FLUX LAYER KK TO THE TOP) SEPARATELY 720 C AS VECTORIZATION IS IMPROVED,AND OZONE CTS TRANSMISSIVITY 721 C MAY BE EXTRACTED HERE. 722 DO 3021 K=1,L 723 DO 3021 I=MYIS,MYIE 724 AVEPHI(I,K)=TOTPHI(I,K+1) 725 3021 CONTINUE 726 C---IN ORDER TO PROPERLY EVALUATE EMISS INTEGRATED OVER THE (LP1) 727 C LAYER, A SPECIAL EVALUATION OF EMISS IS DONE. THIS REQUIRES 728 C A SPECIAL COMPUTATION OF AVEPHI, AND IT IS STORED IN THE 729 C (OTHERWISE VACANT) LP1'TH POSITION 730 C 731 DO 803 I=MYIS,MYIE 732 AVEPHI(I,LP1)=AVEPHI(I,LM1)+EMX1(I) 733 803 CONTINUE 734 C COMPUTE FLUXES FOR K=1 735 CALL E1E290(E1CTS1,E1CTS2,E1FLX,E1CTW1,E1CTW2,EMISS, 736 1 FXO,DT,FXOE2,DTE2,AVEPHI,TEMP,T) 737 DO 302 K=1,L 738 DO 302 I=MYIS,MYIE 739 FAC1(I,K)=BO3RND(2)*TPHIO3(I,K+1)/TOTO3(I,K+1) 740 TO3SPC(I,K)=HAF*(FAC1(I,K)* 741 1 (SQRT(ONE+(FOUR*AO3RND(2)*TOTO3(I,K+1))/FAC1(I,K))-ONE)) 742 C FOR K=1, TO3SP IS USED INSTEAD OF TO31D (THEY ARE EQUAL IN THIS 743 C CASE); TO3SP IS PASSED TO SPA90, WHILE TO31D IS A WORK-ARRAY. 744 TO3SP(I,K)=EXP(HM1EZ*(TO3SPC(I,K)+SKO3R*TOTVO2(I,K+1))) 745 OVER1D(I,K)=EXP(HM1EZ*(SQRT(AB15WD*TOTPHI(I,K+1))+ 746 1 SKC1R*TOTVO2(I,K+1))) 747 C---BECAUSE ALL CONTINUUM TRANSMISSIVITIES ARE OBTAINED FROM THE 748 C 2-D QUANTITY CNTTAU (AND ITS RECIPROCAL TOTEVV) WE STORE BOTH 749 C OF THESE HERE. FOR K=1, CONT1D EQUALS CNTTAU 750 CNTTAU(I,K)=EXP(HM1EZ*TOTVO2(I,K+1)) 751 TOTEVV(I,K)=1./CNTTAU(I,K) 752 302 CONTINUE 753 DO 3022 K=1,L 754 DO 3022 I=MYIS,MYIE 755 CO2SP(I,K+1)=OVER1D(I,K)*CO21(I,1,K+1) 756 3022 CONTINUE 757 DO 3023 K=1,L 758 DO 3023 I=MYIS,MYIE 759 CO21(I,K+1,1)=CO21(I,K+1,1)*OVER1D(I,K) 760 3023 CONTINUE 761 C---RLOG IS THE NBL AMOUNT FOR THE 15 UM BAND CALCULATION 762 DO 1808 I=MYIS,MYIE 763 RLOG(I,1)=OVER1D(I,1)*CO2NBL(I,1) 764 1808 CONTINUE 765 C---THE TERMS WHEN KP=1 FOR ALL K ARE THE PHOTON EXCHANGE WITH Page 13 Source Listing FST88 2025-03-12 18:22 FST88.F 766 C THE TOP OF THE ATMOSPHERE, AND ARE OBTAINED DIFFERENTLY THAN 767 C THE OTHER CALCULATIONS 768 DO 305 K=2,LP1 769 DO 305 I=MYIS,MYIE 770 FLX(I,K)= (TC(I,1)*E1FLX(I,K) 771 1 +SS1(I,1)*CNTTAU(I,K-1) 772 2 +SORC(I,1,13)*TO3SP(I,K-1) 773 3 +CSOUR(I,1)*CO2SP(I,K)) 774 4 *CLDFAC(I,1,K) 775 305 CONTINUE 776 DO 307 I=MYIS,MYIE 777 FLX(I,1)= TC(I,1)*E1FLX(I,1)+SS1(I,1)+SORC(I,1,13) 778 1 +CSOUR(I,1) 779 307 CONTINUE 780 C---THE KP TERMS FOR K=1... 781 DO 303 KP=2,LP1 782 DO 303 I=MYIS,MYIE 783 FLX(I,1)=FLX(I,1)+(OSS(I,KP)*TO3SP(I,KP-1) 784 1 +SS2(I,KP)*CNTTAU(I,KP-1) 785 2 +CSS(I,KP)*CO21(I,KP,1) 786 3 +DTC(I,KP)*EMISS(I,KP-1))*CLDFAC(I,KP,1) 787 303 CONTINUE 788 C SUBROUTINE SPA88 IS CALLED TO OBTAIN EXACT CTS FOR WATER 789 C CO2 AND O3, AND APPROXIMATE CTS CO2 AND O3 CALCULATIONS. 790 C 791 CALL SPA88(EXCTS,CTSO3,GXCTS,SORC,CSOUR, 792 1 CLDFAC,TEMP,PRESS,VAR1,VAR2, 793 2 P,DELP,DELP2,TOTVO2,TO3SP,TO3SPC, 794 3 CO2SP1,CO2SP2,CO2SP) 795 C 796 C THIS SECTION COMPUTES THE EMISSIVITY CTS HEATING RATES FOR 2 797 C EMISSIVITY BANDS: THE 0-160,1200-2200 CM-1 BAND AND THE 800- 798 C 990,1070-1200 CM-1 BAND. THE REMAINING CTS COMTRIBUTIONS ARE 799 C CONTAINED IN CTSO3, COMPUTED IN SPA88. 800 C 801 DO 998 I=MYIS,MYIE 802 VTMP3(I,1)=1. 803 998 CONTINUE 804 DO 999 K=1,L 805 DO 999 I=MYIS,MYIE 806 VTMP3(I,K+1)=CNTTAU(I,K)*CLDFAC(I,K+1,1) 807 999 CONTINUE 808 DO 1001 K=1,L 809 DO 1001 I=MYIS,MYIE 810 CTS(I,K)=RADCON*DELP(I,K)*(TC(I,K)* 811 1 (E1CTW2(I,K)*CLDFAC(I,K+1,1)-E1CTW1(I,K)*CLDFAC(I,K,1)) + 812 2 SS1(I,K)*(VTMP3(I,K+1)-VTMP3(I,K))) 813 1001 CONTINUE 814 C 815 DO 1011 K=1,L 816 DO 1011 I=MYIS,MYIE 817 VTMP3(I,K)=TC(I,K)*(CLDFAC(I,K,1)*(E1CTS1(I,K)-E1CTW1(I,K)) - 818 1 CLDFAC(I,K+1,1)*(E1CTS2(I,K)-E1CTW2(I,K))) 819 1011 CONTINUE 820 DO 1012 I=MYIS,MYIE 821 FLX1E1(I)=TC(I,LP1)*CLDFAC(I,LP1,1)* 822 1 (E1CTS1(I,LP1)-E1CTW1(I,LP1)) Page 14 Source Listing FST88 2025-03-12 18:22 FST88.F 823 c if(mype.eq.13.and.i.eq.40) then 824 c print*,'i,lp1=',i,lp1 825 c print*,'tc(i,lp1)=',tc(i,lp1) 826 c print*,'cldfac(i,lp1,1)=',cldfac(i,lp1,1) 827 c print*,'e1cts1(i,lp1)=',e1cts1(i,lp1) 828 c print*,'e1ctw1(i,lp1)=',e1ctw1(i,lp1) 829 c endif 830 1012 CONTINUE 831 DO 1014 K=1,L 832 DO 1013 I=MYIS,MYIE 833 FLX1E1(I)=FLX1E1(I)+VTMP3(I,K) 834 1013 CONTINUE 835 1014 CONTINUE 836 C 837 C---NOW REPEAT FLUX CALCULATIONS FOR THE K=2..LM1 CASES. 838 C CALCULATIONS FOR FLUX LEVEL L AND LP1 ARE DONE SEPARATELY, AS ALL 839 C EMISSIVITY AND CO2 CALCULATIONS ARE SPECIAL CASES OR NEARBY LAYERS. 840 C 841 DO 321 K=2,LM1 842 KLEN=K 843 C 844 DO 3218 KK=1,LP1-K 845 DO 3218 I=MYIS,MYIE 846 AVEPHI(I,KK+K-1)=TOTPHI(I,KK+K)-TOTPHI(I,K) 847 3218 CONTINUE 848 DO 1803 I=MYIS,MYIE 849 AVEPHI(I,LP1)=AVEPHI(I,LM1)+EMX1(I) 850 1803 CONTINUE 851 C---COMPUTE EMISSIVITY FLUXES (E2) FOR THIS CASE. NOTE THAT 852 C WE HAVE OMITTED THE NEARBY LATER CASE (EMISS(I,K,K)) AS WELL 853 C AS ALL CASES WITH K=L OR LP1. BUT THESE CASES HAVE ALWAYS 854 C BEEN HANDLED AS SPECIAL CASES, SO WE MAY AS WELL COMPUTE 855 C THEIR FLUXES SEPARASTELY. 856 C 857 CALL E290(EMISSB,EMISS,AVEPHI,KLEN,FXOE2,DTE2) 858 DO 322 KK=1,LP1-K 859 DO 322 I=MYIS,MYIE 860 AVMO3(I,KK+K-1)=TOTO3(I,KK+K)-TOTO3(I,K) 861 AVPHO3(I,KK+K-1)=TPHIO3(I,KK+K)-TPHIO3(I,K) 862 AVVO2(I,KK+K-1)=TOTVO2(I,KK+K)-TOTVO2(I,K) 863 CONT1D(I,KK+K-1)=CNTTAU(I,KK+K-1)*TOTEVV(I,K-1) 864 322 CONTINUE 865 C 866 DO 3221 KK=1,LP1-K 867 DO 3221 I=MYIS,MYIE 868 FAC1(I,K+KK-1)=BO3RND(2)*AVPHO3(I,K+KK-1)/AVMO3(I,K+KK-1) 869 VTMP3(I,K+KK-1)=HAF*(FAC1(I,K+KK-1)* 870 1 (SQRT(ONE+(FOUR*AO3RND(2)*AVMO3(I,K+KK-1))/ 871 2 FAC1(I,K+KK-1))-ONE)) 872 TO31D(I,K+KK-1)=EXP(HM1EZ*(VTMP3(I,K+KK-1) 873 1 +SKO3R*AVVO2(I,K+KK-1))) 874 OVER1D(I,K+KK-1)=EXP(HM1EZ*(SQRT(AB15WD*AVEPHI(I,K+KK-1))+ 875 1 SKC1R*AVVO2(I,K+KK-1))) 876 CO21(I,K+KK,K)=OVER1D(I,K+KK-1)*CO21(I,K+KK,K) 877 3221 CONTINUE 878 DO 3223 KP=K+1,LP1 879 DO 3223 I=MYIS,MYIE Page 15 Source Listing FST88 2025-03-12 18:22 FST88.F 880 CO21(I,K,KP)=OVER1D(I,KP-1)*CO21(I,K,KP) 881 3223 CONTINUE 882 C---RLOG IS THE NBL AMOUNT FOR THE 15 UM BAND CALCULATION 883 DO 1804 I=MYIS,MYIE 884 RLOG(I,K)=OVER1D(I,K)*CO2NBL(I,K) 885 1804 CONTINUE 886 C---THE KP TERMS FOR ARBIRRARY K.. 887 DO 3423 KP=K+1,LP1 888 DO 3423 I=MYIS,MYIE 889 FLX(I,K)=FLX(I,K)+(OSS(I,KP)*TO31D(I,KP-1) 890 1 +SS2(I,KP)*CONT1D(I,KP-1) 891 2 +CSS(I,KP)*CO21(I,KP,K) 892 3 +DTC(I,KP)*EMISS(I,KP-1))*CLDFAC(I,KP,K) 893 3423 CONTINUE 894 DO 3425 KP=K+1,LP1 895 DO 3425 I=MYIS,MYIE 896 FLX(I,KP)=FLX(I,KP)+(OSS(I,K)*TO31D(I,KP-1) 897 1 +SS2(I,K)*CONT1D(I,KP-1) 898 2 +CSS(I,K)*CO21(I,K,KP) 899 3 +DTC(I,K)*EMISSB(I,KP-1))*CLDFAC(I,K,KP) 900 3425 CONTINUE 901 321 CONTINUE 902 C 903 C NOW DO K=L CASE. SINCE THE KP LOOP IS LENGTH 1, MANY SIMPLIFI- 904 C CATIONS OCCUR. ALSO, THE CO2 QUANTITIES (AS WELL AS THE EMISS 905 C QUANTITIES) ARE COMPUTED IN THE NBL SEDCTION; THEREFORE, WE WANT 906 C ONLY OVER,TO3 AND CONT1D (OVER(I,L),TO31D(I,L) AND CONT1D(I,L) 907 C ACCORDING TO THE NOTATION. THUS NO CALL IS MADE TO THE E290 908 C SUBROUTINE. 909 C THE THIRD SECTION CALCULATES BOUNDARY LAYER AND NEARBY LAYER 910 C CORRECTIONS TO THE TRANSMISSION FUNCTIONS OBTAINED ABOVE. METHODS 911 C ARE GIVEN IN REF. (4). 912 C THE FOLLOWING RATIOS ARE USED IN VARIOUS NBL CALCULATIONS: 913 C 914 C THE REMAINING CALCULATIONS ARE FOR : 915 C 1) THE (K,K) TERMS, K=2,LM1; 916 C 2) THE (L,L) TERM 917 C 3) THE (L,LP1) TERM 918 C 4) THE (LP1,L) TERM 919 C 5) THE (LP1,LP1) TERM. 920 C EACH IS UNIQUELY HANDLED; DIFFERENT FLUX TERMS ARE COMPUTED 921 C DIFFERENTLY 922 C 923 C 924 C FOURTH SECTION OBTAINS WATER TRANSMISSION FUNCTIONS 925 C USED IN Q(APPROX) CALCULATIONS AND ALSO MAKES NBL CORRECTIONS: 926 C 1) EMISS (I,J) IS THE TRANSMISSION FUNCTION MATRIX OBTAINED 927 C BY CALLING SUBROUTINE E1E288; 928 C 2) "NEARBY LAYER" CORRECTIONS (EMISS(I,I)) ARE OBTAINED 929 C USING SUBROUTINE E3V88; 930 C 3) SPECIAL VALUES AT THE SURFACE (EMISS(L,LP1),EMISS(LP1,L), 931 C EMISS(LP1,LP1)) ARE CALCULATED. 932 C 933 C 934 C OBTAIN ARGUMENTS FOR E1E288 AND E3V88: 935 C 936 DO 821 I=MYIS,MYIE Page 16 Source Listing FST88 2025-03-12 18:22 FST88.F 937 TPL(I,1)=TEMP(I,L) 938 TPL(I,LP1)=HAF*(T(I,LP1)+TEMP(I,L)) 939 TPL(I,LLP1)=HAF*(T(I,L)+TEMP(I,L)) 940 821 CONTINUE 941 DO 823 K=2,L 942 DO 823 I=MYIS,MYIE 943 TPL(I,K)=T(I,K) 944 TPL(I,K+L)=T(I,K) 945 823 CONTINUE 946 C 947 C---E2 FUNCTIONS ARE REQUIRED IN THE NBL CALCULATIONS FOR 2 CASES, 948 C DENOTED (IN OLD CODE) AS (L,LP1) AND (LP1,LP1) 949 DO 833 I=MYIS,MYIE 950 AVEPHI(I,1)=VAR2(I,L) 951 AVEPHI(I,2)=VAR2(I,L)+EMPL(I,L) 952 833 CONTINUE 953 CALL E2SPEC(EMISS,AVEPHI,FXOSP,DTSP) 954 C 955 C CALL E3V88 FOR NBL H2O TRANSMISSIVITIES 956 CALL E3V88(EMD,TPL,EMPL) 957 C 958 C COMPUTE NEARBY LAYER AND SPECIAL-CASE TRANSMISSIVITIES FOR EMISS 959 C USING METHODS FOR H2O GIVEN IN REF. (4) 960 DO 851 K=2,L 961 DO 851 I=MYIS,MYIE 962 EMISDG(I,K)=EMD(I,K+L)+EMD(I,K) 963 851 CONTINUE 964 C 965 C NOTE THAT EMX1/2 (PRESSURE SCALED PATHS) ARE NOW COMPUTED IN 966 C LWR88 967 DO 861 I=MYIS,MYIE 968 EMSPEC(I,1)=(EMD(I,1)*EMPL(I,1)-EMD(I,LP1)*EMPL(I,LP1))/ 969 1 EMX1(I) + QUARTR*(EMISS(I,1)+EMISS(I,2)) 970 EMISDG(I,LP1)=TWO*EMD(I,LP1) 971 EMSPEC(I,2)=TWO*(EMD(I,1)*EMPL(I,1)-EMD(I,LLP1)*EMPL(I,LLP1))/ 972 * EMX2(I) 973 861 CONTINUE 974 DO 331 I=MYIS,MYIE 975 FAC1(I,L)=BO3RND(2)*VAR4(I,L)/VAR3(I,L) 976 VTMP3(I,L)=HAF*(FAC1(I,L)* 977 1 (SQRT(ONE+(FOUR*AO3RND(2)*VAR3(I,L))/FAC1(I,L))-ONE)) 978 TO31D(I,L)=EXP(HM1EZ*(VTMP3(I,L)+SKO3R*CNTVAL(I,L))) 979 OVER1D(I,L)=EXP(HM1EZ*(SQRT(AB15WD*VAR2(I,L))+ 980 1 SKC1R*CNTVAL(I,L))) 981 CONT1D(I,L)=CNTTAU(I,L)*TOTEVV(I,LM1) 982 RLOG(I,L)=OVER1D(I,L)*CO2NBL(I,L) 983 331 CONTINUE 984 DO 618 K=1,L 985 DO 618 I=MYIS,MYIE 986 RLOG(I,K)=LOG(RLOG(I,K)) 987 618 CONTINUE 988 DO 601 K=1,LM1 989 DO 601 I=MYIS,MYIE 990 DELPR1(I,K+1)=DELP(I,K+1)*(PRESS(I,K+1)-P(I,K+1)) 991 ALP(I,LP1+K-1)=-SQRT(DELPR1(I,K+1))*RLOG(I,K+1) 992 601 CONTINUE 993 DO 603 K=1,L Page 17 Source Listing FST88 2025-03-12 18:22 FST88.F 994 DO 603 I=MYIS,MYIE 995 DELPR2(I,K+1)=DELP(I,K)*(P(I,K+1)-PRESS(I,K)) 996 ALP(I,K)=-SQRT(DELPR2(I,K+1))*RLOG(I,K) 997 603 CONTINUE 998 DO 625 I=MYIS,MYIE 999 ALP(I,LL)=-RLOG(I,L) 1000 ALP(I,LLP1)=-RLOG(I,L)*SQRT(DELP(I,L)*(P(I,LP1)-PRESS(I,LM1))) 1001 625 CONTINUE 1002 C THE FIRST COMPUTATION IS FOR THE 15 UM BAND,WITH THE 1003 C FOR THE COMBINED H2O AND CO2 TRANSMISSION FUNCTION. 1004 C 1005 C PERFORM NBL COMPUTATIONS FOR THE 15 UM BAND 1006 C***THE STATEMENT FUNCTION SF IN PREV. VERSIONS IS NOW EXPLICITLY 1007 C EVALUATED. 1008 DO 631 K=1,LLP1 1009 DO 631 I=MYIS,MYIE 1010 C(I,K)=ALP(I,K)*(HMP66667+ALP(I,K)*(QUARTR+ALP(I,K)*HM6666M2)) 1011 631 CONTINUE 1012 DO 641 I=MYIS,MYIE 1013 CO21(I,LP1,LP1)=ONE+C(I,L) 1014 CO21(I,LP1,L)=ONE+(DELP2(I,L)*C(I,LL)-(PRESS(I,L)-P(I,L))* 1015 1 C(I,LLM1))/(P(I,LP1)-PRESS(I,L)) 1016 CO21(I,L,LP1)=ONE+((P(I,LP1)-PRESS(I,LM1))*C(I,LLP1)- 1017 1 (P(I,LP1)-PRESS(I,L))*C(I,L))/(PRESS(I,L)-PRESS(I,LM1)) 1018 641 CONTINUE 1019 DO 643 K=2,L 1020 DO 643 I=MYIS,MYIE 1021 CO21(I,K,K)=ONE+HAF*(C(I,LM1+K)+C(I,K-1)) 1022 643 CONTINUE 1023 C 1024 C COMPUTE NEARBY-LAYER TRANSMISSIVITIES FOR THE O3 BAND AND FOR THE 1025 C ONE-BAND CONTINUUM BAND (TO3 AND EMISS2). THE SF2 FUNCTION IS 1026 C USED. THE METHOD IS THE SAME AS DESCRIBED FOR CO2 IN REF (4). 1027 DO 651 K=1,LM1 1028 DO 651 I=MYIS,MYIE 1029 CSUB(I,K+1)=CNTVAL(I,K+1)*DELPR1(I,K+1) 1030 CSUB(I,LP1+K-1)=CNTVAL(I,K)*DELPR2(I,K+1) 1031 651 CONTINUE 1032 C---THE SF2 FUNCTION IN PREV. VERSIONS IS NOW EXPLICITLY EVALUATED 1033 DO 655 K=1,LLM2 1034 DO 655 I=MYIS,MYIE 1035 CSUB2(I,K+1)=SKO3R*CSUB(I,K+1) 1036 C(I,K+1)=CSUB(I,K+1)*(HMP5+CSUB(I,K+1)* 1037 1 (HP166666-CSUB(I,K+1)*H41666M2)) 1038 C2(I,K+1)=CSUB2(I,K+1)*(HMP5+CSUB2(I,K+1)* 1039 1 (HP166666-CSUB2(I,K+1)*H41666M2)) 1040 655 CONTINUE 1041 DO 661 I=MYIS,MYIE 1042 CONTDG(I,LP1)=1.+C(I,LLM1) 1043 TO3DG(I,LP1)=1.+C2(I,LLM1) 1044 661 CONTINUE 1045 DO 663 K=2,L 1046 DO 663 I=MYIS,MYIE 1047 CONTDG(I,K)=ONE+HAF*(C(I,K)+C(I,LM1+K)) 1048 TO3DG(I,K)=ONE+HAF*(C2(I,K)+C2(I,LM1+K)) 1049 663 CONTINUE 1050 C---NOW OBTAIN FLUXES Page 18 Source Listing FST88 2025-03-12 18:22 FST88.F 1051 C 1052 C FOR THE DIAGONAL TERMS... 1053 DO 871 K=2,LP1 1054 DO 871 I=MYIS,MYIE 1055 FLX(I,K)=FLX(I,K)+(DTC(I,K)*EMISDG(I,K) 1056 1 +SS2(I,K)*CONTDG(I,K) 1057 2 +OSS(I,K)*TO3DG(I,K) 1058 3 +CSS(I,K)*CO21(I,K,K))*CLDFAC(I,K,K) 1059 871 CONTINUE 1060 C FOR THE TWO OFF-DIAGONAL TERMS... 1061 DO 873 I=MYIS,MYIE 1062 FLX(I,L)=FLX(I,L)+(CSS(I,LP1)*CO21(I,LP1,L) 1063 1 +DTC(I,LP1)*EMSPEC(I,2) 1064 2 +OSS(I,LP1)*TO31D(I,L) 1065 3 +SS2(I,LP1)*CONT1D(I,L))*CLDFAC(I,LP1,L) 1066 FLX(I,LP1)=FLX(I,LP1)+(CSS(I,L)*CO21(I,L,LP1) 1067 1 +OSS(I,L)*TO31D(I,L) 1068 2 +SS2(I,L)*CONT1D(I,L) 1069 3 +DTC(I,L)*EMSPEC(I,1))*CLDFAC(I,L,LP1) 1070 873 CONTINUE 1071 C 1072 C FINAL SECTION OBTAINS EMISSIVITY HEATING RATES, 1073 C TOTAL HEATING RATES AND THE FLUX AT THE GROUND 1074 C 1075 C .....CALCULATE THE EMISSIVITY HEATING RATES 1076 DO 1101 K=1,L 1077 DO 1101 I=MYIS,MYIE 1078 HEATEM(I,K)=RADCON*(FLX(I,K+1)-FLX(I,K))*DELP(I,K) 1079 1101 CONTINUE 1080 C .....CALCULATE THE TOTAL HEATING RATES 1081 DO 1103 K=1,L 1082 DO 1103 I=MYIS,MYIE 1083 c if(mype.eq.13.and.i.eq.40) then 1084 c print*,'k=',k 1085 c print*,'cts(i,k)=',cts(i,k) 1086 c print*,'ctso(i,k)=',ctso3(i,k) 1087 c print*,'excts(i,k)=',excts(i,k) 1088 c endif 1089 HEATRA(I,K)=HEATEM(I,K)-CTS(I,K)-CTSO3(I,K)+EXCTS(I,K) 1090 1103 CONTINUE 1091 C .....CALCULATE THE FLUX AT EACH FLUX LEVEL USING THE FLUX AT THE 1092 C TOP (FLX1E1+GXCTS) AND THE INTEGRAL OF THE HEATING RATES (VSUM1) 1093 DO 1111 K=1,L 1094 DO 1111 I=MYIS,MYIE 1095 VSUM1(I,K)=HEATRA(I,K)*DELP2(I,K)*RADCON1 1096 1111 CONTINUE 1097 DO 1115 I=MYIS,MYIE 1098 TOPFLX(I)=FLX1E1(I)+GXCTS(I) 1099 c if(mype.eq.13.and.i.eq.40) then 1100 c print*,'flx1e1(i),gxcts(i)=',flx1e1(i),gxcts(i) 1101 c print*,'topflx(i)=',topflx(i) 1102 c endif 1103 FLXNET(I,1)=TOPFLX(I) 1104 1115 CONTINUE 1105 C---ONLY THE SURFACE VALUE OF FLUX (GRNFLX) IS NEEDED UNLESS 1106 C THE THICK CLOUD SECTION IS INVOKED. 1107 DO 1123 K=2,LP1 Page 19 Source Listing FST88 2025-03-12 18:22 FST88.F 1108 DO 1123 I=MYIS,MYIE 1109 c if(mype.eq.13.and.i.eq.40) then 1110 c print*,'k,k-1,flxnet(i,k-1),vsum1(i,k-1)=', 1111 c * k,k-1,flxnet(i,k-1),vsum1(i,k-1) 1112 c endif 1113 FLXNET(I,K)=FLXNET(I,K-1)+VSUM1(I,K-1) 1114 1123 CONTINUE 1115 DO 1125 I=MYIS,MYIE 1116 GRNFLX(I)=FLXNET(I,LP1) 1117 1125 CONTINUE 1118 C 1119 C THIS IS THE THICK CLOUD SECTION.OPTIONALLY,IF THICK CLOUD 1120 C FLUXES ARE TO BE "CONVECTIVELY ADJUSTED",IE,DF/DP IS CONSTANT, 1121 C FOR CLOUDY PART OF GRID POINT, THE FOLLOWING CODE IS EXECUTED. 1122 C***FIRST,COUNT THE NUMBER OF CLOUDS ALONG THE LAT. ROW. SKIP THE 1123 C ENTIRE THICK CLOUD COMPUTATION OF THERE ARE NO CLOUDS. 1124 ICNT=0 1125 DO 1301 I=MYIS,MYIE 1126 ICNT=ICNT+NCLDS(I) 1127 1301 CONTINUE 1128 IF (ICNT.EQ.0) GO TO 6999 1129 C---FIND THE MAXIMUM NUMBER OF CLOUDS IN THE LATITUDE ROW 1130 KCLDS=NCLDS(1) 1131 DO 2106 I=MYIS,MYIE 1132 KCLDS=MAX(NCLDS(I),KCLDS) 1133 2106 CONTINUE 1134 C 1135 C 1136 C***OBTAIN THE PRESSURES AND FLUXES OF THE TOP AND BOTTOM OF 1137 C THE NC'TH CLOUD (IT IS ASSUMED THAT ALL KTOP AND KBTM'S HAVE 1138 C BEEN DEFINED!). 1139 DO 1361 KK=1,KCLDS 1140 KMIN=LP1 1141 KMAX=0 1142 DO 1362 I=MYIS,MYIE 1143 J1=KTOP(I,KK+1) 1144 C IF (J1.EQ.1) GO TO 1362 1145 J3=KBTM(I,KK+1) 1146 IF (J3.GT.J1) THEN 1147 PTOP(I)=P(I,J1) 1148 PBOT(I)=P(I,J3+1) 1149 FTOP(I)=FLXNET(I,J1) 1150 FBOT(I)=FLXNET(I,J3+1) 1151 C***OBTAIN THE "FLUX DERIVATIVE" DF/DP (DELPTC) 1152 DELPTC(I)=(FTOP(I)-FBOT(I))/(PTOP(I)-PBOT(I)) 1153 KMIN=MIN(KMIN,J1) 1154 KMAX=MAX(KMAX,J3) 1155 ENDIF 1156 1362 CONTINUE 1157 KMIN=KMIN+1 1158 C***CALCULATE THE TOT. FLUX CHG. FROM THE TOP OF THE CLOUD, FOR 1159 C ALL LEVELS. 1160 DO 1365 K=KMIN,KMAX 1161 DO 1363 I=MYIS,MYIE 1162 C IF (KTOP(I,KK+1).EQ.1) GO TO 1363 1163 IF(KTOP(I,KK+1).LT.K .AND. K.LE.KBTM(I,KK+1)) THEN 1164 Z1(I,K)=(P(I,K)-PTOP(I))*DELPTC(I)+FTOP(I) Page 20 Source Listing FST88 2025-03-12 18:22 FST88.F 1165 CORIGINAL FLXNET(I,K)=FLXNET(I,K)*(ONE-CAMT(I,KK+1)) + 1166 CORIGINAL1 Z1(I,K)*CAMT(I,KK+1) 1167 c if(mype.eq.13.and.i.eq.40) then 1168 c print*,'k,z1(i,k)=',k,z1(i,k) 1169 c endif 1170 FLXNET(I,K)=Z1(I,K) 1171 ENDIF 1172 1363 CONTINUE 1173 1365 CONTINUE 1174 1361 CONTINUE 1175 C***USING THIS FLUX CHG. IN THE CLOUDY PART OF THE GRID BOX, OBTAIN 1176 C THE NEW FLUXES, WEIGHTING THE CLEAR AND CLOUDY FLUXES:AGAIN, ONLY 1177 C THE FLUXES IN THICK-CLOUD LEVELS WILL EVENTUALLY BE USED. 1178 C DO 6051 K=1,LP1 1179 C DO 6051 I=MYIS,MYIE 1180 C FLXNET(I,K)=FLXNET(I,K)*(ONE-CAMT(I,NC)) + 1181 C 1 Z1(I,K)*CAMT(I,NC) 1182 C051 CONTINUE 1183 C***MERGE FLXTHK INTO FLXNET FOR APPROPRIATE LEVELS. 1184 C DO 1401 K=1,LP1 1185 C DO 1401 I=MYIS,MYIE 1186 C IF (K.GT.ITOP(I) .AND. K.LE.IBOT(I) 1187 C 1 .AND. (NC-1).LE.NCLDS(I)) THEN 1188 if(mype.eq.13.and.i.eq.40) then 1189 print*,'k,flxthk(i,k)=',k,flxthk(i,k) 1190 endif 1191 C FLXNET(I,K)=FLXTHK(I,K) 1192 C ENDIF 1193 C401 CONTINUE 1194 C 1195 C******END OF CLOUD LOOP***** 1196 6001 CONTINUE 1197 6999 CONTINUE 1198 C***THE FINAL STEP IS TO RECOMPUTE THE HEATING RATES BASED ON THE 1199 C REVISED FLUXES: 1200 DO 6101 K=1,L 1201 DO 6101 I=MYIS,MYIE 1202 c if(mype.eq.13.and.i.eq.40) then 1203 c print*,'i,k=',i,k 1204 c print*,'radcon=',radcon 1205 c print*,'flxnet(i,k+1)=',flxnet(i,k+1) 1206 c print*,'flxnet(i,k)=',flxnet(i,k) 1207 c print*,'delp(i,k)=',delp(i,k) 1208 c endif 1209 HEATRA(I,K)=RADCON*(FLXNET(I,K+1)-FLXNET(I,K))*DELP(I,K) 1210 6101 CONTINUE 1211 C THE THICK CLOUD SECTION ENDS HERE. 1212 RETURN 1213 END Page 21 Source Listing FST88 2025-03-12 18:22 Entry Points FST88.F ENTRY POINTS Name fst88_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1001 Label 813 808,809 101 Label 489 480,481 1011 Label 819 815,816 1012 Label 830 820 1013 Label 834 832 1014 Label 835 831 103 Label 496 490,491 105 Label 505 498 1101 Label 1079 1076,1077 1103 Label 1090 1081,1082 1111 Label 1096 1093,1094 1115 Label 1104 1097 1123 Label 1114 1107,1108 1125 Label 1117 1115 1301 Label 1127 1125 131 Label 659 656,657 1361 Label 1174 1139 1362 Label 1156 1142 1363 Label 1172 1161 1365 Label 1173 1160 143 Label 663 660,661 1803 Label 850 848 1804 Label 885 883 1808 Label 764 762 2106 Label 1133 1131 302 Label 752 737,738 3021 Label 725 722,723 3022 Label 756 753,754 3023 Label 760 757,758 303 Label 787 781,782 305 Label 775 768,769 307 Label 779 776 321 Label 901 841 3218 Label 847 844,845 322 Label 864 858,859 3221 Label 877 866,867 3223 Label 881 878,879 331 Label 983 974 3423 Label 893 887,888 3425 Label 900 894,895 4112 Label 516 513,514 4114 Label 512 508,509 4212 Label 526 523,524 4214 Label 522 518,519 Page 22 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References 4312 Label 536 533,534 4314 Label 532 528,529 4412 Label 546 543,544 4414 Label 542 538,539 4512 Label 556 553,554 4514 Label 552 548,549 4612 Label 566 563,564 4614 Label 562 558,559 4712 Label 576 573,574 4714 Label 572 568,569 4812 Label 586 583,584 4814 Label 582 578,579 4912 Label 596 593,594 4914 Label 592 588,589 5012 Label 606 603,604 5014 Label 602 598,599 5112 Label 616 613,614 5114 Label 612 608,609 5212 Label 626 623,624 5214 Label 622 618,619 5312 Label 636 633,634 5314 Label 632 628,629 5412 Label 646 643,644 5414 Label 642 638,639 6001 Label 1196 601 Label 992 988,989 603 Label 997 993,994 6101 Label 1210 1200,1201 618 Label 987 984,985 625 Label 1001 998 631 Label 1011 1008,1009 641 Label 1018 1012 643 Label 1022 1019,1020 651 Label 1031 1027,1028 655 Label 1040 1033,1034 661 Label 1044 1041 663 Label 1049 1045,1046 6999 Label 1197 1128 803 Label 733 731 821 Label 940 936 823 Label 945 941,942 833 Label 952 949 851 Label 963 960,961 861 Label 973 967 871 Label 1059 1053,1054 873 Label 1070 1061 901 Label 675 669,670 903 Label 682 676,677 998 Label 803 801 999 Label 807 804,805 AINT Func 483 scalar 483,493 ALP Local 435 R(4) 4 2 6461 991,996,999,1000,1010 AVEPHI Local 405 R(4) 4 2 3266 724,732,736,846,849,857,874,950,95 1,953 AVMO3 Local 425 R(4) 4 2 3266 860,868,870 Page 23 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References AVPHO3 Local 425 R(4) 4 2 3266 861,868 AVVO2 Local 422 R(4) 4 2 3266 862,873,875 BANDTA Common 286 5900 BDCOMB Common 355 736 SAVE BDWIDE Common 314 56 SAVE C Local 426 R(4) 4 2 6461 1010,1013,1014,1015,1016,1017,1021 ,1036,1042,1047 C2 Local 426 R(4) 4 2 6461 1038,1043,1048 CAMT Dummy 66 R(4) 4 2 3266 ARG,INOUT CLDFAC Dummy 66 R(4) 4 3 150236 ARG,INOUT 774,786,792,806,811,817,818,821,89 2,899,1058,1065,1069 CNTTAU Local 459 R(4) 4 2 3266 750,751,771,784,806,863,981 CNTVAL Dummy 68 R(4) 4 2 3266 ARG,INOUT 978,980,1029,1030 CO21 Dummy 67 R(4) 4 3 150236 ARG,INOUT 755,759,785,876,880,891,898,1013,1 014,1016,1021,1058,1062,1066 CO2NBL Dummy 67 R(4) 4 2 3195 ARG,INOUT 763,884,982 CO2SP Local 415 R(4) 4 2 3266 755,773,794 CO2SP1 Dummy 67 R(4) 4 2 3266 ARG,INOUT 794 CO2SP2 Dummy 67 R(4) 4 2 3266 ARG,INOUT 794 CONT1D Local 424 R(4) 4 2 3266 863,890,897,981,1065,1068 CONTDG Local 439 R(4) 4 2 3266 1042,1047,1056 CSOUR Local 420 R(4) 4 2 3266 662,679,773,778,791 CSS Local 417 R(4) 4 2 3266 679,785,891,898,1058,1062,1066 CSUB Local 435 R(4) 4 2 6461 1029,1030,1035,1036,1037 CSUB2 Local 436 R(4) 4 2 6461 1035,1038,1039 CTS Local 414 R(4) 4 2 3195 810,1089 CTSO3 Local 413 R(4) 4 2 3195 791,1089 DELP Dummy 65 R(4) 4 2 3195 ARG,INOUT 793,810,990,995,1000,1078,1209 DELP2 Dummy 65 R(4) 4 2 3195 ARG,INOUT 793,1014,1095 DELPR1 Local 438 R(4) 4 2 3266 990,991,1029 DELPR2 Local 438 R(4) 4 2 3266 995,996,1030 DELPTC Local 430 R(4) 4 1 71 1152,1164 DSORC Local 434 R(4) 4 2 3266 511,515,521,525,531,535,541,545,55 1,555,561,565,571,575,581,585,591, 595,601,605,611,615,621,625,631,63 5,641,645 DT Local 453 R(4) 4 2 3266 485,500,504,515,525,535,545,555,56 5,575,585,595,605,615,625,635,645, 736 DTC Local 419 R(4) 4 2 3266 680,786,892,899,1055,1063,1069 DTE2 Local 454 R(4) 4 2 3266 495,500,503,736,857 DTSP Local 455 R(4) 4 2 142 503,504,953 E1CTS1 Local 448 R(4) 4 2 3266 735,817,822 E1CTS2 Local 448 R(4) 4 2 3195 735,818 E1CTW1 Local 449 R(4) 4 2 3266 735,811,817,822 E1CTW2 Local 449 R(4) 4 2 3195 735,811,818 E1E290 Subr 735 735 E1FLX Local 414 R(4) 4 2 3266 735,770,777 E290 Subr 857 857 E2SPEC Subr 953 953 E3V88 Subr 956 956 EMD Local 450 R(4) 4 2 6461 956,962,968,970,971 EMISDG Local 439 R(4) 4 2 3266 962,970,1055 EMISS Local 405 R(4) 4 2 3266 735,786,857,892,953,969 EMISSB Local 406 R(4) 4 2 3266 857,899 Page 24 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References EMPL Dummy 70 R(4) 4 2 6461 ARG,INOUT 951,956,968,971 EMSPEC Local 432 R(4) 4 2 142 968,971,1063,1069 EMX1 Dummy 70 R(4) 4 1 71 ARG,INOUT 732,849,969 EMX2 Dummy 70 R(4) 4 1 71 ARG,INOUT 972 EXCTS Local 413 R(4) 4 2 3195 791,1089 EXP Func 744 scalar 744,745,750,872,874,978,979 FAC1 Local 437 R(4) 4 2 3266 739,740,741,868,869,871,975,976,97 7 FBOT Local 432 R(4) 4 1 71 1150,1152 FLX Local 458 R(4) 4 2 3266 770,777,783,889,896,1055,1062,1066 ,1078 FLX1E1 Local 404 R(4) 4 1 71 821,833,1098 FLXNET Local 441 R(4) 4 2 3266 1103,1113,1116,1149,1150,1170,1209 FLXTHK Local 444 R(4) 4 2 3266 1189 FST88 Subr 64 FTOP Local 431 R(4) 4 1 71 1149,1152,1164 FXO Local 453 R(4) 4 2 3266 484,488,499,502,736 FXOE2 Local 454 R(4) 4 2 3266 494,499,501,736,857 FXOSP Local 455 R(4) 4 2 142 501,502,953 GLB_TABLE Common 200 336 SAVE GRNFLX Dummy 64 R(4) 4 1 71 ARG,INOUT 1116 GXCTS Local 404 R(4) 4 1 71 791,1098 HCON Common 78 872 SAVE HEATEM Local 422 R(4) 4 2 3266 1078,1089 HEATRA Dummy 64 R(4) 4 2 3195 ARG,INOUT 1089,1095,1209 I Local 481 I(4) 4 scalar 481,483,484,485,488,491,493,494,49 5,498,499,500,501,502,503,504,508, 510,511,514,515,518,520,521,524,52 5,528,530,531,534,535,538,540,541, 544,545,548,550,551,554,555,558,56 0,561,564,565,568,570,571,574,575, 578,580,581,584,585,588,590,591,59 4,595,598,600,601,604,605,608,610, 611,614,615,618,620,621,624,625,62 8,630,631,634,635,638,640,641,644, 645,657,658,661,662,670,671,677,67 8,679,680,681,723,724,731,732,738, 739,740,741,744,745,746,750,751,75 4,755,758,759,762,763,769,770,771, 772,773,774,776,777,778,782,783,78 4,785,786,801,802,805,806,809,810, 811,812,816,817,818,820,821,822,83 2,833,845,846,848,849,859,860,861, 862,863,867,868,869,870,871,872,87 3,874,875,876,879,880,883,884,888, 889,890,891,892,895,896,897,898,89 9,936,937,938,939,942,943,944,949, 950,951,961,962,967,968,969,970,97 1,972,974,975,976,977,978,979,980, 981,982,985,986,989,990,991,994,99 5,996,998,999,1000,1009,1010,1012, 1013,1014,1015,1016,1017,1020,1021 ,1028,1029,1030,1034,1035,1036,103 7,1038,1039,1041,1042,1043,1046,10 47,1048,1054,1055,1056,1057,1058,1 Page 25 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References 061,1062,1063,1064,1065,1066,1067, 1068,1069,1077,1078,1082,1089,1094 ,1095,1097,1098,1103,1108,1113,111 5,1116,1125,1126,1131,1132,1142,11 43,1145,1147,1148,1149,1150,1152,1 161,1163,1164,1170,1188,1189,1201, 1209 IBOT Local 427 I(4) 4 1 71 ICNT Local 1124 I(4) 4 scalar 1124,1126,1128 IDIM1 Param 141 I(4) 4 scalar 390,391,392,393,394,395,396,397,39 8,399,400,401,402,403,404,405,406, 408,409,410,411,413,414,415,416,41 7,418,419,420,422,423,424,425,426, 427,428,430,431,432,434,435,436,43 7,438,439,440,441,442,443,444,445, 448,449,450,453,454,455,458,459,46 7,468,469,470 IDIM2 Param 141 I(4) 4 scalar 390,391,392,393,394,395,396,397,39 8,399,400,401,402,403,404,405,406, 408,409,410,411,413,414,415,416,41 7,418,419,420,422,423,424,425,426, 427,428,430,431,432,434,435,436,43 7,438,439,440,441,442,443,444,445, 448,449,450,453,454,455,458,459 IGSTL Param 136 I(4) 4 scalar 141 IGSTR Param 136 I(4) 4 scalar 141 IM Param 124 I(4) 4 scalar 138,141,206,207,208,209,210,214,21 5,220,223 IMAX Param 223 I(4) 4 scalar 223,236,383 INDTC Local 428 I(4) 4 1 71 INLTE Param 235 I(4) 4 scalar 235 INLTEP Param 235 I(4) 4 scalar INPES Param 133 I(4) 4 scalar 138,141,190,191,192,201,202 ITAIL Param 138 I(4) 4 scalar ITOP Local 427 I(4) 4 1 71 IXO Local 442 I(4) 4 2 3266 488,510,511,520,521,530,531,540,54 1,550,551,560,561,570,571,580,581, 590,591,600,601,610,611,620,621,63 0,631,640,641 J1 Local 1143 I(4) 4 scalar 1143,1146,1147,1149,1153 J3 Local 1145 I(4) 4 scalar 1145,1146,1148,1150,1154 JDIM1 Param 142 I(4) 4 scalar JDIM2 Param 142 I(4) 4 scalar JGSTL Param 137 I(4) 4 scalar 142 JGSTR Param 137 I(4) 4 scalar 142 JM Param 124 I(4) 4 scalar 139,142,206,207,208,209,210,214,21 5,220 JNPES Param 133 I(4) 4 scalar 139,142,190,191,192,201,202 JTAIL Param 139 I(4) 4 scalar K Local 480 I(4) 4 scalar 480,483,484,485,488,490,493,494,49 5,509,510,511,513,515,519,520,521, 523,525,529,530,531,533,535,539,54 0,541,543,545,549,550,551,553,555, 559,560,561,563,565,569,570,571,57 3,575,579,580,581,583,585,589,590, Page 26 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References 591,593,595,599,600,601,603,605,60 9,610,611,613,615,619,620,621,623, 625,629,630,631,633,635,639,640,64 1,643,645,656,658,660,662,669,671, 676,678,679,680,681,722,724,737,73 9,740,741,744,745,746,750,751,753, 755,757,759,768,770,771,772,773,77 4,804,806,808,810,811,812,815,817, 818,831,833,841,842,844,846,858,86 0,861,862,863,866,868,869,870,871, 872,873,874,875,876,878,880,884,88 7,889,891,892,894,896,897,898,899, 941,943,944,960,962,984,986,988,99 0,991,993,995,996,1008,1010,1019,1 021,1027,1029,1030,1033,1035,1036, 1037,1038,1039,1045,1047,1048,1053 ,1055,1056,1057,1058,1076,1078,108 1,1089,1093,1095,1107,1113,1160,11 63,1164,1170,1189,1200,1209 KBTM Dummy 66 I(4) 4 2 3266 ARG,INOUT 1145,1163 KCLDS Local 1130 I(4) 4 scalar 1130,1132,1139 KK Local 844 I(4) 4 scalar 844,846,858,860,861,862,863,866,86 8,869,870,871,872,873,874,875,876, 1139,1143,1145,1163 KLEN Local 842 I(4) 4 scalar 842,857 KMAX Local 1141 I(4) 4 scalar 1141,1154,1160 KMIN Local 1140 I(4) 4 scalar 1140,1153,1157,1160 KO2 Param 238 I(4) 4 scalar 239 KO21 Param 239 I(4) 4 scalar KO2M Param 239 I(4) 4 scalar KP Local 781 I(4) 4 scalar 781,783,784,785,786,878,880,887,88 9,890,891,892,894,896,897,898,899 KTOP Dummy 66 I(4) 4 2 3266 ARG,INOUT 1143,1163 L Param 222 I(4) 4 scalar 226,227,228,231,233,391,392,397,39 9,400,402,413,414,415,448,449,458, 490,499,500,676,722,737,753,757,80 4,808,815,831,937,938,939,941,944, 950,951,960,962,975,976,977,978,97 9,980,981,982,984,993,999,1000,101 3,1014,1015,1016,1017,1019,1045,10 62,1064,1065,1066,1067,1068,1069,1 076,1081,1093,1200 LL Param 228 I(4) 4 scalar 228,229,999,1014 LL3P Param 233 I(4) 4 scalar 236 LL3PI Param 236 I(4) 4 scalar LLM1 Param 229 I(4) 4 scalar 1015,1042,1043 LLM2 Param 229 I(4) 4 scalar 1033 LLM3 Param 229 I(4) 4 scalar LLP1 Param 228 I(4) 4 scalar 236,411,426,435,436,450,939,971,10 00,1008,1016 LLP1I Param 236 I(4) 4 scalar LLP2 Param 228 I(4) 4 scalar LLP3 Param 228 I(4) 4 scalar LM Param 124 I(4) 4 scalar 215,222 LM1 Param 227 I(4) 4 scalar 501,502,503,504,732,841,849,981,98 Page 27 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References 8,1000,1016,1017,1021,1027,1047,10 48 LM2 Param 227 I(4) 4 scalar LM3 Param 227 I(4) 4 scalar LOG Func 986 scalar 986 LP1 Param 226 I(4) 4 scalar 230,231,232,236,383,390,391,392,39 3,394,395,396,397,398,401,405,406, 408,409,410,414,415,416,417,418,41 9,420,422,423,424,425,434,437,438, 439,440,441,442,443,444,445,448,44 9,453,454,458,459,480,499,500,509, 513,519,523,529,533,539,543,549,55 3,559,563,569,573,579,583,589,593, 599,603,609,613,619,623,629,633,63 9,643,656,660,669,732,768,781,821, 822,844,849,858,866,878,887,894,93 8,968,970,991,1000,1013,1014,1015, 1016,1017,1030,1042,1043,1053,1062 ,1063,1064,1065,1066,1069,1107,111 6,1140 LP121 Param 232 I(4) 4 scalar LP1I Param 236 I(4) 4 scalar LP1M Param 230 I(4) 4 scalar 230 LP1M1 Param 230 I(4) 4 scalar LP1V Param 231 I(4) 4 scalar 383 LP2 Param 226 I(4) 4 scalar 468,470 LP3 Param 226 I(4) 4 scalar LSM Param 124 I(4) 4 scalar MAPPINGS Common 219 5024 SAVE MAX Func 488 scalar 488,1132,1154 MIN Func 1153 scalar 1153 MPPCOM Common 164 1724 SAVE NB Param 234 I(4) 4 scalar 237 NB1 Param 237 I(4) 4 scalar NBLM Param 225 I(4) 4 scalar NBLW Param 224 I(4) 4 scalar 286,287,288 NBLX Param 224 I(4) 4 scalar NBLY Param 224 I(4) 4 scalar 225,232,355,356,357,386,387,420 NCLDS Dummy 66 I(4) 4 1 71 ARG,INOUT 1126,1130,1132 NCOL Param 223 I(4) 4 scalar NNLTE Param 235 I(4) 4 scalar OSS Local 417 R(4) 4 2 3266 678,783,889,896,1057,1064,1067 OVER1D Local 423 R(4) 4 2 3266 745,755,759,763,874,876,880,884,97 9,982 P Dummy 65 R(4) 4 2 3266 ARG,INOUT 793,990,995,1000,1014,1015,1016,10 17,1147,1148,1164 PBOT Local 430 R(4) 4 1 71 1148,1152 PHYCON Common 72 80 SAVE PRESS Dummy 65 R(4) 4 2 3266 ARG,INOUT 792,990,995,1000,1014,1015,1016,10 17 PTOP Local 430 R(4) 4 1 71 1147,1152,1164 QH2O Dummy 65 R(4) 4 2 3266 ARG,INOUT RLOG Local 458 R(4) 4 2 3195 763,884,982,986,991,996,999,1000 SORC Local 420 R(4) 4 3 48990 515,525,535,545,555,565,575,585,59 5,605,615,625,635,645,658,662,678, Page 28 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Object Declared Type Bytes Dimen Elements Attributes References 772,777,791 SPA88 Subr 791 791 SQRT Func 741 scalar 741,745,870,874,977,979,991,996,10 00 SS1 Local 418 R(4) 4 2 3266 658,681,771,777,812 SS2 Local 418 R(4) 4 2 3266 681,784,890,897,1056,1065,1068 T Dummy 65 R(4) 4 2 3266 ARG,INOUT 493,495,736,938,939,943,944 TABCOM Common 383 133920 SAVE TC Local 419 R(4) 4 2 3266 671,680,770,777,810,817,821 TEMP Dummy 65 R(4) 4 2 3266 ARG,INOUT 483,485,671,736,792,937,938,939 TEMPCOM Common 205 6603768 SAVE TO31D Local 424 R(4) 4 2 3266 872,889,896,978,1064,1067 TO3DG Local 440 R(4) 4 2 3266 1043,1048,1057 TO3SP Local 416 R(4) 4 2 3266 744,772,783,793 TO3SPC Local 415 R(4) 4 2 3195 740,744,793 TOPFLX Dummy 64 R(4) 4 1 71 ARG,INOUT 1098,1103 TOPO Common 213 17478548 SAVE TOTEVV Local 459 R(4) 4 2 3266 751,863,981 TOTO3 Dummy 69 R(4) 4 2 3266 ARG,INOUT 739,741,860 TOTPHI Dummy 69 R(4) 4 2 3266 ARG,INOUT 724,745,846 TOTVO2 Dummy 69 R(4) 4 2 3266 ARG,INOUT 744,746,750,793,862 TPHIO3 Dummy 69 R(4) 4 2 3266 ARG,INOUT 739,861 TPL Local 450 R(4) 4 2 6461 937,938,939,943,944,956 VAR1 Dummy 68 R(4) 4 2 3195 ARG,INOUT 792 VAR2 Dummy 68 R(4) 4 2 3195 ARG,INOUT 792,950,951,979 VAR3 Dummy 68 R(4) 4 2 3195 ARG,INOUT 975,977 VAR4 Dummy 68 R(4) 4 2 3195 ARG,INOUT 975 VSUM1 Local 443 R(4) 4 2 3266 1095,1113 VTMP3 Local 434 R(4) 4 2 3266 483,484,485,493,494,495,510,515,52 0,525,530,535,540,545,550,555,560, 565,570,575,580,585,590,595,600,60 5,610,615,620,625,630,635,640,645, 802,806,812,817,833,869,872,976,97 8 Z1 Local 445 R(4) 4 2 3266 1164,1170 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 745,874,979 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 741,870,977 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 Page 29 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References AWIDE R(4) 4 0 scalar COM 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 BO3CM R(4) 4 716 1 3 COM BO3RND R(4) 4 5880 1 3 COM 739,868,975 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 CSUBP R(4) 4 4 scalar COM DIFFCTR R(4) 4 8 scalar COM DSRCE R(4) 4 132240 2 420 COM 511,521,531,541,551,561,571,581,59 1,601,611,621,631,641 EIGHT R(4) 4 20 scalar COM EM1 R(4) 4 9600 2 5040 COM EM1WDE R(4) 4 29760 2 5040 COM EM3 R(4) 4 110400 2 5040 COM FIFTY R(4) 4 12 scalar COM FIVE R(4) 4 24 scalar COM FOUR R(4) 4 28 scalar COM 741,870,977 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 GP0INV R(4) 4 76 scalar COM GRAVDR R(4) 4 16 scalar COM 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 Page 30 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 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 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 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 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 Page 31 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 H25M31 R(4) 4 552 scalar COM H26E2 R(4) 4 696 scalar COM H26M30 R(4) 4 544 scalar COM 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 H3M4 R(4) 4 860 scalar COM H3P5 R(4) 4 828 scalar COM H3P6 R(4) 4 144 scalar COM H41666M2 R(4) 4 636 scalar COM 1037,1039 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 Page 32 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 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 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 740,869,938,939,976,1021,1047,1048 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 744,745,750,872,874,978,979 Page 33 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References HM2M2 R(4) 4 644 scalar COM HM6666M2 R(4) 4 628 scalar COM 1010 HM8E1 R(4) 4 864 scalar COM HMP5 R(4) 4 640 scalar COM 1036,1038 HMP575 R(4) 4 588 scalar COM HMP66667 R(4) 4 624 scalar COM 1010 HMP805 R(4) 4 840 scalar COM HNINETY R(4) 4 4 scalar COM HP1 R(4) 4 196 scalar COM 483,493 HP118666 R(4) 4 820 scalar COM HP144 R(4) 4 712 scalar COM HP166666 R(4) 4 632 scalar COM 1037,1039 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 IND I(4) 4 0 1 237 COM INDX2 I(4) 4 948 1 2116 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 Page 34 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 KMAXV I(4) 4 9412 1 46 COM KMAXVM I(4) 4 9596 scalar COM L2GI I(4) 4 956 1 239 COM 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 481,491,498,508,514,518,524,528,53 4,538,544,548,554,558,564,568,574, 578,584,588,594,598,604,608,614,61 8,624,628,634,638,644,657,661,670, 677,723,731,738,754,758,762,769,77 6,782,801,805,809,816,820,832,845, 848,859,867,879,883,888,895,936,94 2,949,961,967,974,985,989,994,998, 1009,1012,1020,1028,1034,1041,1046 ,1054,1061,1077,1082,1094,1097,110 8,1115,1125,1131,1142,1161,1201 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 481,491,498,508,514,518,524,528,53 4,538,544,548,554,558,564,568,574, 578,584,588,594,598,604,608,614,61 8,624,628,634,638,644,657,661,670, 677,723,731,738,754,758,762,769,77 6,782,801,805,809,816,820,832,845, 848,859,867,879,883,888,895,936,94 2,949,961,967,974,985,989,994,998, 1009,1012,1020,1028,1034,1041,1046 ,1054,1061,1077,1082,1094,1097,110 8,1115,1125,1131,1142,1161,1201 MYIS1 I(4) 4 44 scalar COM Page 35 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 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 Page 36 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 1188 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 741,870,871,977,1013,1014,1016,102 1,1047,1048 P0 R(4) 4 24 scalar COM P0INV R(4) 4 72 scalar COM P0X2 R(4) 4 36 scalar COM P0XZP2 R(4) 4 28 scalar COM P0XZP8 R(4) 4 32 scalar COM QUARTR R(4) 4 48 scalar COM 969,1010 RADCON R(4) 4 40 scalar COM 810,1078,1209 RADCON1 R(4) 4 64 scalar COM 1095 RATCO2MW R(4) 4 56 scalar COM RATH2OMW R(4) 4 60 scalar COM 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 746,875,980 SKO2D R(4) 4 44 scalar COM SKO3R R(4) 4 52 scalar COM 744,873,978,1035 SOURCE R(4) 4 130560 2 420 COM 510,520,530,540,550,560,570,580,59 0,600,610,620,630,640 TABLE1 R(4) 4 49920 2 5040 COM TABLE2 R(4) 4 70080 2 5040 COM TABLE3 R(4) 4 90240 2 5040 COM 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 Page 37 Source Listing FST88 2025-03-12 18:22 Symbol Table FST88.F Name Type Bytes Offset Dimen Elements Attributes References 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 485,495 THREE R(4) 4 32 scalar COM TTVG R(4) 4 371884 2 92971 COM TWO R(4) 4 36 scalar COM 970,971 ZERO R(4) 4 52 scalar COM Page 38 Source Listing FST88 2025-03-12 18:22 Subprograms/Common Blocks FST88.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BANDTA Common 286 5900 BDCOMB Common 355 736 SAVE BDWIDE Common 314 56 SAVE FST88 Subr 64 GLB_TABLE Common 200 336 SAVE HCON Common 78 872 SAVE MAPPINGS Common 219 5024 SAVE MPPCOM Common 164 1724 SAVE PHYCON Common 72 80 SAVE TABCOM Common 383 133920 SAVE TEMPCOM Common 205 6603768 SAVE TOPO Common 213 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 -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check noshape -check nostack Page 39 Source Listing FST88 2025-03-12 18:22 FST88.F -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 -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg -init noarrays -init nohuge -init noinfinity Page 40 Source Listing FST88 2025-03-12 18:22 FST88.F -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 : FST88.lst no -o COMPILER: Intel(R) Fortran 19.1-1655