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