Page 1 Source Listing SPA88 2014-12-17 20:47 /tmpnwprd/ifortgFBDlL.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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 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 2014-12-17 20:47 Entry Points SPA88.F ENTRY POINTS Name spa88_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1001 Label 627 619,620 101 Label 258 254,255 1015 Label 633 630 1017 Label 638 635 1019 Label 639 634 1021 Label 648 641,642 1053 Label 654 650,651 1061 Label 661 656 1101 Label 677 669,670 1115 Label 683 680 1117 Label 688 685 1119 Label 689 684 1121 Label 698 691,692 1153 Label 704 700,701 1161 Label 711 706 1201 Label 727 719,720 1215 Label 733 730 1217 Label 738 735 1219 Label 739 734 1221 Label 748 741,742 1253 Label 754 750,751 1261 Label 761 756 1301 Label 777 769,770 1315 Label 783 780 1317 Label 788 785 1319 Label 789 784 1321 Label 798 791,792 1353 Label 804 800,801 1361 Label 811 806 1401 Label 827 819,820 1415 Label 833 830 1417 Label 838 835 1419 Label 839 834 1421 Label 848 841,842 1453 Label 854 850,851 1461 Label 861 856 1501 Label 877 869,870 1515 Label 883 880 1517 Label 888 885 1519 Label 889 884 1521 Label 898 891,892 1553 Label 904 900,901 1561 Label 911 906 1601 Label 927 919,920 Page 20 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 1615 Label 933 930 1617 Label 938 935 1619 Label 939 934 1621 Label 948 941,942 1653 Label 954 950,951 1661 Label 961 956 1701 Label 995 990,991 1711 Label 989 985,986 1731 Label 972 969,970 1741 Label 979 976,977 301 Label 281 273,274 315 Label 287 284 317 Label 292 289 319 Label 293 288 321 Label 301 295,296 345 Label 265 261 353 Label 306 303,304 361 Label 313 308 401 Label 330 322,323 415 Label 336 333 417 Label 341 338 419 Label 342 337 421 Label 350 344,345 453 Label 356 352,353 461 Label 363 358 501 Label 379 371,372 515 Label 385 382 517 Label 390 387 519 Label 391 386 521 Label 399 393,394 553 Label 405 401,402 561 Label 412 407 601 Label 428 420,421 615 Label 434 431 617 Label 439 436 619 Label 440 435 621 Label 448 442,443 653 Label 454 450,451 661 Label 461 456 701 Label 477 469,470 715 Label 483 480 717 Label 488 485 719 Label 489 484 721 Label 498 491,492 753 Label 504 500,501 761 Label 511 506 801 Label 527 519,520 815 Label 533 530 817 Label 538 535 819 Label 539 534 821 Label 548 541,542 853 Label 554 550,551 861 Label 561 556 901 Label 577 569,570 915 Label 583 580 Page 21 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 917 Label 588 585 919 Label 589 584 921 Label 598 591,592 953 Label 604 600,601 961 Label 611 606 AB15 Scalar 157 R(4) 4 1 2 COM AB15CM Scalar 226 R(4) 4 1 2 COM AB15WD Scalar 186 R(4) 4 scalar COM ACOMB Scalar 223 R(4) 4 1 15 COM 297,346,395,444,493,543,593,643,69 3,743,793,843,893,943 AG Local 249 R(4) 4 2 2025 277,279,326,328,375,377,424,426,47 3,475,523,525,573,575,623,625,673, 675,723,725,773,775,823,825,873,87 5,923,925 AGG Local 249 R(4) 4 2 2025 278,280,327,329,376,378,425,427,47 4,476,524,526,574,576,624,626,674, 676,724,726,774,776,824,826,874,87 6,924,926 AMOLWT Scalar 26 R(4) 4 scalar COM AO3CM Scalar 226 R(4) 4 1 3 COM AO3RND Scalar 157 R(4) 4 1 3 COM AP Scalar 155 R(4) 4 1 163 COM APCM Scalar 224 R(4) 4 1 15 COM 275,324,373,422,471,521,571,621,67 1,721,771,821,871,921 APWD Scalar 184 R(4) 4 scalar COM ARNDM Scalar 155 R(4) 4 1 163 COM ATP Scalar 156 R(4) 4 1 163 COM ATPCM Scalar 224 R(4) 4 1 15 COM 276,325,374,423,472,522,572,622,67 2,722,772,822,872,922 ATPWD Scalar 184 R(4) 4 scalar COM AWIDE Scalar 183 R(4) 4 scalar COM BANDHI Scalar 157 R(4) 4 1 163 COM BANDLO Scalar 156 R(4) 4 1 163 COM BANDTA Common 155 5900 BCOMB Scalar 223 R(4) 4 1 15 COM 298,347,396,445,494,544,594,644,69 4,744,794,844,894,944 BDCOMB Common 223 736 BDHICM Scalar 225 R(4) 4 1 15 COM BDHIWD Scalar 185 R(4) 4 scalar COM BDLOCM Scalar 225 R(4) 4 1 15 COM BDLOWD Scalar 185 R(4) 4 scalar COM BDWIDE Common 183 56 BETACM Scalar 224 R(4) 4 1 15 COM 496,546,596,646,696,746,796,846,89 6,946 BETAD Scalar 155 R(4) 4 1 163 COM BETAWD Scalar 183 R(4) 4 scalar COM BETINC Scalar 225 R(4) 4 scalar COM BETINW Scalar 185 R(4) 4 scalar COM BO3CM Scalar 226 R(4) 4 1 3 COM BO3RND Scalar 157 R(4) 4 1 3 COM BP Scalar 156 R(4) 4 1 163 COM BPCM Scalar 224 R(4) 4 1 15 COM 275,324,373,422,471,521,571,621,67 1,721,771,821,871,921 BPWD Scalar 184 R(4) 4 scalar COM BRNDM Scalar 155 R(4) 4 1 163 COM Page 22 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References BTP Scalar 156 R(4) 4 1 163 COM BTPCM Scalar 225 R(4) 4 1 15 COM 276,325,374,423,472,522,572,622,67 2,722,772,822,872,922 BTPWD Scalar 184 R(4) 4 scalar COM BWIDE Scalar 183 R(4) 4 scalar COM CLDFAC Dummy 22 R(4) 4 3 95220 ARG,INOUT 300,309,349,359,398,408,447,457,49 7,507,547,557,597,607,647,657,697, 707,747,757,797,807,847,857,897,90 7,947,957,987,988 CO2SP Dummy 24 R(4) 4 2 2070 ARG,INOUT 987 CO2SP1 Dummy 24 R(4) 4 2 2070 ARG,INOUT 696 CO2SP2 Dummy 24 R(4) 4 2 2070 ARG,INOUT 746 CSOUR Dummy 21 R(4) 4 2 2070 ARG,INOUT 993 CSUBP Scalar 26 R(4) 4 scalar COM CTMP Local 244 R(4) 4 2 2070 262,300,305,349,355,398,404,447,45 3,497,503,547,553,597,603,647,653, 697,703,747,753,797,803,847,853,89 7,903,947,953 CTMP2 Local 247 R(4) 4 2 2070 263,987,993 CTMP3 Local 247 R(4) 4 2 2070 264,988,994 CTSO3 Dummy 21 R(4) 4 2 2025 ARG,INOUT 992 DELP Dummy 23 R(4) 4 2 2025 ARG,INOUT 310,360,409,458,508,558,608,658,70 8,758,808,858,908,958,978,992 DELP2 Dummy 23 R(4) 4 2 2025 ARG,INOUT DIFFCTR Scalar 26 R(4) 4 scalar COM EIGHT Scalar 32 R(4) 4 scalar COM EXCTS Dummy 21 R(4) 4 2 2025 ARG,INOUT 305,354,403,452,502,552,602,652,70 2,752,802,852,902,952,971,978 EXP Func 299 scalar 299,348,397,446,495,545,595,645,69 5,745,795,845,895,945 F Local 248 R(4) 4 2 2025 275,277,324,326,373,375,422,424,47 1,473,521,523,571,573,621,623,671, 673,721,723,771,773,821,823,871,87 3,921,923 FAC1 Local 243 R(4) 4 2 2025 297,298,299,346,347,348,395,396,39 7,444,445,446,493,494,495,543,544, 545,593,594,595,643,644,645,693,69 4,695,743,744,745,793,794,795,843, 844,845,893,894,895,943,944,945 FAC2 Local 243 R(4) 4 2 2025 298,299,347,348,396,397,445,446,49 4,495,544,545,594,595,644,645,694, 695,744,745,794,795,844,845,894,89 5,944,945 FF Local 248 R(4) 4 2 2025 276,278,325,327,374,376,423,425,47 2,474,522,524,572,574,622,624,672, 674,722,724,772,774,822,824,872,87 4,922,924 FIFTY Scalar 32 R(4) 4 scalar COM FIVE Scalar 32 R(4) 4 scalar COM FOUR Scalar 33 R(4) 4 scalar COM G Scalar 26 R(4) 4 scalar COM G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GINV Scalar 30 R(4) 4 scalar COM GLB_TABLE Common 37 128 Page 23 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References GP0INV Scalar 30 R(4) 4 scalar COM GRAVDR Scalar 26 R(4) 4 scalar COM GXCTS Dummy 21 R(4) 4 1 45 ARG,INOUT 309,359,408,457,507,557,607,657,70 7,757,807,857,907,957,971 H101M16 Scalar 65 R(4) 4 scalar COM H102M5 Scalar 64 R(4) 4 scalar COM H1036E2 Scalar 63 R(4) 4 scalar COM H114M11 Scalar 47 R(4) 4 scalar COM H1174M7 Scalar 42 R(4) 4 scalar COM H11M10 Scalar 44 R(4) 4 scalar COM H11M11 Scalar 47 R(4) 4 scalar COM H1224E3 Scalar 36 R(4) 4 scalar COM H1226E1 Scalar 58 R(4) 4 scalar COM H128M5 Scalar 41 R(4) 4 scalar COM H129M2 Scalar 62 R(4) 4 scalar COM H12M12 Scalar 51 R(4) 4 scalar COM H12M13 Scalar 52 R(4) 4 scalar COM H12M31 Scalar 54 R(4) 4 scalar COM H135M13 Scalar 52 R(4) 4 scalar COM H1386E2 Scalar 62 R(4) 4 scalar COM H1439M5 Scalar 41 R(4) 4 scalar COM H14M10 Scalar 43 R(4) 4 scalar COM H14M11 Scalar 47 R(4) 4 scalar COM H14M12 Scalar 51 R(4) 4 scalar COM H14M14 Scalar 52 R(4) 4 scalar COM H14M30 Scalar 54 R(4) 4 scalar COM H15E2 Scalar 37 R(4) 4 scalar COM H15M11 Scalar 47 R(4) 4 scalar COM H15M14 Scalar 52 R(4) 4 scalar COM H15M5 Scalar 66 R(4) 4 scalar COM H161E1 Scalar 65 R(4) 4 scalar COM H165E5 Scalar 35 R(4) 4 scalar COM H16E1 Scalar 65 R(4) 4 scalar COM H16M12 Scalar 51 R(4) 4 scalar COM H181E1 Scalar 37 R(4) 4 scalar COM H18E1 Scalar 37 R(4) 4 scalar COM H18E3 Scalar 66 R(4) 4 scalar COM H18M11 Scalar 47 R(4) 4 scalar COM H1E11 Scalar 34 R(4) 4 scalar COM H1E13 Scalar 34 R(4) 4 scalar COM H1E15 Scalar 34 R(4) 4 scalar COM H1E4 Scalar 35 R(4) 4 scalar COM H1E6 Scalar 57 R(4) 4 scalar COM H1E8 Scalar 34 R(4) 4 scalar COM H1M10 Scalar 44 R(4) 4 scalar COM H1M11 Scalar 48 R(4) 4 scalar COM H1M13 Scalar 52 R(4) 4 scalar COM H1M16 Scalar 68 R(4) 4 scalar COM H1M17 Scalar 53 R(4) 4 scalar COM H1M18 Scalar 53 R(4) 4 scalar COM H1M19 Scalar 53 R(4) 4 scalar COM H1M2 Scalar 57 R(4) 4 scalar COM H1M20 Scalar 53 R(4) 4 scalar COM H1M21 Scalar 53 R(4) 4 scalar COM H1M22 Scalar 53 R(4) 4 scalar COM Page 24 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References H1M23 Scalar 53 R(4) 4 scalar COM H1M24 Scalar 54 R(4) 4 scalar COM H1M3 Scalar 40 R(4) 4 scalar COM H1M4 Scalar 40 R(4) 4 scalar COM H1M5 Scalar 41 R(4) 4 scalar COM H1M6 Scalar 42 R(4) 4 scalar COM H1M60 Scalar 55 R(4) 4 scalar COM H1M8 Scalar 43 R(4) 4 scalar COM H1P082 Scalar 62 R(4) 4 scalar COM H1P25892 Scalar 38 R(4) 4 scalar COM H1P4 Scalar 38 R(4) 4 scalar COM H1P41819 Scalar 60 R(4) 4 scalar COM 277,278,326,327,375,376,424,425,47 3,474,523,524,573,574,623,624,673, 674,723,724,773,774,823,824,873,87 4,923,924 H1P4387 Scalar 38 R(4) 4 scalar COM H1P8 Scalar 38 R(4) 4 scalar COM H2075E3 Scalar 36 R(4) 4 scalar COM H20788E3 Scalar 35 R(4) 4 scalar COM H2118M2 Scalar 63 R(4) 4 scalar COM H21M12 Scalar 50 R(4) 4 scalar COM H21M31 Scalar 54 R(4) 4 scalar COM H235M3 Scalar 61 R(4) 4 scalar COM H23E2 Scalar 37 R(4) 4 scalar COM H23M10 Scalar 43 R(4) 4 scalar COM H23M11 Scalar 46 R(4) 4 scalar COM H24E3 Scalar 35 R(4) 4 scalar COM H24M11 Scalar 46 R(4) 4 scalar COM H24M12 Scalar 50 R(4) 4 scalar COM H25452M6 Scalar 42 R(4) 4 scalar COM H257M8 Scalar 43 R(4) 4 scalar COM H25E2 Scalar 66 R(4) 4 scalar COM 256 H25M31 Scalar 54 R(4) 4 scalar COM H26E2 Scalar 60 R(4) 4 scalar COM H26M30 Scalar 54 R(4) 4 scalar COM H28571M2 Scalar 68 R(4) 4 scalar COM H285M4 Scalar 40 R(4) 4 scalar COM H28E1 Scalar 71 R(4) 4 scalar COM H28M11 Scalar 46 R(4) 4 scalar COM H28M12 Scalar 50 R(4) 4 scalar COM H29316E2 Scalar 58 R(4) 4 scalar COM H2945E2 Scalar 36 R(4) 4 scalar COM H29M12 Scalar 50 R(4) 4 scalar COM H2E2 Scalar 60 R(4) 4 scalar COM H2E6 Scalar 57 R(4) 4 scalar COM H2M11 Scalar 47 R(4) 4 scalar COM H2P5 Scalar 38 R(4) 4 scalar COM H2P8 Scalar 37 R(4) 4 scalar COM H2P9 Scalar 37 R(4) 4 scalar COM H3082E2 Scalar 36 R(4) 4 scalar COM H3116E1 Scalar 58 R(4) 4 scalar COM H323M4 Scalar 63 R(4) 4 scalar COM H327M8 Scalar 43 R(4) 4 scalar COM H32M11 Scalar 46 R(4) 4 scalar COM H35E1 Scalar 37 R(4) 4 scalar COM Page 25 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References H35M11 Scalar 46 R(4) 4 scalar COM H36M13 Scalar 51 R(4) 4 scalar COM H37412M5 Scalar 41 R(4) 4 scalar COM H37M11 Scalar 45 R(4) 4 scalar COM H37M12 Scalar 50 R(4) 4 scalar COM H38M12 Scalar 49 R(4) 4 scalar COM H391M7 Scalar 42 R(4) 4 scalar COM H394M5 Scalar 41 R(4) 4 scalar COM H3E2 Scalar 36 R(4) 4 scalar COM H3M11 Scalar 46 R(4) 4 scalar COM H3M12 Scalar 50 R(4) 4 scalar COM H3M14 Scalar 52 R(4) 4 scalar COM H3M3 Scalar 65 R(4) 4 scalar COM H3M4 Scalar 69 R(4) 4 scalar COM H3P5 Scalar 66 R(4) 4 scalar COM H3P6 Scalar 37 R(4) 4 scalar COM H41666M2 Scalar 58 R(4) 4 scalar COM H42M11 Scalar 45 R(4) 4 scalar COM H42M2 Scalar 63 R(4) 4 scalar COM H44194M2 Scalar 60 R(4) 4 scalar COM 275,276,324,325,373,374,422,423,47 1,472,521,522,571,572,621,622,671, 672,721,722,771,772,821,822,871,87 2,921,922 H44871M2 Scalar 40 R(4) 4 scalar COM H44M11 Scalar 45 R(4) 4 scalar COM H44M12 Scalar 49 R(4) 4 scalar COM H451M6 Scalar 64 R(4) 4 scalar COM H45M12 Scalar 49 R(4) 4 scalar COM H45M32 Scalar 55 R(4) 4 scalar COM H46M13 Scalar 51 R(4) 4 scalar COM H488E4 Scalar 35 R(4) 4 scalar COM H48M11 Scalar 45 R(4) 4 scalar COM H4999M6 Scalar 42 R(4) 4 scalar COM H4E5 Scalar 34 R(4) 4 scalar COM H4M12 Scalar 49 R(4) 4 scalar COM H4M33 Scalar 55 R(4) 4 scalar COM H53M11 Scalar 45 R(4) 4 scalar COM H559M3 Scalar 40 R(4) 4 scalar COM H55M32 Scalar 55 R(4) 4 scalar COM H5725E4 Scalar 35 R(4) 4 scalar COM H5E2 Scalar 36 R(4) 4 scalar COM H625M2 Scalar 59 R(4) 4 scalar COM H62M12 Scalar 49 R(4) 4 scalar COM H62M34 Scalar 55 R(4) 4 scalar COM H658M2 Scalar 63 R(4) 4 scalar COM H65M12 Scalar 48 R(4) 4 scalar COM H67390E2 Scalar 64 R(4) 4 scalar COM H6938M5 Scalar 41 R(4) 4 scalar COM H69766E5 Scalar 61 R(4) 4 scalar COM H6M12 Scalar 49 R(4) 4 scalar COM H6P08108 Scalar 67 R(4) 4 scalar COM H71E26 Scalar 34 R(4) 4 scalar COM H72M11 Scalar 45 R(4) 4 scalar COM H74M12 Scalar 48 R(4) 4 scalar COM H75826M4 Scalar 62 R(4) 4 scalar COM Page 26 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References H77M11 Scalar 44 R(4) 4 scalar COM H77M12 Scalar 48 R(4) 4 scalar COM H7M6 Scalar 42 R(4) 4 scalar COM H8121E1 Scalar 60 R(4) 4 scalar COM H82M11 Scalar 44 R(4) 4 scalar COM H83E26 Scalar 34 R(4) 4 scalar COM H83M11 Scalar 44 R(4) 4 scalar COM H8725M8 Scalar 43 R(4) 4 scalar COM H8M11 Scalar 44 R(4) 4 scalar COM H8M13 Scalar 51 R(4) 4 scalar COM H93M12 Scalar 48 R(4) 4 scalar COM H96M12 Scalar 48 R(4) 4 scalar COM H987M4 Scalar 40 R(4) 4 scalar COM H9M32 Scalar 55 R(4) 4 scalar COM H9P94 Scalar 59 R(4) 4 scalar COM HAF Scalar 33 R(4) 4 scalar COM 310,360,409,458,508,558,608,658,70 8,758,808,858,908,958 HCON Common 32 872 HM13EZ Scalar 56 R(4) 4 scalar COM HM1597E1 Scalar 66 R(4) 4 scalar COM HM161E1 Scalar 65 R(4) 4 scalar COM HM1797E1 Scalar 59 R(4) 4 scalar COM HM181E1 Scalar 56 R(4) 4 scalar COM HM19EZ Scalar 56 R(4) 4 scalar COM HM1E1 Scalar 56 R(4) 4 scalar COM HM1E2 Scalar 56 R(4) 4 scalar COM HM1EZ Scalar 60 R(4) 4 scalar COM 299,348,397,446,495,545,595,645,69 5,745,795,845,895,945 HM2M2 Scalar 58 R(4) 4 scalar COM HM6666M2 Scalar 57 R(4) 4 scalar COM HM8E1 Scalar 70 R(4) 4 scalar COM HMP5 Scalar 58 R(4) 4 scalar COM HMP575 Scalar 56 R(4) 4 scalar COM HMP66667 Scalar 57 R(4) 4 scalar COM HMP805 Scalar 67 R(4) 4 scalar COM HNINETY Scalar 32 R(4) 4 scalar COM HP1 Scalar 39 R(4) 4 scalar COM HP118666 Scalar 66 R(4) 4 scalar COM HP144 Scalar 61 R(4) 4 scalar COM HP166666 Scalar 57 R(4) 4 scalar COM HP219 Scalar 61 R(4) 4 scalar COM HP228 Scalar 59 R(4) 4 scalar COM HP26 Scalar 61 R(4) 4 scalar COM HP369 Scalar 39 R(4) 4 scalar COM HP3795 Scalar 64 R(4) 4 scalar COM HP5048 Scalar 64 R(4) 4 scalar COM HP518 Scalar 38 R(4) 4 scalar COM HP526315 Scalar 67 R(4) 4 scalar COM HP6 Scalar 59 R(4) 4 scalar COM HP602409 Scalar 67 R(4) 4 scalar COM HP60241 Scalar 59 R(4) 4 scalar COM HP8 Scalar 38 R(4) 4 scalar COM HP805 Scalar 62 R(4) 4 scalar COM HP816 Scalar 61 R(4) 4 scalar COM HTMG Scalar 52 R(4) 4 3 4183695 COM Page 27 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References HUNDRED Scalar 32 R(4) 4 scalar COM I Local 255 I(4) 4 scalar 255,256,257,261,262,263,264,274,27 5,276,277,278,279,280,284,285,286, 289,290,291,296,297,298,299,300,30 4,305,308,309,310,311,312,323,324, 325,326,327,328,329,333,334,335,33 8,339,340,345,346,347,348,349,353, 354,355,358,359,360,361,362,372,37 3,374,375,376,377,378,382,383,384, 387,388,389,394,395,396,397,398,40 2,403,404,407,408,409,410,411,421, 422,423,424,425,426,427,431,432,43 3,436,437,438,443,444,445,446,447, 451,452,453,456,457,458,459,460,47 0,471,472,473,474,475,476,480,481, 482,485,486,487,492,493,494,495,49 6,497,501,502,503,506,507,508,509, 510,520,521,522,523,524,525,526,53 0,531,532,535,536,537,542,543,544, 545,546,547,551,552,553,556,557,55 8,559,560,570,571,572,573,574,575, 576,580,581,582,585,586,587,592,59 3,594,595,596,597,601,602,603,606, 607,608,609,610,620,621,622,623,62 4,625,626,630,631,632,635,636,637, 642,643,644,645,646,647,651,652,65 3,656,657,658,659,660,670,671,672, 673,674,675,676,680,681,682,685,68 6,687,692,693,694,695,696,697,701, 702,703,706,707,708,709,710,720,72 1,722,723,724,725,726,730,731,732, 735,736,737,742,743,744,745,746,74 7,751,752,753,756,757,758,759,760, 770,771,772,773,774,775,776,780,78 1,782,785,786,787,792,793,794,795, 796,797,801,802,803,806,807,808,80 9,810,820,821,822,823,824,825,826, 830,831,832,835,836,837,842,843,84 4,845,846,847,851,852,853,856,857, 858,859,860,870,871,872,873,874,87 5,876,880,881,882,885,886,887,892, 893,894,895,896,897,901,902,903,90 6,907,908,909,910,920,921,922,923, 924,925,926,930,931,932,935,936,93 7,942,943,944,945,946,947,951,952, 953,956,957,958,959,960,970,971,97 7,978,986,987,988,991,992,993,994 IBAND Scalar 223 I(4) 4 1 40 COM IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDIM1 Param 19 I(4) 4 scalar 228,229,230,231,232,233,234,235,23 6,237,238,239,241,242,243,244,245, 246,247,248,249 IDIM2 Param 19 I(4) 4 scalar 228,229,230,231,232,233,234,235,23 6,237,238,239,241,242,243,244,245, Page 28 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 246,247,248,249 IE_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 ILCOL Scalar 22 I(4) 4 scalar COM ILPAD1 Scalar 23 I(4) 4 scalar COM ILPAD2 Scalar 23 I(4) 4 scalar COM ILPAD3 Scalar 23 I(4) 4 scalar COM ILPAD4 Scalar 23 I(4) 4 scalar COM ILPAD5 Scalar 23 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 16,19,43,44,45,46,47,51,52,57,92 IMAX Param 92 I(4) 4 scalar 92,105 INLTE Param 104 I(4) 4 scalar 104 INLTEP Param 104 I(4) 4 scalar INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INUMQ Scalar 34 I(4) 4 1 100 COM IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IRCOL Scalar 22 I(4) 4 scalar COM IRPAD1 Scalar 24 I(4) 4 scalar COM IRPAD2 Scalar 24 I(4) 4 scalar COM IRPAD3 Scalar 24 I(4) 4 scalar COM IRPAD4 Scalar 24 I(4) 4 scalar COM IRPAD5 Scalar 24 I(4) 4 scalar COM IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM ITAIL Param 16 I(4) 4 scalar ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITROW Scalar 22 I(4) 4 scalar COM JBPAD1 Scalar 25 I(4) 4 scalar COM JBPAD2 Scalar 25 I(4) 4 scalar COM JBPAD3 Scalar 25 I(4) 4 scalar COM JBPAD4 Scalar 25 I(4) 4 scalar COM JBPAD5 Scalar 25 I(4) 4 scalar COM JDIM1 Param 20 I(4) 4 scalar JDIM2 Param 20 I(4) 4 scalar JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20,43,44,45,46,47,51,52,57 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JS_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM JTAIL Param 17 I(4) 4 scalar JTPAD1 Scalar 26 I(4) 4 scalar COM JTPAD2 Scalar 26 I(4) 4 scalar COM JTPAD3 Scalar 26 I(4) 4 scalar COM JTPAD4 Scalar 26 I(4) 4 scalar COM JTPAD5 Scalar 26 I(4) 4 scalar COM K Local 254 I(4) 4 scalar 254,256,257,273,275,276,277,278,27 9,280,288,290,291,295,297,298,299, 300,303,305,322,324,325,326,327,32 8,329,337,339,340,344,346,347,348, Page 29 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References 349,352,354,355,371,373,374,375,37 6,377,378,386,388,389,393,395,396, 397,398,401,403,404,420,422,423,42 4,425,426,427,435,437,438,442,444, 445,446,447,450,452,453,469,471,47 2,473,474,475,476,484,486,487,491, 493,494,495,496,497,500,502,503,51 9,521,522,523,524,525,526,534,536, 537,541,543,544,545,546,547,550,55 2,553,569,571,572,573,574,575,576, 584,586,587,591,593,594,595,596,59 7,600,602,603,619,621,622,623,624, 625,626,634,636,637,641,643,644,64 5,646,647,650,652,653,669,671,672, 673,674,675,676,684,686,687,691,69 3,694,695,696,697,700,702,703,719, 721,722,723,724,725,726,734,736,73 7,741,743,744,745,746,747,750,752, 753,769,771,772,773,774,775,776,78 4,786,787,791,793,794,795,796,797, 800,802,803,819,821,822,823,824,82 5,826,834,836,837,841,843,844,845, 846,847,850,852,853,869,871,872,87 3,874,875,876,884,886,887,891,893, 894,895,896,897,900,902,903,919,92 1,922,923,924,925,926,934,936,937, 941,943,944,945,946,947,950,952,95 3,969,971,976,978,985,987,988,990, 992,993,994 KO2 Param 107 I(4) 4 scalar 108 KO21 Param 108 I(4) 4 scalar KO2M Param 108 I(4) 4 scalar L Param 91 I(4) 4 scalar 95,96,97,100,102,231,232,233,234,2 38,241,242,243,244,245,246,248,249 ,254,273,288,295,303,309,310,311,3 12,322,337,344,352,359,360,361,362 ,371,386,393,401,408,409,410,411,4 20,435,442,450,457,458,459,460,469 ,484,491,500,507,508,509,510,519,5 34,541,550,557,558,559,560,569,584 ,591,600,607,608,609,610,619,634,6 41,650,657,658,659,660,669,684,691 ,700,707,708,709,710,719,734,741,7 50,757,758,759,760,769,784,791,800 ,807,808,809,810,819,834,841,850,8 57,858,859,860,869,884,891,900,907 ,908,909,910,919,934,941,950,957,9 58,959,960,969,976,985,990 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LL Param 97 I(4) 4 scalar 97,98 LL3P Param 102 I(4) 4 scalar 105 LL3PI Param 105 I(4) 4 scalar LLM1 Param 98 I(4) 4 scalar LLM2 Param 98 I(4) 4 scalar Page 30 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References LLM3 Param 98 I(4) 4 scalar LLP1 Param 97 I(4) 4 scalar 105 LLP1I Param 105 I(4) 4 scalar LLP2 Param 97 I(4) 4 scalar LLP3 Param 97 I(4) 4 scalar LM Param 3 I(4) 4 scalar 52,91 LM1 Param 96 I(4) 4 scalar 310,360,409,458,508,558,608,658,70 8,758,808,858,908,958 LM2 Param 96 I(4) 4 scalar LM3 Param 96 I(4) 4 scalar LP1 Param 95 I(4) 4 scalar 99,100,101,105,228,229,230,232,234 ,235,236,237,244,247,309,310,311,3 12,359,360,361,362,408,409,410,411 ,457,458,459,460,507,508,509,510,5 57,558,559,560,607,608,609,610,657 ,658,659,660,707,708,709,710,757,7 58,759,760,807,808,809,810,857,858 ,859,860,907,908,909,910,957,958,9 59,960 LP121 Param 101 I(4) 4 scalar LP1I Param 105 I(4) 4 scalar LP1M Param 99 I(4) 4 scalar 99 LP1M1 Param 99 I(4) 4 scalar LP1V Param 100 I(4) 4 scalar LP2 Param 95 I(4) 4 scalar LP3 Param 95 I(4) 4 scalar LSM Param 3 I(4) 4 scalar MAPPINGS Common 56 5024 MPI_COMM_COMP Scalar 30 I(4) 4 scalar COM MPI_COMM_INTER Scalar 31 I(4) 4 scalar COM MPI_COMM_INTER_ARRAY Scalar 32 I(4) 4 1 100 COM MPPCOM Common 1 1464 MYIE Scalar 4 I(4) 4 scalar COM 255,261,274,284,289,296,304,308,32 3,333,338,345,353,358,372,382,387, 394,402,407,421,431,436,443,451,45 6,470,480,485,492,501,506,520,530, 535,542,551,556,570,580,585,592,60 1,606,620,630,635,642,651,656,670, 680,685,692,701,706,720,730,735,74 2,751,756,770,780,785,792,801,806, 820,830,835,842,851,856,870,880,88 5,892,901,906,920,930,935,942,951, 956,970,977,986,991 MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM MYIE1_P3 Scalar 8 I(4) 4 scalar COM MYIE1_P4 Scalar 8 I(4) 4 scalar COM MYIE2 Scalar 4 I(4) 4 scalar COM MYIE2_P1 Scalar 9 I(4) 4 scalar COM MYIE_P1 Scalar 7 I(4) 4 scalar COM MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM MYIE_P4 Scalar 7 I(4) 4 scalar COM MYIE_P5 Scalar 7 I(4) 4 scalar COM Page 31 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References MYIS Scalar 4 I(4) 4 scalar COM 255,261,274,284,289,296,304,308,32 3,333,338,345,353,358,372,382,387, 394,402,407,421,431,436,443,451,45 6,470,480,485,492,501,506,520,530, 535,542,551,556,570,580,585,592,60 1,606,620,630,635,642,651,656,670, 680,685,692,701,706,720,730,735,74 2,751,756,770,780,785,792,801,806, 820,830,835,842,851,856,870,880,88 5,892,901,906,920,930,935,942,951, 956,970,977,986,991 MYIS1 Scalar 4 I(4) 4 scalar COM MYIS1_P1 Scalar 6 I(4) 4 scalar COM MYIS1_P2 Scalar 6 I(4) 4 scalar COM MYIS1_P3 Scalar 6 I(4) 4 scalar COM MYIS1_P4 Scalar 6 I(4) 4 scalar COM MYIS2 Scalar 4 I(4) 4 scalar COM MYIS_P1 Scalar 5 I(4) 4 scalar COM MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM MYIS_P4 Scalar 5 I(4) 4 scalar COM MYIS_P5 Scalar 5 I(4) 4 scalar COM MYJE Scalar 15 I(4) 4 scalar COM MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM MYJE2_P1 Scalar 18 I(4) 4 scalar COM MYJE2_P2 Scalar 18 I(4) 4 scalar COM MYJE2_P3 Scalar 18 I(4) 4 scalar COM MYJE2_P4 Scalar 18 I(4) 4 scalar COM MYJE3 Scalar 15 I(4) 4 scalar COM MYJE3_P4 Scalar 18 I(4) 4 scalar COM MYJE4 Scalar 15 I(4) 4 scalar COM MYJE4_P1 Scalar 19 I(4) 4 scalar COM MYJE4_P4 Scalar 19 I(4) 4 scalar COM MYJE5 Scalar 15 I(4) 4 scalar COM MYJE5_P1 Scalar 19 I(4) 4 scalar COM MYJE5_P2 Scalar 19 I(4) 4 scalar COM MYJE_P1 Scalar 16 I(4) 4 scalar COM MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM MYJE_P4 Scalar 16 I(4) 4 scalar COM MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM MYJS1_P3 Scalar 12 I(4) 4 scalar COM MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM MYJS2_P1 Scalar 13 I(4) 4 scalar COM MYJS2_P2 Scalar 13 I(4) 4 scalar COM Page 32 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM MYJS3 Scalar 10 I(4) 4 scalar COM MYJS3_P4 Scalar 13 I(4) 4 scalar COM MYJS4 Scalar 10 I(4) 4 scalar COM MYJS4_P1 Scalar 14 I(4) 4 scalar COM MYJS4_P4 Scalar 14 I(4) 4 scalar COM MYJS5 Scalar 10 I(4) 4 scalar COM MYJS5_P1 Scalar 14 I(4) 4 scalar COM MYJS5_P2 Scalar 14 I(4) 4 scalar COM MYJS_P1 Scalar 11 I(4) 4 scalar COM MYJS_P2 Scalar 11 I(4) 4 scalar COM MYJS_P3 Scalar 11 I(4) 4 scalar COM MYJS_P4 Scalar 11 I(4) 4 scalar COM MYJS_P5 Scalar 11 I(4) 4 scalar COM MYPE Scalar 2 I(4) 4 scalar COM MY_E Scalar 20 I(4) 4 scalar COM MY_IE_GLB Scalar 2 I(4) 4 scalar COM MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM MY_JS_LOC Scalar 3 I(4) 4 scalar COM MY_N Scalar 20 I(4) 4 scalar COM MY_NE Scalar 21 I(4) 4 scalar COM MY_NEB Scalar 21 I(4) 4 1 8 COM MY_NW Scalar 21 I(4) 4 scalar COM MY_S Scalar 20 I(4) 4 scalar COM MY_SE Scalar 21 I(4) 4 scalar COM MY_SW Scalar 21 I(4) 4 scalar COM MY_W Scalar 20 I(4) 4 scalar COM NB Param 103 I(4) 4 scalar 106 NB1 Param 106 I(4) 4 scalar NBLM Param 94 I(4) 4 scalar NBLW Param 93 I(4) 4 scalar 155,156,157 NBLX Param 93 I(4) 4 scalar NBLY Param 93 I(4) 4 scalar 94,101,223,224,225,228 NCOL Param 92 I(4) 4 scalar NNLTE Param 104 I(4) 4 scalar NPES Scalar 2 I(4) 4 scalar COM O3DIFCTR Scalar 26 R(4) 4 scalar COM ONE Scalar 33 R(4) 4 scalar COM 262,277,278,326,327,375,376,424,42 5,473,474,495,523,524,545,573,574, 595,623,624,645,673,674,695,723,72 4,745,773,774,795,823,824,845,873, 874,895,923,924,945 P Dummy 23 R(4) 4 2 2070 ARG,INOUT 310,311,360,361,409,410,458,459,50 8,509,558,559,608,609,658,659,708, 709,758,759,808,809,858,859,908,90 9,958,959 P0 Scalar 26 R(4) 4 scalar COM P0INV Scalar 30 R(4) 4 scalar COM P0X2 Scalar 27 R(4) 4 scalar COM Page 33 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References P0XZP2 Scalar 27 R(4) 4 scalar COM P0XZP8 Scalar 27 R(4) 4 scalar COM PHITMP Local 241 R(4) 4 2 2025 279,285,290,328,334,339,377,383,38 8,426,432,437,475,481,486,525,531, 536,575,581,586,625,631,636,675,68 1,686,725,731,736,775,781,786,825, 831,836,875,881,886,925,931,936 PHYCON Common 26 80 PRESS Dummy 22 R(4) 4 2 2070 ARG,INOUT 310,311,360,361,409,410,458,459,50 8,509,558,559,608,609,658,659,708, 709,758,759,808,809,858,859,908,90 9,958,959 PSITMP Local 241 R(4) 4 2 2025 280,286,291,329,335,340,378,384,38 9,427,433,438,476,482,487,526,532, 537,576,582,587,626,632,637,676,68 2,687,726,732,737,776,782,787,826, 832,837,876,882,887,926,932,937 QUARTR Scalar 33 R(4) 4 scalar COM RADCON Scalar 27 R(4) 4 scalar COM 978,992 RADCON1 Scalar 29 R(4) 4 scalar COM RATCO2MW Scalar 28 R(4) 4 scalar COM RATH2OMW Scalar 28 R(4) 4 scalar COM RGAS Scalar 27 R(4) 4 scalar COM RGASSP Scalar 27 R(4) 4 scalar COM SECPDA Scalar 27 R(4) 4 scalar COM SIXTY Scalar 32 R(4) 4 scalar COM SKC1R Scalar 186 R(4) 4 scalar COM SKO2D Scalar 186 R(4) 4 scalar COM 496,546,596,646,696,746,796,846,89 6,946 SKO3R Scalar 186 R(4) 4 scalar COM SORC Dummy 21 R(4) 4 3 31050 ARG,INOUT 305,309,312,354,359,362,403,408,41 1,452,457,460,502,507,510,552,557, 560,602,607,610,652,657,660,702,70 7,710,752,757,760,802,807,810,852, 857,860,902,907,910,952,957,960,99 4 SPA88 Subr 21 SQRT Func 299 scalar 299,348,397,446,495,545,595,645,69 5,745,795,845,895,945 TEMP Dummy 22 R(4) 4 2 2070 ARG,INOUT 256 TEMP1 Scalar 43 R(4) 4 2 91719 COM TEMP10 Scalar 45 R(4) 4 2 91719 COM TEMP11 Scalar 45 R(4) 4 2 91719 COM TEMP12 Scalar 45 R(4) 4 2 91719 COM TEMP13 Scalar 46 R(4) 4 2 91719 COM TEMP14 Scalar 46 R(4) 4 2 91719 COM TEMP15 Scalar 46 R(4) 4 2 91719 COM TEMP16 Scalar 46 R(4) 4 2 91719 COM TEMP2 Scalar 43 R(4) 4 2 91719 COM TEMP2X Scalar 51 R(4) 4 2 92971 COM TEMP3 Scalar 43 R(4) 4 2 91719 COM TEMP4 Scalar 43 R(4) 4 2 91719 COM TEMP5 Scalar 44 R(4) 4 2 91719 COM TEMP6 Scalar 44 R(4) 4 2 91719 COM TEMP7 Scalar 44 R(4) 4 2 91719 COM Page 34 Source Listing SPA88 2014-12-17 20:47 Symbol Table SPA88.F Name Object Declared Type Bytes Dimen Elements Attributes References TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 TEN Scalar 32 R(4) 4 scalar COM THREE Scalar 33 R(4) 4 scalar COM TO3SP Dummy 23 R(4) 4 2 2070 ARG,INOUT 988 TO3SPC Dummy 23 R(4) 4 2 2025 ARG,INOUT 896 TOPM Local 246 R(4) 4 2 2025 285,290,297,298,334,339,346,347,38 3,388,395,396,432,437,444,445,481, 486,493,494,531,536,543,544,581,58 6,593,594,631,636,643,644,681,686, 693,694,731,736,743,744,781,786,79 3,794,831,836,843,844,881,886,893, 894,931,936,943,944 TOPO Common 50 17478548 TOPPHI Local 246 R(4) 4 2 2025 286,291,298,335,340,347,384,389,39 6,433,438,445,482,487,494,532,537, 544,582,587,594,632,637,644,682,68 7,694,732,737,744,782,787,794,832, 837,844,882,887,894,932,937,944 TOTVO2 Dummy 23 R(4) 4 2 2070 ARG,INOUT 496,546,596,646,696,746,796,846,89 6,946 TT Local 242 R(4) 4 2 2025 299,300,309,310,311,348,349,359,36 0,361,397,398,408,409,410,446,447, 457,458,459,495,497,507,508,509,54 5,547,557,558,559,595,597,607,608, 609,645,647,657,658,659,695,697,70 7,708,709,745,747,757,758,759,795, 797,807,808,809,845,847,857,858,85 9,895,897,907,908,909,945,947,957, 958,959 TTVG Scalar 51 R(4) 4 2 92971 COM TWO Scalar 33 R(4) 4 scalar COM 311,361,410,459,509,559,609,659,70 9,759,809,859,909,959 VAR1 Dummy 22 R(4) 4 2 2025 ARG,INOUT 279,328,377,426,475,525,575,625,67 5,725,775,825,875,925 VAR2 Dummy 22 R(4) 4 2 2025 ARG,INOUT 280,329,378,427,476,526,576,626,67 6,726,776,826,876,926 X Local 244 R(4) 4 2 2025 256,257,275,276,324,325,373,374,42 2,423,471,472,521,522,571,572,621, 622,671,672,721,722,771,772,821,82 2,871,872,921,922 Y Local 245 R(4) 4 2 2025 257,275,276,324,325,373,374,422,42 3,471,472,521,522,571,572,621,622, 671,672,721,722,771,772,821,822,87 1,872,921,922 ZERO Scalar 33 R(4) 4 scalar COM Page 35 Source Listing SPA88 2014-12-17 20:47 Subprograms/Common Blocks SPA88.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BANDTA Common 155 5900 BDCOMB Common 223 736 BDWIDE Common 183 56 GLB_TABLE Common 37 128 HCON Common 32 872 MAPPINGS Common 56 5024 MPPCOM Common 1 1464 PHYCON Common 26 80 SPA88 Subr 21 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ Page 36 Source Listing SPA88 2014-12-17 20:47 SPA88.F -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : SPA88.lst -o filename : none Page 37 Source Listing SPA88 2014-12-17 20:47 SPA88.F COMPILER: Intel(R) Fortran 12.1-2100