Page 1 Source Listing SWR93 2014-12-17 20:47 /tmpnwprd/ifortISLGob.i 1 # 1 "SWR93.F" 2 SUBROUTINE SWR93(FSWC,HSWC,UFSWC,DFSWC,FSWL,HSWL,UFSWL, 3 & DFSWL, 4 1 PRESS,COSZRO,TAUDAR,RH2O,RRCO2,SSOLAR,QO3, 5 2 NCLDS,KTOPSW,KBTMSW,CAMT,CRR,CTT, 6 A ALVB,ALNB,ALVD,ALND,GDFVB,GDFNB,GDFVD,GDFND) 7 C===> ********************************************************* 8 C --- SWR91SIB --- MODIFIED FROM SWR89-BAND12....YUTAI HOU 9 C 10 C -SW- RADIATION CODE............................ 11 C INPUTS:PRESS,COSZRO,TAUDAR,RH2O,RRCO2,SSOLAR,QO3,NCLDS, 12 C KTOPSW,KBTMSW,CIRAB,CIRRF,CUVRF,CAMT, 13 C ALVB,ALVD,ALNB,ALND; 14 C OUTPUT:FSWC,HSWC,UFSWC,DFSWC,FSWL,HSWL,UFSWL,DFSWL, 15 C GDFVB,GDFVD,GDFNB,GDFND. 16 C --- SWR91SIB --- MODIFIED BY K. CAMPANA..06 MAR 92 17 C INCLUDE HPCON,PARMC CHANGED TO HCON,RDPARM 18 C 6 Q..... VARIABLE NAMES RESTORED TO ORIGINAL 7,8 CHAR 19 C CHANGE O3DIFF,DIFFCC TO O3DIFCTR,DIFFCTR 20 C --- SWR91SIB --- MODIFIED BY Y. HOU FEB 93 21 C INPUTS 12 BANDS CLD REFLECTTANCE AND TRANSMITTANCE 22 C CRR,CTT TO REPLACE CIRAB,CIRRF,CUVRF 23 C===> ********************************************************* 24 C 25 C 26 INCLUDE "parmeta" 49 INCLUDE "HCON.comm" 97 INCLUDE "rdparm" 129 INCLUDE "mpp.h" 130 # 1 "./sp.h" 1 131 # 4 132 133 # 30 "SWR93.F" 2 134 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE: 135 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL 136 C NB IS A SHORTWAVE PARAMETER; OTHER QUANTITIES ARE DERIVED 137 C FROM THE ABOVE PARAMETERS. 138 C --- VARIABLES AS IN ARGUMENT LIST 197 D I M E N S I O N 198 1 FSWC (IDIM1:IDIM2,LP1), HSWC (IDIM1:IDIM2,LP1), 199 2 CRR (IDIM1:IDIM2,NB,LP1) 200 3, FSWL (IDIM1:IDIM2,LP1), HSWL (IDIM1:IDIM2,LP1), 201 4 CTT (IDIM1:IDIM2,NB,LP1) 202 5, UFSWC (IDIM1:IDIM2,LP1), DFSWC (IDIM1:IDIM2,LP1) 203 6, UFSWL (IDIM1:IDIM2,LP1), DFSWL (IDIM1:IDIM2,LP1) 204 7, PRESS (IDIM1:IDIM2,LP1), RH2O (IDIM1:IDIM2,L), 205 8 QO3 (IDIM1:IDIM2,L) 206 9, CAMT (IDIM1:IDIM2,LP1), KTOPSW(IDIM1:IDIM2,LP1), 207 o KBTMSW(IDIM1:IDIM2,LP1) 208 1, COSZRO(IDIM1:IDIM2), TAUDAR(IDIM1:IDIM2), 209 2 NCLDS (IDIM1:IDIM2) 210 3, ALVB (IDIM1:IDIM2), ALNB (IDIM1:IDIM2), 211 4 ALVD (IDIM1:IDIM2), ALND (IDIM1:IDIM2) 212 5, GDFVB (IDIM1:IDIM2), GDFNB (IDIM1:IDIM2), GDFVD (IDIM1:IDIM2), 213 6 GDFND (IDIM1:IDIM2) 214 C --- LOCAL VARIABLES 215 D I M E N S I O N Page 2 Source Listing SWR93 2014-12-17 20:47 SWR93.F 216 1 PP (IDIM1:IDIM2,LP1), DP (IDIM1:IDIM2,LP1), 217 2 PR2 (IDIM1:IDIM2,LP1) 218 3, DU (IDIM1:IDIM2,LP1), DUCO2 (IDIM1:IDIM2,LP1), 219 4 DUO3 (IDIM1:IDIM2,LP1) 220 5, FF (IDIM1:IDIM2,LP1), FFCO2 (IDIM1:IDIM2,LP1), 221 6 FFO3 (IDIM1:IDIM2,LP1) 222 7, RRAY (IDIM1:IDIM2), DFNTOP(IDIM1:IDIM2,NB), 223 8 SECZ (IDIM1:IDIM2) 224 9, REFL (IDIM1:IDIM2), TMP1 (IDIM1:IDIM2), 225 o REFL2 (IDIM1:IDIM2) 226 1, CCMAX (IDIM1:IDIM2), XAMT (IDIM1:IDIM2,LP1) 227 D I M E N S I O N 228 1 UD (IDIM1:IDIM2,LP1), UR (IDIM1:IDIM2,LP1) 229 2, UCO2 (IDIM1:IDIM2,LLP2), UDCO2 (IDIM1:IDIM2,LP1), 230 3 URCO2 (IDIM1:IDIM2,LP1) 231 4, UO3 (IDIM1:IDIM2,LLP2), UDO3 (IDIM1:IDIM2,LP1), 232 5 URO3 (IDIM1:IDIM2,LP1) 233 6, TCO2 (IDIM1:IDIM2,LLP2), TDCO2 (IDIM1:IDIM2,LP1), 234 7 TUCO2 (IDIM1:IDIM2,LP1) 235 8, TO3 (IDIM1:IDIM2,LLP2), TDO3 (IDIM1:IDIM2,LP1), 236 9 TUO3 (IDIM1:IDIM2,LP1) 237 D I M E N S I O N 238 1 DFN (IDIM1:IDIM2,LP1), UFN (IDIM1:IDIM2,LP1), 239 2 CR (IDIM1:IDIM2,LP1) 240 3, TTD (IDIM1:IDIM2,LP1), TTU (IDIM1:IDIM2,LP1), 241 4 CT (IDIM1:IDIM2,LP1) 242 5, PPTOP (IDIM1:IDIM2,LP1), DPCLD (IDIM1:IDIM2,LP1) 243 C --- EQUIVALENCED LOCAL VARIABLES 244 D I M E N S I O N 245 1 TTUB1 (IDIM1:IDIM2,LP1), TUCL1 (IDIM1:IDIM2,LP1) 246 2, TTDB1 (IDIM1:IDIM2,LP1), TDCL1 (IDIM1:IDIM2,LP1), 247 3 TDCL2 (IDIM1:IDIM2,LP1) 248 4, UFNTRN(IDIM1:IDIM2,LP1), UFNCLU(IDIM1:IDIM2,LP1), 249 5 TCLU (IDIM1:IDIM2,LP1) 250 6, DFNTRN(IDIM1:IDIM2,LP1), DFNCLU(IDIM1:IDIM2,LP1), 251 7 TCLD (IDIM1:IDIM2,LP1) 252 8, ALFA (IDIM1:IDIM2,LP1), ALFAU (IDIM1:IDIM2,LP1) 253 E Q U I V A L E N C E 254 1 (UDO3,UO3(IDIM1,1),DFNCLU), (URO3,UO3(IDIM1,LP2), UFNCLU) 255 2, (UDCO2,UCO2(IDIM1,1),TCLD), (URCO2,UCO2(IDIM1,LP2), TCLU) 256 3, (TDO3 ,TO3(IDIM1,1),DFNTRN),(TUO3,TO3(IDIM1,LP2), UFNTRN) 257 4, (TDCO2,TCO2(IDIM1,1) ),(TUCO2,TCO2(IDIM1,LP2) ) 258 5, (FF , ALFA ), (FFCO2 , ALFAU ), (FFO3 , TTDB1 ) 259 6, (DU , TTUB1), (DUCO2 , TUCL1 ), (DUO3 , TDCL1 ) 260 7, (PR2 , TDCL2) 261 C 262 C---COMMON FOR LOCAL DATA VARIABLES--- 263 COMMON /SWRSAV/ ABCFF(NB),PWTS(NB),CFCO2,CFO3,REFLO3,RRAYAV 264 save /SWRSAV/ 265 C D A T A 266 C 1 ABCFF / 2*4.0E-5, 0.002, 0.035, 0.377, 1.95, 9.40, 44.6, 267 C 1 190.0, 989.0, 2706.0, 39011.0 / 268 C 2, PWTS / 0.5000, 0.121416, 0.0698, 0.1558, 0.0631, 0.0362, 269 C 2 0.0243, 0.0158, 0.0087, 0.001467, 0.002342, 0.001075 / 270 C 3, CFCO2, CFO3, REFLO3, RRAYAV / 508.96, 466.64, 1.9, 0.144 / 271 C 1 ABCFF / 2*4.0E-5, .002, .035, .377, 1.95, 9.40, 44.6, 190. / 272 C 2, PWTS /.5000,.1470,.698,.1443,.0584,.0335,.0225,.0158,.0087/ Page 3 Source Listing SWR93 2014-12-17 20:47 SWR93.F 273 C 3, CFCO2, CFO3, REFLO3, RRAYAV / 508.96, 466.64, 1.9, 0.144 / 274 C 275 C CALCULATE SECANT OF ZENITH ANGLE (SECZ),FLUX PRESSURES(PP), 276 C LAYER WIDTH (DP) AND PRESSURE SCALING FACTOR (PR2). 277 DO 100 I=MYIS,MYIE 278 SECZ(I) = H35E1/SQRT(H1224E3*COSZRO(I)*COSZRO(I)+ONE) 279 PP(I,1) = ZERO 280 PP(I,LP1) = PRESS(I,LP1) 281 TMP1(I) = ONE/PRESS(I,LP1) 282 100 CONTINUE 283 DO 110 K=1,LM1 284 DO 110 I=MYIS,MYIE 285 PP(I,K+1) = HAF*(PRESS(I,K+1)+PRESS(I,K)) 286 110 CONTINUE 287 DO 120 K=1,L 288 DO 120 I=MYIS,MYIE 289 DP (I,K) = PP(I,K+1)-PP(I,K) 290 PR2(I,K) = HAF*(PP(I,K)+PP(I,K+1)) 291 120 CONTINUE 292 DO 130 K=1,L 293 DO 130 I=MYIS,MYIE 294 PR2(I,K) = PR2(I,K)*TMP1(I) 295 130 CONTINUE 296 C CALCULATE ENTERING FLUX AT THE TOP FOR EACH BAND(IN CGS UNITS) 297 DO 140 N=1,NB 298 DO 140 IP=MYIS,MYIE 299 DFNTOP(IP,N) = SSOLAR*H69766E5*COSZRO(IP)*TAUDAR(IP)*PWTS(N) 300 140 CONTINUE 301 C EXECUTE THE LACIS-HANSEN REFLECTIVITY PARAMETERIZATION 302 C FOR THE VISIBLE BAND 303 DO 150 I=MYIS,MYIE 304 RRAY(I) = HP219/(ONE+HP816*COSZRO(I)) 305 REFL(I) = RRAY(I) + (ONE-RRAY(I))*(ONE-RRAYAV)*ALVB(I)/ 306 1 (ONE-ALVD(I)*RRAYAV) 307 150 CONTINUE 308 DO 155 I=MYIS,MYIE 309 RRAY(I) = 0.104/(ONE+4.8*COSZRO(I)) 310 REFL2(I)= RRAY(I) + (ONE-RRAY(I))*(ONE-0.093)*ALVB(I)/ 311 1 (ONE-ALVD(I)*0.093) 312 155 CONTINUE 313 C CALCULATE PRESSURE-WEIGHTED OPTICAL PATHS FOR EACH LAYER 314 C IN UNITS OF CM-ATM. PRESSURE WEIGHTING IS USING PR2. 315 C DU= VALUE FOR H2O;DUCO2 FOR CO2;DUO3 FOR O3. 316 DO 160 K=1,L 317 DO 160 I=MYIS,MYIE 318 DU (I,K) = GINV*RH2O(I,K)*DP(I,K)*PR2(I,K) 319 DUCO2(I,K) = (RRCO2*GINV*CFCO2)*DP(I,K)*PR2(I,K) 320 DUO3 (I,K) = (GINV*CFO3)*QO3(I,K)*DP(I,K) 321 160 CONTINUE 322 C 323 C CALCULATE CLEAR SKY SW FLUX 324 C 325 C OBTAIN THE OPTICAL PATH FROM THE TOP OF THE ATMOSPHERE TO THE 326 C FLUX PRESSURE. ANGULAR FACTORS ARE NOW INCLUDED. UD=DOWNWARD 327 C PATH FOR H2O,WIGTH UR THE UPWARD PATH FOR H2O. CORRESPONDING 328 C QUANTITIES FOR CO2,O3 ARE UDCO2/URCO2 AND UDO3/URO3. 329 DO 200 IP=MYIS,MYIE Page 4 Source Listing SWR93 2014-12-17 20:47 SWR93.F 330 UD (IP,1) = ZERO 331 UDCO2(IP,1) = ZERO 332 UDO3 (IP,1) = ZERO 333 200 CONTINUE 334 DO 210 K=2,LP1 335 DO 210 I=MYIS,MYIE 336 UD (I,K) = UD (I,K-1)+DU (I,K-1)*SECZ(I) 337 UDCO2(I,K) = UDCO2(I,K-1)+DUCO2(I,K-1)*SECZ(I) 338 UDO3 (I,K) = UDO3 (I,K-1)+DUO3 (I,K-1)*SECZ(I) 339 210 CONTINUE 340 DO 220 IP=MYIS,MYIE 341 UR (IP,LP1) = UD (IP,LP1) 342 URCO2(IP,LP1) = UDCO2(IP,LP1) 343 URO3 (IP,LP1) = UDO3 (IP,LP1) 344 220 CONTINUE 345 DO 230 K=L,1,-1 346 DO 230 IP=MYIS,MYIE 347 UR (IP,K) = UR (IP,K+1)+DU (IP,K)*DIFFCTR 348 URCO2(IP,K) = URCO2(IP,K+1)+DUCO2(IP,K)*DIFFCTR 349 URO3 (IP,K) = URO3 (IP,K+1)+DUO3 (IP,K)*O3DIFCTR 350 230 CONTINUE 351 C CALCULATE CO2 ABSORPTIONS . THEY WILL BE USED IN NEAR INFRARED 352 C BANDS.SINCE THE ABSORPTION AMOUNT IS GIVEN (IN THE FORMULA USED 353 C BELOW, DERIVED FROM SASAMORI) IN TERMS OF THE TOTAL SOLAR FLUX, 354 C AND THE ABSORPTION IS ONLY INCLUDED IN THE NEAR IR (50 PERCENT 355 C OF THE SOLAR SPECTRUM), THE ABSORPTIONS ARE MULTIPLIED BY 2. 356 C SINCE CODE ACTUALLY REQUIRES TRANSMISSIONS, THESE ARE THE 357 C VALUES ACTUALLY STORED IN TCO2. 358 DO 240 K=1,LL 359 DO 240 I=MYIS,MYIE 360 TCO2(I,K+1)=ONE-TWO*(H235M3*EXP(HP26*LOG(UCO2(I,K+1)+H129M2)) 361 1 -H75826M4) 362 240 CONTINUE 363 C NOW CALCULATE OZONE ABSORPTIONS. THESE WILL BE USED IN 364 C THE VISIBLE BAND.JUST AS IN THE CO2 CASE, SINCE THIS BAND IS 365 C 50 PERCENT OF THE SOLAR SPECTRUM,THE ABSORPTIONS ARE MULTIPLIED 366 C BY 2. THE TRANSMISSIONS ARE STORED IN TO3. 367 HTEMP = H1036E2*H1036E2*H1036E2 368 DO 250 K=1,LL 369 DO 250 I=MYIS,MYIE 370 TO3(I,K+1)=ONE-TWO*UO3(I,K+1)* 371 1 (H1P082*EXP(HMP805*LOG(ONE+H1386E2*UO3(I,K+1)))+ 372 2 H658M2/(ONE+HTEMP*UO3(I,K+1)*UO3(I,K+1)*UO3(I,K+1))+ 373 3 H2118M2/(ONE+UO3(I,K+1)*(H42M2+H323M4*UO3(I,K+1)))) 374 250 CONTINUE 375 C START FREQUENCY LOOP (ON N) HERE 376 C 377 C--- BAND 1 (VISIBLE) INCLUDES O3 AND H2O ABSORPTION 378 DO 260 K=1,L 379 DO 260 I=MYIS,MYIE 380 TTD(I,K+1) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UD(I,K+1))) 381 TTU(I,K) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UR(I,K))) 382 DFN(I,K+1) = TTD(I,K+1)*TDO3(I,K+1) 383 UFN(I,K) = TTU(I,K)*TUO3(I,K) 384 260 CONTINUE 385 DO 270 I=MYIS,MYIE 386 DFN(I,1) = ONE Page 5 Source Listing SWR93 2014-12-17 20:47 SWR93.F 387 UFN(I,LP1) = DFN(I,LP1) 388 270 CONTINUE 389 C SCALE VISIBLE BAND FLUXES BY SOLAR FLUX AT THE TOP OF THE 390 C ATMOSPHERE (DFNTOP(I,1)) 391 C DFSW/UFSW WILL BE THE FLUXES, SUMMED OVER ALL BANDS 392 DO 280 K=1,LP1 393 DO 280 I=MYIS,MYIE 394 DFSWL(I,K) = DFN(I,K)*DFNTOP(I,1) 395 UFSWL(I,K) = REFL(I)*UFN(I,K)*DFNTOP(I,1) 396 280 CONTINUE 397 DO 285 I=MYIS,MYIE 398 GDFVB(I) = DFSWL(I,LP1)*EXP(-0.15746*SECZ(I)) 399 GDFVD(I) = ((ONE-REFL2(I))*DFSWL(I,LP1) - 400 1 (ONE-ALVB(I)) *GDFVB(I)) / (ONE-ALVD(I)) 401 GDFNB(I) = ZERO 402 GDFND(I) = ZERO 403 285 CONTINUE 404 C---NOW OBTAIN FLUXES FOR THE NEAR IR BANDS. THE METHODS ARE THE SAME 405 C AS FOR THE VISIBLE BAND, EXCEPT THAT THE REFLECTION AND 406 C TRANSMISSION COEFFICIENTS (OBTAINED BELOW) ARE DIFFERENT, AS 407 C RAYLEIGH SCATTERING NEED NOT BE CONSIDERED. 408 DO 350 N=2,NB 409 IF (N.EQ.2) THEN 410 C THE WATER VAPOR TRANSMISSION FUNCTION FOR BAND 2 IS EQUAL TO 411 C THAT OF BAND 1 (SAVED AS TTD,TTU) 412 C--- BAND 2-9 (NEAR-IR) INCLUDES O3, CO2 AND H2O ABSORPTION 413 DO 290 K=1,L 414 DO 290 I=MYIS,MYIE 415 DFN(I,K+1) = TTD(I,K+1)*TDCO2(I,K+1) 416 UFN(I,K) = TTU(I,K)*TUCO2(I,K) 417 290 CONTINUE 418 ELSE 419 C CALCULATE WATER VAPOR TRANSMISSION FUNCTIONS FOR NEAR INFRARED 420 C BANDS. INCLUDE CO2 TRANSMISSION (TDCO2/TUCO2), WHICH 421 C IS THE SAME FOR ALL INFRARED BANDS. 422 DO 300 K=1,L 423 DO 300 I=MYIS,MYIE 424 DFN(I,K+1)=EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UD(I,K+1))) 425 1 *TDCO2(I,K+1) 426 UFN(I,K)=EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UR(I,K))) 427 1 *TUCO2(I,K) 428 300 CONTINUE 429 ENDIF 430 C---AT THIS POINT,INCLUDE DFN(1),UFN(LP1), NOTING THAT DFN(1)=1 FOR 431 C ALL BANDS, AND THAT UFN(LP1)=DFN(LP1) FOR ALL BANDS. 432 DO 310 I=MYIS,MYIE 433 DFN(I,1) = ONE 434 UFN(I,LP1) = DFN(I,LP1) 435 310 CONTINUE 436 C SCALE THE PREVIOUSLY COMPUTED FLUXES BY THE FLUX AT THE TOP 437 C AND SUM OVER BANDS 438 DO 320 K=1,LP1 439 DO 320 I=MYIS,MYIE 440 DFSWL(I,K) = DFSWL(I,K) + DFN(I,K)*DFNTOP(I,N) 441 UFSWL(I,K) = UFSWL(I,K) + ALNB(I)*UFN(I,K)*DFNTOP(I,N) 442 320 CONTINUE 443 DO 330 I=MYIS,MYIE Page 6 Source Listing SWR93 2014-12-17 20:47 SWR93.F 444 GDFNB(I) = GDFNB(I) + DFN(I,LP1)*DFNTOP(I,N) 445 330 CONTINUE 446 350 CONTINUE 447 DO 360 K=1,LP1 448 DO 360 I=MYIS,MYIE 449 FSWL(I,K) = UFSWL(I,K)-DFSWL(I,K) 450 360 CONTINUE 451 DO 370 K=1,L 452 DO 370 I=MYIS,MYIE 453 HSWL(I,K)=RADCON*(FSWL(I,K+1)-FSWL(I,K))/DP(I,K) 454 370 CONTINUE 455 C 456 C---END OF FREQUENCY LOOP (OVER N) 457 C 458 C CALCULATE CLOUDY SKY SW FLUX 459 C 460 KCLDS=NCLDS(MYIS) 461 DO 400 I=MYIS1,MYIE 462 KCLDS=MAX(NCLDS(I),KCLDS) 463 400 CONTINUE 464 DO 410 K=1,LP1 465 DO 410 I=MYIS,MYIE 466 DFSWC(I,K) = DFSWL(I,K) 467 UFSWC(I,K) = UFSWL(I,K) 468 FSWC (I,K) = FSWL (I,K) 469 410 CONTINUE 470 DO 420 K=1,L 471 DO 420 I=MYIS,MYIE 472 HSWC(I,K) = HSWL(I,K) 473 420 CONTINUE 474 C******************************************************************* 475 IF (KCLDS .EQ. 0) RETURN 476 C******************************************************************* 477 DO 430 K=1,LP1 478 DO 430 I=MYIS,MYIE 479 XAMT(I,K) = CAMT(I,K) 480 430 CONTINUE 481 DO 470 I=MYIS,MYIE 482 NNCLDS = NCLDS(I) 483 CCMAX(I) = ZERO 484 IF (NNCLDS .LE. 0) GO TO 470 485 CCMAX(I) = ONE 486 DO 450 K=1,NNCLDS 487 CCMAX(I) = CCMAX(I) * (ONE - CAMT(I,K+1)) 488 450 CONTINUE 489 CCMAX(I) = ONE - CCMAX(I) 490 IF (CCMAX(I) .GT. ZERO) THEN 491 DO 460 K=1,NNCLDS 492 XAMT(I,K+1) = CAMT(I,K+1)/CCMAX(I) 493 460 CONTINUE 494 END IF 495 470 CONTINUE 496 DO 480 K=1,LP1 497 DO 480 I=MYIS,MYIE 498 FF (I,K) = DIFFCTR 499 FFCO2(I,K) = DIFFCTR 500 FFO3 (I,K) = O3DIFCTR Page 7 Source Listing SWR93 2014-12-17 20:47 SWR93.F 501 480 CONTINUE 502 DO 490 IP=MYIS,MYIE 503 JTOP = KTOPSW(IP,NCLDS(IP)+1) 504 DO 490 K=1,JTOP 505 FF (IP,K) = SECZ(IP) 506 FFCO2(IP,K) = SECZ(IP) 507 FFO3 (IP,K) = SECZ(IP) 508 490 CONTINUE 509 DO 500 I=MYIS,MYIE 510 RRAY(I) = HP219/(ONE+HP816*COSZRO(I)) 511 REFL(I) = RRAY(I) + (ONE-RRAY(I))*(ONE-RRAYAV)*ALVD(I)/ 512 1 (ONE-ALVD(I)*RRAYAV) 513 500 CONTINUE 514 DO 510 IP=MYIS,MYIE 515 UD (IP,1) = ZERO 516 UDCO2(IP,1) = ZERO 517 UDO3 (IP,1) = ZERO 518 510 CONTINUE 519 DO 520 K=2,LP1 520 DO 520 I=MYIS,MYIE 521 UD (I,K) = UD (I,K-1)+DU (I,K-1)*FF (I,K) 522 UDCO2(I,K) = UDCO2(I,K-1)+DUCO2(I,K-1)*FFCO2(I,K) 523 UDO3 (I,K) = UDO3 (I,K-1)+DUO3 (I,K-1)*FFO3 (I,K) 524 520 CONTINUE 525 DO 530 IP=MYIS,MYIE 526 UR (IP,LP1) = UD (IP,LP1) 527 URCO2(IP,LP1) = UDCO2(IP,LP1) 528 URO3 (IP,LP1) = UDO3 (IP,LP1) 529 530 CONTINUE 530 DO 540 K=L,1,-1 531 DO 540 IP=MYIS,MYIE 532 UR (IP,K) = UR (IP,K+1)+DU (IP,K)*DIFFCTR 533 URCO2(IP,K) = URCO2(IP,K+1)+DUCO2(IP,K)*DIFFCTR 534 URO3 (IP,K) = URO3 (IP,K+1)+DUO3 (IP,K)*O3DIFCTR 535 540 CONTINUE 536 DO 550 K=1,LL 537 DO 550 I=MYIS,MYIE 538 TCO2(I,K+1)=ONE-TWO*(H235M3*EXP(HP26*LOG(UCO2(I,K+1)+H129M2)) 539 1 -H75826M4) 540 550 CONTINUE 541 DO 560 K=1,LL 542 DO 560 I=MYIS,MYIE 543 TO3(I,K+1)=ONE-TWO*UO3(I,K+1)* 544 1 (H1P082*EXP(HMP805*LOG(ONE+H1386E2*UO3(I,K+1)))+ 545 2 H658M2/(ONE+HTEMP*UO3(I,K+1)*UO3(I,K+1)*UO3(I,K+1))+ 546 3 H2118M2/(ONE+UO3(I,K+1)*(H42M2+H323M4*UO3(I,K+1)))) 547 560 CONTINUE 548 C******************************************************************** 549 C---THE FIRST CLOUD IS THE GROUND; ITS PROPERTIES ARE GIVEN 550 C BY REFL (THE TRANSMISSION (0) IS IRRELEVANT FOR NOW!). 551 C******************************************************************** 552 DO 570 I=MYIS,MYIE 553 CR(I,1) = REFL(I) 554 570 CONTINUE 555 C***OBTAIN CLOUD REFLECTION AND TRANSMISSION COEFFICIENTS FOR 556 C REMAINING CLOUDS (IF ANY) IN THE VISIBLE BAND 557 C---THE MAXIMUM NO OF CLOUDS IN THE ROW (KCLDS) IS USED. THIS CREATES Page 8 Source Listing SWR93 2014-12-17 20:47 SWR93.F 558 C EXTRA WORK (MAY BE REMOVED IN A SUBSEQUENT UPDATE). 559 DO 581 I=MYIS,MYIE 560 KCLDS=NCLDS(I) 561 IF(KCLDS.EQ.0) GO TO 581 562 DO 580 KK=2,KCLDS+1 563 CR(I,KK) = CRR(I,1,KK)*XAMT(I,KK) 564 CT(I,KK) = ONE - (ONE-CTT(I,1,KK))*XAMT(I,KK) 565 580 CONTINUE 566 581 CONTINUE 567 C---OBTAIN THE PRESSURE AT THE TOP,BOTTOM AND THE THICKNESS OF 568 C "THICK" CLOUDS (THOSE AT LEAST 2 LAYERS THICK). THIS IS USED 569 C LATER IS OBTAINING FLUXES INSIDE THE THICK CLOUDS, FOR ALL 570 C FREQUENCY BANDS. 571 DO 591 I=MYIS,MYIE 572 KCLDS=NCLDS(I) 573 IF(KCLDS.EQ.0) GO TO 591 574 DO 590 KK=1,KCLDS 575 IF ((KBTMSW(I,KK+1)-1).GT.KTOPSW(I,KK+1)) THEN 576 PPTOP(I,KK)=PP(I,KTOPSW(I,KK+1)) 577 DPCLD(I,KK)=ONE/(PPTOP(I,KK)-PP(I,KBTMSW(I,KK+1))) 578 ENDIF 579 590 CONTINUE 580 591 CONTINUE 581 DO 600 K=1,L 582 DO 600 I=MYIS,MYIE 583 TTDB1(I,K+1) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UD(I,K+1))) 584 TTUB1(I,K) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UR(I,K))) 585 TTD (I,K+1) = TTDB1(I,K+1)*TDO3(I,K+1) 586 TTU (I,K) = TTUB1(I,K)*TUO3(I,K) 587 600 CONTINUE 588 DO 610 I=MYIS,MYIE 589 TTD(I,1) = ONE 590 TTU(I,LP1) = TTD(I,LP1) 591 610 CONTINUE 592 C***FOR EXECUTION OF THE CLOUD LOOP, IT IS NECESSARY TO SEPARATE OUT 593 C TRANSMISSION FCTNS AT THE TOP AND BOTTOM OF THE CLOUDS, FOR 594 C EACH BAND N. THE REQUIRED QUANTITIES ARE: 595 C TTD(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1: 596 C TTU(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1: 597 C TTD(I,KBTMSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1: 598 C AND INVERSES OF THE FIRST TWO. THE ABOVE QUANTITIES ARE 599 C STORED IN TDCL1,TUCL1,TDCL2, AND DFNTRN,UFNTRN, RESPECTIVELY, 600 C AS THEY HAVE MULTIPLE USE IN THE PGM. 601 C---FOR FIRST CLOUD LAYER (GROUND) TDCL1,TUCL1 ARE KNOWN: 602 DO 620 I=MYIS,MYIE 603 TDCL1 (I,1) = TTD(I,LP1) 604 TUCL1 (I,1) = TTU(I,LP1) 605 TDCL2 (I,1) = TDCL1(I,1) 606 DFNTRN(I,1) = ONE/TDCL1(I,1) 607 UFNTRN(I,1) = DFNTRN(I,1) 608 620 CONTINUE 609 DO 631 I=MYIS,MYIE 610 KCLDS=NCLDS(I) 611 IF(KCLDS.EQ.0) GO TO 631 612 DO 630 KK=2,KCLDS+1 613 TDCL1(I,KK) = TTD(I,KTOPSW(I,KK)) 614 TUCL1(I,KK) = TTU(I,KTOPSW(I,KK)) Page 9 Source Listing SWR93 2014-12-17 20:47 SWR93.F 615 TDCL2(I,KK) = TTD(I,KBTMSW(I,KK)) 616 630 CONTINUE 617 631 CONTINUE 618 C---COMPUTE INVERSES 619 DO 641 I=MYIS,MYIE 620 KCLDS=NCLDS(I) 621 IF(KCLDS.EQ.0) GO TO 641 622 DO 640 KK=2,KCLDS 623 DFNTRN(I,KK) = ONE/TDCL1(I,KK) 624 UFNTRN(I,KK) = ONE/TUCL1(I,KK) 625 640 CONTINUE 626 641 CONTINUE 627 C---COMPUTE THE TRANSMISSIVITY FROM THE TOP OF CLOUD (K+1) TO THE 628 C TOP OF CLOUD (K). THE CLOUD TRANSMISSION (CT) IS INCLUDED. THIS 629 C QUANTITY IS CALLED TCLU (INDEX K). ALSO, OBTAIN THE TRANSMISSIVITY 630 C FROM THE BOTTOM OF CLOUD (K+1) TO THE TOP OF CLOUD (K)(A PATH 631 C ENTIRELY OUTSIDE CLOUDS). THIS QUANTITY IS CALLED TCLD (INDEX K). 632 DO 651 I=MYIS,MYIE 633 KCLDS=NCLDS(I) 634 IF(KCLDS.EQ.0) GO TO 651 635 DO 650 KK=1,KCLDS 636 TCLU(I,KK) = TDCL1(I,KK)*DFNTRN(I,KK+1)*CT(I,KK+1) 637 TCLD(I,KK) = TDCL1(I,KK)/TDCL2(I,KK+1) 638 650 CONTINUE 639 651 CONTINUE 640 C***THE FOLLOWING IS THE RECURSION RELATION FOR ALFA: THE REFLECTION 641 C COEFFICIENT FOR A SYSTEM INCLUDING THE CLOUD IN QUESTION AND THE 642 C FLUX COMING OUT OF THE CLOUD SYSTEM INCLUDING ALL CLOUDS BELOW 643 C THE CLOUD IN QUESTION. 644 C---ALFAU IS ALFA WITHOUT THE REFLECTION OF THE CLOUD IN QUESTION 645 DO 660 I=MYIS,MYIE 646 KCLDS=NCLDS(I) 647 IF(KCLDS.EQ.0) GO TO 660 648 ALFA (I,1)=CR(I,1) 649 ALFAU(I,1)=ZERO 650 660 CONTINUE 651 C---AGAIN,EXCESSIVE CALCULATIONS-MAY BE CHANGED LATER! 652 DO 671 I=MYIS,MYIE 653 KCLDS=NCLDS(I) 654 IF(KCLDS.EQ.0) GO TO 671 655 DO 670 KK=2,KCLDS+1 656 ALFAU(I,KK)= TCLU(I,KK-1)*TCLU(I,KK-1)*ALFA(I,KK-1)/ 657 1 (ONE - TCLD(I,KK-1)*TCLD(I,KK-1)*ALFA(I,KK-1)*CR(I,KK)) 658 ALFA (I,KK)= ALFAU(I,KK)+CR(I,KK) 659 670 CONTINUE 660 671 CONTINUE 661 C CALCULATE UFN AT CLOUD TOPS AND DFN AT CLOUD BOTTOMS 662 C---NOTE THAT UFNCLU(I,KCLDS+1) GIVES THE UPWARD FLUX AT THE TOP 663 C OF THE HIGHEST REAL CLOUD (IF NCLDS(I)=KCLDS). IT GIVES THE FLUX 664 C AT THE TOP OF THE ATMOSPHERE IF NCLDS(I) < KCLDS. IN THE FIRST 665 C CASE, TDCL1 EQUALS THE TRANSMISSION FCTN TO THE TOP OF THE 666 C HIGHEST CLOUD, AS WE WANT. IN THE SECOND CASE, TDCL1=1, SO UFNCLU 667 C EQUALS ALFA. THIS IS ALSO CORRECT. 668 DO 680 I=MYIS,MYIE 669 KCLDS=NCLDS(I) 670 IF(KCLDS.EQ.0) GO TO 680 671 UFNCLU(I,KCLDS+1) = ALFA(I,KCLDS+1)*TDCL1(I,KCLDS+1) Page 10 Source Listing SWR93 2014-12-17 20:47 SWR93.F 672 DFNCLU(I,KCLDS+1) = TDCL1(I,KCLDS+1) 673 680 CONTINUE 674 C---THIS CALCULATION IS THE REVERSE OF THE RECURSION RELATION USED 675 C ABOVE 676 DO 691 I=MYIS,MYIE 677 KCLDS=NCLDS(I) 678 IF(KCLDS.EQ.0) GO TO 691 679 DO 690 KK=KCLDS,1,-1 680 UFNCLU(I,KK) = UFNCLU(I,KK+1)*ALFAU(I,KK+1)/(ALFA(I,KK+1)* 681 1 TCLU(I,KK)) 682 DFNCLU(I,KK) = UFNCLU(I,KK)/ALFA(I,KK) 683 690 CONTINUE 684 691 CONTINUE 685 DO 701 I=MYIS,MYIE 686 KCLDS=NCLDS(I) 687 IF(KCLDS.EQ.0) GO TO 701 688 DO 700 KK=1,KCLDS+1 689 UFNTRN(I,KK) = UFNCLU(I,KK)*UFNTRN(I,KK) 690 DFNTRN(I,KK) = DFNCLU(I,KK)*DFNTRN(I,KK) 691 700 CONTINUE 692 701 CONTINUE 693 C---CASE OF KK=1( FROM THE GROUND TO THE BOTTOM OF THE LOWEST CLOUD) 694 DO 720 I=MYIS,MYIE 695 KCLDS=NCLDS(I) 696 IF(KCLDS.EQ.0) GO TO 720 697 J2=KBTMSW(I,2) 698 DO 710 K=J2,LP1 699 UFN(I,K) = UFNTRN(I,1)*TTU(I,K) 700 DFN(I,K) = DFNTRN(I,1)*TTD(I,K) 701 710 CONTINUE 702 720 CONTINUE 703 C---REMAINING LEVELS (IF ANY!) 704 DO 760 I=MYIS,MYIE 705 KCLDS=NCLDS(I) 706 IF(KCLDS.EQ.0) GO TO 760 707 DO 755 KK=2,KCLDS+1 708 J1=KTOPSW(I,KK) 709 J2=KBTMSW(I,KK+1) 710 IF (J1.EQ.1) GO TO 755 711 DO 730 K=J2,J1 712 UFN(I,K) = UFNTRN(I,KK)*TTU(I,K) 713 DFN(I,K) = DFNTRN(I,KK)*TTD(I,K) 714 730 CONTINUE 715 C---FOR THE THICK CLOUDS, THE FLUX DIVERGENCE THROUGH THE CLOUD 716 C LAYER IS ASSUMED TO BE CONSTANT. THE FLUX DERIVATIVE IS GIVEN BY 717 C TEMPF (FOR THE UPWARD FLUX) AND TEMPG (FOR THE DOWNWARD FLUX). 718 J3=KBTMSW(I,KK) 719 IF ((J3-J1).GT.1) THEN 720 TEMPF = (UFNCLU(I,KK)-UFN(I,J3))*DPCLD(I,KK-1) 721 TEMPG = (DFNCLU(I,KK)-DFN(I,J3))*DPCLD(I,KK-1) 722 DO 740 K=J1+1,J3-1 723 UFN(I,K) = UFNCLU(I,KK)+TEMPF*(PP(I,K)-PPTOP(I,KK-1)) 724 DFN(I,K) = DFNCLU(I,KK)+TEMPG*(PP(I,K)-PPTOP(I,KK-1)) 725 740 CONTINUE 726 ENDIF 727 755 CONTINUE 728 760 CONTINUE Page 11 Source Listing SWR93 2014-12-17 20:47 SWR93.F 729 DO 770 I=MYIS,MYIE 730 KCLDS=NCLDS(I) 731 IF(KCLDS.EQ.0) GO TO 770 732 DO 771 K=1,LP1 733 DFSWC(I,K) = DFN(I,K)*DFNTOP(I,1) 734 UFSWC(I,K) = UFN(I,K)*DFNTOP(I,1) 735 771 CONTINUE 736 770 CONTINUE 737 DO 780 I=MYIS,MYIE 738 KCLDS=NCLDS(I) 739 IF(KCLDS.EQ.0) GO TO 780 740 TMP1(I) = ONE - CCMAX(I) 741 GDFVB(I) = TMP1(I)*GDFVB(I) 742 GDFNB(I) = TMP1(I)*GDFNB(I) 743 GDFVD(I) = TMP1(I)*GDFVD(I) + CCMAX(I)*DFSWC(I,LP1) 744 780 CONTINUE 745 C---NOW OBTAIN FLUXES FOR THE NEAR IR BANDS. THE METHODS ARE THE SAME 746 C AS FOR THE VISIBLE BAND, EXCEPT THAT THE REFLECTION AND 747 C TRANSMISSION COEFFICIENTS ARE DIFFERENT, AS 748 C RAYLEIGH SCATTERING NEED NOT BE CONSIDERED. 749 C 750 DO 1000 N=2,NB 751 CYH93 752 DO 791 I=MYIS,MYIE 753 KCLDS=NCLDS(I) 754 IF(KCLDS.EQ.0) GO TO 791 755 DO 790 K=1,KCLDS+1 756 CR(I,K) = CRR(I,N,K)*XAMT(I,K) 757 CT(I,K) = ONE - (ONE-CTT(I,N,K))*XAMT(I,K) 758 790 CONTINUE 759 791 CONTINUE 760 CYH93 761 IF (N.EQ.2) THEN 762 C THE WATER VAPOR TRANSMISSION FUNCTION FOR BAND 2 IS EQUAL TO 763 C THAT OF BAND 1 (SAVED AS TTDB1,TTUB1) 764 DO 800 I=MYIS,MYIE 765 KCLDS=NCLDS(I) 766 IF(KCLDS.EQ.0) GO TO 800 767 DO 801 KK=2,LP1 768 TTD(I,KK) = TTDB1(I,KK)*TDCO2(I,KK) 769 801 CONTINUE 770 DO 802 KK=1,L 771 TTU(I,KK) = TTUB1(I,KK)*TUCO2(I,KK) 772 802 CONTINUE 773 800 CONTINUE 774 ELSE 775 DO 810 I=MYIS,MYIE 776 KCLDS=NCLDS(I) 777 IF(KCLDS.EQ.0) GO TO 810 778 DO 811 KK=2,LP1 779 TTD(I,KK) = EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UD(I,KK))) 780 1 * TDCO2(I,KK) 781 811 CONTINUE 782 DO 812 KK=1,L 783 TTU(I,KK) = EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UR(I,KK))) 784 1 * TUCO2(I,KK) 785 812 CONTINUE Page 12 Source Listing SWR93 2014-12-17 20:47 SWR93.F 786 810 CONTINUE 787 ENDIF 788 C---AT THIS POINT,INCLUDE TTD(1),TTU(LP1), NOTING THAT TTD(1)=1 FOR 789 C ALL BANDS, AND THAT TTU(LP1)=TTD(LP1) FOR ALL BANDS. 790 DO 820 I=MYIS,MYIE 791 KCLDS=NCLDS(I) 792 IF(KCLDS.EQ.0) GO TO 820 793 TTU(I,LP1) = TTD(I,LP1) 794 TTD(I,1) = ONE 795 820 CONTINUE 796 C***FOR EXECUTION OF THE CLOUD LOOP, IT IS NECESSARY TO SEPARATE OUT 797 C TRANSMISSION FCTNS AT THE TOP AND BOTTOM OF THE CLOUDS, FOR 798 C EACH BAND N. THE REQUIRED QUANTITIES ARE: 799 C TTD(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1: 800 C TTD(I,KBTMSW(I,K),N) K RUNS FROM 2 TO NCLDS(I)+1: 801 C TTU(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1: 802 C AND INVERSES OF THE ABOVE. THE ABOVE QUANTITIES ARE STORED 803 C IN TDCL1,TDCL2,TUCL1,AND DFNTRN,UFNTRN,RESPECTIVELY, AS 804 C THEY HAVE MULTIPLE USE IN THE PGM. 805 C---FOR FIRST CLOUD LAYER (GROUND) TDCL1,TUCL1 ARE KNOWN: 806 DO 830 I=MYIS,MYIE 807 KCLDS=NCLDS(I) 808 IF(KCLDS.EQ.0) GO TO 830 809 TDCL1 (I,1) = TTD(I,LP1) 810 TUCL1 (I,1) = TTU(I,LP1) 811 TDCL2 (I,1) = TDCL1(I,1) 812 DFNTRN(I,1) = ONE/TDCL1(I,1) 813 UFNTRN(I,1) = DFNTRN(I,1) 814 830 CONTINUE 815 DO 841 I=MYIS,MYIE 816 KCLDS=NCLDS(I) 817 IF(KCLDS.EQ.0) GO TO 841 818 DO 840 KK=2,KCLDS+1 819 TDCL1(I,KK) = TTD(I,KTOPSW(I,KK)) 820 TUCL1(I,KK) = TTU(I,KTOPSW(I,KK)) 821 TDCL2(I,KK) = TTD(I,KBTMSW(I,KK)) 822 840 CONTINUE 823 841 CONTINUE 824 DO 851 I=MYIS,MYIE 825 KCLDS=NCLDS(I) 826 IF(KCLDS.EQ.0) GO TO 851 827 DO 850 KK=2,KCLDS+1 828 DFNTRN(I,KK) = ONE/TDCL1(I,KK) 829 UFNTRN(I,KK) = ONE/TUCL1(I,KK) 830 850 CONTINUE 831 851 CONTINUE 832 DO 861 I=MYIS,MYIE 833 KCLDS=NCLDS(I) 834 IF(KCLDS.EQ.0) GO TO 861 835 DO 860 KK=1,KCLDS 836 TCLU(I,KK) = TDCL1(I,KK)*DFNTRN(I,KK+1)*CT(I,KK+1) 837 TCLD(I,KK) = TDCL1(I,KK)/TDCL2(I,KK+1) 838 860 CONTINUE 839 861 CONTINUE 840 C***THE FOLLOWING IS THE RECURSION RELATION FOR ALFA: THE REFLECTION 841 C COEFFICIENT FOR A SYSTEM INCLUDING THE CLOUD IN QUESTION AND THE 842 C FLUX COMING OUT OF THE CLOUD SYSTEM INCLUDING ALL CLOUDS BELOW Page 13 Source Listing SWR93 2014-12-17 20:47 SWR93.F 843 C THE CLOUD IN QUESTION. 844 DO 870 I=MYIS,MYIE 845 KCLDS=NCLDS(I) 846 IF(KCLDS.EQ.0) GO TO 870 847 ALFA (I,1) = CR(I,1) 848 ALFAU(I,1) = ZERO 849 870 CONTINUE 850 C---AGAIN,EXCESSIVE CALCULATIONS-MAY BE CHANGED LATER! 851 DO 881 I=MYIS,MYIE 852 KCLDS=NCLDS(I) 853 IF(KCLDS.EQ.0) GO TO 881 854 DO 880 KK=2,KCLDS+1 855 ALFAU(I,KK) = TCLU(I,KK-1)*TCLU(I,KK-1)*ALFA(I,KK-1)/(ONE - 856 1 TCLD(I,KK-1)*TCLD(I,KK-1)*ALFA(I,KK-1)*CR(I,KK)) 857 ALFA (I,KK) = ALFAU(I,KK)+CR(I,KK) 858 880 CONTINUE 859 881 CONTINUE 860 C CALCULATE UFN AT CLOUD TOPS AND DFN AT CLOUD BOTTOMS 861 C---NOTE THAT UFNCLU(I,KCLDS+1) GIVES THE UPWARD FLUX AT THE TOP 862 C OF THE HIGHEST REAL CLOUD (IF NCLDS(I)=KCLDS). IT GIVES THE FLUX 863 C AT THE TOP OF THE ATMOSPHERE IF NCLDS(I) < KCLDS. IT THE FIRST 864 C CASE, TDCL1 EQUALS THE TRANSMISSION FCTN TO THE TOP OF THE 865 C HIGHEST CLOUD, AS WE WANT. IN THE SECOND CASE, TDCL1=1, SO UFNCLU 866 C EQUALS ALFA. THIS IS ALSO CORRECT. 867 DO 890 I=MYIS,MYIE 868 KCLDS=NCLDS(I) 869 IF(KCLDS.EQ.0) GO TO 890 870 UFNCLU(I,KCLDS+1) = ALFA(I,KCLDS+1)*TDCL1(I,KCLDS+1) 871 DFNCLU(I,KCLDS+1) = TDCL1(I,KCLDS+1) 872 890 CONTINUE 873 DO 901 I=MYIS,MYIE 874 KCLDS=NCLDS(I) 875 IF(KCLDS.EQ.0) GO TO 901 876 DO 900 KK=KCLDS,1,-1 877 ! 878 !--- Ferrier, 6/17/02: Emergency change to eliminate problematic 879 ! features of unrealistically small cloud amounts 880 ! 881 DENOM=ALFA(I,KK+1)*TCLU(I,KK) 882 IF (DENOM .GT. 1.E-15) THEN 883 UFNCLU(I,KK)=UFNCLU(I,KK+1)*ALFAU(I,KK+1)/DENOM 884 ELSE 885 c print * 886 c &, ' xnum=',UFNCLU(I,KK+1)*ALFAU(I,KK+1) 887 c &, ' xden=',DENOM 888 c &, ' ALFA(I,KK+1)=',ALFA(I,KK+1) 889 c &, ' TCLU(I,KK)=',TCLU(I,KK) 890 c &, ' UFNCLU(I,KK+1)=',UFNCLU(I,KK+1) 891 c &, ' ALFAU(I,KK+1)=',ALFAU(I,KK+1) 892 c &, ' UFNCLU(I,KK)=',UFNCLU(I,KK) 893 c &, ' ALFA(I,KK)=',ALFA(I,KK) 894 c &, ' I=',I 895 c &, ' KK+1=',KK+1 896 UFNCLU(I,KK)=0. 897 ENDIF 898 DFNCLU(I,KK) = UFNCLU(I,KK)/ALFA(I,KK) 899 900 CONTINUE Page 14 Source Listing SWR93 2014-12-17 20:47 SWR93.F 900 901 CONTINUE 901 C NOW OBTAIN DFN AND UFN FOR LEVELS BETWEEN THE CLOUDS 902 DO 911 I=MYIS,MYIE 903 KCLDS=NCLDS(I) 904 IF(KCLDS.EQ.0) GO TO 911 905 DO 910 KK=1,KCLDS+1 906 UFNTRN(I,KK) = UFNCLU(I,KK)*UFNTRN(I,KK) 907 DFNTRN(I,KK) = DFNCLU(I,KK)*DFNTRN(I,KK) 908 910 CONTINUE 909 911 CONTINUE 910 DO 930 I=MYIS,MYIE 911 KCLDS=NCLDS(I) 912 IF(KCLDS.EQ.0) GO TO 930 913 J2=KBTMSW(I,2) 914 DO 920 K=J2,LP1 915 UFN(I,K) = UFNTRN(I,1)*TTU(I,K) 916 DFN(I,K) = DFNTRN(I,1)*TTD(I,K) 917 920 CONTINUE 918 930 CONTINUE 919 DO 970 I=MYIS,MYIE 920 KCLDS=NCLDS(I) 921 IF(KCLDS.EQ.0) GO TO 970 922 DO 965 KK=2,KCLDS+1 923 J1 = KTOPSW(I,KK) 924 J2 = KBTMSW(I,KK+1) 925 IF (J1.EQ.1) GO TO 965 926 DO 940 K=J2,J1 927 UFN(I,K) = UFNTRN(I,KK)*TTU(I,K) 928 DFN(I,K) = DFNTRN(I,KK)*TTD(I,K) 929 940 CONTINUE 930 J3 = KBTMSW(I,KK) 931 IF ((J3-J1).GT.1) THEN 932 TEMPF = (UFNCLU(I,KK)-UFN(I,J3))*DPCLD(I,KK-1) 933 TEMPG = (DFNCLU(I,KK)-DFN(I,J3))*DPCLD(I,KK-1) 934 DO 950 K=J1+1,J3-1 935 UFN(I,K) = UFNCLU(I,KK)+TEMPF*(PP(I,K)-PPTOP(I,KK-1)) 936 DFN(I,K) = DFNCLU(I,KK)+TEMPG*(PP(I,K)-PPTOP(I,KK-1)) 937 950 CONTINUE 938 ENDIF 939 965 CONTINUE 940 970 CONTINUE 941 DO 980 I=MYIS,MYIE 942 KCLDS=NCLDS(I) 943 IF(KCLDS.EQ.0) GO TO 980 944 DO 981 K=1,LP1 945 DFSWC(I,K) = DFSWC(I,K) + DFN(I,K)*DFNTOP(I,N) 946 UFSWC(I,K) = UFSWC(I,K) + UFN(I,K)*DFNTOP(I,N) 947 981 CONTINUE 948 980 CONTINUE 949 DO 990 I=MYIS,MYIE 950 KCLDS=NCLDS(I) 951 IF(KCLDS.EQ.0) GO TO 990 952 GDFND(I) = GDFND(I) + CCMAX(I)*DFN(I,LP1)*DFNTOP(I,N) 953 990 CONTINUE 954 1000 CONTINUE 955 DO 1100 I=MYIS,MYIE 956 KCLDS=NCLDS(I) Page 15 Source Listing SWR93 2014-12-17 20:47 SWR93.F 957 IF(KCLDS.EQ.0) GO TO 1100 958 DO 1101 K=1,LP1 959 DFSWC(I,K) = TMP1(I)*DFSWL(I,K) + CCMAX(I)*DFSWC(I,K) 960 UFSWC(I,K) = TMP1(I)*UFSWL(I,K) + CCMAX(I)*UFSWC(I,K) 961 1101 CONTINUE 962 1100 CONTINUE 963 DO 1200 I=MYIS,MYIE 964 KCLDS=NCLDS(I) 965 IF(KCLDS.EQ.0) GO TO 1200 966 DO 1201 KK=1,LP1 967 FSWC(I,KK) = UFSWC(I,KK)-DFSWC(I,KK) 968 1201 CONTINUE 969 1200 CONTINUE 970 DO 1250 I=MYIS,MYIE 971 KCLDS=NCLDS(I) 972 IF(KCLDS.EQ.0) GO TO 1250 973 DO 1251 KK=1, L 974 HSWC(I,KK) = RADCON*(FSWC(I,KK+1)-FSWC(I,KK))/DP(I,KK) 975 1251 CONTINUE 976 1250 CONTINUE 977 RETURN 978 END ENTRY POINTS Name swr93_ Page 16 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 120 115 1000 Label 792 588 110 Label 124 121,122 1100 Label 800 793,795 1101 Label 799 796 120 Label 129 125,126 1200 Label 807 801,803 1201 Label 806 804 1250 Label 814 808,810 1251 Label 813 811 130 Label 133 130,131 140 Label 138 135,136 150 Label 145 141 155 Label 150 146 160 Label 159 154,155 200 Label 171 167 210 Label 177 172,173 220 Label 182 178 230 Label 188 183,184 240 Label 200 196,197 250 Label 212 206,207 260 Label 222 216,217 270 Label 226 223 280 Label 234 230,231 285 Label 241 235 290 Label 255 251,252 300 Label 266 260,261 310 Label 273 270 320 Label 280 276,277 330 Label 283 281 350 Label 284 246 360 Label 288 285,286 370 Label 292 289,290 400 Label 301 299 410 Label 307 302,303 420 Label 311 308,309 430 Label 318 315,316 450 Label 326 324 460 Label 331 329 470 Label 333 319,322 480 Label 339 334,335 490 Label 346 340,342 500 Label 351 347 510 Label 356 352 520 Label 362 357,358 530 Label 367 363 540 Label 373 368,369 550 Label 378 374,375 560 Label 385 379,380 570 Label 392 390 580 Label 403 400 Page 17 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 581 Label 404 397,399 590 Label 417 412 591 Label 418 409,411 600 Label 425 419,420 610 Label 429 426 620 Label 446 440 630 Label 454 450 631 Label 455 447,449 640 Label 463 460 641 Label 464 457,459 650 Label 476 473 651 Label 477 470,472 660 Label 488 483,485 670 Label 497 493 671 Label 498 490,492 680 Label 511 506,508 690 Label 521 517 691 Label 522 514,516 700 Label 529 526 701 Label 530 523,525 710 Label 539 536 720 Label 540 532,534 730 Label 552 549 740 Label 563 560 755 Label 565 545,548 760 Label 566 542,544 770 Label 574 567,569 771 Label 573 570 780 Label 582 575,577 790 Label 596 593 791 Label 597 590,592 800 Label 611 602,604 801 Label 607 605 802 Label 610 608 810 Label 624 613,615 811 Label 619 616 812 Label 623 620 820 Label 633 628,630 830 Label 652 644,646 840 Label 660 656 841 Label 661 653,655 850 Label 668 665 851 Label 669 662,664 860 Label 676 673 861 Label 677 670,672 870 Label 687 682,684 880 Label 696 692 881 Label 697 689,691 890 Label 710 705,707 900 Label 737 714 901 Label 738 711,713 910 Label 746 743 911 Label 747 740,742 920 Label 755 752 930 Label 756 748,750 Page 18 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 940 Label 767 764 950 Label 775 772 965 Label 777 760,763 970 Label 778 757,759 980 Label 786 779,781 981 Label 785 782 990 Label 791 787,789 ABCFF Scalar 101 R(4) 4 1 12 COM 218,219,262,264,421,422,617,621 ALFA Local 90 R(4) 4 2 2070 486,494,495,496,509,518,520,685,69 3,694,695,708,719,736 ALFAU Local 90 R(4) 4 2 2070 487,494,496,518,686,693,695,721 ALNB Dummy 5 R(4) 4 1 45 ARG,INOUT 279 ALND Dummy 5 R(4) 4 1 45 ARG,INOUT ALVB Dummy 5 R(4) 4 1 45 ARG,INOUT 143,148,238 ALVD Dummy 5 R(4) 4 1 45 ARG,INOUT 144,149,238,349,350 AMOLWT Scalar 1 R(4) 4 scalar COM CAMT Dummy 4 R(4) 4 2 2070 ARG,INOUT 317,325,330 CCMAX Local 64 R(4) 4 1 45 321,323,325,327,328,330,578,581,79 0,797,798 CFCO2 Scalar 101 R(4) 4 scalar COM 157 CFO3 Scalar 101 R(4) 4 scalar COM 158 COSZRO Dummy 3 R(4) 4 1 45 ARG,INOUT 116,137,142,147,348 CR Local 77 R(4) 4 2 2070 391,401,486,495,496,594,685,694,69 5 CRR Dummy 4 R(4) 4 3 24840 ARG,INOUT 401,594 CSUBP Scalar 1 R(4) 4 scalar COM CT Local 79 R(4) 4 2 2070 402,474,595,674 CTT Dummy 4 R(4) 4 3 24840 ARG,INOUT 402,595 DENOM Local 719 R(4) 4 scalar 719,720,721 DFN Local 76 R(4) 4 2 2070 220,224,225,232,253,262,271,272,27 8,282,538,551,559,562,571,754,766, 771,774,783,790 DFNCLU Local 88 R(4) 4 2 2070 510,520,528,559,562,709,736,745,77 1,774 DFNTOP Local 60 R(4) 4 2 540 137,232,233,278,279,282,571,572,78 3,784,790 DFNTRN Local 88 R(4) 4 2 2070 444,445,461,474,528,538,551,650,65 1,666,674,745,754,766 DFSWC Dummy 1 R(4) 4 2 2070 ARG,INOUT 304,571,581,783,797,805 DFSWL Dummy 2 R(4) 4 2 2070 ARG,INOUT 232,236,237,278,287,304,797 DIFFCTR Scalar 1 R(4) 4 scalar COM 185,186,336,337,370,371 DP Local 54 R(4) 4 2 2070 127,156,157,158,291,812 DPCLD Local 80 R(4) 4 2 2070 415,558,559,770,771 DU Local 56 R(4) 4 2 2070 156,174,185,359,370 DUCO2 Local 56 R(4) 4 2 2070 157,175,186,360,371 DUO3 Local 57 R(4) 4 2 2070 158,176,187,361,372 EIGHT Scalar 7 R(4) 4 scalar COM EXP Func 198 scalar 198,209,218,219,236,262,264,376,38 2,421,422,617,621 FF Local 58 R(4) 4 2 2070 336,343,359 FFCO2 Local 58 R(4) 4 2 2070 337,344,360 FFO3 Local 59 R(4) 4 2 2070 338,345,361 FIFTY Scalar 7 R(4) 4 scalar COM 218,219,262,264,421,422,617,621 FIVE Scalar 7 R(4) 4 scalar COM FOUR Scalar 8 R(4) 4 scalar COM Page 19 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References FSWC Dummy 1 R(4) 4 2 2070 ARG,INOUT 306,805,812 FSWL Dummy 1 R(4) 4 2 2070 ARG,INOUT 287,291,306 G Scalar 1 R(4) 4 scalar COM G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GDFNB Dummy 5 R(4) 4 1 45 ARG,INOUT 239,282,580 GDFND Dummy 5 R(4) 4 1 45 ARG,INOUT 240,790 GDFVB Dummy 5 R(4) 4 1 45 ARG,INOUT 236,238,579 GDFVD Dummy 5 R(4) 4 1 45 ARG,INOUT 237,581 GINV Scalar 5 R(4) 4 scalar COM 156,157,158 GLB_TABLE Common 37 128 GP0INV Scalar 5 R(4) 4 scalar COM GRAVDR Scalar 1 R(4) 4 scalar COM H101M16 Scalar 40 R(4) 4 scalar COM H102M5 Scalar 39 R(4) 4 scalar COM H1036E2 Scalar 38 R(4) 4 scalar COM 205 H114M11 Scalar 22 R(4) 4 scalar COM H1174M7 Scalar 17 R(4) 4 scalar COM H11M10 Scalar 19 R(4) 4 scalar COM H11M11 Scalar 22 R(4) 4 scalar COM H1224E3 Scalar 11 R(4) 4 scalar COM 116 H1226E1 Scalar 33 R(4) 4 scalar COM H128M5 Scalar 16 R(4) 4 scalar COM H129M2 Scalar 37 R(4) 4 scalar COM 198,376 H12M12 Scalar 26 R(4) 4 scalar COM H12M13 Scalar 27 R(4) 4 scalar COM H12M31 Scalar 29 R(4) 4 scalar COM H135M13 Scalar 27 R(4) 4 scalar COM H1386E2 Scalar 37 R(4) 4 scalar COM 209,382 H1439M5 Scalar 16 R(4) 4 scalar COM H14M10 Scalar 18 R(4) 4 scalar COM H14M11 Scalar 22 R(4) 4 scalar COM H14M12 Scalar 26 R(4) 4 scalar COM H14M14 Scalar 27 R(4) 4 scalar COM H14M30 Scalar 29 R(4) 4 scalar COM H15E2 Scalar 12 R(4) 4 scalar COM H15M11 Scalar 22 R(4) 4 scalar COM H15M14 Scalar 27 R(4) 4 scalar COM H15M5 Scalar 41 R(4) 4 scalar COM H161E1 Scalar 40 R(4) 4 scalar COM H165E5 Scalar 10 R(4) 4 scalar COM H16E1 Scalar 40 R(4) 4 scalar COM H16M12 Scalar 26 R(4) 4 scalar COM H181E1 Scalar 12 R(4) 4 scalar COM H18E1 Scalar 12 R(4) 4 scalar COM H18E3 Scalar 41 R(4) 4 scalar COM H18M11 Scalar 22 R(4) 4 scalar COM H1E11 Scalar 9 R(4) 4 scalar COM H1E13 Scalar 9 R(4) 4 scalar COM H1E15 Scalar 9 R(4) 4 scalar COM H1E4 Scalar 10 R(4) 4 scalar COM H1E6 Scalar 32 R(4) 4 scalar COM H1E8 Scalar 9 R(4) 4 scalar COM H1M10 Scalar 19 R(4) 4 scalar COM H1M11 Scalar 23 R(4) 4 scalar COM Page 20 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References H1M13 Scalar 27 R(4) 4 scalar COM H1M16 Scalar 43 R(4) 4 scalar COM H1M17 Scalar 28 R(4) 4 scalar COM H1M18 Scalar 28 R(4) 4 scalar COM H1M19 Scalar 28 R(4) 4 scalar COM H1M2 Scalar 32 R(4) 4 scalar COM H1M20 Scalar 28 R(4) 4 scalar COM H1M21 Scalar 28 R(4) 4 scalar COM H1M22 Scalar 28 R(4) 4 scalar COM H1M23 Scalar 28 R(4) 4 scalar COM H1M24 Scalar 29 R(4) 4 scalar COM H1M3 Scalar 15 R(4) 4 scalar COM H1M4 Scalar 15 R(4) 4 scalar COM H1M5 Scalar 16 R(4) 4 scalar COM H1M6 Scalar 17 R(4) 4 scalar COM H1M60 Scalar 30 R(4) 4 scalar COM H1M8 Scalar 18 R(4) 4 scalar COM H1P082 Scalar 37 R(4) 4 scalar COM 209,382 H1P25892 Scalar 13 R(4) 4 scalar COM H1P4 Scalar 13 R(4) 4 scalar COM H1P41819 Scalar 35 R(4) 4 scalar COM H1P4387 Scalar 13 R(4) 4 scalar COM H1P8 Scalar 13 R(4) 4 scalar COM H2075E3 Scalar 11 R(4) 4 scalar COM H20788E3 Scalar 10 R(4) 4 scalar COM H2118M2 Scalar 38 R(4) 4 scalar COM 211,384 H21M12 Scalar 25 R(4) 4 scalar COM H21M31 Scalar 29 R(4) 4 scalar COM H235M3 Scalar 36 R(4) 4 scalar COM 198,376 H23E2 Scalar 12 R(4) 4 scalar COM H23M10 Scalar 18 R(4) 4 scalar COM H23M11 Scalar 21 R(4) 4 scalar COM H24E3 Scalar 10 R(4) 4 scalar COM H24M11 Scalar 21 R(4) 4 scalar COM H24M12 Scalar 25 R(4) 4 scalar COM H25452M6 Scalar 17 R(4) 4 scalar COM H257M8 Scalar 18 R(4) 4 scalar COM H25E2 Scalar 41 R(4) 4 scalar COM H25M31 Scalar 29 R(4) 4 scalar COM H26E2 Scalar 35 R(4) 4 scalar COM H26M30 Scalar 29 R(4) 4 scalar COM H28571M2 Scalar 43 R(4) 4 scalar COM H285M4 Scalar 15 R(4) 4 scalar COM H28E1 Scalar 46 R(4) 4 scalar COM H28M11 Scalar 21 R(4) 4 scalar COM H28M12 Scalar 25 R(4) 4 scalar COM H29316E2 Scalar 33 R(4) 4 scalar COM H2945E2 Scalar 11 R(4) 4 scalar COM H29M12 Scalar 25 R(4) 4 scalar COM H2E2 Scalar 35 R(4) 4 scalar COM H2E6 Scalar 32 R(4) 4 scalar COM H2M11 Scalar 22 R(4) 4 scalar COM H2P5 Scalar 13 R(4) 4 scalar COM H2P8 Scalar 12 R(4) 4 scalar COM H2P9 Scalar 12 R(4) 4 scalar COM Page 21 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References H3082E2 Scalar 11 R(4) 4 scalar COM H3116E1 Scalar 33 R(4) 4 scalar COM H323M4 Scalar 38 R(4) 4 scalar COM 211,384 H327M8 Scalar 18 R(4) 4 scalar COM H32M11 Scalar 21 R(4) 4 scalar COM H35E1 Scalar 12 R(4) 4 scalar COM 116 H35M11 Scalar 21 R(4) 4 scalar COM H36M13 Scalar 26 R(4) 4 scalar COM H37412M5 Scalar 16 R(4) 4 scalar COM H37M11 Scalar 20 R(4) 4 scalar COM H37M12 Scalar 25 R(4) 4 scalar COM H38M12 Scalar 24 R(4) 4 scalar COM H391M7 Scalar 17 R(4) 4 scalar COM H394M5 Scalar 16 R(4) 4 scalar COM H3E2 Scalar 11 R(4) 4 scalar COM H3M11 Scalar 21 R(4) 4 scalar COM H3M12 Scalar 25 R(4) 4 scalar COM H3M14 Scalar 27 R(4) 4 scalar COM H3M3 Scalar 40 R(4) 4 scalar COM H3M4 Scalar 44 R(4) 4 scalar COM H3P5 Scalar 41 R(4) 4 scalar COM H3P6 Scalar 12 R(4) 4 scalar COM H41666M2 Scalar 33 R(4) 4 scalar COM H42M11 Scalar 20 R(4) 4 scalar COM H42M2 Scalar 38 R(4) 4 scalar COM 211,384 H44194M2 Scalar 35 R(4) 4 scalar COM H44871M2 Scalar 15 R(4) 4 scalar COM H44M11 Scalar 20 R(4) 4 scalar COM H44M12 Scalar 24 R(4) 4 scalar COM H451M6 Scalar 39 R(4) 4 scalar COM H45M12 Scalar 24 R(4) 4 scalar COM H45M32 Scalar 30 R(4) 4 scalar COM H46M13 Scalar 26 R(4) 4 scalar COM H488E4 Scalar 10 R(4) 4 scalar COM H48M11 Scalar 20 R(4) 4 scalar COM H4999M6 Scalar 17 R(4) 4 scalar COM H4E5 Scalar 9 R(4) 4 scalar COM H4M12 Scalar 24 R(4) 4 scalar COM H4M33 Scalar 30 R(4) 4 scalar COM H53M11 Scalar 20 R(4) 4 scalar COM H559M3 Scalar 15 R(4) 4 scalar COM H55M32 Scalar 30 R(4) 4 scalar COM H5725E4 Scalar 10 R(4) 4 scalar COM H5E2 Scalar 11 R(4) 4 scalar COM H625M2 Scalar 34 R(4) 4 scalar COM H62M12 Scalar 24 R(4) 4 scalar COM H62M34 Scalar 30 R(4) 4 scalar COM H658M2 Scalar 38 R(4) 4 scalar COM 210,383 H65M12 Scalar 23 R(4) 4 scalar COM H67390E2 Scalar 39 R(4) 4 scalar COM H6938M5 Scalar 16 R(4) 4 scalar COM H69766E5 Scalar 36 R(4) 4 scalar COM 137 H6M12 Scalar 24 R(4) 4 scalar COM H6P08108 Scalar 42 R(4) 4 scalar COM H71E26 Scalar 9 R(4) 4 scalar COM Page 22 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References H72M11 Scalar 20 R(4) 4 scalar COM H74M12 Scalar 23 R(4) 4 scalar COM H75826M4 Scalar 37 R(4) 4 scalar COM 199,377 H77M11 Scalar 19 R(4) 4 scalar COM H77M12 Scalar 23 R(4) 4 scalar COM H7M6 Scalar 17 R(4) 4 scalar COM H8121E1 Scalar 35 R(4) 4 scalar COM H82M11 Scalar 19 R(4) 4 scalar COM H83E26 Scalar 9 R(4) 4 scalar COM H83M11 Scalar 19 R(4) 4 scalar COM H8725M8 Scalar 18 R(4) 4 scalar COM H8M11 Scalar 19 R(4) 4 scalar COM H8M13 Scalar 26 R(4) 4 scalar COM H93M12 Scalar 23 R(4) 4 scalar COM H96M12 Scalar 23 R(4) 4 scalar COM H987M4 Scalar 15 R(4) 4 scalar COM H9M32 Scalar 30 R(4) 4 scalar COM H9P94 Scalar 34 R(4) 4 scalar COM HAF Scalar 8 R(4) 4 scalar COM 123,128 HCON Common 7 872 HM13EZ Scalar 31 R(4) 4 scalar COM HM1597E1 Scalar 41 R(4) 4 scalar COM HM161E1 Scalar 40 R(4) 4 scalar COM HM1797E1 Scalar 34 R(4) 4 scalar COM HM181E1 Scalar 31 R(4) 4 scalar COM HM19EZ Scalar 31 R(4) 4 scalar COM HM1E1 Scalar 31 R(4) 4 scalar COM HM1E2 Scalar 31 R(4) 4 scalar COM HM1EZ Scalar 35 R(4) 4 scalar COM 218,219,262,264,421,422,617,621 HM2M2 Scalar 33 R(4) 4 scalar COM HM6666M2 Scalar 32 R(4) 4 scalar COM HM8E1 Scalar 45 R(4) 4 scalar COM HMP5 Scalar 33 R(4) 4 scalar COM HMP575 Scalar 31 R(4) 4 scalar COM HMP66667 Scalar 32 R(4) 4 scalar COM HMP805 Scalar 42 R(4) 4 scalar COM 209,382 HNINETY Scalar 7 R(4) 4 scalar COM HP1 Scalar 14 R(4) 4 scalar COM HP118666 Scalar 41 R(4) 4 scalar COM HP144 Scalar 36 R(4) 4 scalar COM HP166666 Scalar 32 R(4) 4 scalar COM HP219 Scalar 36 R(4) 4 scalar COM 142,348 HP228 Scalar 34 R(4) 4 scalar COM HP26 Scalar 36 R(4) 4 scalar COM 198,376 HP369 Scalar 14 R(4) 4 scalar COM HP3795 Scalar 39 R(4) 4 scalar COM HP5048 Scalar 39 R(4) 4 scalar COM HP518 Scalar 13 R(4) 4 scalar COM HP526315 Scalar 42 R(4) 4 scalar COM HP6 Scalar 34 R(4) 4 scalar COM HP602409 Scalar 42 R(4) 4 scalar COM HP60241 Scalar 34 R(4) 4 scalar COM HP8 Scalar 13 R(4) 4 scalar COM HP805 Scalar 37 R(4) 4 scalar COM HP816 Scalar 36 R(4) 4 scalar COM 142,348 Page 23 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References HSWC Dummy 1 R(4) 4 2 2070 ARG,INOUT 310,812 HSWL Dummy 1 R(4) 4 2 2070 ARG,INOUT 291,310 HTEMP Local 205 R(4) 4 scalar 205,210,383 HTMG Scalar 52 R(4) 4 3 4183695 COM HUNDRED Scalar 7 R(4) 4 scalar COM I Local 115 I(4) 4 scalar 115,116,117,118,119,122,123,126,12 7,128,131,132,141,142,143,144,146, 147,148,149,155,156,157,158,173,17 4,175,176,197,198,207,208,209,210, 211,217,218,219,220,221,223,224,22 5,231,232,233,235,236,237,238,239, 240,252,253,254,261,262,263,264,26 5,270,271,272,277,278,279,281,282, 286,287,290,291,299,300,303,304,30 5,306,309,310,316,317,319,320,321, 323,325,327,328,330,335,336,337,33 8,347,348,349,350,358,359,360,361, 375,376,380,381,382,383,384,390,39 1,397,398,401,402,409,410,413,414, 415,420,421,422,423,424,426,427,42 8,440,441,442,443,444,445,447,448, 451,452,453,457,458,461,462,470,47 1,474,475,483,484,486,487,490,491, 494,495,496,506,507,509,510,514,51 5,518,519,520,523,524,527,528,532, 533,535,537,538,542,543,546,547,55 0,551,556,558,559,561,562,567,568, 571,572,575,576,578,579,580,581,59 0,591,594,595,602,603,606,609,613, 614,617,618,621,622,628,629,631,63 2,644,645,647,648,649,650,651,653, 654,657,658,659,662,663,666,667,67 0,671,674,675,682,683,685,686,689, 690,693,694,695,705,706,708,709,71 1,712,719,721,734,736,740,741,744, 745,748,749,751,753,754,757,758,76 1,762,765,766,768,770,771,773,774, 779,780,783,784,787,788,790,793,79 4,797,798,801,802,805,808,809,812 IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDIM1 Param 19 I(4) 4 scalar 36,37,38,39,40,41,42,43,44,45,46,4 7,48,49,50,51,54,55,56,57,58,59,60 ,61,62,63,64,66,67,68,69,70,71,72, 73,74,76,77,78,79,80,83,84,85,86,8 7,88,89,90,92,93,94,95 IDIM2 Param 19 I(4) 4 scalar 36,37,38,39,40,41,42,43,44,45,46,4 7,48,49,50,51,54,55,56,57,58,59,60 ,61,62,63,64,66,67,68,69,70,71,72, 73,74,76,77,78,79,80,83,84,85,86,8 7,88,89,90 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 Page 24 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 15,16,19,43,44,45,46,47,51,52,57 IMAX Param 15 I(4) 4 scalar 28 INLTE Param 27 I(4) 4 scalar 27 INLTEP Param 27 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 IP Local 136 I(4) 4 scalar 136,137,167,168,169,170,178,179,18 0,181,184,185,186,187,340,341,343, 344,345,352,353,354,355,363,364,36 5,366,369,370,371,372 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 J1 Local 546 I(4) 4 scalar 546,548,549,557,560,761,763,764,76 9,772 J2 Local 535 I(4) 4 scalar 535,536,547,549,751,752,762,764 J3 Local 556 I(4) 4 scalar 556,557,558,559,560,768,769,770,77 1,772 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 JTOP Local 341 I(4) 4 scalar 341,342 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 Page 25 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References JTPAD5 Scalar 26 I(4) 4 scalar COM K Local 121 I(4) 4 scalar 121,123,125,127,128,130,132,154,15 6,157,158,172,174,175,176,183,185, 186,187,196,198,206,208,209,210,21 1,216,218,219,220,221,230,232,233, 251,253,254,260,262,263,264,265,27 6,278,279,285,287,289,291,302,304, 305,306,308,310,315,317,324,325,32 9,330,334,336,337,338,342,343,344, 345,357,359,360,361,368,370,371,37 2,374,376,379,381,382,383,384,419, 421,422,423,424,536,537,538,549,55 0,551,560,561,562,570,571,572,593, 594,595,752,753,754,764,765,766,77 2,773,774,782,783,784,796,797,798 KBTMSW Dummy 4 I(4) 4 2 2070 ARG,INOUT 413,415,453,535,547,556,659,751,76 2,768 KCLDS Local 298 I(4) 4 scalar 298,300,313,398,399,400,410,411,41 2,448,449,450,458,459,460,471,472, 473,484,485,491,492,493,507,508,50 9,510,515,516,517,524,525,526,533, 534,543,544,545,568,569,576,577,59 1,592,593,603,604,614,615,629,630, 645,646,654,655,656,663,664,665,67 1,672,673,683,684,690,691,692,706, 707,708,709,712,713,714,741,742,74 3,749,750,758,759,760,780,781,788, 789,794,795,802,803,809,810 KK Local 400 I(4) 4 scalar 400,401,402,412,413,414,415,450,45 1,452,453,460,461,462,473,474,475, 493,494,495,496,517,518,519,520,52 6,527,528,545,546,547,550,551,556, 558,559,561,562,605,606,608,609,61 6,617,618,620,621,622,656,657,658, 659,665,666,667,673,674,675,692,69 3,694,695,714,719,721,734,736,743, 744,745,760,761,762,765,766,768,77 0,771,773,774,804,805,811,812 KO2 Param 30 I(4) 4 scalar 31 KO21 Param 31 I(4) 4 scalar KO2M Param 31 I(4) 4 scalar KTOPSW Dummy 4 I(4) 4 2 2070 ARG,INOUT 341,413,414,451,452,546,657,658,76 1 L Param 14 I(4) 4 scalar 18,19,20,23,25,42,43,125,130,154,1 83,216,251,260,289,308,368,419,608 ,620,811 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LL Param 20 I(4) 4 scalar 20,21,196,206,374,379 LL3P Param 25 I(4) 4 scalar 28 LL3PI Param 28 I(4) 4 scalar LLM1 Param 21 I(4) 4 scalar LLM2 Param 21 I(4) 4 scalar LLM3 Param 21 I(4) 4 scalar LLP1 Param 20 I(4) 4 scalar 28 Page 26 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References LLP1I Param 28 I(4) 4 scalar LLP2 Param 20 I(4) 4 scalar 67,69,71,73 LLP3 Param 20 I(4) 4 scalar LM Param 3 I(4) 4 scalar 14,52 LM1 Param 19 I(4) 4 scalar 121 LM2 Param 19 I(4) 4 scalar LM3 Param 19 I(4) 4 scalar LOG Func 198 scalar 198,209,376,382 LP1 Param 18 I(4) 4 scalar 22,23,24,28,36,37,38,39,40,41,42,4 4,45,54,55,56,57,58,59,64,66,67,68 ,69,70,71,72,73,74,76,77,78,79,80, 83,84,85,86,87,88,89,90,118,119,17 2,179,180,181,225,230,236,237,272, 276,282,285,302,315,334,357,364,36 5,366,428,441,442,536,570,581,605, 616,631,647,648,752,782,790,796,80 4 LP121 Param 24 I(4) 4 scalar LP1I Param 28 I(4) 4 scalar LP1M Param 22 I(4) 4 scalar 22 LP1M1 Param 22 I(4) 4 scalar LP1V Param 23 I(4) 4 scalar LP2 Param 18 I(4) 4 scalar 92,93,94,95 LP3 Param 18 I(4) 4 scalar LSM Param 3 I(4) 4 scalar MAPPINGS Common 56 5024 MAX Func 300 scalar 300 MIN Func 218 scalar 218,219,262,264,421,422,617,621 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 115,122,126,131,136,141,146,155,16 7,173,178,184,197,207,217,223,231, 235,252,261,270,277,281,286,290,29 9,303,309,316,319,335,340,347,352, 358,363,369,375,380,390,397,409,42 0,426,440,447,457,470,483,490,506, 514,523,532,542,567,575,590,602,61 3,628,644,653,662,670,682,689,705, 711,740,748,757,779,787,793,801,80 8 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 MYIS Scalar 4 I(4) 4 scalar COM 115,122,126,131,136,141,146,155,16 Page 27 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 7,173,178,184,197,207,217,223,231, 235,252,261,270,277,281,286,290,29 8,303,309,316,319,335,340,347,352, 358,363,369,375,380,390,397,409,42 0,426,440,447,457,470,483,490,506, 514,523,532,542,567,575,590,602,61 3,628,644,653,662,670,682,689,705, 711,740,748,757,779,787,793,801,80 8 MYIS1 Scalar 4 I(4) 4 scalar COM 299 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 MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM Page 28 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 N Local 135 I(4) 4 scalar 135,137,246,247,262,264,278,279,28 2,588,594,595,599,617,621,783,784, 790 NB Param 26 I(4) 4 scalar 29,37,39,60,101,135,246,588 NB1 Param 29 I(4) 4 scalar NBLM Param 17 I(4) 4 scalar NBLW Param 16 I(4) 4 scalar NBLX Param 16 I(4) 4 scalar NBLY Param 16 I(4) 4 scalar 17,24 NCLDS Dummy 4 I(4) 4 1 45 ARG,INOUT 298,300,320,341,398,410,448,458,47 1,484,491,507,515,524,533,543,568, 576,591,603,614,629,645,654,663,67 1,683,690,706,712,741,749,758,780, 788,794,802,809 NNCLDS Local 320 I(4) 4 scalar 320,322,324,329 NNLTE Param 27 I(4) 4 scalar NPES Scalar 2 I(4) 4 scalar COM O3DIFCTR Scalar 1 R(4) 4 scalar COM 187,338,372 ONE Scalar 8 R(4) 4 scalar COM 116,119,142,143,144,147,148,149,19 8,208,209,210,211,224,237,238,271, 323,325,327,348,349,350,376,381,38 2,383,384,402,415,427,444,461,462, 495,578,595,632,650,666,667,693 P0 Scalar 1 R(4) 4 scalar COM Page 29 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References P0INV Scalar 5 R(4) 4 scalar COM P0X2 Scalar 2 R(4) 4 scalar COM P0XZP2 Scalar 2 R(4) 4 scalar COM P0XZP8 Scalar 2 R(4) 4 scalar COM PHYCON Common 1 80 PP Local 54 R(4) 4 2 2070 117,118,123,127,128,414,415,561,56 2,773,774 PPTOP Local 80 R(4) 4 2 2070 414,415,561,562,773,774 PR2 Local 55 R(4) 4 2 2070 128,132,156,157 PRESS Dummy 3 R(4) 4 2 2070 ARG,INOUT 118,119,123 PWTS Scalar 101 R(4) 4 1 12 COM 137 QO3 Dummy 3 R(4) 4 2 2025 ARG,INOUT 158 QUARTR Scalar 8 R(4) 4 scalar COM RADCON Scalar 2 R(4) 4 scalar COM 291,812 RADCON1 Scalar 4 R(4) 4 scalar COM RATCO2MW Scalar 3 R(4) 4 scalar COM RATH2OMW Scalar 3 R(4) 4 scalar COM REFL Local 62 R(4) 4 1 45 143,233,349,391 REFL2 Local 63 R(4) 4 1 45 148,237 REFLO3 Scalar 101 R(4) 4 scalar COM RGAS Scalar 2 R(4) 4 scalar COM RGASSP Scalar 2 R(4) 4 scalar COM RH2O Dummy 3 R(4) 4 2 2025 ARG,INOUT 156 RRAY Local 60 R(4) 4 1 45 142,143,147,148,348,349 RRAYAV Scalar 101 R(4) 4 scalar COM 143,144,349,350 RRCO2 Dummy 3 R(4) 4 scalar ARG,INOUT 157 SECPDA Scalar 2 R(4) 4 scalar COM SECZ Local 61 R(4) 4 1 45 116,174,175,176,236,343,344,345 SIXTY Scalar 7 R(4) 4 scalar COM SQRT Func 116 scalar 116 SSOLAR Dummy 3 R(4) 4 scalar ARG,INOUT 137 SWR93 Subr 1 SWRSAV Common 101 112 TAUDAR Dummy 3 R(4) 4 1 45 ARG,INOUT 137 TCLD Local 89 R(4) 4 2 2070 475,495,675,694 TCLU Local 87 R(4) 4 2 2070 474,494,519,674,693,719 TCO2 Local 71 R(4) 4 2 4140 198,376 TDCL1 Local 84 R(4) 4 2 2070 441,443,444,451,461,474,475,509,51 0,647,649,650,657,666,674,675,708, 709 TDCL2 Local 85 R(4) 4 2 2070 443,453,475,649,659,675 TDCO2 Local 71 R(4) 4 2 2070 253,263,606,618 TDO3 Local 73 R(4) 4 2 2070 220,423 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 Page 30 Source Listing SWR93 2014-12-17 20:47 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 TEMPF Local 558 R(4) 4 scalar 558,561,770,773 TEMPG Local 559 R(4) 4 scalar 559,562,771,774 TEN Scalar 7 R(4) 4 scalar COM THREE Scalar 8 R(4) 4 scalar COM TMP1 Local 62 R(4) 4 1 45 119,132,578,579,580,581,797,798 TO3 Local 73 R(4) 4 2 4140 208,381 TOPO Common 50 17478548 TTD Local 78 R(4) 4 2 2070 218,220,253,423,427,428,441,451,45 3,538,551,606,617,631,632,647,657, 659,754,766 TTDB1 Local 84 R(4) 4 2 2070 421,423,606 TTU Local 78 R(4) 4 2 2070 219,221,254,424,428,442,452,537,55 0,609,621,631,648,658,753,765 TTUB1 Local 83 R(4) 4 2 2070 422,424,609 TTVG Scalar 51 R(4) 4 2 92971 COM TUCL1 Local 83 R(4) 4 2 2070 442,452,462,648,658,667 TUCO2 Local 72 R(4) 4 2 2070 254,265,609,622 TUO3 Local 74 R(4) 4 2 2070 221,424 TWO Scalar 8 R(4) 4 scalar COM 198,208,376,381 UCO2 Local 67 R(4) 4 2 4140 198,376 UD Local 66 R(4) 4 2 2070 168,174,179,218,262,353,359,364,42 1,617 UDCO2 Local 67 R(4) 4 2 2070 169,175,180,354,360,365 UDO3 Local 69 R(4) 4 2 2070 170,176,181,355,361,366 UFN Local 76 R(4) 4 2 2070 221,225,233,254,264,272,279,537,55 0,558,561,572,753,765,770,773,784 UFNCLU Local 86 R(4) 4 2 2070 509,518,520,527,558,561,708,721,73 4,736,744,770,773 UFNTRN Local 86 R(4) 4 2 2070 445,462,527,537,550,651,667,744,75 3,765 UFSWC Dummy 1 R(4) 4 2 2070 ARG,INOUT 305,572,784,798,805 UFSWL Dummy 1 R(4) 4 2 2070 ARG,INOUT 233,279,287,305,798 UO3 Local 69 R(4) 4 2 4140 208,209,210,211,381,382,383,384 UR Local 66 R(4) 4 2 2070 179,185,219,264,364,370,422,621 URCO2 Local 68 R(4) 4 2 2070 180,186,365,371 URO3 Local 70 R(4) 4 2 2070 181,187,366,372 XAMT Local 64 R(4) 4 2 2070 317,330,401,402,594,595 ZERO Scalar 8 R(4) 4 scalar COM 117,168,169,170,239,240,321,328,35 3,354,355,487,686 Page 31 Source Listing SWR93 2014-12-17 20:47 Subprograms/Common Blocks SWR93.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 37 128 HCON Common 7 872 MAPPINGS Common 56 5024 MPPCOM Common 1 1464 PHYCON Common 1 80 SWR93 Subr 1 SWRSAV Common 101 112 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__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ Page 32 Source Listing SWR93 2014-12-17 20:47 SWR93.F -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 : SWR93.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100