Page 1 Source Listing SPA88 2025-03-12 18:23 /tmp/ifortyBmJlu.i 1 # 1 "SPA88.F" 2 C SUBROUTINE SPA88 COMPUTES EXACT CTS HEATING RATES AND FLUXES AND 3 C CORRESPONDING CTS EMISSIVITY QUANTITIES FOR H2O,CO2 AND O3. 4 C INPUTS: (COMMON BLOCKS) 5 C ACOMB,BCOMB,APCM,BPCM BDCOMB 6 C ATPCM,BTPCM,BETACM BDCOMB 7 C BETINW BDWIDE 8 C TEMP,PRESS RADISW 9 C VAR1,VAR2,P,DELP,DELP2 KDACOM 10 C TOTVO2,TO3SP,TO3SPC TFCOM 11 C CO2SP1,CO2SP2,CO2SP TFCOM 12 C CLDFAC CLDCOM 13 C SKO2D TABCOM 14 C SORC,CSOUR SRCCOM 15 C OUTPUTS: 16 C EXCTS,CTSO3 TFCOM 17 C GXCTS RDFLUX 18 C CALLED BY: 19 C FST88 20 C CALLS: 21 C 22 SUBROUTINE SPA88(EXCTS,CTSO3,GXCTS,SORC,CSOUR, 23 1 CLDFAC,TEMP,PRESS,VAR1,VAR2, 24 2 P,DELP,DELP2,TOTVO2,TO3SP,TO3SPC, 25 3 CO2SP1,CO2SP2,CO2SP) 26 C 27 COMMON/PHYCON/AMOLWT,CSUBP,DIFFCTR,G,GRAVDR,O3DIFCTR,P0, 28 * P0XZP2,P0XZP8,P0X2,RADCON,RGAS,RGASSP,SECPDA 29 COMMON/PHYCON/RATCO2MW,RATH2OMW 30 COMMON/PHYCON/RADCON1 31 COMMON/PHYCON/GINV,P0INV,GP0INV 32 save /PHYCON/ 33 COMMON/HCON/HUNDRED,HNINETY,SIXTY,FIFTY,TEN,EIGHT,FIVE, 34 * FOUR,THREE,TWO,ONE,HAF,QUARTR,ZERO 35 COMMON/HCON/H83E26,H71E26,H1E15,H1E13,H1E11,H1E8,H4E5, 36 * H165E5,H5725E4,H488E4,H1E4,H24E3,H20788E3, 37 * H2075E3,H1224E3,H5E2,H3082E2,H3E2,H2945E2, 38 * H23E2,H15E2,H35E1,H3P6,H181E1,H18E1,H2P9,H2P8, 39 * H2P5,H1P8,H1P4387,H1P4,H1P25892,HP8,HP518, 40 * HP369,HP1 41 COMMON/HCON/H44871M2,H559M3,H1M3,H987M4,H285M4,H1M4, 42 * H6938M5,H394M5,H37412M5,H1439M5,H128M5,H1M5, 43 * H7M6,H4999M6,H25452M6,H1M6,H391M7,H1174M7, 44 * H8725M8,H327M8,H257M8,H1M8,H23M10,H14M10, 45 * H11M10,H1M10,H83M11,H82M11,H8M11,H77M11, 46 * H72M11,H53M11,H48M11,H44M11,H42M11,H37M11, 47 * H35M11,H32M11,H3M11,H28M11,H24M11,H23M11, 48 * H2M11,H18M11,H15M11,H14M11,H114M11,H11M11, 49 * H1M11,H96M12,H93M12,H77M12,H74M12,H65M12, 50 * H62M12,H6M12,H45M12,H44M12,H4M12,H38M12, 51 * H37M12,H3M12,H29M12,H28M12,H24M12,H21M12, 52 * H16M12,H14M12,H12M12,H8M13,H46M13,H36M13, 53 * H135M13,H12M13,H1M13,H3M14,H15M14,H14M14, 54 * H1M17,H1M18,H1M19,H1M20,H1M21,H1M22,H1M23, 55 * H1M24,H26M30,H14M30,H25M31,H21M31,H12M31, 56 * H9M32,H55M32,H45M32,H4M33,H62M34,H1M60 57 COMMON/HCON/HMP575,HM13EZ,HM19EZ,HM1E1,HM181E1,HM1E2 Page 2 Source Listing SPA88 2025-03-12 18:23 SPA88.F 58 COMMON/HCON/H1E6,H2E6,H1M2,HMP66667,HM6666M2,HP166666, 59 * H41666M2,HMP5,HM2M2,H29316E2,H1226E1,H3116E1, 60 * H9P94,HP6,H625M2,HP228,HP60241,HM1797E1, 61 * H8121E1,H2E2,HM1EZ,H26E2,H44194M2,H1P41819 62 COMMON/HCON/HP219,HP144,HP816,H69766E5,H235M3,HP26, 63 * H129M2,H75826M4,H1P082,HP805,H1386E2, 64 * H658M2,H1036E2,H2118M2,H42M2,H323M4, 65 * H67390E2,HP3795,HP5048,H102M5,H451M6 66 COMMON/HCON/H16E1,HM161E1,H161E1,H3M3,H101M16, 67 * HM1597E1,H25E2,HP118666,H15M5,H3P5,H18E3, 68 * H6P08108,HMP805,HP602409,HP526315, 69 * H28571M2,H1M16 70 COMMON/HCON/H3M4 71 COMMON/HCON/HM8E1 72 COMMON/HCON/H28E1 73 save /HCON/ 74 C----------------------------------------------------------------------- 75 INCLUDE "parmeta" 99 INCLUDE "mpp.h" 100 # 1 "./sp.h" 1 101 # 4 102 103 # 77 "SPA88.F" 2 104 C----------------------------------------------------------------------- 105 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 106 C IMAX = NO. POINTS ALONG THE LAT. CIRCLE USED IN CALCS. 107 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 108 C***NOTE: THE USER NORMALLY WILL MODIFY ONLY THE IMAX AND L PARAMETERS 109 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 110 C BANDTA FOR DEFINITION 111 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 112 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE 113 C BDCOMB FOR DEFINITION 114 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 115 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 116 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 117 C FROM THE ABOVE PARAMETERS. 176 PARAMETER (L=LM) 177 PARAMETER (IMAX=IM,NCOL=IMAX) 178 PARAMETER (NBLW=163,NBLX=47,NBLY=15) 179 PARAMETER (NBLM=NBLY-1) 180 PARAMETER (LP1=L+1,LP2=L+2,LP3=L+3) 181 PARAMETER (LM1=L-1,LM2=L-2,LM3=L-3) 182 PARAMETER (LL=2*L,LLP1=LL+1,LLP2=LL+2,LLP3=LL+3) 183 PARAMETER (LLM1=LL-1,LLM2=LL-2,LLM3=LL-3) 184 PARAMETER (LP1M=LP1*LP1,LP1M1=LP1M-1) 185 PARAMETER (LP1V=LP1*(1+2*L/2)) 186 PARAMETER (LP121=LP1*NBLY) 187 PARAMETER (LL3P=3*L+2) 188 PARAMETER (NB=12) 189 PARAMETER (INLTE=3,INLTEP=INLTE+1,NNLTE=56) 190 PARAMETER (LP1I=IMAX*LP1,LLP1I=IMAX*LLP1,LL3PI=IMAX*LL3P) 191 PARAMETER (NB1=NB-1) 192 PARAMETER (KO2=12) 193 PARAMETER (KO21=KO2+1,KO2M=KO2-1) 194 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 195 C IMAX = NO. POINTS SENT TO RADFS Page 3 Source Listing SPA88 2025-03-12 18:23 SPA88.F 196 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 197 C***NOTE: THE USER NORMALLY WILL MODIFY ONLY THE IMAX AND L PARAMETERS 198 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE 199 C BANDTA FOR DEFINITION 200 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS 201 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE 202 C BDCOMB FOR DEFINITION 203 C INLTE = NO. LEVELS USED FOR NLTE CALCS. 204 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS. 205 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED 206 C FROM THE ABOVE PARAMETERS. 207 C COMMON BLOCK BANDTA CONTAINS RANDOM BAND PARAMETERS FOR THE LW 208 C CALCULATIONS USING 10 CM-1 WIDE BANDS.THE 15 UM CO2 COMPLEX 209 C IS 2 BANDS,560-670 AND 670-800 CM-1. OZONE COEFFICIENTS ARE 210 C IN 3 BANDS,670-800 (14.1 UM),990-1070 AND 1070-1200 (9.6 UM). 211 C THE (NBLW) BANDS NOW INCLUDE: 212 C 56 BANDS, 10 CM-1 WIDE 0 - 560 CM-1 213 C 2 BANDS, 15 UM COMPLEX 560 - 670 CM-1 214 C 670 - 800 CM-1 215 C 3 "CONTINUUM" BANDS 800 - 900 CM-1 216 C 900 - 990 CM-1 217 C 1070 - 1200 CM-1 218 C 1 BAND FOR 9.6 UM BAND 990 - 1070 CM-1 219 C 100 BANDS, 10 CM-1 WIDE 1200 - 2200 CM-1 220 C 1 BAND FOR 4.3 UM SRC 2270 - 2380 CM-1 221 C THUS NBLW PRESENTLY EQUALS 163 222 C ALL BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 223 C 224 C ARNDM = RANDOM "A" PARAMETER FOR (NBLW) BANDS 225 C BRNDM = RANDOM "B" PARAMETER FOR (NBLW) BANDS 226 C BETAD = CONTINUUM COEFFICIENTS FOR (NBLW) BANDS 227 C AP,BP = CAPPHI COEFFICIENTS FOR (NBLW) BANDS 228 C ATP,BTP = CAPPSI COEFFICIENTS FOR (NBLW) BANDS 229 C BANDLO = LOWEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 230 C BANDHI = HIGHEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS 231 C AO3RND = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 232 C BANDS 233 C BO3RND = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 234 C BANDS 235 C AB15 = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS 236 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 237 C DATA FOR ARNDM,BRNDM,AP,BP,ATP,BTP,AO3RND,BO3RND ARE OBTAINED BY 238 C USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS ARE FROM 239 C ROBERTS (1976). 240 COMMON / BANDTA / ARNDM(NBLW),BRNDM(NBLW),BETAD(NBLW),AP(NBLW), 241 1 BP(NBLW),ATP(NBLW),BTP(NBLW),BANDLO(NBLW), 242 2 BANDHI(NBLW),AO3RND(3),BO3RND(3),AB15(2) 243 C 244 C COMMON BLOCK BDWIDE CONTAINS RANDOM BAND PARAMETERS FOR SPECIFIC 245 C WIDE BANDS. AT PRESENT,THE INFORMATION CONSISTS OF 1) RANDOM 246 C MODEL PARAMETERS FOR THE 15 UM BAND,560-800 CM-1; 2) THE 247 C CONTINUUM COEFFICIENT FOR THE 800-990,1070-1200 CM-1 BAND 248 C SPECIFICALLY: 249 C AWIDE = RANDOM "A" PARAMETER FOR BAND 250 C BWIDE = RANDOM "B" PARAMETER FOR BAND 251 C BETAWD = CONTINUUM COEFFICIENTS FOR BAND 252 C APWD,BPWD = CAPPHI COEFFICIENTS FOR BAND Page 4 Source Listing SPA88 2025-03-12 18:23 SPA88.F 253 C ATPWD,BTPWD = CAPPSI COEFFICIENTS FOR BAND 254 C BDLOWD = LOWEST FREQUENCY IN EACH FREQ BAND 255 C BDHIWD = HIGHEST FREQUENCY IN EACH FREQ BAND 256 C AB15WD = THE PRODUCT ARNDM*BRNDM FOR THE ONE BAND 257 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 258 C BETINW = CONT.COEFFICIENT FOR A SPECIFIED WIDE 259 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 260 C SKO2D = 1./BETINW, USED IN SPA88 FOR CONT. COEFFS 261 C SKC1R = BETAWD/BETINW, USED FOR CONT. COEFF. FOR 262 C 15 UM BAND IN FST88 263 C SKO3R = RATIO OF CONT. COEFF. FOR 9.9 UM BAND TO 264 C BETINW, USED FOR 9.6 UM CONT COEFF IN FST88 265 C DATA FOR AWIDE,BWIDE,APWD,BPWD,ATPWD,BTPWD,AO3WD,BO3WD ARE 266 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 267 C ARE FROM ROBERTS (1976). 268 COMMON / BDWIDE / AWIDE,BWIDE,BETAWD, 269 1 APWD,BPWD,ATPWD,BTPWD, 270 2 BDLOWD,BDHIWD,BETINW, 271 3 AB15WD,SKO2D,SKC1R,SKO3R 272 C 273 C COMMON BLOCK BDCOMB CONTAINS RANDOM BAND PARAMETERS FOR THE LW 274 C CALCULATIONS USING COMBINED WIDE FREQUENCY BANDS BETWEEN 160 AND 275 C 1200 CM-1,AS WELL AS THE 2270-2380 BAND FOR SOURCE CALC. 276 C BANDS 1-8: COMBINED WIDE FREQUENCY BANDS FOR 160-560 CM-1 277 C BANDS 9-14: FREQUENCY BANDS,AS IN BANDTA (NARROW BANDS) 278 C FOR 560-1200 CM-1 279 C BAND 15: FREQUENCY BAND 2270-2380 CM-1,USED FOR SOURCE 280 C CALCULATION ONLY 281 C THUS NBLY PRESENTLY EQUALS 15 282 C 283 C BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER 284 C ACOMB = RANDOM "A" PARAMETER FOR (NBLY) BANDS 285 C BCOMB = RANDOM "B" PARAMETER FOR (NBLY) BANDS 286 C BETACM = CONTINUUM COEFFICIENTS FOR (NBLY) BANDS 287 C APCM,BPCM = CAPPHI COEFFICIENTS FOR (NBLY) BANDS 288 C ATPCM,BTPCM = CAPPSI COEFFICIENTS FOR (NBLY) BANDS 289 C BDLOCM = LOWEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 290 C BDHICM = HIGHEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS 291 C AO3CM = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE 292 C BANDS 293 C BO3CM = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE 294 C BANDS 295 C AB15CM = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS 296 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2 297 C BETINC = CONT.COEFFICIENT FOR A SPECIFIED WIDE 298 C FREQ.BAND (800-990 AND 1070-1200 CM-1). 299 C IBAND = INDEX NO OF THE 40 WIDE BANDS USED IN 300 C COMBINED WIDE BAND CALCULATIONS. IN OTHER 301 C WORDS,INDEX TELLING WHICH OF THE 40 WIDE 302 C BANDS BETWEEN 160-560 CM-1 ARE INCLUDED IN 303 C EACH OF THE FIRST 8 COMBINED WIDE BANDS 304 C DATA FOR ACOMB,BCOMB,APCM,BPCM,ATPCM,BTPCM,AO3CM,BO3CM ARE 305 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS 306 C ARE FROM ROBERTS (1976). IBAND INDEX VALUES ARE OBTAINED BY 307 C EXPERIMENTATION. 308 COMMON / BDCOMB / IBAND(40),ACOMB(NBLY),BCOMB(NBLY), 309 1 BETACM(NBLY),APCM(NBLY),BPCM(NBLY),ATPCM(NBLY), Page 5 Source Listing SPA88 2025-03-12 18:23 SPA88.F 310 2 BTPCM(NBLY),BDLOCM(NBLY),BDHICM(NBLY),BETINC, 311 3 AO3CM(3),BO3CM(3),AB15CM(2) 312 C 313 DIMENSION SORC(IDIM1:IDIM2,LP1,NBLY),CSOUR(IDIM1:IDIM2,LP1) 314 DIMENSION CLDFAC(IDIM1:IDIM2,LP1,LP1) 315 DIMENSION TEMP(IDIM1:IDIM2,LP1),PRESS(IDIM1:IDIM2,LP1) 316 DIMENSION VAR1(IDIM1:IDIM2,L),VAR2(IDIM1:IDIM2,L) 317 DIMENSION P(IDIM1:IDIM2,LP1),DELP(IDIM1:IDIM2,L), 318 1 DELP2(IDIM1:IDIM2,L) 319 DIMENSION TOTVO2(IDIM1:IDIM2,LP1),TO3SPC(IDIM1:IDIM2,L), 320 1 TO3SP(IDIM1:IDIM2,LP1) 321 DIMENSION CO2SP1(IDIM1:IDIM2,LP1),CO2SP2(IDIM1:IDIM2,LP1), 322 1 CO2SP(IDIM1:IDIM2,LP1) 323 DIMENSION EXCTS(IDIM1:IDIM2,L),CTSO3(IDIM1:IDIM2,L), 324 1 GXCTS(IDIM1:IDIM2) 325 C 326 DIMENSION PHITMP(IDIM1:IDIM2,L),PSITMP(IDIM1:IDIM2,L), 327 1 TT(IDIM1:IDIM2,L), 328 2 FAC1(IDIM1:IDIM2,L),FAC2(IDIM1:IDIM2,L), 329 3 CTMP(IDIM1:IDIM2,LP1),X(IDIM1:IDIM2,L), 330 4 Y(IDIM1:IDIM2,L), 331 5 TOPM(IDIM1:IDIM2,L),TOPPHI(IDIM1:IDIM2,L), 332 6 CTMP3(IDIM1:IDIM2,LP1),CTMP2(IDIM1:IDIM2,LP1) 333 DIMENSION F(IDIM1:IDIM2,L),FF(IDIM1:IDIM2,L), 334 1 AG(IDIM1:IDIM2,L),AGG(IDIM1:IDIM2,L) 335 C 336 EQUIVALENCE (F,AG,PHITMP) 337 EQUIVALENCE (FF,AGG,PSITMP) 338 C---COMPUTE TEMPERATURE QUANTITIES FOR USE IN PROGRAM 339 DO 101 K=1,L 340 DO 101 I=MYIS,MYIE 341 X(I,K)=TEMP(I,K)-H25E2 342 Y(I,K)=X(I,K)*X(I,K) 343 101 CONTINUE 344 C---INITIALIZE CTMP(I,1),CTMP2(I,1),CTMP3(I,1) TO UNITY; THESE ARE 345 C TRANSMISSION FCTNS AT THE TOP. 346 DO 345 I=MYIS,MYIE 347 CTMP(I,1)=ONE 348 CTMP2(I,1)=1. 349 CTMP3(I,1)=1. 350 345 CONTINUE 351 C***BEGIN LOOP ON FREQUENCY BANDS (1)*** 352 C 353 C---CALCULATION FOR BAND 1 (COMBINED BAND 1) 354 C 355 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 356 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 357 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 358 DO 301 K=1,L 359 DO 301 I=MYIS,MYIE 360 F(I,K)=H44194M2*(APCM(1)*X(I,K)+BPCM(1)*Y(I,K)) 361 FF(I,K)=H44194M2*(ATPCM(1)*X(I,K)+BTPCM(1)*Y(I,K)) 362 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 363 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 364 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 365 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 366 301 CONTINUE Page 6 Source Listing SPA88 2025-03-12 18:23 SPA88.F 367 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 368 C P(K) (TOPM,TOPPHI) 369 DO 315 I=MYIS,MYIE 370 TOPM(I,1)=PHITMP(I,1) 371 TOPPHI(I,1)=PSITMP(I,1) 372 315 CONTINUE 373 DO 319 K=2,L 374 DO 317 I=MYIS,MYIE 375 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 376 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 377 317 CONTINUE 378 319 CONTINUE 379 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 380 DO 321 K=1,L 381 DO 321 I=MYIS,MYIE 382 FAC1(I,K)=ACOMB(1)*TOPM(I,K) 383 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(1)*TOPPHI(I,K)) 384 TT(I,K)=EXP(HM1EZ*FAC1(I,K)/SQRT(1.+FAC2(I,K))) 385 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 386 321 CONTINUE 387 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 388 DO 353 K=1,L 389 DO 353 I=MYIS,MYIE 390 EXCTS(I,K)=SORC(I,K,1)*(CTMP(I,K+1)-CTMP(I,K)) 391 353 CONTINUE 392 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 393 DO 361 I=MYIS,MYIE 394 GXCTS(I)=CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,1)+ 395 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 396 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 397 3 (SORC(I,LP1,1)-SORC(I,L,1))) 398 361 CONTINUE 399 C 400 C 401 C-----CALCULATION FOR BAND 2 (COMBINED BAND 2) 402 C 403 C 404 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 405 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 406 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 407 DO 401 K=1,L 408 DO 401 I=MYIS,MYIE 409 F(I,K)=H44194M2*(APCM(2)*X(I,K)+BPCM(2)*Y(I,K)) 410 FF(I,K)=H44194M2*(ATPCM(2)*X(I,K)+BTPCM(2)*Y(I,K)) 411 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 412 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 413 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 414 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 415 401 CONTINUE 416 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 417 C P(K) (TOPM,TOPPHI) 418 DO 415 I=MYIS,MYIE 419 TOPM(I,1)=PHITMP(I,1) 420 TOPPHI(I,1)=PSITMP(I,1) 421 415 CONTINUE 422 DO 419 K=2,L 423 DO 417 I=MYIS,MYIE Page 7 Source Listing SPA88 2025-03-12 18:23 SPA88.F 424 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 425 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 426 417 CONTINUE 427 419 CONTINUE 428 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 429 DO 421 K=1,L 430 DO 421 I=MYIS,MYIE 431 FAC1(I,K)=ACOMB(2)*TOPM(I,K) 432 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(2)*TOPPHI(I,K)) 433 TT(I,K)=EXP(HM1EZ*FAC1(I,K)/SQRT(1.+FAC2(I,K))) 434 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 435 421 CONTINUE 436 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 437 DO 453 K=1,L 438 DO 453 I=MYIS,MYIE 439 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,2)* 440 1 (CTMP(I,K+1)-CTMP(I,K)) 441 453 CONTINUE 442 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 443 DO 461 I=MYIS,MYIE 444 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,2)+ 445 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 446 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 447 3 (SORC(I,LP1,2)-SORC(I,L,2))) 448 461 CONTINUE 449 C 450 C-----CALCULATION FOR BAND 3 (COMBINED BAND 3) 451 C 452 C 453 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 454 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 455 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 456 DO 501 K=1,L 457 DO 501 I=MYIS,MYIE 458 F(I,K)=H44194M2*(APCM(3)*X(I,K)+BPCM(3)*Y(I,K)) 459 FF(I,K)=H44194M2*(ATPCM(3)*X(I,K)+BTPCM(3)*Y(I,K)) 460 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 461 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 462 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 463 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 464 501 CONTINUE 465 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 466 C P(K) (TOPM,TOPPHI) 467 DO 515 I=MYIS,MYIE 468 TOPM(I,1)=PHITMP(I,1) 469 TOPPHI(I,1)=PSITMP(I,1) 470 515 CONTINUE 471 DO 519 K=2,L 472 DO 517 I=MYIS,MYIE 473 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 474 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 475 517 CONTINUE 476 519 CONTINUE 477 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 478 DO 521 K=1,L 479 DO 521 I=MYIS,MYIE 480 FAC1(I,K)=ACOMB(3)*TOPM(I,K) Page 8 Source Listing SPA88 2025-03-12 18:23 SPA88.F 481 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(3)*TOPPHI(I,K)) 482 TT(I,K)=EXP(HM1EZ*FAC1(I,K)/SQRT(1.+FAC2(I,K))) 483 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 484 521 CONTINUE 485 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 486 DO 553 K=1,L 487 DO 553 I=MYIS,MYIE 488 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,3)* 489 1 (CTMP(I,K+1)-CTMP(I,K)) 490 553 CONTINUE 491 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 492 DO 561 I=MYIS,MYIE 493 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,3)+ 494 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 495 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 496 3 (SORC(I,LP1,3)-SORC(I,L,3))) 497 561 CONTINUE 498 C 499 C-----CALCULATION FOR BAND 4 (COMBINED BAND 4) 500 C 501 C 502 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 503 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 504 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 505 DO 601 K=1,L 506 DO 601 I=MYIS,MYIE 507 F(I,K)=H44194M2*(APCM(4)*X(I,K)+BPCM(4)*Y(I,K)) 508 FF(I,K)=H44194M2*(ATPCM(4)*X(I,K)+BTPCM(4)*Y(I,K)) 509 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 510 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 511 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 512 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 513 601 CONTINUE 514 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 515 C P(K) (TOPM,TOPPHI) 516 DO 615 I=MYIS,MYIE 517 TOPM(I,1)=PHITMP(I,1) 518 TOPPHI(I,1)=PSITMP(I,1) 519 615 CONTINUE 520 DO 619 K=2,L 521 DO 617 I=MYIS,MYIE 522 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 523 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 524 617 CONTINUE 525 619 CONTINUE 526 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 527 DO 621 K=1,L 528 DO 621 I=MYIS,MYIE 529 FAC1(I,K)=ACOMB(4)*TOPM(I,K) 530 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(4)*TOPPHI(I,K)) 531 TT(I,K)=EXP(HM1EZ*FAC1(I,K)/SQRT(1.+FAC2(I,K))) 532 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 533 621 CONTINUE 534 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 535 DO 653 K=1,L 536 DO 653 I=MYIS,MYIE 537 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,4)* Page 9 Source Listing SPA88 2025-03-12 18:23 SPA88.F 538 1 (CTMP(I,K+1)-CTMP(I,K)) 539 653 CONTINUE 540 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 541 DO 661 I=MYIS,MYIE 542 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,4)+ 543 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 544 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 545 3 (SORC(I,LP1,4)-SORC(I,L,4))) 546 661 CONTINUE 547 C 548 C-----CALCULATION FOR BAND 5 (COMBINED BAND 5) 549 C 550 C 551 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 552 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 553 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 554 DO 701 K=1,L 555 DO 701 I=MYIS,MYIE 556 F(I,K)=H44194M2*(APCM(5)*X(I,K)+BPCM(5)*Y(I,K)) 557 FF(I,K)=H44194M2*(ATPCM(5)*X(I,K)+BTPCM(5)*Y(I,K)) 558 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 559 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 560 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 561 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 562 701 CONTINUE 563 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 564 C P(K) (TOPM,TOPPHI) 565 DO 715 I=MYIS,MYIE 566 TOPM(I,1)=PHITMP(I,1) 567 TOPPHI(I,1)=PSITMP(I,1) 568 715 CONTINUE 569 DO 719 K=2,L 570 DO 717 I=MYIS,MYIE 571 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 572 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 573 717 CONTINUE 574 719 CONTINUE 575 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 576 DO 721 K=1,L 577 DO 721 I=MYIS,MYIE 578 FAC1(I,K)=ACOMB(5)*TOPM(I,K) 579 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(5)*TOPPHI(I,K)) 580 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 581 1 BETACM(5)*TOTVO2(I,K+1)*SKO2D)) 582 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 583 721 CONTINUE 584 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 585 DO 753 K=1,L 586 DO 753 I=MYIS,MYIE 587 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,5)* 588 1 (CTMP(I,K+1)-CTMP(I,K)) 589 753 CONTINUE 590 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 591 DO 761 I=MYIS,MYIE 592 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,5)+ 593 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 594 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * Page 10 Source Listing SPA88 2025-03-12 18:23 SPA88.F 595 3 (SORC(I,LP1,5)-SORC(I,L,5))) 596 761 CONTINUE 597 C 598 C-----CALCULATION FOR BAND 6 (COMBINED BAND 6) 599 C 600 C 601 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 602 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 603 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 604 DO 801 K=1,L 605 DO 801 I=MYIS,MYIE 606 F(I,K)=H44194M2*(APCM(6)*X(I,K)+BPCM(6)*Y(I,K)) 607 FF(I,K)=H44194M2*(ATPCM(6)*X(I,K)+BTPCM(6)*Y(I,K)) 608 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 609 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 610 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 611 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 612 801 CONTINUE 613 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 614 C P(K) (TOPM,TOPPHI) 615 DO 815 I=MYIS,MYIE 616 TOPM(I,1)=PHITMP(I,1) 617 TOPPHI(I,1)=PSITMP(I,1) 618 815 CONTINUE 619 DO 819 K=2,L 620 DO 817 I=MYIS,MYIE 621 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 622 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 623 817 CONTINUE 624 819 CONTINUE 625 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 626 DO 821 K=1,L 627 DO 821 I=MYIS,MYIE 628 FAC1(I,K)=ACOMB(6)*TOPM(I,K) 629 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(6)*TOPPHI(I,K)) 630 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 631 1 BETACM(6)*TOTVO2(I,K+1)*SKO2D)) 632 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 633 821 CONTINUE 634 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 635 DO 853 K=1,L 636 DO 853 I=MYIS,MYIE 637 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,6)* 638 1 (CTMP(I,K+1)-CTMP(I,K)) 639 853 CONTINUE 640 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 641 DO 861 I=MYIS,MYIE 642 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,6)+ 643 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 644 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 645 3 (SORC(I,LP1,6)-SORC(I,L,6))) 646 861 CONTINUE 647 C 648 C-----CALCULATION FOR BAND 7 (COMBINED BAND 7) 649 C 650 C 651 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY Page 11 Source Listing SPA88 2025-03-12 18:23 SPA88.F 652 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 653 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 654 DO 901 K=1,L 655 DO 901 I=MYIS,MYIE 656 F(I,K)=H44194M2*(APCM(7)*X(I,K)+BPCM(7)*Y(I,K)) 657 FF(I,K)=H44194M2*(ATPCM(7)*X(I,K)+BTPCM(7)*Y(I,K)) 658 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 659 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 660 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 661 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 662 901 CONTINUE 663 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 664 C P(K) (TOPM,TOPPHI) 665 DO 915 I=MYIS,MYIE 666 TOPM(I,1)=PHITMP(I,1) 667 TOPPHI(I,1)=PSITMP(I,1) 668 915 CONTINUE 669 DO 919 K=2,L 670 DO 917 I=MYIS,MYIE 671 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 672 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 673 917 CONTINUE 674 919 CONTINUE 675 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 676 DO 921 K=1,L 677 DO 921 I=MYIS,MYIE 678 FAC1(I,K)=ACOMB(7)*TOPM(I,K) 679 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(7)*TOPPHI(I,K)) 680 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 681 1 BETACM(7)*TOTVO2(I,K+1)*SKO2D)) 682 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 683 921 CONTINUE 684 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 685 DO 953 K=1,L 686 DO 953 I=MYIS,MYIE 687 EXCTS(I,K)=EXCTS(I,K)+SORC(I,k,7)* 688 1 (CTMP(I,K+1)-CTMP(I,K)) 689 953 CONTINUE 690 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 691 DO 961 I=MYIS,MYIE 692 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,7)+ 693 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 694 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 695 3 (SORC(I,LP1,7)-SORC(I,L,7))) 696 961 CONTINUE 697 C 698 C-----CALCULATION FOR BAND 8 (COMBINED BAND 8) 699 C 700 C 701 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 702 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 703 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 704 DO 1001 K=1,L 705 DO 1001 I=MYIS,MYIE 706 F(I,K)=H44194M2*(APCM(8)*X(I,K)+BPCM(8)*Y(I,K)) 707 FF(I,K)=H44194M2*(ATPCM(8)*X(I,K)+BTPCM(8)*Y(I,K)) 708 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE Page 12 Source Listing SPA88 2025-03-12 18:23 SPA88.F 709 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 710 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 711 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 712 1001 CONTINUE 713 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 714 C P(K) (TOPM,TOPPHI) 715 DO 1015 I=MYIS,MYIE 716 TOPM(I,1)=PHITMP(I,1) 717 TOPPHI(I,1)=PSITMP(I,1) 718 1015 CONTINUE 719 DO 1019 K=2,L 720 DO 1017 I=MYIS,MYIE 721 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 722 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 723 1017 CONTINUE 724 1019 CONTINUE 725 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 726 DO 1021 K=1,L 727 DO 1021 I=MYIS,MYIE 728 FAC1(I,K)=ACOMB(8)*TOPM(I,K) 729 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(8)*TOPPHI(I,K)) 730 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 731 1 BETACM(8)*TOTVO2(I,K+1)*SKO2D)) 732 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 733 1021 CONTINUE 734 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 735 DO 1053 K=1,L 736 DO 1053 I=MYIS,MYIE 737 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,8)* 738 1 (CTMP(I,K+1)-CTMP(I,K)) 739 1053 CONTINUE 740 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 741 DO 1061 I=MYIS,MYIE 742 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,8)+ 743 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 744 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 745 3 (SORC(I,LP1,8)-SORC(I,L,8))) 746 1061 CONTINUE 747 C 748 C-----CALCULATION FOR BAND 9 ( 560-670 CM-1; INCLUDES CO2) 749 C 750 C 751 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 752 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 753 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 754 DO 1101 K=1,L 755 DO 1101 I=MYIS,MYIE 756 F(I,K)=H44194M2*(APCM(9)*X(I,K)+BPCM(9)*Y(I,K)) 757 FF(I,K)=H44194M2*(ATPCM(9)*X(I,K)+BTPCM(9)*Y(I,K)) 758 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 759 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 760 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 761 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 762 1101 CONTINUE 763 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 764 C P(K) (TOPM,TOPPHI) 765 DO 1115 I=MYIS,MYIE Page 13 Source Listing SPA88 2025-03-12 18:23 SPA88.F 766 TOPM(I,1)=PHITMP(I,1) 767 TOPPHI(I,1)=PSITMP(I,1) 768 1115 CONTINUE 769 DO 1119 K=2,L 770 DO 1117 I=MYIS,MYIE 771 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 772 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 773 1117 CONTINUE 774 1119 CONTINUE 775 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 776 DO 1121 K=1,L 777 DO 1121 I=MYIS,MYIE 778 FAC1(I,K)=ACOMB(9)*TOPM(I,K) 779 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(9)*TOPPHI(I,K)) 780 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 781 1 BETACM(9)*TOTVO2(I,K+1)*SKO2D))*CO2SP1(I,K+1) 782 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 783 1121 CONTINUE 784 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 785 DO 1153 K=1,L 786 DO 1153 I=MYIS,MYIE 787 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,9)* 788 1 (CTMP(I,K+1)-CTMP(I,K)) 789 1153 CONTINUE 790 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 791 DO 1161 I=MYIS,MYIE 792 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,9)+ 793 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 794 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 795 3 (SORC(I,LP1,9)-SORC(I,L,9))) 796 1161 CONTINUE 797 C 798 C-----CALCULATION FOR BAND 10 (670-800 CM-1; INCLUDES CO2) 799 C 800 C 801 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 802 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 803 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 804 DO 1201 K=1,L 805 DO 1201 I=MYIS,MYIE 806 F(I,K)=H44194M2*(APCM(10)*X(I,K)+BPCM(10)*Y(I,K)) 807 FF(I,K)=H44194M2*(ATPCM(10)*X(I,K)+BTPCM(10)*Y(I,K)) 808 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 809 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 810 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 811 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 812 1201 CONTINUE 813 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 814 C P(K) (TOPM,TOPPHI) 815 DO 1215 I=MYIS,MYIE 816 TOPM(I,1)=PHITMP(I,1) 817 TOPPHI(I,1)=PSITMP(I,1) 818 1215 CONTINUE 819 DO 1219 K=2,L 820 DO 1217 I=MYIS,MYIE 821 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 822 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) Page 14 Source Listing SPA88 2025-03-12 18:23 SPA88.F 823 1217 CONTINUE 824 1219 CONTINUE 825 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 826 DO 1221 K=1,L 827 DO 1221 I=MYIS,MYIE 828 FAC1(I,K)=ACOMB(10)*TOPM(I,K) 829 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(10)*TOPPHI(I,K)) 830 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 831 1 BETACM(10)*TOTVO2(I,K+1)*SKO2D))*CO2SP2(I,K+1) 832 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 833 1221 CONTINUE 834 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 835 DO 1253 K=1,L 836 DO 1253 I=MYIS,MYIE 837 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,10)* 838 1 (CTMP(I,K+1)-CTMP(I,K)) 839 1253 CONTINUE 840 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 841 DO 1261 I=MYIS,MYIE 842 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,10)+ 843 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 844 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 845 3 (SORC(I,LP1,10)-SORC(I,L,10))) 846 1261 CONTINUE 847 C 848 C-----CALCULATION FOR BAND 11 (800-900 CM-1) 849 C 850 C 851 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 852 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 853 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 854 DO 1301 K=1,L 855 DO 1301 I=MYIS,MYIE 856 F(I,K)=H44194M2*(APCM(11)*X(I,K)+BPCM(11)*Y(I,K)) 857 FF(I,K)=H44194M2*(ATPCM(11)*X(I,K)+BTPCM(11)*Y(I,K)) 858 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 859 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 860 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 861 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 862 1301 CONTINUE 863 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 864 C P(K) (TOPM,TOPPHI) 865 DO 1315 I=MYIS,MYIE 866 TOPM(I,1)=PHITMP(I,1) 867 TOPPHI(I,1)=PSITMP(I,1) 868 1315 CONTINUE 869 DO 1319 K=2,L 870 DO 1317 I=MYIS,MYIE 871 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 872 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 873 1317 CONTINUE 874 1319 CONTINUE 875 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 876 DO 1321 K=1,L 877 DO 1321 I=MYIS,MYIE 878 FAC1(I,K)=ACOMB(11)*TOPM(I,K) 879 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(11)*TOPPHI(I,K)) Page 15 Source Listing SPA88 2025-03-12 18:23 SPA88.F 880 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 881 1 BETACM(11)*TOTVO2(I,K+1)*SKO2D)) 882 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 883 1321 CONTINUE 884 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 885 DO 1353 K=1,L 886 DO 1353 I=MYIS,MYIE 887 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,11)* 888 1 (CTMP(I,K+1)-CTMP(I,K)) 889 1353 CONTINUE 890 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 891 DO 1361 I=MYIS,MYIE 892 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,11)+ 893 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 894 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 895 3 (SORC(I,LP1,11)-SORC(I,L,11))) 896 1361 CONTINUE 897 C 898 C-----CALCULATION FOR BAND 12 (900-990 CM-1) 899 C 900 C 901 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 902 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 903 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 904 DO 1401 K=1,L 905 DO 1401 I=MYIS,MYIE 906 F(I,K)=H44194M2*(APCM(12)*X(I,K)+BPCM(12)*Y(I,K)) 907 FF(I,K)=H44194M2*(ATPCM(12)*X(I,K)+BTPCM(12)*Y(I,K)) 908 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 909 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 910 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 911 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 912 1401 CONTINUE 913 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 914 C P(K) (TOPM,TOPPHI) 915 DO 1415 I=MYIS,MYIE 916 TOPM(I,1)=PHITMP(I,1) 917 TOPPHI(I,1)=PSITMP(I,1) 918 1415 CONTINUE 919 DO 1419 K=2,L 920 DO 1417 I=MYIS,MYIE 921 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 922 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 923 1417 CONTINUE 924 1419 CONTINUE 925 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 926 DO 1421 K=1,L 927 DO 1421 I=MYIS,MYIE 928 FAC1(I,K)=ACOMB(12)*TOPM(I,K) 929 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(12)*TOPPHI(I,K)) 930 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 931 1 BETACM(12)*TOTVO2(I,K+1)*SKO2D)) 932 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 933 1421 CONTINUE 934 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 935 DO 1453 K=1,L 936 DO 1453 I=MYIS,MYIE Page 16 Source Listing SPA88 2025-03-12 18:23 SPA88.F 937 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,12)* 938 1 (CTMP(I,K+1)-CTMP(I,K)) 939 1453 CONTINUE 940 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 941 DO 1461 I=MYIS,MYIE 942 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,12)+ 943 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 944 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 945 3 (SORC(I,LP1,12)-SORC(I,L,12))) 946 1461 CONTINUE 947 C 948 C-----CALCULATION FOR BAND 13 (990-1070 CM-1; INCLUDES O3)) 949 C 950 C 951 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 952 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 953 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 954 DO 1501 K=1,L 955 DO 1501 I=MYIS,MYIE 956 F(I,K)=H44194M2*(APCM(13)*X(I,K)+BPCM(13)*Y(I,K)) 957 FF(I,K)=H44194M2*(ATPCM(13)*X(I,K)+BTPCM(13)*Y(I,K)) 958 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 959 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 960 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 961 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 962 1501 CONTINUE 963 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 964 C P(K) (TOPM,TOPPHI) 965 DO 1515 I=MYIS,MYIE 966 TOPM(I,1)=PHITMP(I,1) 967 TOPPHI(I,1)=PSITMP(I,1) 968 1515 CONTINUE 969 DO 1519 K=2,L 970 DO 1517 I=MYIS,MYIE 971 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 972 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 973 1517 CONTINUE 974 1519 CONTINUE 975 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 976 DO 1521 K=1,L 977 DO 1521 I=MYIS,MYIE 978 FAC1(I,K)=ACOMB(13)*TOPM(I,K) 979 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(13)*TOPPHI(I,K)) 980 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 981 1 BETACM(13)*TOTVO2(I,K+1)*SKO2D+TO3SPC(I,K))) 982 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 983 1521 CONTINUE 984 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 985 DO 1553 K=1,L 986 DO 1553 I=MYIS,MYIE 987 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,13)* 988 1 (CTMP(I,K+1)-CTMP(I,K)) 989 1553 CONTINUE 990 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 991 DO 1561 I=MYIS,MYIE 992 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,13)+ 993 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + Page 17 Source Listing SPA88 2025-03-12 18:23 SPA88.F 994 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 995 3 (SORC(I,LP1,13)-SORC(I,L,13))) 996 1561 CONTINUE 997 C 998 C-----CALCULATION FOR BAND 14 (1070-1200 CM-1) 999 C 1000 C 1001 C---OBTAIN TEMPERATURE CORRECTION (CAPPHI,CAPPSI),THEN MULTIPLY 1002 C BY OPTICAL PATH (VAR1,VAR2) TO COMPUTE TEMPERATURE-CORRECTED 1003 C OPTICAL PATH AND MEAN PRESSURE FOR A LAYER (PHITMP,PSITMP) 1004 DO 1601 K=1,L 1005 DO 1601 I=MYIS,MYIE 1006 F(I,K)=H44194M2*(APCM(14)*X(I,K)+BPCM(14)*Y(I,K)) 1007 FF(I,K)=H44194M2*(ATPCM(14)*X(I,K)+BTPCM(14)*Y(I,K)) 1008 AG(I,K)=(H1P41819+F(I,K))*F(I,K)+ONE 1009 AGG(I,K)=(H1P41819+FF(I,K))*FF(I,K)+ONE 1010 PHITMP(I,K)=VAR1(I,K)*(((( AG(I,K)*AG(I,K))**2)**2)**2) 1011 PSITMP(I,K)=VAR2(I,K)*(((( AGG(I,K)*AGG(I,K))**2)**2)**2) 1012 1601 CONTINUE 1013 C---OBTAIN OPTICAL PATH,MEAN PRESSURE FROM THE TOP TO THE PRESSURE 1014 C P(K) (TOPM,TOPPHI) 1015 DO 1615 I=MYIS,MYIE 1016 TOPM(I,1)=PHITMP(I,1) 1017 TOPPHI(I,1)=PSITMP(I,1) 1018 1615 CONTINUE 1019 DO 1619 K=2,L 1020 DO 1617 I=MYIS,MYIE 1021 TOPM(I,K)=TOPM(I,K-1)+PHITMP(I,K) 1022 TOPPHI(I,K)=TOPPHI(I,K-1)+PSITMP(I,K) 1023 1617 CONTINUE 1024 1619 CONTINUE 1025 C---TT IS THE CLOUD-FREE CTS TRANSMISSION FUNCTION 1026 DO 1621 K=1,L 1027 DO 1621 I=MYIS,MYIE 1028 FAC1(I,K)=ACOMB(14)*TOPM(I,K) 1029 FAC2(I,K)=FAC1(I,K)*TOPM(I,K)/(BCOMB(14)*TOPPHI(I,K)) 1030 TT(I,K)=EXP(HM1EZ*(FAC1(I,K)/SQRT(ONE+FAC2(I,K))+ 1031 1 BETACM(14)*TOTVO2(I,K+1)*SKO2D)) 1032 CTMP(I,K+1)=TT(I,K)*CLDFAC(I,K+1,1) 1033 1621 CONTINUE 1034 C---EXCTS IS THE CTS COOLING RATE ACCUMULATED OVER FREQUENCY BANDS 1035 DO 1653 K=1,L 1036 DO 1653 I=MYIS,MYIE 1037 EXCTS(I,K)=EXCTS(I,K)+SORC(I,K,14)* 1038 1 (CTMP(I,K+1)-CTMP(I,K)) 1039 1653 CONTINUE 1040 C---GXCTS IS THE EXACT CTS TOP FLUX ACCUMULATED OVER FREQUENCY BANDS 1041 DO 1661 I=MYIS,MYIE 1042 GXCTS(I)=GXCTS(I)+CLDFAC(I,LP1,1)*(TT(I,L)*SORC(I,L,14)+ 1043 1 (HAF*DELP(I,L)*(TT(I,LM1)*(P(I,LP1)-PRESS(I,L)) + 1044 2 TT(I,L)*(P(I,LP1)+PRESS(I,L)-TWO*P(I,L)))) * 1045 3 (SORC(I,LP1,14)-SORC(I,L,14))) 1046 1661 CONTINUE 1047 C 1048 C 1049 C OBTAIN CTS FLUX AT THE TOP BY INTEGRATION OF HEATING RATES AND 1050 C USING CTS FLUX AT THE BOTTOM (CURRENT VALUE OF GXCTS). NOTE Page 18 Source Listing SPA88 2025-03-12 18:23 SPA88.F 1051 C THAT THE PRESSURE QUANTITIES AND CONVERSION FACTORS HAVE NOT 1052 C BEEN INCLUDED EITHER IN EXCTS OR IN GXCTS. THESE CANCEL OUT, THUS 1053 C REDUCING COMPUTATIONS! 1054 DO 1731 K=1,L 1055 DO 1731 I=MYIS,MYIE 1056 GXCTS(I)=GXCTS(I)-EXCTS(I,K) 1057 1731 CONTINUE 1058 C 1059 C NOW SCALE THE COOLING RATE (EXCTS) BY INCLUDING THE PRESSURE 1060 C FACTOR (DELP) AND THE CONVERSION FACTOR (RADCON) 1061 DO 1741 K=1,L 1062 DO 1741 I=MYIS,MYIE 1063 EXCTS(I,K)=EXCTS(I,K)*RADCON*DELP(I,K) 1064 1741 CONTINUE 1065 C---THIS IS THE END OF THE EXACT CTS COMPUTATIONS; AT THIS POINT 1066 C EXCTS HAS ITS APPROPRIATE VALUE. 1067 C 1068 C*** COMPUTE APPROXIMATE CTS HEATING RATES FOR 15UM AND 9.6 UM BANDS 1069 C (CTSO3) 1070 DO 1711 K=1,L 1071 DO 1711 I=MYIS,MYIE 1072 CTMP2(I,K+1)=CO2SP(I,K+1)*CLDFAC(I,K+1,1) 1073 CTMP3(I,K+1)=TO3SP(I,K)*CLDFAC(I,K+1,1) 1074 1711 CONTINUE 1075 DO 1701 K=1,L 1076 DO 1701 I=MYIS,MYIE 1077 CTSO3(I,K)=RADCON*DELP(I,K)* 1078 1 (CSOUR(I,K)*(CTMP2(I,K+1)-CTMP2(I,K)) + 1079 2 SORC(I,K,13)*(CTMP3(I,K+1)-CTMP3(I,K))) 1080 1701 CONTINUE 1081 RETURN 1082 END Page 19 Source Listing SPA88 2025-03-12 18:23 Entry Points SPA88.F ENTRY POINTS Name spa88_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1001 Label 712 704,705 101 Label 343 339,340 1015 Label 718 715 1017 Label 723 720 1019 Label 724 719 1021 Label 733 726,727 1053 Label 739 735,736 1061 Label 746 741 1101 Label 762 754,755 1115 Label 768 765 1117 Label 773 770 1119 Label 774 769 1121 Label 783 776,777 1153 Label 789 785,786 1161 Label 796 791 1201 Label 812 804,805 1215 Label 818 815 1217 Label 823 820 1219 Label 824 819 1221 Label 833 826,827 1253 Label 839 835,836 1261 Label 846 841 1301 Label 862 854,855 1315 Label 868 865 1317 Label 873 870 1319 Label 874 869 1321 Label 883 876,877 1353 Label 889 885,886 1361 Label 896 891 1401 Label 912 904,905 1415 Label 918 915 1417 Label 923 920 1419 Label 924 919 1421 Label 933 926,927 1453 Label 939 935,936 1461 Label 946 941 1501 Label 962 954,955 1515 Label 968 965 1517 Label 973 970 1519 Label 974 969 1521 Label 983 976,977 1553 Label 989 985,986 1561 Label 996 991 1601 Label 1012 1004,1005 Page 20 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 1615 Label 1018 1015 1617 Label 1023 1020 1619 Label 1024 1019 1621 Label 1033 1026,1027 1653 Label 1039 1035,1036 1661 Label 1046 1041 1701 Label 1080 1075,1076 1711 Label 1074 1070,1071 1731 Label 1057 1054,1055 1741 Label 1064 1061,1062 301 Label 366 358,359 315 Label 372 369 317 Label 377 374 319 Label 378 373 321 Label 386 380,381 345 Label 350 346 353 Label 391 388,389 361 Label 398 393 401 Label 415 407,408 415 Label 421 418 417 Label 426 423 419 Label 427 422 421 Label 435 429,430 453 Label 441 437,438 461 Label 448 443 501 Label 464 456,457 515 Label 470 467 517 Label 475 472 519 Label 476 471 521 Label 484 478,479 553 Label 490 486,487 561 Label 497 492 601 Label 513 505,506 615 Label 519 516 617 Label 524 521 619 Label 525 520 621 Label 533 527,528 653 Label 539 535,536 661 Label 546 541 701 Label 562 554,555 715 Label 568 565 717 Label 573 570 719 Label 574 569 721 Label 583 576,577 753 Label 589 585,586 761 Label 596 591 801 Label 612 604,605 815 Label 618 615 817 Label 623 620 819 Label 624 619 821 Label 633 626,627 853 Label 639 635,636 861 Label 646 641 901 Label 662 654,655 915 Label 668 665 Page 21 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 917 Label 673 670 919 Label 674 669 921 Label 683 676,677 953 Label 689 685,686 961 Label 696 691 AG Local 334 R(4) 4 2 3195 362,364,411,413,460,462,509,511,55 8,560,608,610,658,660,708,710,758, 760,808,810,858,860,908,910,958,96 0,1008,1010 AGG Local 334 R(4) 4 2 3195 363,365,412,414,461,463,510,512,55 9,561,609,611,659,661,709,711,759, 761,809,811,859,861,909,911,959,96 1,1009,1011 BANDTA Common 240 5900 BDCOMB Common 308 736 BDWIDE Common 268 56 CLDFAC Dummy 23 R(4) 4 3 150236 ARG,INOUT 385,394,434,444,483,493,532,542,58 2,592,632,642,682,692,732,742,782, 792,832,842,882,892,932,942,982,99 2,1032,1042,1072,1073 CO2SP Dummy 25 R(4) 4 2 3266 ARG,INOUT 1072 CO2SP1 Dummy 25 R(4) 4 2 3266 ARG,INOUT 781 CO2SP2 Dummy 25 R(4) 4 2 3266 ARG,INOUT 831 CSOUR Dummy 22 R(4) 4 2 3266 ARG,INOUT 1078 CTMP Local 329 R(4) 4 2 3266 347,385,390,434,440,483,489,532,53 8,582,588,632,638,682,688,732,738, 782,788,832,838,882,888,932,938,98 2,988,1032,1038 CTMP2 Local 332 R(4) 4 2 3266 348,1072,1078 CTMP3 Local 332 R(4) 4 2 3266 349,1073,1079 CTSO3 Dummy 22 R(4) 4 2 3195 ARG,INOUT 1077 DELP Dummy 24 R(4) 4 2 3195 ARG,INOUT 395,445,494,543,593,643,693,743,79 3,843,893,943,993,1043,1063,1077 DELP2 Dummy 24 R(4) 4 2 3195 ARG,INOUT EXCTS Dummy 22 R(4) 4 2 3195 ARG,INOUT 390,439,488,537,587,637,687,737,78 7,837,887,937,987,1037,1056,1063 EXP Func 384 scalar 384,433,482,531,580,630,680,730,78 0,830,880,930,980,1030 F Local 333 R(4) 4 2 3195 360,362,409,411,458,460,507,509,55 6,558,606,608,656,658,706,708,756, 758,806,808,856,858,906,908,956,95 8,1006,1008 FAC1 Local 328 R(4) 4 2 3195 382,383,384,431,432,433,480,481,48 2,529,530,531,578,579,580,628,629, 630,678,679,680,728,729,730,778,77 9,780,828,829,830,878,879,880,928, 929,930,978,979,980,1028,1029,1030 FAC2 Local 328 R(4) 4 2 3195 383,384,432,433,481,482,530,531,57 9,580,629,630,679,680,729,730,779, 780,829,830,879,880,929,930,979,98 0,1029,1030 FF Local 333 R(4) 4 2 3195 361,363,410,412,459,461,508,510,55 7,559,607,609,657,659,707,709,757, 759,807,809,857,859,907,909,957,95 9,1007,1009 Page 22 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 154 336 SAVE GXCTS Dummy 22 R(4) 4 1 71 ARG,INOUT 394,444,493,542,592,642,692,742,79 2,842,892,942,992,1042,1056 HCON Common 33 872 SAVE I Local 340 I(4) 4 scalar 340,341,342,346,347,348,349,359,36 0,361,362,363,364,365,369,370,371, 374,375,376,381,382,383,384,385,38 9,390,393,394,395,396,397,408,409, 410,411,412,413,414,418,419,420,42 3,424,425,430,431,432,433,434,438, 439,440,443,444,445,446,447,457,45 8,459,460,461,462,463,467,468,469, 472,473,474,479,480,481,482,483,48 7,488,489,492,493,494,495,496,506, 507,508,509,510,511,512,516,517,51 8,521,522,523,528,529,530,531,532, 536,537,538,541,542,543,544,545,55 5,556,557,558,559,560,561,565,566, 567,570,571,572,577,578,579,580,58 1,582,586,587,588,591,592,593,594, 595,605,606,607,608,609,610,611,61 5,616,617,620,621,622,627,628,629, 630,631,632,636,637,638,641,642,64 3,644,645,655,656,657,658,659,660, 661,665,666,667,670,671,672,677,67 8,679,680,681,682,686,687,688,691, 692,693,694,695,705,706,707,708,70 9,710,711,715,716,717,720,721,722, 727,728,729,730,731,732,736,737,73 8,741,742,743,744,745,755,756,757, 758,759,760,761,765,766,767,770,77 1,772,777,778,779,780,781,782,786, 787,788,791,792,793,794,795,805,80 6,807,808,809,810,811,815,816,817, 820,821,822,827,828,829,830,831,83 2,836,837,838,841,842,843,844,845, 855,856,857,858,859,860,861,865,86 6,867,870,871,872,877,878,879,880, 881,882,886,887,888,891,892,893,89 4,895,905,906,907,908,909,910,911, 915,916,917,920,921,922,927,928,92 9,930,931,932,936,937,938,941,942, 943,944,945,955,956,957,958,959,96 0,961,965,966,967,970,971,972,977, 978,979,980,981,982,986,987,988,99 1,992,993,994,995,1005,1006,1007,1 008,1009,1010,1011,1015,1016,1017, 1020,1021,1022,1027,1028,1029,1030 ,1031,1032,1036,1037,1038,1041,104 2,1043,1044,1045,1055,1056,1062,10 63,1071,1072,1073,1076,1077,1078,1 079 IDIM1 Param 95 I(4) 4 scalar 313,314,315,316,317,318,319,320,32 1,322,323,324,326,327,328,329,330, 331,332,333,334 Page 23 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References IDIM2 Param 95 I(4) 4 scalar 313,314,315,316,317,318,319,320,32 1,322,323,324,326,327,328,329,330, 331,332,333,334 IGSTL Param 90 I(4) 4 scalar 95 IGSTR Param 90 I(4) 4 scalar 95 IM Param 78 I(4) 4 scalar 92,95,160,161,162,163,164,168,169, 174,177 IMAX Param 177 I(4) 4 scalar 177,190 INLTE Param 189 I(4) 4 scalar 189 INLTEP Param 189 I(4) 4 scalar INPES Param 87 I(4) 4 scalar 92,95,144,145,146,155,156 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,160,161,162,163,164,168,169, 174 JNPES Param 87 I(4) 4 scalar 93,96,144,145,146,155,156 JTAIL Param 93 I(4) 4 scalar K Local 339 I(4) 4 scalar 339,341,342,358,360,361,362,363,36 4,365,373,375,376,380,382,383,384, 385,388,390,407,409,410,411,412,41 3,414,422,424,425,429,431,432,433, 434,437,439,440,456,458,459,460,46 1,462,463,471,473,474,478,480,481, 482,483,486,488,489,505,507,508,50 9,510,511,512,520,522,523,527,529, 530,531,532,535,537,538,554,556,55 7,558,559,560,561,569,571,572,576, 578,579,580,581,582,585,587,588,60 4,606,607,608,609,610,611,619,621, 622,626,628,629,630,631,632,635,63 7,638,654,656,657,658,659,660,661, 669,671,672,676,678,679,680,681,68 2,685,687,688,704,706,707,708,709, 710,711,719,721,722,726,728,729,73 0,731,732,735,737,738,754,756,757, 758,759,760,761,769,771,772,776,77 8,779,780,781,782,785,787,788,804, 806,807,808,809,810,811,819,821,82 2,826,828,829,830,831,832,835,837, 838,854,856,857,858,859,860,861,86 9,871,872,876,878,879,880,881,882, 885,887,888,904,906,907,908,909,91 0,911,919,921,922,926,928,929,930, 931,932,935,937,938,954,956,957,95 8,959,960,961,969,971,972,976,978, 979,980,981,982,985,987,988,1004,1 006,1007,1008,1009,1010,1011,1019, 1021,1022,1026,1028,1029,1030,1031 ,1032,1035,1037,1038,1054,1056,106 1,1063,1070,1072,1073,1075,1077,10 78,1079 KO2 Param 192 I(4) 4 scalar 193 Page 24 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References KO21 Param 193 I(4) 4 scalar KO2M Param 193 I(4) 4 scalar L Param 176 I(4) 4 scalar 180,181,182,185,187,316,317,318,31 9,323,326,327,328,329,330,331,333, 334,339,358,373,380,388,394,395,39 6,397,407,422,429,437,444,445,446, 447,456,471,478,486,493,494,495,49 6,505,520,527,535,542,543,544,545, 554,569,576,585,592,593,594,595,60 4,619,626,635,642,643,644,645,654, 669,676,685,692,693,694,695,704,71 9,726,735,742,743,744,745,754,769, 776,785,792,793,794,795,804,819,82 6,835,842,843,844,845,854,869,876, 885,892,893,894,895,904,919,926,93 5,942,943,944,945,954,969,976,985, 992,993,994,995,1004,1019,1026,103 5,1042,1043,1044,1045,1054,1061,10 70,1075 LL Param 182 I(4) 4 scalar 182,183 LL3P Param 187 I(4) 4 scalar 190 LL3PI Param 190 I(4) 4 scalar LLM1 Param 183 I(4) 4 scalar LLM2 Param 183 I(4) 4 scalar LLM3 Param 183 I(4) 4 scalar LLP1 Param 182 I(4) 4 scalar 190 LLP1I Param 190 I(4) 4 scalar LLP2 Param 182 I(4) 4 scalar LLP3 Param 182 I(4) 4 scalar LM Param 78 I(4) 4 scalar 169,176 LM1 Param 181 I(4) 4 scalar 395,445,494,543,593,643,693,743,79 3,843,893,943,993,1043 LM2 Param 181 I(4) 4 scalar LM3 Param 181 I(4) 4 scalar LP1 Param 180 I(4) 4 scalar 184,185,186,190,313,314,315,317,31 9,320,321,322,329,332,394,395,396, 397,444,445,446,447,493,494,495,49 6,542,543,544,545,592,593,594,595, 642,643,644,645,692,693,694,695,74 2,743,744,745,792,793,794,795,842, 843,844,845,892,893,894,895,942,94 3,944,945,992,993,994,995,1042,104 3,1044,1045 LP121 Param 186 I(4) 4 scalar LP1I Param 190 I(4) 4 scalar LP1M Param 184 I(4) 4 scalar 184 LP1M1 Param 184 I(4) 4 scalar LP1V Param 185 I(4) 4 scalar LP2 Param 180 I(4) 4 scalar LP3 Param 180 I(4) 4 scalar LSM Param 78 I(4) 4 scalar MAPPINGS Common 173 5024 SAVE MPPCOM Common 118 1724 SAVE NB Param 188 I(4) 4 scalar 191 NB1 Param 191 I(4) 4 scalar Page 25 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References NBLM Param 179 I(4) 4 scalar NBLW Param 178 I(4) 4 scalar 240,241,242 NBLX Param 178 I(4) 4 scalar NBLY Param 178 I(4) 4 scalar 179,186,308,309,310,313 NCOL Param 177 I(4) 4 scalar NNLTE Param 189 I(4) 4 scalar P Dummy 24 R(4) 4 2 3266 ARG,INOUT 395,396,445,446,494,495,543,544,59 3,594,643,644,693,694,743,744,793, 794,843,844,893,894,943,944,993,99 4,1043,1044 PHITMP Local 326 R(4) 4 2 3195 364,370,375,413,419,424,462,468,47 3,511,517,522,560,566,571,610,616, 621,660,666,671,710,716,721,760,76 6,771,810,816,821,860,866,871,910, 916,921,960,966,971,1010,1016,1021 PHYCON Common 27 80 SAVE PRESS Dummy 23 R(4) 4 2 3266 ARG,INOUT 395,396,445,446,494,495,543,544,59 3,594,643,644,693,694,743,744,793, 794,843,844,893,894,943,944,993,99 4,1043,1044 PSITMP Local 326 R(4) 4 2 3195 365,371,376,414,420,425,463,469,47 4,512,518,523,561,567,572,611,617, 622,661,667,672,711,717,722,761,76 7,772,811,817,822,861,867,872,911, 917,922,961,967,972,1011,1017,1022 SORC Dummy 22 R(4) 4 3 48990 ARG,INOUT 390,394,397,439,444,447,488,493,49 6,537,542,545,587,592,595,637,642, 645,687,692,695,737,742,745,787,79 2,795,837,842,845,887,892,895,937, 942,945,987,992,995,1037,1042,1045 ,1079 SPA88 Subr 22 SQRT Func 384 scalar 384,433,482,531,580,630,680,730,78 0,830,880,930,980,1030 TEMP Dummy 23 R(4) 4 2 3266 ARG,INOUT 341 TEMPCOM Common 159 6603768 SAVE TO3SP Dummy 24 R(4) 4 2 3266 ARG,INOUT 1073 TO3SPC Dummy 24 R(4) 4 2 3195 ARG,INOUT 981 TOPM Local 331 R(4) 4 2 3195 370,375,382,383,419,424,431,432,46 8,473,480,481,517,522,529,530,566, 571,578,579,616,621,628,629,666,67 1,678,679,716,721,728,729,766,771, 778,779,816,821,828,829,866,871,87 8,879,916,921,928,929,966,971,978, 979,1016,1021,1028,1029 TOPO Common 167 17478548 SAVE TOPPHI Local 331 R(4) 4 2 3195 371,376,383,420,425,432,469,474,48 1,518,523,530,567,572,579,617,622, 629,667,672,679,717,722,729,767,77 2,779,817,822,829,867,872,879,917, 922,929,967,972,979,1017,1022,1029 TOTVO2 Dummy 24 R(4) 4 2 3266 ARG,INOUT 581,631,681,731,781,831,881,931,98 1,1031 TT Local 327 R(4) 4 2 3195 384,385,394,395,396,433,434,444,44 5,446,482,483,493,494,495,531,532, Page 26 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 542,543,544,580,582,592,593,594,63 0,632,642,643,644,680,682,692,693, 694,730,732,742,743,744,780,782,79 2,793,794,830,832,842,843,844,880, 882,892,893,894,930,932,942,943,94 4,980,982,992,993,994,1030,1032,10 42,1043,1044 VAR1 Dummy 23 R(4) 4 2 3195 ARG,INOUT 364,413,462,511,560,610,660,710,76 0,810,860,910,960,1010 VAR2 Dummy 23 R(4) 4 2 3195 ARG,INOUT 365,414,463,512,561,611,661,711,76 1,811,861,911,961,1011 X Local 329 R(4) 4 2 3195 341,342,360,361,409,410,458,459,50 7,508,556,557,606,607,656,657,706, 707,756,757,806,807,856,857,906,90 7,956,957,1006,1007 Y Local 330 R(4) 4 2 3195 342,360,361,409,410,458,459,507,50 8,556,557,606,607,656,657,706,707, 756,757,806,807,856,857,906,907,95 6,957,1006,1007 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 382,431,480,529,578,628,678,728,77 8,828,878,928,978,1028 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 360,409,458,507,556,606,656,706,75 6,806,856,906,956,1006 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 361,410,459,508,557,607,657,707,75 7,807,857,907,957,1007 ATPWD R(4) 4 20 scalar COM 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 383,432,481,530,579,629,679,729,77 9,829,879,929,979,1029 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 581,631,681,731,781,831,881,931,98 1,1031 BETAD R(4) 4 1304 1 163 COM BETAWD R(4) 4 8 scalar COM Page 27 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 BP R(4) 4 2608 1 163 COM BPCM R(4) 4 400 1 15 COM 360,409,458,507,556,606,656,706,75 6,806,856,906,956,1006 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 361,410,459,508,557,607,657,707,75 7,807,857,907,957,1007 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 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 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 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 Page 28 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 362,363,411,412,460,461,509,510,55 8,559,608,609,658,659,708,709,758, 759,808,809,858,859,908,909,958,95 9,1008,1009 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 341 H25M31 R(4) 4 552 scalar COM Page 29 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 H42M11 R(4) 4 336 scalar COM H42M2 R(4) 4 764 scalar COM H44194M2 R(4) 4 700 scalar COM 360,361,409,410,458,459,507,508,55 6,557,606,607,656,657,706,707,756, 757,806,807,856,857,906,907,956,95 7,1006,1007 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 Page 30 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 395,445,494,543,593,643,693,743,79 3,843,893,943,993,1043 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 384,433,482,531,580,630,680,730,78 0,830,880,930,980,1030 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 Page 31 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 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 Page 32 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 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 340,346,359,369,374,381,389,393,40 8,418,423,430,438,443,457,467,472, 479,487,492,506,516,521,528,536,54 1,555,565,570,577,586,591,605,615, 620,627,636,641,655,665,670,677,68 6,691,705,715,720,727,736,741,755, 765,770,777,786,791,805,815,820,82 7,836,841,855,865,870,877,886,891, 905,915,920,927,936,941,955,965,97 0,977,986,991,1005,1015,1020,1027, 1036,1041,1055,1062,1071,1076 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 340,346,359,369,374,381,389,393,40 8,418,423,430,438,443,457,467,472, 479,487,492,506,516,521,528,536,54 1,555,565,570,577,586,591,605,615, 620,627,636,641,655,665,670,677,68 6,691,705,715,720,727,736,741,755, 765,770,777,786,791,805,815,820,82 7,836,841,855,865,870,877,886,891, 905,915,920,927,936,941,955,965,97 0,977,986,991,1005,1015,1020,1027, 1036,1041,1055,1062,1071,1076 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 Page 33 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 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 Page 34 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References 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 347,362,363,411,412,460,461,509,51 0,558,559,580,608,609,630,658,659, 680,708,709,730,758,759,780,808,80 9,830,858,859,880,908,909,930,958, 959,980,1008,1009,1030 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 RADCON R(4) 4 40 scalar COM 1063,1077 RADCON1 R(4) 4 64 scalar COM 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 SKO2D R(4) 4 44 scalar COM 581,631,681,731,781,831,881,931,98 1,1031 SKO3R R(4) 4 52 scalar 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 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 Page 35 Source Listing SPA88 2025-03-12 18:23 Symbol Table SPA88.F Name Type Bytes Offset Dimen Elements Attributes References TTVG R(4) 4 371884 2 92971 COM TWO R(4) 4 36 scalar COM 396,446,495,544,594,644,694,744,79 4,844,894,944,994,1044 ZERO R(4) 4 52 scalar COM Page 36 Source Listing SPA88 2025-03-12 18:23 Subprograms/Common Blocks SPA88.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BANDTA Common 240 5900 BDCOMB Common 308 736 BDWIDE Common 268 56 GLB_TABLE Common 154 336 SAVE HCON Common 33 872 SAVE MAPPINGS Common 173 5024 SAVE MPPCOM Common 118 1724 SAVE PHYCON Common 27 80 SAVE SPA88 Subr 22 TEMPCOM Common 159 6603768 SAVE TOPO Common 167 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 -check nouninitialized -check noudio_iostat Page 37 Source Listing SPA88 2025-03-12 18:23 SPA88.F -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 -init nominus_huge -init nominus_infinity Page 38 Source Listing SPA88 2025-03-12 18:23 SPA88.F -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 : SPA88.lst no -o COMPILER: Intel(R) Fortran 19.1-1655