Page 1 Source Listing SWR93 2025-03-12 18:22 /tmp/ifortvVA69S.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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 Symbol Table SWR93.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 282 277 1000 Label 954 750 110 Label 286 283,284 1100 Label 962 955,957 1101 Label 961 958 120 Label 291 287,288 1200 Label 969 963,965 1201 Label 968 966 1250 Label 976 970,972 1251 Label 975 973 130 Label 295 292,293 140 Label 300 297,298 150 Label 307 303 155 Label 312 308 160 Label 321 316,317 200 Label 333 329 210 Label 339 334,335 220 Label 344 340 230 Label 350 345,346 240 Label 362 358,359 250 Label 374 368,369 260 Label 384 378,379 270 Label 388 385 280 Label 396 392,393 285 Label 403 397 290 Label 417 413,414 300 Label 428 422,423 310 Label 435 432 320 Label 442 438,439 330 Label 445 443 350 Label 446 408 360 Label 450 447,448 370 Label 454 451,452 400 Label 463 461 410 Label 469 464,465 420 Label 473 470,471 430 Label 480 477,478 450 Label 488 486 460 Label 493 491 470 Label 495 481,484 480 Label 501 496,497 490 Label 508 502,504 500 Label 513 509 510 Label 518 514 520 Label 524 519,520 530 Label 529 525 540 Label 535 530,531 550 Label 540 536,537 560 Label 547 541,542 570 Label 554 552 580 Label 565 562 Page 17 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 581 Label 566 559,561 590 Label 579 574 591 Label 580 571,573 600 Label 587 581,582 610 Label 591 588 620 Label 608 602 630 Label 616 612 631 Label 617 609,611 640 Label 625 622 641 Label 626 619,621 650 Label 638 635 651 Label 639 632,634 660 Label 650 645,647 670 Label 659 655 671 Label 660 652,654 680 Label 673 668,670 690 Label 683 679 691 Label 684 676,678 700 Label 691 688 701 Label 692 685,687 710 Label 701 698 720 Label 702 694,696 730 Label 714 711 740 Label 725 722 755 Label 727 707,710 760 Label 728 704,706 770 Label 736 729,731 771 Label 735 732 780 Label 744 737,739 790 Label 758 755 791 Label 759 752,754 800 Label 773 764,766 801 Label 769 767 802 Label 772 770 810 Label 786 775,777 811 Label 781 778 812 Label 785 782 820 Label 795 790,792 830 Label 814 806,808 840 Label 822 818 841 Label 823 815,817 850 Label 830 827 851 Label 831 824,826 860 Label 838 835 861 Label 839 832,834 870 Label 849 844,846 880 Label 858 854 881 Label 859 851,853 890 Label 872 867,869 900 Label 899 876 901 Label 900 873,875 910 Label 908 905 911 Label 909 902,904 920 Label 917 914 930 Label 918 910,912 Page 18 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 940 Label 929 926 950 Label 937 934 965 Label 939 922,925 970 Label 940 919,921 980 Label 948 941,943 981 Label 947 944 990 Label 953 949,951 ALFA Local 252 R(4) 4 2 2070 648,656,657,658,671,680,682,847,85 5,856,857,870,881,898 ALFAU Local 252 R(4) 4 2 2070 649,656,658,680,848,855,857,883 ALNB Dummy 6 R(4) 4 1 45 ARG,INOUT 441 ALND Dummy 6 R(4) 4 1 45 ARG,INOUT ALVB Dummy 6 R(4) 4 1 45 ARG,INOUT 305,310,400 ALVD Dummy 6 R(4) 4 1 45 ARG,INOUT 306,311,400,511,512 CAMT Dummy 5 R(4) 4 2 2070 ARG,INOUT 479,487,492 CCMAX Local 226 R(4) 4 1 45 483,485,487,489,490,492,740,743,95 2,959,960 COSZRO Dummy 4 R(4) 4 1 45 ARG,INOUT 278,299,304,309,510 CR Local 239 R(4) 4 2 2070 553,563,648,657,658,756,847,856,85 7 CRR Dummy 5 R(4) 4 3 24840 ARG,INOUT 563,756 CT Local 241 R(4) 4 2 2070 564,636,757,836 CTT Dummy 5 R(4) 4 3 24840 ARG,INOUT 564,757 DENOM Local 881 R(4) 4 scalar 881,882,883 DFN Local 238 R(4) 4 2 2070 382,386,387,394,415,424,433,434,44 0,444,700,713,721,724,733,916,928, 933,936,945,952 DFNCLU Local 250 R(4) 4 2 2070 672,682,690,721,724,871,898,907,93 3,936 DFNTOP Local 222 R(4) 4 2 540 299,394,395,440,441,444,733,734,94 5,946,952 DFNTRN Local 250 R(4) 4 2 2070 606,607,623,636,690,700,713,812,81 3,828,836,907,916,928 DFSWC Dummy 2 R(4) 4 2 2070 ARG,INOUT 466,733,743,945,959,967 DFSWL Dummy 3 R(4) 4 2 2070 ARG,INOUT 394,398,399,440,449,466,959 DP Local 216 R(4) 4 2 2070 289,318,319,320,453,974 DPCLD Local 242 R(4) 4 2 2070 577,720,721,932,933 DU Local 218 R(4) 4 2 2070 318,336,347,521,532 DUCO2 Local 218 R(4) 4 2 2070 319,337,348,522,533 DUO3 Local 219 R(4) 4 2 2070 320,338,349,523,534 EXP Func 360 scalar 360,371,380,381,398,424,426,538,54 4,583,584,779,783 FF Local 220 R(4) 4 2 2070 498,505,521 FFCO2 Local 220 R(4) 4 2 2070 499,506,522 FFO3 Local 221 R(4) 4 2 2070 500,507,523 FSWC Dummy 2 R(4) 4 2 2070 ARG,INOUT 468,967,974 FSWL Dummy 2 R(4) 4 2 2070 ARG,INOUT 449,453,468 GDFNB Dummy 6 R(4) 4 1 45 ARG,INOUT 401,444,742 GDFND Dummy 6 R(4) 4 1 45 ARG,INOUT 402,952 GDFVB Dummy 6 R(4) 4 1 45 ARG,INOUT 398,400,741 GDFVD Dummy 6 R(4) 4 1 45 ARG,INOUT 399,743 GLB_TABLE Common 175 128 SAVE HCON Common 56 872 SAVE HSWC Dummy 2 R(4) 4 2 2070 ARG,INOUT 472,974 HSWL Dummy 2 R(4) 4 2 2070 ARG,INOUT 453,472 Page 19 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References HTEMP Local 367 R(4) 4 scalar 367,372,545 I Local 277 I(4) 4 scalar 277,278,279,280,281,284,285,288,28 9,290,293,294,303,304,305,306,308, 309,310,311,317,318,319,320,335,33 6,337,338,359,360,369,370,371,372, 373,379,380,381,382,383,385,386,38 7,393,394,395,397,398,399,400,401, 402,414,415,416,423,424,425,426,42 7,432,433,434,439,440,441,443,444, 448,449,452,453,461,462,465,466,46 7,468,471,472,478,479,481,482,483, 485,487,489,490,492,497,498,499,50 0,509,510,511,512,520,521,522,523, 537,538,542,543,544,545,546,552,55 3,559,560,563,564,571,572,575,576, 577,582,583,584,585,586,588,589,59 0,602,603,604,605,606,607,609,610, 613,614,615,619,620,623,624,632,63 3,636,637,645,646,648,649,652,653, 656,657,658,668,669,671,672,676,67 7,680,681,682,685,686,689,690,694, 695,697,699,700,704,705,708,709,71 2,713,718,720,721,723,724,729,730, 733,734,737,738,740,741,742,743,75 2,753,756,757,764,765,768,771,775, 776,779,780,783,784,790,791,793,79 4,806,807,809,810,811,812,813,815, 816,819,820,821,824,825,828,829,83 2,833,836,837,844,845,847,848,851, 852,855,856,857,867,868,870,871,87 3,874,881,883,896,898,902,903,906, 907,910,911,913,915,916,919,920,92 3,924,927,928,930,932,933,935,936, 941,942,945,946,949,950,952,955,95 6,959,960,963,964,967,970,971,974 IDIM1 Param 45 I(4) 4 scalar 198,199,200,201,202,203,204,205,20 6,207,208,209,210,211,212,213,216, 217,218,219,220,221,222,223,224,22 5,226,228,229,230,231,232,233,234, 235,236,238,239,240,241,242,245,24 6,247,248,249,250,251,252,254,255, 256,257 IDIM2 Param 45 I(4) 4 scalar 198,199,200,201,202,203,204,205,20 6,207,208,209,210,211,212,213,216, 217,218,219,220,221,222,223,224,22 5,226,228,229,230,231,232,233,234, 235,236,238,239,240,241,242,245,24 6,247,248,249,250,251,252 IGSTL Param 40 I(4) 4 scalar 45 IGSTR Param 40 I(4) 4 scalar 45 IM Param 29 I(4) 4 scalar 42,45,112,181,182,183,184,185,189, 190,195 IMAX Param 112 I(4) 4 scalar 125 INLTE Param 124 I(4) 4 scalar 124 INLTEP Param 124 I(4) 4 scalar Page 20 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References INPES Param 37 I(4) 4 scalar 42,45,165,166,167,176,177 IP Local 298 I(4) 4 scalar 298,299,329,330,331,332,340,341,34 2,343,346,347,348,349,502,503,505, 506,507,514,515,516,517,525,526,52 7,528,531,532,533,534 ITAIL Param 42 I(4) 4 scalar J1 Local 708 I(4) 4 scalar 708,710,711,719,722,923,925,926,93 1,934 J2 Local 697 I(4) 4 scalar 697,698,709,711,913,914,924,926 J3 Local 718 I(4) 4 scalar 718,719,720,721,722,930,931,932,93 3,934 JDIM1 Param 46 I(4) 4 scalar JDIM2 Param 46 I(4) 4 scalar JGSTL Param 41 I(4) 4 scalar 46 JGSTR Param 41 I(4) 4 scalar 46 JM Param 29 I(4) 4 scalar 43,46,181,182,183,184,185,189,190, 195 JNPES Param 37 I(4) 4 scalar 43,46,165,166,167,176,177 JTAIL Param 43 I(4) 4 scalar JTOP Local 503 I(4) 4 scalar 503,504 K Local 283 I(4) 4 scalar 283,285,287,289,290,292,294,316,31 8,319,320,334,336,337,338,345,347, 348,349,358,360,368,370,371,372,37 3,378,380,381,382,383,392,394,395, 413,415,416,422,424,425,426,427,43 8,440,441,447,449,451,453,464,466, 467,468,470,472,477,479,486,487,49 1,492,496,498,499,500,504,505,506, 507,519,521,522,523,530,532,533,53 4,536,538,541,543,544,545,546,581, 583,584,585,586,698,699,700,711,71 2,713,722,723,724,732,733,734,755, 756,757,914,915,916,926,927,928,93 4,935,936,944,945,946,958,959,960 KBTMSW Dummy 5 I(4) 4 2 2070 ARG,INOUT 575,577,615,697,709,718,821,913,92 4,930 KCLDS Local 460 I(4) 4 scalar 460,462,475,560,561,562,572,573,57 4,610,611,612,620,621,622,633,634, 635,646,647,653,654,655,669,670,67 1,672,677,678,679,686,687,688,695, 696,705,706,707,730,731,738,739,75 3,754,755,765,766,776,777,791,792, 807,808,816,817,818,825,826,827,83 3,834,835,845,846,852,853,854,868, 869,870,871,874,875,876,903,904,90 5,911,912,920,921,922,942,943,950, 951,956,957,964,965,971,972 KK Local 562 I(4) 4 scalar 562,563,564,574,575,576,577,612,61 3,614,615,622,623,624,635,636,637, 655,656,657,658,679,680,681,682,68 8,689,690,707,708,709,712,713,718, 720,721,723,724,767,768,770,771,77 8,779,780,782,783,784,818,819,820, 821,827,828,829,835,836,837,854,85 5,856,857,876,881,883,896,898,905, Page 21 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References 906,907,922,923,924,927,928,930,93 2,933,935,936,966,967,973,974 KO2 Param 127 I(4) 4 scalar 128 KO21 Param 128 I(4) 4 scalar KO2M Param 128 I(4) 4 scalar KTOPSW Dummy 5 I(4) 4 2 2070 ARG,INOUT 503,575,576,613,614,708,819,820,92 3 L Param 111 I(4) 4 scalar 115,116,117,120,122,204,205,287,29 2,316,345,378,413,422,451,470,530, 581,770,782,973 LL Param 117 I(4) 4 scalar 117,118,358,368,536,541 LL3P Param 122 I(4) 4 scalar 125 LL3PI Param 125 I(4) 4 scalar LLM1 Param 118 I(4) 4 scalar LLM2 Param 118 I(4) 4 scalar LLM3 Param 118 I(4) 4 scalar LLP1 Param 117 I(4) 4 scalar 125 LLP1I Param 125 I(4) 4 scalar LLP2 Param 117 I(4) 4 scalar 229,231,233,235 LLP3 Param 117 I(4) 4 scalar LM Param 29 I(4) 4 scalar 111,190 LM1 Param 116 I(4) 4 scalar 283 LM2 Param 116 I(4) 4 scalar LM3 Param 116 I(4) 4 scalar LOG Func 360 scalar 360,371,538,544 LP1 Param 115 I(4) 4 scalar 119,120,121,125,198,199,200,201,20 2,203,204,206,207,216,217,218,219, 220,221,226,228,229,230,231,232,23 3,234,235,236,238,239,240,241,242, 245,246,247,248,249,250,251,252,28 0,281,334,341,342,343,387,392,398, 399,434,438,444,447,464,477,496,51 9,526,527,528,590,603,604,698,732, 743,767,778,793,809,810,914,944,95 2,958,966 LP121 Param 121 I(4) 4 scalar LP1I Param 125 I(4) 4 scalar LP1M Param 119 I(4) 4 scalar 119 LP1M1 Param 119 I(4) 4 scalar LP1V Param 120 I(4) 4 scalar LP2 Param 115 I(4) 4 scalar 254,255,256,257 LP3 Param 115 I(4) 4 scalar LSM Param 29 I(4) 4 scalar MAPPINGS Common 194 5024 SAVE MAX Func 462 scalar 462 MIN Func 380 scalar 380,381,424,426,583,584,779,783 MPPCOM Common 139 1464 SAVE N Local 297 I(4) 4 scalar 297,299,408,409,424,426,440,441,44 4,750,756,757,761,779,783,945,946, 952 NB Param 123 I(4) 4 scalar 126,199,201,222,263,297,408,750 NB1 Param 126 I(4) 4 scalar NBLM Param 114 I(4) 4 scalar NBLW Param 113 I(4) 4 scalar NBLX Param 113 I(4) 4 scalar Page 22 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References NBLY Param 113 I(4) 4 scalar 114,121 NCLDS Dummy 5 I(4) 4 1 45 ARG,INOUT 460,462,482,503,560,572,610,620,63 3,646,653,669,677,686,695,705,730, 738,753,765,776,791,807,816,825,83 3,845,852,868,874,903,911,920,942, 950,956,964,971 NNCLDS Local 482 I(4) 4 scalar 482,484,486,491 NNLTE Param 124 I(4) 4 scalar PHYCON Common 50 80 SAVE PP Local 216 R(4) 4 2 2070 279,280,285,289,290,576,577,723,72 4,935,936 PPTOP Local 242 R(4) 4 2 2070 576,577,723,724,935,936 PR2 Local 217 R(4) 4 2 2070 290,294,318,319 PRESS Dummy 4 R(4) 4 2 2070 ARG,INOUT 280,281,285 QO3 Dummy 4 R(4) 4 2 2025 ARG,INOUT 320 REFL Local 224 R(4) 4 1 45 305,395,511,553 REFL2 Local 225 R(4) 4 1 45 310,399 RH2O Dummy 4 R(4) 4 2 2025 ARG,INOUT 318 RRAY Local 222 R(4) 4 1 45 304,305,309,310,510,511 RRCO2 Dummy 4 R(4) 4 scalar ARG,INOUT 319 SECZ Local 223 R(4) 4 1 45 278,336,337,338,398,505,506,507 SQRT Func 278 scalar 278 SSOLAR Dummy 4 R(4) 4 scalar ARG,INOUT 299 SWR93 Subr 2 SWRSAV Common 263 112 SAVE TAUDAR Dummy 4 R(4) 4 1 45 ARG,INOUT 299 TCLD Local 251 R(4) 4 2 2070 637,657,837,856 TCLU Local 249 R(4) 4 2 2070 636,656,681,836,855,881 TCO2 Local 233 R(4) 4 2 4140 360,538 TDCL1 Local 246 R(4) 4 2 2070 603,605,606,613,623,636,637,671,67 2,809,811,812,819,828,836,837,870, 871 TDCL2 Local 247 R(4) 4 2 2070 605,615,637,811,821,837 TDCO2 Local 233 R(4) 4 2 2070 415,425,768,780 TDO3 Local 235 R(4) 4 2 2070 382,585 TEMPCOM Common 180 6603768 SAVE TEMPF Local 720 R(4) 4 scalar 720,723,932,935 TEMPG Local 721 R(4) 4 scalar 721,724,933,936 TMP1 Local 224 R(4) 4 1 45 281,294,740,741,742,743,959,960 TO3 Local 235 R(4) 4 2 4140 370,543 TOPO Common 188 17478548 SAVE TTD Local 240 R(4) 4 2 2070 380,382,415,585,589,590,603,613,61 5,700,713,768,779,793,794,809,819, 821,916,928 TTDB1 Local 246 R(4) 4 2 2070 583,585,768 TTU Local 240 R(4) 4 2 2070 381,383,416,586,590,604,614,699,71 2,771,783,793,810,820,915,927 TTUB1 Local 245 R(4) 4 2 2070 584,586,771 TUCL1 Local 245 R(4) 4 2 2070 604,614,624,810,820,829 TUCO2 Local 234 R(4) 4 2 2070 416,427,771,784 TUO3 Local 236 R(4) 4 2 2070 383,586 UCO2 Local 229 R(4) 4 2 4140 360,538 UD Local 228 R(4) 4 2 2070 330,336,341,380,424,515,521,526,58 3,779 UDCO2 Local 229 R(4) 4 2 2070 331,337,342,516,522,527 Page 23 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Object Declared Type Bytes Dimen Elements Attributes References UDO3 Local 231 R(4) 4 2 2070 332,338,343,517,523,528 UFN Local 238 R(4) 4 2 2070 383,387,395,416,426,434,441,699,71 2,720,723,734,915,927,932,935,946 UFNCLU Local 248 R(4) 4 2 2070 671,680,682,689,720,723,870,883,89 6,898,906,932,935 UFNTRN Local 248 R(4) 4 2 2070 607,624,689,699,712,813,829,906,91 5,927 UFSWC Dummy 2 R(4) 4 2 2070 ARG,INOUT 467,734,946,960,967 UFSWL Dummy 2 R(4) 4 2 2070 ARG,INOUT 395,441,449,467,960 UO3 Local 231 R(4) 4 2 4140 370,371,372,373,543,544,545,546 UR Local 228 R(4) 4 2 2070 341,347,381,426,526,532,584,783 URCO2 Local 230 R(4) 4 2 2070 342,348,527,533 URO3 Local 232 R(4) 4 2 2070 343,349,528,534 XAMT Local 226 R(4) 4 2 2070 479,492,563,564,756,757 TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References ABCFF R(4) 4 0 1 12 COM 380,381,424,426,583,584,779,783 AMOLWT R(4) 4 0 scalar COM CFCO2 R(4) 4 96 scalar COM 319 CFO3 R(4) 4 100 scalar COM 320 CSUBP R(4) 4 4 scalar COM DIFFCTR R(4) 4 8 scalar COM 347,348,498,499,532,533 EIGHT R(4) 4 20 scalar COM FIFTY R(4) 4 12 scalar COM 380,381,424,426,583,584,779,783 FIVE R(4) 4 24 scalar COM FOUR R(4) 4 28 scalar COM G R(4) 4 12 scalar COM G2LI I(4) 4 0 1 239 COM G2LJ I(4) 4 1912 1 389 COM GINV R(4) 4 68 scalar COM 318,319,320 GP0INV R(4) 4 76 scalar COM GRAVDR R(4) 4 16 scalar COM H101M16 R(4) 4 808 scalar COM H102M5 R(4) 4 784 scalar COM H1036E2 R(4) 4 756 scalar COM 367 H114M11 R(4) 4 384 scalar COM H1174M7 R(4) 4 268 scalar COM H11M10 R(4) 4 296 scalar COM H11M11 R(4) 4 388 scalar COM H1224E3 R(4) 4 112 scalar COM 278 H1226E1 R(4) 4 652 scalar COM H128M5 R(4) 4 240 scalar COM H129M2 R(4) 4 732 scalar COM 360,538 H12M12 R(4) 4 472 scalar COM H12M13 R(4) 4 492 scalar COM H12M31 R(4) 4 560 scalar COM H135M13 R(4) 4 488 scalar COM H1386E2 R(4) 4 748 scalar COM 371,544 H1439M5 R(4) 4 236 scalar COM H14M10 R(4) 4 292 scalar COM H14M11 R(4) 4 380 scalar COM Page 24 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References H14M12 R(4) 4 468 scalar COM H14M14 R(4) 4 508 scalar COM H14M30 R(4) 4 548 scalar COM H15E2 R(4) 4 136 scalar COM H15M11 R(4) 4 376 scalar COM H15M14 R(4) 4 504 scalar COM H15M5 R(4) 4 824 scalar COM H161E1 R(4) 4 800 scalar COM H165E5 R(4) 4 84 scalar COM H16E1 R(4) 4 792 scalar COM H16M12 R(4) 4 464 scalar COM H181E1 R(4) 4 148 scalar COM H18E1 R(4) 4 152 scalar COM H18E3 R(4) 4 832 scalar COM H18M11 R(4) 4 372 scalar COM H1E11 R(4) 4 72 scalar COM H1E13 R(4) 4 68 scalar COM H1E15 R(4) 4 64 scalar COM H1E4 R(4) 4 96 scalar COM H1E6 R(4) 4 612 scalar COM H1E8 R(4) 4 76 scalar COM H1M10 R(4) 4 300 scalar COM H1M11 R(4) 4 392 scalar COM H1M13 R(4) 4 496 scalar COM H1M16 R(4) 4 856 scalar COM H1M17 R(4) 4 512 scalar COM H1M18 R(4) 4 516 scalar COM H1M19 R(4) 4 520 scalar COM H1M2 R(4) 4 620 scalar COM H1M20 R(4) 4 524 scalar COM H1M21 R(4) 4 528 scalar COM H1M22 R(4) 4 532 scalar COM H1M23 R(4) 4 536 scalar COM H1M24 R(4) 4 540 scalar COM H1M3 R(4) 4 208 scalar COM H1M4 R(4) 4 220 scalar COM H1M5 R(4) 4 244 scalar COM H1M6 R(4) 4 260 scalar COM H1M60 R(4) 4 584 scalar COM H1M8 R(4) 4 284 scalar COM H1P082 R(4) 4 740 scalar COM 371,544 H1P25892 R(4) 4 180 scalar COM H1P4 R(4) 4 176 scalar COM H1P41819 R(4) 4 704 scalar COM H1P4387 R(4) 4 172 scalar COM H1P8 R(4) 4 168 scalar COM H2075E3 R(4) 4 108 scalar COM H20788E3 R(4) 4 104 scalar COM H2118M2 R(4) 4 760 scalar COM 373,546 H21M12 R(4) 4 460 scalar COM H21M31 R(4) 4 556 scalar COM H235M3 R(4) 4 724 scalar COM 360,538 H23E2 R(4) 4 132 scalar COM H23M10 R(4) 4 288 scalar COM H23M11 R(4) 4 364 scalar COM Page 25 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References H24E3 R(4) 4 100 scalar COM H24M11 R(4) 4 360 scalar COM H24M12 R(4) 4 456 scalar COM H25452M6 R(4) 4 256 scalar COM H257M8 R(4) 4 280 scalar COM H25E2 R(4) 4 816 scalar COM H25M31 R(4) 4 552 scalar COM H26E2 R(4) 4 696 scalar COM H26M30 R(4) 4 544 scalar COM H28571M2 R(4) 4 852 scalar COM H285M4 R(4) 4 216 scalar COM H28E1 R(4) 4 868 scalar COM H28M11 R(4) 4 356 scalar COM H28M12 R(4) 4 452 scalar COM H29316E2 R(4) 4 648 scalar COM H2945E2 R(4) 4 128 scalar COM H29M12 R(4) 4 448 scalar COM H2E2 R(4) 4 688 scalar COM H2E6 R(4) 4 616 scalar COM H2M11 R(4) 4 368 scalar COM H2P5 R(4) 4 164 scalar COM H2P8 R(4) 4 160 scalar COM H2P9 R(4) 4 156 scalar COM H3082E2 R(4) 4 120 scalar COM H3116E1 R(4) 4 656 scalar COM H323M4 R(4) 4 768 scalar COM 373,546 H327M8 R(4) 4 276 scalar COM H32M11 R(4) 4 348 scalar COM H35E1 R(4) 4 140 scalar COM 278 H35M11 R(4) 4 344 scalar COM H36M13 R(4) 4 484 scalar COM H37412M5 R(4) 4 232 scalar COM H37M11 R(4) 4 340 scalar COM H37M12 R(4) 4 440 scalar COM H38M12 R(4) 4 436 scalar COM H391M7 R(4) 4 264 scalar COM H394M5 R(4) 4 228 scalar COM H3E2 R(4) 4 124 scalar COM H3M11 R(4) 4 352 scalar COM H3M12 R(4) 4 444 scalar COM H3M14 R(4) 4 500 scalar COM H3M3 R(4) 4 804 scalar COM H3M4 R(4) 4 860 scalar COM H3P5 R(4) 4 828 scalar COM H3P6 R(4) 4 144 scalar COM H41666M2 R(4) 4 636 scalar COM H42M11 R(4) 4 336 scalar COM H42M2 R(4) 4 764 scalar COM 373,546 H44194M2 R(4) 4 700 scalar COM H44871M2 R(4) 4 200 scalar COM H44M11 R(4) 4 332 scalar COM H44M12 R(4) 4 428 scalar COM H451M6 R(4) 4 788 scalar COM H45M12 R(4) 4 424 scalar COM H45M32 R(4) 4 572 scalar COM Page 26 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References H46M13 R(4) 4 480 scalar COM H488E4 R(4) 4 92 scalar COM H48M11 R(4) 4 328 scalar COM H4999M6 R(4) 4 252 scalar COM H4E5 R(4) 4 80 scalar COM H4M12 R(4) 4 432 scalar COM H4M33 R(4) 4 576 scalar COM H53M11 R(4) 4 324 scalar COM H559M3 R(4) 4 204 scalar COM H55M32 R(4) 4 568 scalar COM H5725E4 R(4) 4 88 scalar COM H5E2 R(4) 4 116 scalar COM H625M2 R(4) 4 668 scalar COM H62M12 R(4) 4 416 scalar COM H62M34 R(4) 4 580 scalar COM H658M2 R(4) 4 752 scalar COM 372,545 H65M12 R(4) 4 412 scalar COM H67390E2 R(4) 4 772 scalar COM H6938M5 R(4) 4 224 scalar COM H69766E5 R(4) 4 720 scalar COM 299 H6M12 R(4) 4 420 scalar COM H6P08108 R(4) 4 836 scalar COM H71E26 R(4) 4 60 scalar COM H72M11 R(4) 4 320 scalar COM H74M12 R(4) 4 408 scalar COM H75826M4 R(4) 4 736 scalar COM 361,539 H77M11 R(4) 4 316 scalar COM H77M12 R(4) 4 404 scalar COM H7M6 R(4) 4 248 scalar COM H8121E1 R(4) 4 684 scalar COM H82M11 R(4) 4 308 scalar COM H83E26 R(4) 4 56 scalar COM H83M11 R(4) 4 304 scalar COM H8725M8 R(4) 4 272 scalar COM H8M11 R(4) 4 312 scalar COM H8M13 R(4) 4 476 scalar COM H93M12 R(4) 4 400 scalar COM H96M12 R(4) 4 396 scalar COM H987M4 R(4) 4 212 scalar COM H9M32 R(4) 4 564 scalar COM H9P94 R(4) 4 660 scalar COM HAF R(4) 4 44 scalar COM 285,290 HM13EZ R(4) 4 592 scalar COM HM1597E1 R(4) 4 812 scalar COM HM161E1 R(4) 4 796 scalar COM HM1797E1 R(4) 4 680 scalar COM HM181E1 R(4) 4 604 scalar COM HM19EZ R(4) 4 596 scalar COM HM1E1 R(4) 4 600 scalar COM HM1E2 R(4) 4 608 scalar COM HM1EZ R(4) 4 692 scalar COM 380,381,424,426,583,584,779,783 HM2M2 R(4) 4 644 scalar COM HM6666M2 R(4) 4 628 scalar COM HM8E1 R(4) 4 864 scalar COM HMP5 R(4) 4 640 scalar COM Page 27 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References HMP575 R(4) 4 588 scalar COM HMP66667 R(4) 4 624 scalar COM HMP805 R(4) 4 840 scalar COM 371,544 HNINETY R(4) 4 4 scalar COM HP1 R(4) 4 196 scalar COM HP118666 R(4) 4 820 scalar COM HP144 R(4) 4 712 scalar COM HP166666 R(4) 4 632 scalar COM HP219 R(4) 4 708 scalar COM 304,510 HP228 R(4) 4 672 scalar COM HP26 R(4) 4 728 scalar COM 360,538 HP369 R(4) 4 192 scalar COM HP3795 R(4) 4 776 scalar COM HP5048 R(4) 4 780 scalar COM HP518 R(4) 4 188 scalar COM HP526315 R(4) 4 848 scalar COM HP6 R(4) 4 664 scalar COM HP602409 R(4) 4 844 scalar COM HP60241 R(4) 4 676 scalar COM HP8 R(4) 4 184 scalar COM HP805 R(4) 4 744 scalar COM HP816 R(4) 4 716 scalar COM 304,510 HTMG R(4) 4 743768 3 4183695 COM HUNDRED R(4) 4 0 scalar COM IBROW I(4) 4 404 scalar COM ICHUNKTAB I(4) 4 620 1 8 COM IE_GLB_TABLE I(4) 4 32 1 8 COM IE_LOC_TABLE I(4) 4 556 1 8 COM ILCOL I(4) 4 396 scalar COM ILPAD1 I(4) 4 412 scalar COM ILPAD2 I(4) 4 416 scalar COM ILPAD3 I(4) 4 420 scalar COM ILPAD4 I(4) 4 424 scalar COM ILPAD5 I(4) 4 428 scalar COM INUMQ I(4) 4 1064 1 100 COM IQUILT_GROUP I(4) 4 1060 scalar COM IRCOL I(4) 4 400 scalar COM IRPAD1 I(4) 4 432 scalar COM IRPAD2 I(4) 4 436 scalar COM IRPAD3 I(4) 4 440 scalar COM IRPAD4 I(4) 4 444 scalar COM IRPAD5 I(4) 4 448 scalar COM IS_GLB_TABLE I(4) 4 0 1 8 COM IS_LOC_TABLE I(4) 4 492 1 8 COM ITEMP I(4) 4 5870016 2 91719 COM ITEMP2 I(4) 4 6236892 2 91719 COM ITROW I(4) 4 408 scalar COM JBPAD1 I(4) 4 452 scalar COM JBPAD2 I(4) 4 456 scalar COM JBPAD3 I(4) 4 460 scalar COM JBPAD4 I(4) 4 464 scalar COM JBPAD5 I(4) 4 468 scalar COM JE_GLB_TABLE I(4) 4 96 1 8 COM JE_LOC_TABLE I(4) 4 588 1 8 COM JS_GLB_TABLE I(4) 4 64 1 8 COM Page 28 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References JS_LOC_TABLE I(4) 4 524 1 8 COM JTPAD1 I(4) 4 472 scalar COM JTPAD2 I(4) 4 476 scalar COM JTPAD3 I(4) 4 480 scalar COM JTPAD4 I(4) 4 484 scalar COM JTPAD5 I(4) 4 488 scalar COM L2GI I(4) 4 956 1 239 COM L2GJ I(4) 4 3468 1 389 COM MPI_COMM_COMP I(4) 4 652 scalar COM MPI_COMM_INTER I(4) 4 656 scalar COM MPI_COMM_INTER_ARRAY I(4) 4 660 1 100 COM MYIE I(4) 4 52 scalar COM 277,284,288,293,298,303,308,317,32 9,335,340,346,359,369,379,385,393, 397,414,423,432,439,443,448,452,46 1,465,471,478,481,497,502,509,514, 520,525,531,537,542,552,559,571,58 2,588,602,609,619,632,645,652,668, 676,685,694,704,729,737,752,764,77 5,790,806,815,824,832,844,851,867, 873,902,910,919,941,949,955,963,97 0 MYIE1 I(4) 4 56 scalar COM MYIE1_P1 I(4) 4 120 scalar COM MYIE1_P2 I(4) 4 124 scalar COM MYIE1_P3 I(4) 4 128 scalar COM MYIE1_P4 I(4) 4 132 scalar COM MYIE2 I(4) 4 60 scalar COM MYIE2_P1 I(4) 4 136 scalar COM MYIE_P1 I(4) 4 100 scalar COM MYIE_P2 I(4) 4 104 scalar COM MYIE_P3 I(4) 4 108 scalar COM MYIE_P4 I(4) 4 112 scalar COM MYIE_P5 I(4) 4 116 scalar COM MYIS I(4) 4 40 scalar COM 277,284,288,293,298,303,308,317,32 9,335,340,346,359,369,379,385,393, 397,414,423,432,439,443,448,452,46 0,465,471,478,481,497,502,509,514, 520,525,531,537,542,552,559,571,58 2,588,602,609,619,632,645,652,668, 676,685,694,704,729,737,752,764,77 5,790,806,815,824,832,844,851,867, 873,902,910,919,941,949,955,963,97 0 MYIS1 I(4) 4 44 scalar COM 461 MYIS1_P1 I(4) 4 84 scalar COM MYIS1_P2 I(4) 4 88 scalar COM MYIS1_P3 I(4) 4 92 scalar COM MYIS1_P4 I(4) 4 96 scalar COM MYIS2 I(4) 4 48 scalar COM MYIS_P1 I(4) 4 64 scalar COM MYIS_P2 I(4) 4 68 scalar COM MYIS_P3 I(4) 4 72 scalar COM MYIS_P4 I(4) 4 76 scalar COM MYIS_P5 I(4) 4 80 scalar COM MYJE I(4) 4 236 scalar COM Page 29 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References MYJE1 I(4) 4 240 scalar COM MYJE1_P1 I(4) 4 280 scalar COM MYJE1_P2 I(4) 4 284 scalar COM MYJE1_P3 I(4) 4 288 scalar COM MYJE1_P4 I(4) 4 292 scalar COM MYJE2 I(4) 4 244 scalar COM MYJE2_P1 I(4) 4 296 scalar COM MYJE2_P2 I(4) 4 300 scalar COM MYJE2_P3 I(4) 4 304 scalar COM MYJE2_P4 I(4) 4 308 scalar COM MYJE3 I(4) 4 248 scalar COM MYJE3_P4 I(4) 4 312 scalar COM MYJE4 I(4) 4 252 scalar COM MYJE4_P1 I(4) 4 316 scalar COM MYJE4_P4 I(4) 4 320 scalar COM MYJE5 I(4) 4 256 scalar COM MYJE5_P1 I(4) 4 324 scalar COM MYJE5_P2 I(4) 4 328 scalar COM MYJE_P1 I(4) 4 260 scalar COM MYJE_P2 I(4) 4 264 scalar COM MYJE_P3 I(4) 4 268 scalar COM MYJE_P4 I(4) 4 272 scalar COM MYJE_P5 I(4) 4 276 scalar COM MYJS I(4) 4 140 scalar COM MYJS1 I(4) 4 144 scalar COM MYJS1_P1 I(4) 4 184 scalar COM MYJS1_P2 I(4) 4 188 scalar COM MYJS1_P3 I(4) 4 192 scalar COM MYJS1_P4 I(4) 4 196 scalar COM MYJS2 I(4) 4 148 scalar COM MYJS2_P1 I(4) 4 200 scalar COM MYJS2_P2 I(4) 4 204 scalar COM MYJS2_P3 I(4) 4 208 scalar COM MYJS2_P4 I(4) 4 212 scalar COM MYJS3 I(4) 4 152 scalar COM MYJS3_P4 I(4) 4 216 scalar COM MYJS4 I(4) 4 156 scalar COM MYJS4_P1 I(4) 4 220 scalar COM MYJS4_P4 I(4) 4 224 scalar COM MYJS5 I(4) 4 160 scalar COM MYJS5_P1 I(4) 4 228 scalar COM MYJS5_P2 I(4) 4 232 scalar COM MYJS_P1 I(4) 4 164 scalar COM MYJS_P2 I(4) 4 168 scalar COM MYJS_P3 I(4) 4 172 scalar COM MYJS_P4 I(4) 4 176 scalar COM MYJS_P5 I(4) 4 180 scalar COM MYPE I(4) 4 0 scalar COM MY_E I(4) 4 336 scalar COM MY_IE_GLB I(4) 4 12 scalar COM MY_IE_LOC I(4) 4 28 scalar COM MY_IS_GLB I(4) 4 8 scalar COM MY_IS_LOC I(4) 4 24 scalar COM MY_JE_GLB I(4) 4 20 scalar COM MY_JE_LOC I(4) 4 36 scalar COM Page 30 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References MY_JS_GLB I(4) 4 16 scalar COM MY_JS_LOC I(4) 4 32 scalar COM MY_N I(4) 4 332 scalar COM MY_NE I(4) 4 348 scalar COM MY_NEB I(4) 4 364 1 8 COM MY_NW I(4) 4 360 scalar COM MY_S I(4) 4 340 scalar COM MY_SE I(4) 4 352 scalar COM MY_SW I(4) 4 356 scalar COM MY_W I(4) 4 344 scalar COM NPES I(4) 4 4 scalar COM O3DIFCTR R(4) 4 20 scalar COM 349,500,534 ONE R(4) 4 40 scalar COM 278,281,304,305,306,309,310,311,36 0,370,371,372,373,386,399,400,433, 485,487,489,510,511,512,538,543,54 4,545,546,564,577,589,606,623,624, 657,740,757,794,812,828,829,855 P0 R(4) 4 24 scalar COM P0INV R(4) 4 72 scalar COM P0X2 R(4) 4 36 scalar COM P0XZP2 R(4) 4 28 scalar COM P0XZP8 R(4) 4 32 scalar COM PWTS R(4) 4 48 1 12 COM 299 QUARTR R(4) 4 48 scalar COM RADCON R(4) 4 40 scalar COM 453,974 RADCON1 R(4) 4 64 scalar COM RATCO2MW R(4) 4 56 scalar COM RATH2OMW R(4) 4 60 scalar COM REFLO3 R(4) 4 104 scalar COM RGAS R(4) 4 44 scalar COM RGASSP R(4) 4 48 scalar COM RRAYAV R(4) 4 108 scalar COM 305,306,511,512 SECPDA R(4) 4 52 scalar COM SIXTY R(4) 4 8 scalar COM TEMP1 R(4) 4 0 2 91719 COM TEMP10 R(4) 4 3301884 2 91719 COM TEMP11 R(4) 4 3668760 2 91719 COM TEMP12 R(4) 4 4035636 2 91719 COM TEMP13 R(4) 4 4402512 2 91719 COM TEMP14 R(4) 4 4769388 2 91719 COM TEMP15 R(4) 4 5136264 2 91719 COM TEMP16 R(4) 4 5503140 2 91719 COM TEMP2 R(4) 4 366876 2 91719 COM TEMP2X R(4) 4 0 2 92971 COM TEMP3 R(4) 4 733752 2 91719 COM TEMP4 R(4) 4 1100628 2 91719 COM TEMP5 R(4) 4 1467504 2 91719 COM TEMP6 R(4) 4 1834380 2 91719 COM TEMP7 R(4) 4 2201256 2 91719 COM TEMP8 R(4) 4 2568132 2 91719 COM TEMP9 R(4) 4 2935008 2 91719 COM TEN R(4) 4 16 scalar COM THREE R(4) 4 32 scalar COM TTVG R(4) 4 371884 2 92971 COM TWO R(4) 4 36 scalar COM 360,370,538,543 Page 31 Source Listing SWR93 2025-03-12 18:22 Symbol Table SWR93.F Name Type Bytes Offset Dimen Elements Attributes References ZERO R(4) 4 52 scalar COM 279,330,331,332,401,402,483,490,51 5,516,517,649,848 Page 32 Source Listing SWR93 2025-03-12 18:22 Subprograms/Common Blocks SWR93.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 175 128 SAVE HCON Common 56 872 SAVE MAPPINGS Common 194 5024 SAVE MPPCOM Common 139 1464 SAVE PHYCON Common 50 80 SAVE SWR93 Subr 2 SWRSAV Common 263 112 SAVE TEMPCOM Common 180 6603768 SAVE TOPO Common 188 17478548 SAVE COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -align norec32byte -align norec64byte -align noarray8byte -align noarray16byte -align noarray32byte -align noarray64byte -align noarray128byte -align noarray256byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume nobuffered_stdout -assume byterecl -assume nocontiguous_assumed_shape -assume nocontiguous_pointer -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_complex_align -assume old_unit_star -assume old_inquire_recl -assume old_ldout_format -assume old_ldout_zero -assume noold_logical_assign -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume noprotect_allocates -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume std_minus0_rounding -assume nostd_mod_proc_name -assume std_value -assume realloc_lhs -assume underscore -assume no2underscores -assume norecursion no -auto -auto_scalar no -bintext -ccdefault default -check noarg_temp_created -check noassume -check nobounds -check nocontiguous -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check noshape -check nostack -check nouninitialized -check noudio_iostat -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference Page 33 Source Listing SWR93 2025-03-12 18:22 SWR93.F -D __INTEL_COMPILER=1910 -D __INTEL_COMPILER_UPDATE=3 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D __amd64 -D __amd64__ -D __INTEL_COMPILER_BUILD_DATE=20200925 -D __INTEL_OFFLOAD -D __MMX__ -D __SSE__ -D __SSE_MATH__ -D __SSE2__ -D __SSE2_MATH__ -D __SSE3__ -D __SSSE3__ -D __SSE4_1__ -D __SSE4_2__ -D __POPCNT__ -D __PCLMUL__ -D __AES__ -D __AVX__ -D __F16C__ -D __AVX_I__ -D __RDRND__ -D __FMA__ -D __FP_FAST_FMA -D __FP_FAST_FMAF -D __BMI__ -D __LZCNT__ -D __AVX2__ -D __haswell -D __haswell__ -D __tune_haswell__ -D __core_avx2 -D __core_avx2__ -D __tune_core_avx2__ -D __CRAY_X86_ROME -D __CRAYXT_COMPUTE_LINUX_TARGET -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg -init noarrays -init nohuge -init noinfinity -init nominus_huge -init nominus_infinity -init nominus_tiny -init nonan -init nosnan -init notiny Page 34 Source Listing SWR93 2025-03-12 18:22 SWR93.F -init nozero no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude no -o -offload-build=host -openmp-simd -O2 no -pad_source -real_size 32 no -recursive -reentrancy threaded -vec=simd -show nofullpath -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w nodeclarations -w noexternals -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w noshape -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage no -wrap-margins -includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, .f90,./.f90,/opt/cray/pe/mpich/8.1.12/ofi/intel/19.0/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.f90,/usr/lib64/gcc/x86_64-suse-linux/7/include/.f90, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.f90,/usr/include/.f90,/usr/include/.f90,/usr/include/.f90 -list filename : SWR93.lst no -o COMPILER: Intel(R) Fortran 19.1-1655