Page 1 Source Listing SWR93
2025-03-12 18:23 /tmp/ifortLi223W.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"
50 INCLUDE "HCON.comm"
98 INCLUDE "rdparm"
130 INCLUDE "mpp.h"
131 # 1 "./sp.h" 1
132 # 4
133
134 # 30 "SWR93.F" 2
135 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE:
136 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL
137 C NB IS A SHORTWAVE PARAMETER; OTHER QUANTITIES ARE DERIVED
138 C FROM THE ABOVE PARAMETERS.
139 C --- VARIABLES AS IN ARGUMENT LIST
198 D I M E N S I O N
199 1 FSWC (IDIM1:IDIM2,LP1), HSWC (IDIM1:IDIM2,LP1),
200 2 CRR (IDIM1:IDIM2,NB,LP1)
201 3, FSWL (IDIM1:IDIM2,LP1), HSWL (IDIM1:IDIM2,LP1),
202 4 CTT (IDIM1:IDIM2,NB,LP1)
203 5, UFSWC (IDIM1:IDIM2,LP1), DFSWC (IDIM1:IDIM2,LP1)
204 6, UFSWL (IDIM1:IDIM2,LP1), DFSWL (IDIM1:IDIM2,LP1)
205 7, PRESS (IDIM1:IDIM2,LP1), RH2O (IDIM1:IDIM2,L),
206 8 QO3 (IDIM1:IDIM2,L)
207 9, CAMT (IDIM1:IDIM2,LP1), KTOPSW(IDIM1:IDIM2,LP1),
208 o KBTMSW(IDIM1:IDIM2,LP1)
209 1, COSZRO(IDIM1:IDIM2), TAUDAR(IDIM1:IDIM2),
210 2 NCLDS (IDIM1:IDIM2)
211 3, ALVB (IDIM1:IDIM2), ALNB (IDIM1:IDIM2),
212 4 ALVD (IDIM1:IDIM2), ALND (IDIM1:IDIM2)
213 5, GDFVB (IDIM1:IDIM2), GDFNB (IDIM1:IDIM2), GDFVD (IDIM1:IDIM2),
214 6 GDFND (IDIM1:IDIM2)
215 C --- LOCAL VARIABLES
216 D I M E N S I O N
Page 2 Source Listing SWR93
2025-03-12 18:23 SWR93.F
217 1 PP (IDIM1:IDIM2,LP1), DP (IDIM1:IDIM2,LP1),
218 2 PR2 (IDIM1:IDIM2,LP1)
219 3, DU (IDIM1:IDIM2,LP1), DUCO2 (IDIM1:IDIM2,LP1),
220 4 DUO3 (IDIM1:IDIM2,LP1)
221 5, FF (IDIM1:IDIM2,LP1), FFCO2 (IDIM1:IDIM2,LP1),
222 6 FFO3 (IDIM1:IDIM2,LP1)
223 7, RRAY (IDIM1:IDIM2), DFNTOP(IDIM1:IDIM2,NB),
224 8 SECZ (IDIM1:IDIM2)
225 9, REFL (IDIM1:IDIM2), TMP1 (IDIM1:IDIM2),
226 o REFL2 (IDIM1:IDIM2)
227 1, CCMAX (IDIM1:IDIM2), XAMT (IDIM1:IDIM2,LP1)
228 D I M E N S I O N
229 1 UD (IDIM1:IDIM2,LP1), UR (IDIM1:IDIM2,LP1)
230 2, UCO2 (IDIM1:IDIM2,LLP2), UDCO2 (IDIM1:IDIM2,LP1),
231 3 URCO2 (IDIM1:IDIM2,LP1)
232 4, UO3 (IDIM1:IDIM2,LLP2), UDO3 (IDIM1:IDIM2,LP1),
233 5 URO3 (IDIM1:IDIM2,LP1)
234 6, TCO2 (IDIM1:IDIM2,LLP2), TDCO2 (IDIM1:IDIM2,LP1),
235 7 TUCO2 (IDIM1:IDIM2,LP1)
236 8, TO3 (IDIM1:IDIM2,LLP2), TDO3 (IDIM1:IDIM2,LP1),
237 9 TUO3 (IDIM1:IDIM2,LP1)
238 D I M E N S I O N
239 1 DFN (IDIM1:IDIM2,LP1), UFN (IDIM1:IDIM2,LP1),
240 2 CR (IDIM1:IDIM2,LP1)
241 3, TTD (IDIM1:IDIM2,LP1), TTU (IDIM1:IDIM2,LP1),
242 4 CT (IDIM1:IDIM2,LP1)
243 5, PPTOP (IDIM1:IDIM2,LP1), DPCLD (IDIM1:IDIM2,LP1)
244 C --- EQUIVALENCED LOCAL VARIABLES
245 D I M E N S I O N
246 1 TTUB1 (IDIM1:IDIM2,LP1), TUCL1 (IDIM1:IDIM2,LP1)
247 2, TTDB1 (IDIM1:IDIM2,LP1), TDCL1 (IDIM1:IDIM2,LP1),
248 3 TDCL2 (IDIM1:IDIM2,LP1)
249 4, UFNTRN(IDIM1:IDIM2,LP1), UFNCLU(IDIM1:IDIM2,LP1),
250 5 TCLU (IDIM1:IDIM2,LP1)
251 6, DFNTRN(IDIM1:IDIM2,LP1), DFNCLU(IDIM1:IDIM2,LP1),
252 7 TCLD (IDIM1:IDIM2,LP1)
253 8, ALFA (IDIM1:IDIM2,LP1), ALFAU (IDIM1:IDIM2,LP1)
254 E Q U I V A L E N C E
255 1 (UDO3,UO3(IDIM1,1),DFNCLU), (URO3,UO3(IDIM1,LP2), UFNCLU)
256 2, (UDCO2,UCO2(IDIM1,1),TCLD), (URCO2,UCO2(IDIM1,LP2), TCLU)
257 3, (TDO3 ,TO3(IDIM1,1),DFNTRN),(TUO3,TO3(IDIM1,LP2), UFNTRN)
258 4, (TDCO2,TCO2(IDIM1,1) ),(TUCO2,TCO2(IDIM1,LP2) )
259 5, (FF , ALFA ), (FFCO2 , ALFAU ), (FFO3 , TTDB1 )
260 6, (DU , TTUB1), (DUCO2 , TUCL1 ), (DUO3 , TDCL1 )
261 7, (PR2 , TDCL2)
262 C
263 C---COMMON FOR LOCAL DATA VARIABLES---
264 COMMON /SWRSAV/ ABCFF(NB),PWTS(NB),CFCO2,CFO3,REFLO3,RRAYAV
265 save /SWRSAV/
266 C D A T A
267 C 1 ABCFF / 2*4.0E-5, 0.002, 0.035, 0.377, 1.95, 9.40, 44.6,
268 C 1 190.0, 989.0, 2706.0, 39011.0 /
269 C 2, PWTS / 0.5000, 0.121416, 0.0698, 0.1558, 0.0631, 0.0362,
270 C 2 0.0243, 0.0158, 0.0087, 0.001467, 0.002342, 0.001075 /
271 C 3, CFCO2, CFO3, REFLO3, RRAYAV / 508.96, 466.64, 1.9, 0.144 /
272 C 1 ABCFF / 2*4.0E-5, .002, .035, .377, 1.95, 9.40, 44.6, 190. /
273 C 2, PWTS /.5000,.1470,.698,.1443,.0584,.0335,.0225,.0158,.0087/
Page 3 Source Listing SWR93
2025-03-12 18:23 SWR93.F
274 C 3, CFCO2, CFO3, REFLO3, RRAYAV / 508.96, 466.64, 1.9, 0.144 /
275 C
276 C CALCULATE SECANT OF ZENITH ANGLE (SECZ),FLUX PRESSURES(PP),
277 C LAYER WIDTH (DP) AND PRESSURE SCALING FACTOR (PR2).
278 DO 100 I=MYIS,MYIE
279 SECZ(I) = H35E1/SQRT(H1224E3*COSZRO(I)*COSZRO(I)+ONE)
280 PP(I,1) = ZERO
281 PP(I,LP1) = PRESS(I,LP1)
282 TMP1(I) = ONE/PRESS(I,LP1)
283 100 CONTINUE
284 DO 110 K=1,LM1
285 DO 110 I=MYIS,MYIE
286 PP(I,K+1) = HAF*(PRESS(I,K+1)+PRESS(I,K))
287 110 CONTINUE
288 DO 120 K=1,L
289 DO 120 I=MYIS,MYIE
290 DP (I,K) = PP(I,K+1)-PP(I,K)
291 PR2(I,K) = HAF*(PP(I,K)+PP(I,K+1))
292 120 CONTINUE
293 DO 130 K=1,L
294 DO 130 I=MYIS,MYIE
295 PR2(I,K) = PR2(I,K)*TMP1(I)
296 130 CONTINUE
297 C CALCULATE ENTERING FLUX AT THE TOP FOR EACH BAND(IN CGS UNITS)
298 DO 140 N=1,NB
299 DO 140 IP=MYIS,MYIE
300 DFNTOP(IP,N) = SSOLAR*H69766E5*COSZRO(IP)*TAUDAR(IP)*PWTS(N)
301 140 CONTINUE
302 C EXECUTE THE LACIS-HANSEN REFLECTIVITY PARAMETERIZATION
303 C FOR THE VISIBLE BAND
304 DO 150 I=MYIS,MYIE
305 RRAY(I) = HP219/(ONE+HP816*COSZRO(I))
306 REFL(I) = RRAY(I) + (ONE-RRAY(I))*(ONE-RRAYAV)*ALVB(I)/
307 1 (ONE-ALVD(I)*RRAYAV)
308 150 CONTINUE
309 DO 155 I=MYIS,MYIE
310 RRAY(I) = 0.104/(ONE+4.8*COSZRO(I))
311 REFL2(I)= RRAY(I) + (ONE-RRAY(I))*(ONE-0.093)*ALVB(I)/
312 1 (ONE-ALVD(I)*0.093)
313 155 CONTINUE
314 C CALCULATE PRESSURE-WEIGHTED OPTICAL PATHS FOR EACH LAYER
315 C IN UNITS OF CM-ATM. PRESSURE WEIGHTING IS USING PR2.
316 C DU= VALUE FOR H2O;DUCO2 FOR CO2;DUO3 FOR O3.
317 DO 160 K=1,L
318 DO 160 I=MYIS,MYIE
319 DU (I,K) = GINV*RH2O(I,K)*DP(I,K)*PR2(I,K)
320 DUCO2(I,K) = (RRCO2*GINV*CFCO2)*DP(I,K)*PR2(I,K)
321 DUO3 (I,K) = (GINV*CFO3)*QO3(I,K)*DP(I,K)
322 160 CONTINUE
323 C
324 C CALCULATE CLEAR SKY SW FLUX
325 C
326 C OBTAIN THE OPTICAL PATH FROM THE TOP OF THE ATMOSPHERE TO THE
327 C FLUX PRESSURE. ANGULAR FACTORS ARE NOW INCLUDED. UD=DOWNWARD
328 C PATH FOR H2O,WIGTH UR THE UPWARD PATH FOR H2O. CORRESPONDING
329 C QUANTITIES FOR CO2,O3 ARE UDCO2/URCO2 AND UDO3/URO3.
330 DO 200 IP=MYIS,MYIE
Page 4 Source Listing SWR93
2025-03-12 18:23 SWR93.F
331 UD (IP,1) = ZERO
332 UDCO2(IP,1) = ZERO
333 UDO3 (IP,1) = ZERO
334 200 CONTINUE
335 DO 210 K=2,LP1
336 DO 210 I=MYIS,MYIE
337 UD (I,K) = UD (I,K-1)+DU (I,K-1)*SECZ(I)
338 UDCO2(I,K) = UDCO2(I,K-1)+DUCO2(I,K-1)*SECZ(I)
339 UDO3 (I,K) = UDO3 (I,K-1)+DUO3 (I,K-1)*SECZ(I)
340 210 CONTINUE
341 DO 220 IP=MYIS,MYIE
342 UR (IP,LP1) = UD (IP,LP1)
343 URCO2(IP,LP1) = UDCO2(IP,LP1)
344 URO3 (IP,LP1) = UDO3 (IP,LP1)
345 220 CONTINUE
346 DO 230 K=L,1,-1
347 DO 230 IP=MYIS,MYIE
348 UR (IP,K) = UR (IP,K+1)+DU (IP,K)*DIFFCTR
349 URCO2(IP,K) = URCO2(IP,K+1)+DUCO2(IP,K)*DIFFCTR
350 URO3 (IP,K) = URO3 (IP,K+1)+DUO3 (IP,K)*O3DIFCTR
351 230 CONTINUE
352 C CALCULATE CO2 ABSORPTIONS . THEY WILL BE USED IN NEAR INFRARED
353 C BANDS.SINCE THE ABSORPTION AMOUNT IS GIVEN (IN THE FORMULA USED
354 C BELOW, DERIVED FROM SASAMORI) IN TERMS OF THE TOTAL SOLAR FLUX,
355 C AND THE ABSORPTION IS ONLY INCLUDED IN THE NEAR IR (50 PERCENT
356 C OF THE SOLAR SPECTRUM), THE ABSORPTIONS ARE MULTIPLIED BY 2.
357 C SINCE CODE ACTUALLY REQUIRES TRANSMISSIONS, THESE ARE THE
358 C VALUES ACTUALLY STORED IN TCO2.
359 DO 240 K=1,LL
360 DO 240 I=MYIS,MYIE
361 TCO2(I,K+1)=ONE-TWO*(H235M3*EXP(HP26*LOG(UCO2(I,K+1)+H129M2))
362 1 -H75826M4)
363 240 CONTINUE
364 C NOW CALCULATE OZONE ABSORPTIONS. THESE WILL BE USED IN
365 C THE VISIBLE BAND.JUST AS IN THE CO2 CASE, SINCE THIS BAND IS
366 C 50 PERCENT OF THE SOLAR SPECTRUM,THE ABSORPTIONS ARE MULTIPLIED
367 C BY 2. THE TRANSMISSIONS ARE STORED IN TO3.
368 HTEMP = H1036E2*H1036E2*H1036E2
369 DO 250 K=1,LL
370 DO 250 I=MYIS,MYIE
371 TO3(I,K+1)=ONE-TWO*UO3(I,K+1)*
372 1 (H1P082*EXP(HMP805*LOG(ONE+H1386E2*UO3(I,K+1)))+
373 2 H658M2/(ONE+HTEMP*UO3(I,K+1)*UO3(I,K+1)*UO3(I,K+1))+
374 3 H2118M2/(ONE+UO3(I,K+1)*(H42M2+H323M4*UO3(I,K+1))))
375 250 CONTINUE
376 C START FREQUENCY LOOP (ON N) HERE
377 C
378 C--- BAND 1 (VISIBLE) INCLUDES O3 AND H2O ABSORPTION
379 DO 260 K=1,L
380 DO 260 I=MYIS,MYIE
381 TTD(I,K+1) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UD(I,K+1)))
382 TTU(I,K) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UR(I,K)))
383 DFN(I,K+1) = TTD(I,K+1)*TDO3(I,K+1)
384 UFN(I,K) = TTU(I,K)*TUO3(I,K)
385 260 CONTINUE
386 DO 270 I=MYIS,MYIE
387 DFN(I,1) = ONE
Page 5 Source Listing SWR93
2025-03-12 18:23 SWR93.F
388 UFN(I,LP1) = DFN(I,LP1)
389 270 CONTINUE
390 C SCALE VISIBLE BAND FLUXES BY SOLAR FLUX AT THE TOP OF THE
391 C ATMOSPHERE (DFNTOP(I,1))
392 C DFSW/UFSW WILL BE THE FLUXES, SUMMED OVER ALL BANDS
393 DO 280 K=1,LP1
394 DO 280 I=MYIS,MYIE
395 DFSWL(I,K) = DFN(I,K)*DFNTOP(I,1)
396 UFSWL(I,K) = REFL(I)*UFN(I,K)*DFNTOP(I,1)
397 280 CONTINUE
398 DO 285 I=MYIS,MYIE
399 GDFVB(I) = DFSWL(I,LP1)*EXP(-0.15746*SECZ(I))
400 GDFVD(I) = ((ONE-REFL2(I))*DFSWL(I,LP1) -
401 1 (ONE-ALVB(I)) *GDFVB(I)) / (ONE-ALVD(I))
402 GDFNB(I) = ZERO
403 GDFND(I) = ZERO
404 285 CONTINUE
405 C---NOW OBTAIN FLUXES FOR THE NEAR IR BANDS. THE METHODS ARE THE SAME
406 C AS FOR THE VISIBLE BAND, EXCEPT THAT THE REFLECTION AND
407 C TRANSMISSION COEFFICIENTS (OBTAINED BELOW) ARE DIFFERENT, AS
408 C RAYLEIGH SCATTERING NEED NOT BE CONSIDERED.
409 DO 350 N=2,NB
410 IF (N.EQ.2) THEN
411 C THE WATER VAPOR TRANSMISSION FUNCTION FOR BAND 2 IS EQUAL TO
412 C THAT OF BAND 1 (SAVED AS TTD,TTU)
413 C--- BAND 2-9 (NEAR-IR) INCLUDES O3, CO2 AND H2O ABSORPTION
414 DO 290 K=1,L
415 DO 290 I=MYIS,MYIE
416 DFN(I,K+1) = TTD(I,K+1)*TDCO2(I,K+1)
417 UFN(I,K) = TTU(I,K)*TUCO2(I,K)
418 290 CONTINUE
419 ELSE
420 C CALCULATE WATER VAPOR TRANSMISSION FUNCTIONS FOR NEAR INFRARED
421 C BANDS. INCLUDE CO2 TRANSMISSION (TDCO2/TUCO2), WHICH
422 C IS THE SAME FOR ALL INFRARED BANDS.
423 DO 300 K=1,L
424 DO 300 I=MYIS,MYIE
425 DFN(I,K+1)=EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UD(I,K+1)))
426 1 *TDCO2(I,K+1)
427 UFN(I,K)=EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UR(I,K)))
428 1 *TUCO2(I,K)
429 300 CONTINUE
430 ENDIF
431 C---AT THIS POINT,INCLUDE DFN(1),UFN(LP1), NOTING THAT DFN(1)=1 FOR
432 C ALL BANDS, AND THAT UFN(LP1)=DFN(LP1) FOR ALL BANDS.
433 DO 310 I=MYIS,MYIE
434 DFN(I,1) = ONE
435 UFN(I,LP1) = DFN(I,LP1)
436 310 CONTINUE
437 C SCALE THE PREVIOUSLY COMPUTED FLUXES BY THE FLUX AT THE TOP
438 C AND SUM OVER BANDS
439 DO 320 K=1,LP1
440 DO 320 I=MYIS,MYIE
441 DFSWL(I,K) = DFSWL(I,K) + DFN(I,K)*DFNTOP(I,N)
442 UFSWL(I,K) = UFSWL(I,K) + ALNB(I)*UFN(I,K)*DFNTOP(I,N)
443 320 CONTINUE
444 DO 330 I=MYIS,MYIE
Page 6 Source Listing SWR93
2025-03-12 18:23 SWR93.F
445 GDFNB(I) = GDFNB(I) + DFN(I,LP1)*DFNTOP(I,N)
446 330 CONTINUE
447 350 CONTINUE
448 DO 360 K=1,LP1
449 DO 360 I=MYIS,MYIE
450 FSWL(I,K) = UFSWL(I,K)-DFSWL(I,K)
451 360 CONTINUE
452 DO 370 K=1,L
453 DO 370 I=MYIS,MYIE
454 HSWL(I,K)=RADCON*(FSWL(I,K+1)-FSWL(I,K))/DP(I,K)
455 370 CONTINUE
456 C
457 C---END OF FREQUENCY LOOP (OVER N)
458 C
459 C CALCULATE CLOUDY SKY SW FLUX
460 C
461 KCLDS=NCLDS(MYIS)
462 DO 400 I=MYIS1,MYIE
463 KCLDS=MAX(NCLDS(I),KCLDS)
464 400 CONTINUE
465 DO 410 K=1,LP1
466 DO 410 I=MYIS,MYIE
467 DFSWC(I,K) = DFSWL(I,K)
468 UFSWC(I,K) = UFSWL(I,K)
469 FSWC (I,K) = FSWL (I,K)
470 410 CONTINUE
471 DO 420 K=1,L
472 DO 420 I=MYIS,MYIE
473 HSWC(I,K) = HSWL(I,K)
474 420 CONTINUE
475 C*******************************************************************
476 IF (KCLDS .EQ. 0) RETURN
477 C*******************************************************************
478 DO 430 K=1,LP1
479 DO 430 I=MYIS,MYIE
480 XAMT(I,K) = CAMT(I,K)
481 430 CONTINUE
482 DO 470 I=MYIS,MYIE
483 NNCLDS = NCLDS(I)
484 CCMAX(I) = ZERO
485 IF (NNCLDS .LE. 0) GO TO 470
486 CCMAX(I) = ONE
487 DO 450 K=1,NNCLDS
488 CCMAX(I) = CCMAX(I) * (ONE - CAMT(I,K+1))
489 450 CONTINUE
490 CCMAX(I) = ONE - CCMAX(I)
491 IF (CCMAX(I) .GT. ZERO) THEN
492 DO 460 K=1,NNCLDS
493 XAMT(I,K+1) = CAMT(I,K+1)/CCMAX(I)
494 460 CONTINUE
495 END IF
496 470 CONTINUE
497 DO 480 K=1,LP1
498 DO 480 I=MYIS,MYIE
499 FF (I,K) = DIFFCTR
500 FFCO2(I,K) = DIFFCTR
501 FFO3 (I,K) = O3DIFCTR
Page 7 Source Listing SWR93
2025-03-12 18:23 SWR93.F
502 480 CONTINUE
503 DO 490 IP=MYIS,MYIE
504 JTOP = KTOPSW(IP,NCLDS(IP)+1)
505 DO 490 K=1,JTOP
506 FF (IP,K) = SECZ(IP)
507 FFCO2(IP,K) = SECZ(IP)
508 FFO3 (IP,K) = SECZ(IP)
509 490 CONTINUE
510 DO 500 I=MYIS,MYIE
511 RRAY(I) = HP219/(ONE+HP816*COSZRO(I))
512 REFL(I) = RRAY(I) + (ONE-RRAY(I))*(ONE-RRAYAV)*ALVD(I)/
513 1 (ONE-ALVD(I)*RRAYAV)
514 500 CONTINUE
515 DO 510 IP=MYIS,MYIE
516 UD (IP,1) = ZERO
517 UDCO2(IP,1) = ZERO
518 UDO3 (IP,1) = ZERO
519 510 CONTINUE
520 DO 520 K=2,LP1
521 DO 520 I=MYIS,MYIE
522 UD (I,K) = UD (I,K-1)+DU (I,K-1)*FF (I,K)
523 UDCO2(I,K) = UDCO2(I,K-1)+DUCO2(I,K-1)*FFCO2(I,K)
524 UDO3 (I,K) = UDO3 (I,K-1)+DUO3 (I,K-1)*FFO3 (I,K)
525 520 CONTINUE
526 DO 530 IP=MYIS,MYIE
527 UR (IP,LP1) = UD (IP,LP1)
528 URCO2(IP,LP1) = UDCO2(IP,LP1)
529 URO3 (IP,LP1) = UDO3 (IP,LP1)
530 530 CONTINUE
531 DO 540 K=L,1,-1
532 DO 540 IP=MYIS,MYIE
533 UR (IP,K) = UR (IP,K+1)+DU (IP,K)*DIFFCTR
534 URCO2(IP,K) = URCO2(IP,K+1)+DUCO2(IP,K)*DIFFCTR
535 URO3 (IP,K) = URO3 (IP,K+1)+DUO3 (IP,K)*O3DIFCTR
536 540 CONTINUE
537 DO 550 K=1,LL
538 DO 550 I=MYIS,MYIE
539 TCO2(I,K+1)=ONE-TWO*(H235M3*EXP(HP26*LOG(UCO2(I,K+1)+H129M2))
540 1 -H75826M4)
541 550 CONTINUE
542 DO 560 K=1,LL
543 DO 560 I=MYIS,MYIE
544 TO3(I,K+1)=ONE-TWO*UO3(I,K+1)*
545 1 (H1P082*EXP(HMP805*LOG(ONE+H1386E2*UO3(I,K+1)))+
546 2 H658M2/(ONE+HTEMP*UO3(I,K+1)*UO3(I,K+1)*UO3(I,K+1))+
547 3 H2118M2/(ONE+UO3(I,K+1)*(H42M2+H323M4*UO3(I,K+1))))
548 560 CONTINUE
549 C********************************************************************
550 C---THE FIRST CLOUD IS THE GROUND; ITS PROPERTIES ARE GIVEN
551 C BY REFL (THE TRANSMISSION (0) IS IRRELEVANT FOR NOW!).
552 C********************************************************************
553 DO 570 I=MYIS,MYIE
554 CR(I,1) = REFL(I)
555 570 CONTINUE
556 C***OBTAIN CLOUD REFLECTION AND TRANSMISSION COEFFICIENTS FOR
557 C REMAINING CLOUDS (IF ANY) IN THE VISIBLE BAND
558 C---THE MAXIMUM NO OF CLOUDS IN THE ROW (KCLDS) IS USED. THIS CREATES
Page 8 Source Listing SWR93
2025-03-12 18:23 SWR93.F
559 C EXTRA WORK (MAY BE REMOVED IN A SUBSEQUENT UPDATE).
560 DO 581 I=MYIS,MYIE
561 KCLDS=NCLDS(I)
562 IF(KCLDS.EQ.0) GO TO 581
563 DO 580 KK=2,KCLDS+1
564 CR(I,KK) = CRR(I,1,KK)*XAMT(I,KK)
565 CT(I,KK) = ONE - (ONE-CTT(I,1,KK))*XAMT(I,KK)
566 580 CONTINUE
567 581 CONTINUE
568 C---OBTAIN THE PRESSURE AT THE TOP,BOTTOM AND THE THICKNESS OF
569 C "THICK" CLOUDS (THOSE AT LEAST 2 LAYERS THICK). THIS IS USED
570 C LATER IS OBTAINING FLUXES INSIDE THE THICK CLOUDS, FOR ALL
571 C FREQUENCY BANDS.
572 DO 591 I=MYIS,MYIE
573 KCLDS=NCLDS(I)
574 IF(KCLDS.EQ.0) GO TO 591
575 DO 590 KK=1,KCLDS
576 IF ((KBTMSW(I,KK+1)-1).GT.KTOPSW(I,KK+1)) THEN
577 PPTOP(I,KK)=PP(I,KTOPSW(I,KK+1))
578 DPCLD(I,KK)=ONE/(PPTOP(I,KK)-PP(I,KBTMSW(I,KK+1)))
579 ENDIF
580 590 CONTINUE
581 591 CONTINUE
582 DO 600 K=1,L
583 DO 600 I=MYIS,MYIE
584 TTDB1(I,K+1) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UD(I,K+1)))
585 TTUB1(I,K) = EXP(HM1EZ*MIN(FIFTY,ABCFF(1)*UR(I,K)))
586 TTD (I,K+1) = TTDB1(I,K+1)*TDO3(I,K+1)
587 TTU (I,K) = TTUB1(I,K)*TUO3(I,K)
588 600 CONTINUE
589 DO 610 I=MYIS,MYIE
590 TTD(I,1) = ONE
591 TTU(I,LP1) = TTD(I,LP1)
592 610 CONTINUE
593 C***FOR EXECUTION OF THE CLOUD LOOP, IT IS NECESSARY TO SEPARATE OUT
594 C TRANSMISSION FCTNS AT THE TOP AND BOTTOM OF THE CLOUDS, FOR
595 C EACH BAND N. THE REQUIRED QUANTITIES ARE:
596 C TTD(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1:
597 C TTU(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1:
598 C TTD(I,KBTMSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1:
599 C AND INVERSES OF THE FIRST TWO. THE ABOVE QUANTITIES ARE
600 C STORED IN TDCL1,TUCL1,TDCL2, AND DFNTRN,UFNTRN, RESPECTIVELY,
601 C AS THEY HAVE MULTIPLE USE IN THE PGM.
602 C---FOR FIRST CLOUD LAYER (GROUND) TDCL1,TUCL1 ARE KNOWN:
603 DO 620 I=MYIS,MYIE
604 TDCL1 (I,1) = TTD(I,LP1)
605 TUCL1 (I,1) = TTU(I,LP1)
606 TDCL2 (I,1) = TDCL1(I,1)
607 DFNTRN(I,1) = ONE/TDCL1(I,1)
608 UFNTRN(I,1) = DFNTRN(I,1)
609 620 CONTINUE
610 DO 631 I=MYIS,MYIE
611 KCLDS=NCLDS(I)
612 IF(KCLDS.EQ.0) GO TO 631
613 DO 630 KK=2,KCLDS+1
614 TDCL1(I,KK) = TTD(I,KTOPSW(I,KK))
615 TUCL1(I,KK) = TTU(I,KTOPSW(I,KK))
Page 9 Source Listing SWR93
2025-03-12 18:23 SWR93.F
616 TDCL2(I,KK) = TTD(I,KBTMSW(I,KK))
617 630 CONTINUE
618 631 CONTINUE
619 C---COMPUTE INVERSES
620 DO 641 I=MYIS,MYIE
621 KCLDS=NCLDS(I)
622 IF(KCLDS.EQ.0) GO TO 641
623 DO 640 KK=2,KCLDS
624 DFNTRN(I,KK) = ONE/TDCL1(I,KK)
625 UFNTRN(I,KK) = ONE/TUCL1(I,KK)
626 640 CONTINUE
627 641 CONTINUE
628 C---COMPUTE THE TRANSMISSIVITY FROM THE TOP OF CLOUD (K+1) TO THE
629 C TOP OF CLOUD (K). THE CLOUD TRANSMISSION (CT) IS INCLUDED. THIS
630 C QUANTITY IS CALLED TCLU (INDEX K). ALSO, OBTAIN THE TRANSMISSIVITY
631 C FROM THE BOTTOM OF CLOUD (K+1) TO THE TOP OF CLOUD (K)(A PATH
632 C ENTIRELY OUTSIDE CLOUDS). THIS QUANTITY IS CALLED TCLD (INDEX K).
633 DO 651 I=MYIS,MYIE
634 KCLDS=NCLDS(I)
635 IF(KCLDS.EQ.0) GO TO 651
636 DO 650 KK=1,KCLDS
637 TCLU(I,KK) = TDCL1(I,KK)*DFNTRN(I,KK+1)*CT(I,KK+1)
638 TCLD(I,KK) = TDCL1(I,KK)/TDCL2(I,KK+1)
639 650 CONTINUE
640 651 CONTINUE
641 C***THE FOLLOWING IS THE RECURSION RELATION FOR ALFA: THE REFLECTION
642 C COEFFICIENT FOR A SYSTEM INCLUDING THE CLOUD IN QUESTION AND THE
643 C FLUX COMING OUT OF THE CLOUD SYSTEM INCLUDING ALL CLOUDS BELOW
644 C THE CLOUD IN QUESTION.
645 C---ALFAU IS ALFA WITHOUT THE REFLECTION OF THE CLOUD IN QUESTION
646 DO 660 I=MYIS,MYIE
647 KCLDS=NCLDS(I)
648 IF(KCLDS.EQ.0) GO TO 660
649 ALFA (I,1)=CR(I,1)
650 ALFAU(I,1)=ZERO
651 660 CONTINUE
652 C---AGAIN,EXCESSIVE CALCULATIONS-MAY BE CHANGED LATER!
653 DO 671 I=MYIS,MYIE
654 KCLDS=NCLDS(I)
655 IF(KCLDS.EQ.0) GO TO 671
656 DO 670 KK=2,KCLDS+1
657 ALFAU(I,KK)= TCLU(I,KK-1)*TCLU(I,KK-1)*ALFA(I,KK-1)/
658 1 (ONE - TCLD(I,KK-1)*TCLD(I,KK-1)*ALFA(I,KK-1)*CR(I,KK))
659 ALFA (I,KK)= ALFAU(I,KK)+CR(I,KK)
660 670 CONTINUE
661 671 CONTINUE
662 C CALCULATE UFN AT CLOUD TOPS AND DFN AT CLOUD BOTTOMS
663 C---NOTE THAT UFNCLU(I,KCLDS+1) GIVES THE UPWARD FLUX AT THE TOP
664 C OF THE HIGHEST REAL CLOUD (IF NCLDS(I)=KCLDS). IT GIVES THE FLUX
665 C AT THE TOP OF THE ATMOSPHERE IF NCLDS(I) < KCLDS. IN THE FIRST
666 C CASE, TDCL1 EQUALS THE TRANSMISSION FCTN TO THE TOP OF THE
667 C HIGHEST CLOUD, AS WE WANT. IN THE SECOND CASE, TDCL1=1, SO UFNCLU
668 C EQUALS ALFA. THIS IS ALSO CORRECT.
669 DO 680 I=MYIS,MYIE
670 KCLDS=NCLDS(I)
671 IF(KCLDS.EQ.0) GO TO 680
672 UFNCLU(I,KCLDS+1) = ALFA(I,KCLDS+1)*TDCL1(I,KCLDS+1)
Page 10 Source Listing SWR93
2025-03-12 18:23 SWR93.F
673 DFNCLU(I,KCLDS+1) = TDCL1(I,KCLDS+1)
674 680 CONTINUE
675 C---THIS CALCULATION IS THE REVERSE OF THE RECURSION RELATION USED
676 C ABOVE
677 DO 691 I=MYIS,MYIE
678 KCLDS=NCLDS(I)
679 IF(KCLDS.EQ.0) GO TO 691
680 DO 690 KK=KCLDS,1,-1
681 UFNCLU(I,KK) = UFNCLU(I,KK+1)*ALFAU(I,KK+1)/(ALFA(I,KK+1)*
682 1 TCLU(I,KK))
683 DFNCLU(I,KK) = UFNCLU(I,KK)/ALFA(I,KK)
684 690 CONTINUE
685 691 CONTINUE
686 DO 701 I=MYIS,MYIE
687 KCLDS=NCLDS(I)
688 IF(KCLDS.EQ.0) GO TO 701
689 DO 700 KK=1,KCLDS+1
690 UFNTRN(I,KK) = UFNCLU(I,KK)*UFNTRN(I,KK)
691 DFNTRN(I,KK) = DFNCLU(I,KK)*DFNTRN(I,KK)
692 700 CONTINUE
693 701 CONTINUE
694 C---CASE OF KK=1( FROM THE GROUND TO THE BOTTOM OF THE LOWEST CLOUD)
695 DO 720 I=MYIS,MYIE
696 KCLDS=NCLDS(I)
697 IF(KCLDS.EQ.0) GO TO 720
698 J2=KBTMSW(I,2)
699 DO 710 K=J2,LP1
700 UFN(I,K) = UFNTRN(I,1)*TTU(I,K)
701 DFN(I,K) = DFNTRN(I,1)*TTD(I,K)
702 710 CONTINUE
703 720 CONTINUE
704 C---REMAINING LEVELS (IF ANY!)
705 DO 760 I=MYIS,MYIE
706 KCLDS=NCLDS(I)
707 IF(KCLDS.EQ.0) GO TO 760
708 DO 755 KK=2,KCLDS+1
709 J1=KTOPSW(I,KK)
710 J2=KBTMSW(I,KK+1)
711 IF (J1.EQ.1) GO TO 755
712 DO 730 K=J2,J1
713 UFN(I,K) = UFNTRN(I,KK)*TTU(I,K)
714 DFN(I,K) = DFNTRN(I,KK)*TTD(I,K)
715 730 CONTINUE
716 C---FOR THE THICK CLOUDS, THE FLUX DIVERGENCE THROUGH THE CLOUD
717 C LAYER IS ASSUMED TO BE CONSTANT. THE FLUX DERIVATIVE IS GIVEN BY
718 C TEMPF (FOR THE UPWARD FLUX) AND TEMPG (FOR THE DOWNWARD FLUX).
719 J3=KBTMSW(I,KK)
720 IF ((J3-J1).GT.1) THEN
721 TEMPF = (UFNCLU(I,KK)-UFN(I,J3))*DPCLD(I,KK-1)
722 TEMPG = (DFNCLU(I,KK)-DFN(I,J3))*DPCLD(I,KK-1)
723 DO 740 K=J1+1,J3-1
724 UFN(I,K) = UFNCLU(I,KK)+TEMPF*(PP(I,K)-PPTOP(I,KK-1))
725 DFN(I,K) = DFNCLU(I,KK)+TEMPG*(PP(I,K)-PPTOP(I,KK-1))
726 740 CONTINUE
727 ENDIF
728 755 CONTINUE
729 760 CONTINUE
Page 11 Source Listing SWR93
2025-03-12 18:23 SWR93.F
730 DO 770 I=MYIS,MYIE
731 KCLDS=NCLDS(I)
732 IF(KCLDS.EQ.0) GO TO 770
733 DO 771 K=1,LP1
734 DFSWC(I,K) = DFN(I,K)*DFNTOP(I,1)
735 UFSWC(I,K) = UFN(I,K)*DFNTOP(I,1)
736 771 CONTINUE
737 770 CONTINUE
738 DO 780 I=MYIS,MYIE
739 KCLDS=NCLDS(I)
740 IF(KCLDS.EQ.0) GO TO 780
741 TMP1(I) = ONE - CCMAX(I)
742 GDFVB(I) = TMP1(I)*GDFVB(I)
743 GDFNB(I) = TMP1(I)*GDFNB(I)
744 GDFVD(I) = TMP1(I)*GDFVD(I) + CCMAX(I)*DFSWC(I,LP1)
745 780 CONTINUE
746 C---NOW OBTAIN FLUXES FOR THE NEAR IR BANDS. THE METHODS ARE THE SAME
747 C AS FOR THE VISIBLE BAND, EXCEPT THAT THE REFLECTION AND
748 C TRANSMISSION COEFFICIENTS ARE DIFFERENT, AS
749 C RAYLEIGH SCATTERING NEED NOT BE CONSIDERED.
750 C
751 DO 1000 N=2,NB
752 CYH93
753 DO 791 I=MYIS,MYIE
754 KCLDS=NCLDS(I)
755 IF(KCLDS.EQ.0) GO TO 791
756 DO 790 K=1,KCLDS+1
757 CR(I,K) = CRR(I,N,K)*XAMT(I,K)
758 CT(I,K) = ONE - (ONE-CTT(I,N,K))*XAMT(I,K)
759 790 CONTINUE
760 791 CONTINUE
761 CYH93
762 IF (N.EQ.2) THEN
763 C THE WATER VAPOR TRANSMISSION FUNCTION FOR BAND 2 IS EQUAL TO
764 C THAT OF BAND 1 (SAVED AS TTDB1,TTUB1)
765 DO 800 I=MYIS,MYIE
766 KCLDS=NCLDS(I)
767 IF(KCLDS.EQ.0) GO TO 800
768 DO 801 KK=2,LP1
769 TTD(I,KK) = TTDB1(I,KK)*TDCO2(I,KK)
770 801 CONTINUE
771 DO 802 KK=1,L
772 TTU(I,KK) = TTUB1(I,KK)*TUCO2(I,KK)
773 802 CONTINUE
774 800 CONTINUE
775 ELSE
776 DO 810 I=MYIS,MYIE
777 KCLDS=NCLDS(I)
778 IF(KCLDS.EQ.0) GO TO 810
779 DO 811 KK=2,LP1
780 TTD(I,KK) = EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UD(I,KK)))
781 1 * TDCO2(I,KK)
782 811 CONTINUE
783 DO 812 KK=1,L
784 TTU(I,KK) = EXP(HM1EZ*MIN(FIFTY,ABCFF(N)*UR(I,KK)))
785 1 * TUCO2(I,KK)
786 812 CONTINUE
Page 12 Source Listing SWR93
2025-03-12 18:23 SWR93.F
787 810 CONTINUE
788 ENDIF
789 C---AT THIS POINT,INCLUDE TTD(1),TTU(LP1), NOTING THAT TTD(1)=1 FOR
790 C ALL BANDS, AND THAT TTU(LP1)=TTD(LP1) FOR ALL BANDS.
791 DO 820 I=MYIS,MYIE
792 KCLDS=NCLDS(I)
793 IF(KCLDS.EQ.0) GO TO 820
794 TTU(I,LP1) = TTD(I,LP1)
795 TTD(I,1) = ONE
796 820 CONTINUE
797 C***FOR EXECUTION OF THE CLOUD LOOP, IT IS NECESSARY TO SEPARATE OUT
798 C TRANSMISSION FCTNS AT THE TOP AND BOTTOM OF THE CLOUDS, FOR
799 C EACH BAND N. THE REQUIRED QUANTITIES ARE:
800 C TTD(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1:
801 C TTD(I,KBTMSW(I,K),N) K RUNS FROM 2 TO NCLDS(I)+1:
802 C TTU(I,KTOPSW(I,K),N) K RUNS FROM 1 TO NCLDS(I)+1:
803 C AND INVERSES OF THE ABOVE. THE ABOVE QUANTITIES ARE STORED
804 C IN TDCL1,TDCL2,TUCL1,AND DFNTRN,UFNTRN,RESPECTIVELY, AS
805 C THEY HAVE MULTIPLE USE IN THE PGM.
806 C---FOR FIRST CLOUD LAYER (GROUND) TDCL1,TUCL1 ARE KNOWN:
807 DO 830 I=MYIS,MYIE
808 KCLDS=NCLDS(I)
809 IF(KCLDS.EQ.0) GO TO 830
810 TDCL1 (I,1) = TTD(I,LP1)
811 TUCL1 (I,1) = TTU(I,LP1)
812 TDCL2 (I,1) = TDCL1(I,1)
813 DFNTRN(I,1) = ONE/TDCL1(I,1)
814 UFNTRN(I,1) = DFNTRN(I,1)
815 830 CONTINUE
816 DO 841 I=MYIS,MYIE
817 KCLDS=NCLDS(I)
818 IF(KCLDS.EQ.0) GO TO 841
819 DO 840 KK=2,KCLDS+1
820 TDCL1(I,KK) = TTD(I,KTOPSW(I,KK))
821 TUCL1(I,KK) = TTU(I,KTOPSW(I,KK))
822 TDCL2(I,KK) = TTD(I,KBTMSW(I,KK))
823 840 CONTINUE
824 841 CONTINUE
825 DO 851 I=MYIS,MYIE
826 KCLDS=NCLDS(I)
827 IF(KCLDS.EQ.0) GO TO 851
828 DO 850 KK=2,KCLDS+1
829 DFNTRN(I,KK) = ONE/TDCL1(I,KK)
830 UFNTRN(I,KK) = ONE/TUCL1(I,KK)
831 850 CONTINUE
832 851 CONTINUE
833 DO 861 I=MYIS,MYIE
834 KCLDS=NCLDS(I)
835 IF(KCLDS.EQ.0) GO TO 861
836 DO 860 KK=1,KCLDS
837 TCLU(I,KK) = TDCL1(I,KK)*DFNTRN(I,KK+1)*CT(I,KK+1)
838 TCLD(I,KK) = TDCL1(I,KK)/TDCL2(I,KK+1)
839 860 CONTINUE
840 861 CONTINUE
841 C***THE FOLLOWING IS THE RECURSION RELATION FOR ALFA: THE REFLECTION
842 C COEFFICIENT FOR A SYSTEM INCLUDING THE CLOUD IN QUESTION AND THE
843 C FLUX COMING OUT OF THE CLOUD SYSTEM INCLUDING ALL CLOUDS BELOW
Page 13 Source Listing SWR93
2025-03-12 18:23 SWR93.F
844 C THE CLOUD IN QUESTION.
845 DO 870 I=MYIS,MYIE
846 KCLDS=NCLDS(I)
847 IF(KCLDS.EQ.0) GO TO 870
848 ALFA (I,1) = CR(I,1)
849 ALFAU(I,1) = ZERO
850 870 CONTINUE
851 C---AGAIN,EXCESSIVE CALCULATIONS-MAY BE CHANGED LATER!
852 DO 881 I=MYIS,MYIE
853 KCLDS=NCLDS(I)
854 IF(KCLDS.EQ.0) GO TO 881
855 DO 880 KK=2,KCLDS+1
856 ALFAU(I,KK) = TCLU(I,KK-1)*TCLU(I,KK-1)*ALFA(I,KK-1)/(ONE -
857 1 TCLD(I,KK-1)*TCLD(I,KK-1)*ALFA(I,KK-1)*CR(I,KK))
858 ALFA (I,KK) = ALFAU(I,KK)+CR(I,KK)
859 880 CONTINUE
860 881 CONTINUE
861 C CALCULATE UFN AT CLOUD TOPS AND DFN AT CLOUD BOTTOMS
862 C---NOTE THAT UFNCLU(I,KCLDS+1) GIVES THE UPWARD FLUX AT THE TOP
863 C OF THE HIGHEST REAL CLOUD (IF NCLDS(I)=KCLDS). IT GIVES THE FLUX
864 C AT THE TOP OF THE ATMOSPHERE IF NCLDS(I) < KCLDS. IT THE FIRST
865 C CASE, TDCL1 EQUALS THE TRANSMISSION FCTN TO THE TOP OF THE
866 C HIGHEST CLOUD, AS WE WANT. IN THE SECOND CASE, TDCL1=1, SO UFNCLU
867 C EQUALS ALFA. THIS IS ALSO CORRECT.
868 DO 890 I=MYIS,MYIE
869 KCLDS=NCLDS(I)
870 IF(KCLDS.EQ.0) GO TO 890
871 UFNCLU(I,KCLDS+1) = ALFA(I,KCLDS+1)*TDCL1(I,KCLDS+1)
872 DFNCLU(I,KCLDS+1) = TDCL1(I,KCLDS+1)
873 890 CONTINUE
874 DO 901 I=MYIS,MYIE
875 KCLDS=NCLDS(I)
876 IF(KCLDS.EQ.0) GO TO 901
877 DO 900 KK=KCLDS,1,-1
878 !
879 !--- Ferrier, 6/17/02: Emergency change to eliminate problematic
880 ! features of unrealistically small cloud amounts
881 !
882 DENOM=ALFA(I,KK+1)*TCLU(I,KK)
883 IF (DENOM .GT. 1.E-15) THEN
884 UFNCLU(I,KK)=UFNCLU(I,KK+1)*ALFAU(I,KK+1)/DENOM
885 ELSE
886 c print *
887 c &, ' xnum=',UFNCLU(I,KK+1)*ALFAU(I,KK+1)
888 c &, ' xden=',DENOM
889 c &, ' ALFA(I,KK+1)=',ALFA(I,KK+1)
890 c &, ' TCLU(I,KK)=',TCLU(I,KK)
891 c &, ' UFNCLU(I,KK+1)=',UFNCLU(I,KK+1)
892 c &, ' ALFAU(I,KK+1)=',ALFAU(I,KK+1)
893 c &, ' UFNCLU(I,KK)=',UFNCLU(I,KK)
894 c &, ' ALFA(I,KK)=',ALFA(I,KK)
895 c &, ' I=',I
896 c &, ' KK+1=',KK+1
897 UFNCLU(I,KK)=0.
898 ENDIF
899 DFNCLU(I,KK) = UFNCLU(I,KK)/ALFA(I,KK)
900 900 CONTINUE
Page 14 Source Listing SWR93
2025-03-12 18:23 SWR93.F
901 901 CONTINUE
902 C NOW OBTAIN DFN AND UFN FOR LEVELS BETWEEN THE CLOUDS
903 DO 911 I=MYIS,MYIE
904 KCLDS=NCLDS(I)
905 IF(KCLDS.EQ.0) GO TO 911
906 DO 910 KK=1,KCLDS+1
907 UFNTRN(I,KK) = UFNCLU(I,KK)*UFNTRN(I,KK)
908 DFNTRN(I,KK) = DFNCLU(I,KK)*DFNTRN(I,KK)
909 910 CONTINUE
910 911 CONTINUE
911 DO 930 I=MYIS,MYIE
912 KCLDS=NCLDS(I)
913 IF(KCLDS.EQ.0) GO TO 930
914 J2=KBTMSW(I,2)
915 DO 920 K=J2,LP1
916 UFN(I,K) = UFNTRN(I,1)*TTU(I,K)
917 DFN(I,K) = DFNTRN(I,1)*TTD(I,K)
918 920 CONTINUE
919 930 CONTINUE
920 DO 970 I=MYIS,MYIE
921 KCLDS=NCLDS(I)
922 IF(KCLDS.EQ.0) GO TO 970
923 DO 965 KK=2,KCLDS+1
924 J1 = KTOPSW(I,KK)
925 J2 = KBTMSW(I,KK+1)
926 IF (J1.EQ.1) GO TO 965
927 DO 940 K=J2,J1
928 UFN(I,K) = UFNTRN(I,KK)*TTU(I,K)
929 DFN(I,K) = DFNTRN(I,KK)*TTD(I,K)
930 940 CONTINUE
931 J3 = KBTMSW(I,KK)
932 IF ((J3-J1).GT.1) THEN
933 TEMPF = (UFNCLU(I,KK)-UFN(I,J3))*DPCLD(I,KK-1)
934 TEMPG = (DFNCLU(I,KK)-DFN(I,J3))*DPCLD(I,KK-1)
935 DO 950 K=J1+1,J3-1
936 UFN(I,K) = UFNCLU(I,KK)+TEMPF*(PP(I,K)-PPTOP(I,KK-1))
937 DFN(I,K) = DFNCLU(I,KK)+TEMPG*(PP(I,K)-PPTOP(I,KK-1))
938 950 CONTINUE
939 ENDIF
940 965 CONTINUE
941 970 CONTINUE
942 DO 980 I=MYIS,MYIE
943 KCLDS=NCLDS(I)
944 IF(KCLDS.EQ.0) GO TO 980
945 DO 981 K=1,LP1
946 DFSWC(I,K) = DFSWC(I,K) + DFN(I,K)*DFNTOP(I,N)
947 UFSWC(I,K) = UFSWC(I,K) + UFN(I,K)*DFNTOP(I,N)
948 981 CONTINUE
949 980 CONTINUE
950 DO 990 I=MYIS,MYIE
951 KCLDS=NCLDS(I)
952 IF(KCLDS.EQ.0) GO TO 990
953 GDFND(I) = GDFND(I) + CCMAX(I)*DFN(I,LP1)*DFNTOP(I,N)
954 990 CONTINUE
955 1000 CONTINUE
956 DO 1100 I=MYIS,MYIE
957 KCLDS=NCLDS(I)
Page 15 Source Listing SWR93
2025-03-12 18:23 SWR93.F
958 IF(KCLDS.EQ.0) GO TO 1100
959 DO 1101 K=1,LP1
960 DFSWC(I,K) = TMP1(I)*DFSWL(I,K) + CCMAX(I)*DFSWC(I,K)
961 UFSWC(I,K) = TMP1(I)*UFSWL(I,K) + CCMAX(I)*UFSWC(I,K)
962 1101 CONTINUE
963 1100 CONTINUE
964 DO 1200 I=MYIS,MYIE
965 KCLDS=NCLDS(I)
966 IF(KCLDS.EQ.0) GO TO 1200
967 DO 1201 KK=1,LP1
968 FSWC(I,KK) = UFSWC(I,KK)-DFSWC(I,KK)
969 1201 CONTINUE
970 1200 CONTINUE
971 DO 1250 I=MYIS,MYIE
972 KCLDS=NCLDS(I)
973 IF(KCLDS.EQ.0) GO TO 1250
974 DO 1251 KK=1, L
975 HSWC(I,KK) = RADCON*(FSWC(I,KK+1)-FSWC(I,KK))/DP(I,KK)
976 1251 CONTINUE
977 1250 CONTINUE
978 RETURN
979 END
ENTRY POINTS
Name
swr93_
Page 16 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
SYMBOL CROSS REFERENCE
Name Object Declared Type Bytes Dimen Elements Attributes References
100 Label 283 278
1000 Label 955 751
110 Label 287 284,285
1100 Label 963 956,958
1101 Label 962 959
120 Label 292 288,289
1200 Label 970 964,966
1201 Label 969 967
1250 Label 977 971,973
1251 Label 976 974
130 Label 296 293,294
140 Label 301 298,299
150 Label 308 304
155 Label 313 309
160 Label 322 317,318
200 Label 334 330
210 Label 340 335,336
220 Label 345 341
230 Label 351 346,347
240 Label 363 359,360
250 Label 375 369,370
260 Label 385 379,380
270 Label 389 386
280 Label 397 393,394
285 Label 404 398
290 Label 418 414,415
300 Label 429 423,424
310 Label 436 433
320 Label 443 439,440
330 Label 446 444
350 Label 447 409
360 Label 451 448,449
370 Label 455 452,453
400 Label 464 462
410 Label 470 465,466
420 Label 474 471,472
430 Label 481 478,479
450 Label 489 487
460 Label 494 492
470 Label 496 482,485
480 Label 502 497,498
490 Label 509 503,505
500 Label 514 510
510 Label 519 515
520 Label 525 520,521
530 Label 530 526
540 Label 536 531,532
550 Label 541 537,538
560 Label 548 542,543
570 Label 555 553
580 Label 566 563
Page 17 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
581 Label 567 560,562
590 Label 580 575
591 Label 581 572,574
600 Label 588 582,583
610 Label 592 589
620 Label 609 603
630 Label 617 613
631 Label 618 610,612
640 Label 626 623
641 Label 627 620,622
650 Label 639 636
651 Label 640 633,635
660 Label 651 646,648
670 Label 660 656
671 Label 661 653,655
680 Label 674 669,671
690 Label 684 680
691 Label 685 677,679
700 Label 692 689
701 Label 693 686,688
710 Label 702 699
720 Label 703 695,697
730 Label 715 712
740 Label 726 723
755 Label 728 708,711
760 Label 729 705,707
770 Label 737 730,732
771 Label 736 733
780 Label 745 738,740
790 Label 759 756
791 Label 760 753,755
800 Label 774 765,767
801 Label 770 768
802 Label 773 771
810 Label 787 776,778
811 Label 782 779
812 Label 786 783
820 Label 796 791,793
830 Label 815 807,809
840 Label 823 819
841 Label 824 816,818
850 Label 831 828
851 Label 832 825,827
860 Label 839 836
861 Label 840 833,835
870 Label 850 845,847
880 Label 859 855
881 Label 860 852,854
890 Label 873 868,870
900 Label 900 877
901 Label 901 874,876
910 Label 909 906
911 Label 910 903,905
920 Label 918 915
930 Label 919 911,913
Page 18 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
940 Label 930 927
950 Label 938 935
965 Label 940 923,926
970 Label 941 920,922
980 Label 949 942,944
981 Label 948 945
990 Label 954 950,952
ALFA Local 253 R(4) 4 2 3266 649,657,658,659,672,681,683,848,85
6,857,858,871,882,899
ALFAU Local 253 R(4) 4 2 3266 650,657,659,681,849,856,858,884
ALNB Dummy 6 R(4) 4 1 71 ARG,INOUT 442
ALND Dummy 6 R(4) 4 1 71 ARG,INOUT
ALVB Dummy 6 R(4) 4 1 71 ARG,INOUT 306,311,401
ALVD Dummy 6 R(4) 4 1 71 ARG,INOUT 307,312,401,512,513
CAMT Dummy 5 R(4) 4 2 3266 ARG,INOUT 480,488,493
CCMAX Local 227 R(4) 4 1 71 484,486,488,490,491,493,741,744,95
3,960,961
COSZRO Dummy 4 R(4) 4 1 71 ARG,INOUT 279,300,305,310,511
CR Local 240 R(4) 4 2 3266 554,564,649,658,659,757,848,857,85
8
CRR Dummy 5 R(4) 4 3 39192 ARG,INOUT 564,757
CT Local 242 R(4) 4 2 3266 565,637,758,837
CTT Dummy 5 R(4) 4 3 39192 ARG,INOUT 565,758
DENOM Local 882 R(4) 4 scalar 882,883,884
DFN Local 239 R(4) 4 2 3266 383,387,388,395,416,425,434,435,44
1,445,701,714,722,725,734,917,929,
934,937,946,953
DFNCLU Local 251 R(4) 4 2 3266 673,683,691,722,725,872,899,908,93
4,937
DFNTOP Local 223 R(4) 4 2 852 300,395,396,441,442,445,734,735,94
6,947,953
DFNTRN Local 251 R(4) 4 2 3266 607,608,624,637,691,701,714,813,81
4,829,837,908,917,929
DFSWC Dummy 2 R(4) 4 2 3266 ARG,INOUT 467,734,744,946,960,968
DFSWL Dummy 3 R(4) 4 2 3266 ARG,INOUT 395,399,400,441,450,467,960
DP Local 217 R(4) 4 2 3266 290,319,320,321,454,975
DPCLD Local 243 R(4) 4 2 3266 578,721,722,933,934
DU Local 219 R(4) 4 2 3266 319,337,348,522,533
DUCO2 Local 219 R(4) 4 2 3266 320,338,349,523,534
DUO3 Local 220 R(4) 4 2 3266 321,339,350,524,535
EXP Func 361 scalar 361,372,381,382,399,425,427,539,54
5,584,585,780,784
FF Local 221 R(4) 4 2 3266 499,506,522
FFCO2 Local 221 R(4) 4 2 3266 500,507,523
FFO3 Local 222 R(4) 4 2 3266 501,508,524
FSWC Dummy 2 R(4) 4 2 3266 ARG,INOUT 469,968,975
FSWL Dummy 2 R(4) 4 2 3266 ARG,INOUT 450,454,469
GDFNB Dummy 6 R(4) 4 1 71 ARG,INOUT 402,445,743
GDFND Dummy 6 R(4) 4 1 71 ARG,INOUT 403,953
GDFVB Dummy 6 R(4) 4 1 71 ARG,INOUT 399,401,742
GDFVD Dummy 6 R(4) 4 1 71 ARG,INOUT 400,744
GLB_TABLE Common 176 336 SAVE
HCON Common 57 872 SAVE
HSWC Dummy 2 R(4) 4 2 3266 ARG,INOUT 473,975
HSWL Dummy 2 R(4) 4 2 3266 ARG,INOUT 454,473
Page 19 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
HTEMP Local 368 R(4) 4 scalar 368,373,546
I Local 278 I(4) 4 scalar 278,279,280,281,282,285,286,289,29
0,291,294,295,304,305,306,307,309,
310,311,312,318,319,320,321,336,33
7,338,339,360,361,370,371,372,373,
374,380,381,382,383,384,386,387,38
8,394,395,396,398,399,400,401,402,
403,415,416,417,424,425,426,427,42
8,433,434,435,440,441,442,444,445,
449,450,453,454,462,463,466,467,46
8,469,472,473,479,480,482,483,484,
486,488,490,491,493,498,499,500,50
1,510,511,512,513,521,522,523,524,
538,539,543,544,545,546,547,553,55
4,560,561,564,565,572,573,576,577,
578,583,584,585,586,587,589,590,59
1,603,604,605,606,607,608,610,611,
614,615,616,620,621,624,625,633,63
4,637,638,646,647,649,650,653,654,
657,658,659,669,670,672,673,677,67
8,681,682,683,686,687,690,691,695,
696,698,700,701,705,706,709,710,71
3,714,719,721,722,724,725,730,731,
734,735,738,739,741,742,743,744,75
3,754,757,758,765,766,769,772,776,
777,780,781,784,785,791,792,794,79
5,807,808,810,811,812,813,814,816,
817,820,821,822,825,826,829,830,83
3,834,837,838,845,846,848,849,852,
853,856,857,858,868,869,871,872,87
4,875,882,884,897,899,903,904,907,
908,911,912,914,916,917,920,921,92
4,925,928,929,931,933,934,936,937,
942,943,946,947,950,951,953,956,95
7,960,961,964,965,968,971,972,975
IDIM1 Param 46 I(4) 4 scalar 199,200,201,202,203,204,205,206,20
7,208,209,210,211,212,213,214,217,
218,219,220,221,222,223,224,225,22
6,227,229,230,231,232,233,234,235,
236,237,239,240,241,242,243,246,24
7,248,249,250,251,252,253,255,256,
257,258
IDIM2 Param 46 I(4) 4 scalar 199,200,201,202,203,204,205,206,20
7,208,209,210,211,212,213,214,217,
218,219,220,221,222,223,224,225,22
6,227,229,230,231,232,233,234,235,
236,237,239,240,241,242,243,246,24
7,248,249,250,251,252,253
IGSTL Param 41 I(4) 4 scalar 46
IGSTR Param 41 I(4) 4 scalar 46
IM Param 29 I(4) 4 scalar 43,46,113,182,183,184,185,186,190,
191,196
IMAX Param 113 I(4) 4 scalar 126
INLTE Param 125 I(4) 4 scalar 125
INLTEP Param 125 I(4) 4 scalar
Page 20 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
INPES Param 38 I(4) 4 scalar 43,46,166,167,168,177,178
IP Local 299 I(4) 4 scalar 299,300,330,331,332,333,341,342,34
3,344,347,348,349,350,503,504,506,
507,508,515,516,517,518,526,527,52
8,529,532,533,534,535
ITAIL Param 43 I(4) 4 scalar
J1 Local 709 I(4) 4 scalar 709,711,712,720,723,924,926,927,93
2,935
J2 Local 698 I(4) 4 scalar 698,699,710,712,914,915,925,927
J3 Local 719 I(4) 4 scalar 719,720,721,722,723,931,932,933,93
4,935
JDIM1 Param 47 I(4) 4 scalar
JDIM2 Param 47 I(4) 4 scalar
JGSTL Param 42 I(4) 4 scalar 47
JGSTR Param 42 I(4) 4 scalar 47
JM Param 29 I(4) 4 scalar 44,47,182,183,184,185,186,190,191,
196
JNPES Param 38 I(4) 4 scalar 44,47,166,167,168,177,178
JTAIL Param 44 I(4) 4 scalar
JTOP Local 504 I(4) 4 scalar 504,505
K Local 284 I(4) 4 scalar 284,286,288,290,291,293,295,317,31
9,320,321,335,337,338,339,346,348,
349,350,359,361,369,371,372,373,37
4,379,381,382,383,384,393,395,396,
414,416,417,423,425,426,427,428,43
9,441,442,448,450,452,454,465,467,
468,469,471,473,478,480,487,488,49
2,493,497,499,500,501,505,506,507,
508,520,522,523,524,531,533,534,53
5,537,539,542,544,545,546,547,582,
584,585,586,587,699,700,701,712,71
3,714,723,724,725,733,734,735,756,
757,758,915,916,917,927,928,929,93
5,936,937,945,946,947,959,960,961
KBTMSW Dummy 5 I(4) 4 2 3266 ARG,INOUT 576,578,616,698,710,719,822,914,92
5,931
KCLDS Local 461 I(4) 4 scalar 461,463,476,561,562,563,573,574,57
5,611,612,613,621,622,623,634,635,
636,647,648,654,655,656,670,671,67
2,673,678,679,680,687,688,689,696,
697,706,707,708,731,732,739,740,75
4,755,756,766,767,777,778,792,793,
808,809,817,818,819,826,827,828,83
4,835,836,846,847,853,854,855,869,
870,871,872,875,876,877,904,905,90
6,912,913,921,922,923,943,944,951,
952,957,958,965,966,972,973
KK Local 563 I(4) 4 scalar 563,564,565,575,576,577,578,613,61
4,615,616,623,624,625,636,637,638,
656,657,658,659,680,681,682,683,68
9,690,691,708,709,710,713,714,719,
721,722,724,725,768,769,771,772,77
9,780,781,783,784,785,819,820,821,
822,828,829,830,836,837,838,855,85
6,857,858,877,882,884,897,899,906,
Page 21 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
907,908,923,924,925,928,929,931,93
3,934,936,937,967,968,974,975
KO2 Param 128 I(4) 4 scalar 129
KO21 Param 129 I(4) 4 scalar
KO2M Param 129 I(4) 4 scalar
KTOPSW Dummy 5 I(4) 4 2 3266 ARG,INOUT 504,576,577,614,615,709,820,821,92
4
L Param 112 I(4) 4 scalar 116,117,118,121,123,205,206,288,29
3,317,346,379,414,423,452,471,531,
582,771,783,974
LL Param 118 I(4) 4 scalar 118,119,359,369,537,542
LL3P Param 123 I(4) 4 scalar 126
LL3PI Param 126 I(4) 4 scalar
LLM1 Param 119 I(4) 4 scalar
LLM2 Param 119 I(4) 4 scalar
LLM3 Param 119 I(4) 4 scalar
LLP1 Param 118 I(4) 4 scalar 126
LLP1I Param 126 I(4) 4 scalar
LLP2 Param 118 I(4) 4 scalar 230,232,234,236
LLP3 Param 118 I(4) 4 scalar
LM Param 29 I(4) 4 scalar 112,191
LM1 Param 117 I(4) 4 scalar 284
LM2 Param 117 I(4) 4 scalar
LM3 Param 117 I(4) 4 scalar
LOG Func 361 scalar 361,372,539,545
LP1 Param 116 I(4) 4 scalar 120,121,122,126,199,200,201,202,20
3,204,205,207,208,217,218,219,220,
221,222,227,229,230,231,232,233,23
4,235,236,237,239,240,241,242,243,
246,247,248,249,250,251,252,253,28
1,282,335,342,343,344,388,393,399,
400,435,439,445,448,465,478,497,52
0,527,528,529,591,604,605,699,733,
744,768,779,794,810,811,915,945,95
3,959,967
LP121 Param 122 I(4) 4 scalar
LP1I Param 126 I(4) 4 scalar
LP1M Param 120 I(4) 4 scalar 120
LP1M1 Param 120 I(4) 4 scalar
LP1V Param 121 I(4) 4 scalar
LP2 Param 116 I(4) 4 scalar 255,256,257,258
LP3 Param 116 I(4) 4 scalar
LSM Param 29 I(4) 4 scalar
MAPPINGS Common 195 5024 SAVE
MAX Func 463 scalar 463
MIN Func 381 scalar 381,382,425,427,584,585,780,784
MPPCOM Common 140 1724 SAVE
N Local 298 I(4) 4 scalar 298,300,409,410,425,427,441,442,44
5,751,757,758,762,780,784,946,947,
953
NB Param 124 I(4) 4 scalar 127,200,202,223,264,298,409,751
NB1 Param 127 I(4) 4 scalar
NBLM Param 115 I(4) 4 scalar
NBLW Param 114 I(4) 4 scalar
NBLX Param 114 I(4) 4 scalar
Page 22 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
NBLY Param 114 I(4) 4 scalar 115,122
NCLDS Dummy 5 I(4) 4 1 71 ARG,INOUT 461,463,483,504,561,573,611,621,63
4,647,654,670,678,687,696,706,731,
739,754,766,777,792,808,817,826,83
4,846,853,869,875,904,912,921,943,
951,957,965,972
NNCLDS Local 483 I(4) 4 scalar 483,485,487,492
NNLTE Param 125 I(4) 4 scalar
PHYCON Common 51 80 SAVE
PP Local 217 R(4) 4 2 3266 280,281,286,290,291,577,578,724,72
5,936,937
PPTOP Local 243 R(4) 4 2 3266 577,578,724,725,936,937
PR2 Local 218 R(4) 4 2 3266 291,295,319,320
PRESS Dummy 4 R(4) 4 2 3266 ARG,INOUT 281,282,286
QO3 Dummy 4 R(4) 4 2 3195 ARG,INOUT 321
REFL Local 225 R(4) 4 1 71 306,396,512,554
REFL2 Local 226 R(4) 4 1 71 311,400
RH2O Dummy 4 R(4) 4 2 3195 ARG,INOUT 319
RRAY Local 223 R(4) 4 1 71 305,306,310,311,511,512
RRCO2 Dummy 4 R(4) 4 scalar ARG,INOUT 320
SECZ Local 224 R(4) 4 1 71 279,337,338,339,399,506,507,508
SQRT Func 279 scalar 279
SSOLAR Dummy 4 R(4) 4 scalar ARG,INOUT 300
SWR93 Subr 2
SWRSAV Common 264 112 SAVE
TAUDAR Dummy 4 R(4) 4 1 71 ARG,INOUT 300
TCLD Local 252 R(4) 4 2 3266 638,658,838,857
TCLU Local 250 R(4) 4 2 3266 637,657,682,837,856,882
TCO2 Local 234 R(4) 4 2 6532 361,539
TDCL1 Local 247 R(4) 4 2 3266 604,606,607,614,624,637,638,672,67
3,810,812,813,820,829,837,838,871,
872
TDCL2 Local 248 R(4) 4 2 3266 606,616,638,812,822,838
TDCO2 Local 234 R(4) 4 2 3266 416,426,769,781
TDO3 Local 236 R(4) 4 2 3266 383,586
TEMPCOM Common 181 6603768 SAVE
TEMPF Local 721 R(4) 4 scalar 721,724,933,936
TEMPG Local 722 R(4) 4 scalar 722,725,934,937
TMP1 Local 225 R(4) 4 1 71 282,295,741,742,743,744,960,961
TO3 Local 236 R(4) 4 2 6532 371,544
TOPO Common 189 17478548 SAVE
TTD Local 241 R(4) 4 2 3266 381,383,416,586,590,591,604,614,61
6,701,714,769,780,794,795,810,820,
822,917,929
TTDB1 Local 247 R(4) 4 2 3266 584,586,769
TTU Local 241 R(4) 4 2 3266 382,384,417,587,591,605,615,700,71
3,772,784,794,811,821,916,928
TTUB1 Local 246 R(4) 4 2 3266 585,587,772
TUCL1 Local 246 R(4) 4 2 3266 605,615,625,811,821,830
TUCO2 Local 235 R(4) 4 2 3266 417,428,772,785
TUO3 Local 237 R(4) 4 2 3266 384,587
UCO2 Local 230 R(4) 4 2 6532 361,539
UD Local 229 R(4) 4 2 3266 331,337,342,381,425,516,522,527,58
4,780
UDCO2 Local 230 R(4) 4 2 3266 332,338,343,517,523,528
Page 23 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Object Declared Type Bytes Dimen Elements Attributes References
UDO3 Local 232 R(4) 4 2 3266 333,339,344,518,524,529
UFN Local 239 R(4) 4 2 3266 384,388,396,417,427,435,442,700,71
3,721,724,735,916,928,933,936,947
UFNCLU Local 249 R(4) 4 2 3266 672,681,683,690,721,724,871,884,89
7,899,907,933,936
UFNTRN Local 249 R(4) 4 2 3266 608,625,690,700,713,814,830,907,91
6,928
UFSWC Dummy 2 R(4) 4 2 3266 ARG,INOUT 468,735,947,961,968
UFSWL Dummy 2 R(4) 4 2 3266 ARG,INOUT 396,442,450,468,961
UO3 Local 232 R(4) 4 2 6532 371,372,373,374,544,545,546,547
UR Local 229 R(4) 4 2 3266 342,348,382,427,527,533,585,784
URCO2 Local 231 R(4) 4 2 3266 343,349,528,534
URO3 Local 233 R(4) 4 2 3266 344,350,529,535
XAMT Local 227 R(4) 4 2 3266 480,493,564,565,757,758
TYPE COMPONENTS/COMMON VARIABLES
Name Type Bytes Offset Dimen Elements Attributes References
ABCFF R(4) 4 0 1 12 COM 381,382,425,427,584,585,780,784
AMOLWT R(4) 4 0 scalar COM
CFCO2 R(4) 4 96 scalar COM 320
CFO3 R(4) 4 100 scalar COM 321
CSUBP R(4) 4 4 scalar COM
DIFFCTR R(4) 4 8 scalar COM 348,349,499,500,533,534
EIGHT R(4) 4 20 scalar COM
FIFTY R(4) 4 12 scalar COM 381,382,425,427,584,585,780,784
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 319,320,321
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 368
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 279
H1226E1 R(4) 4 652 scalar COM
H128M5 R(4) 4 240 scalar COM
H129M2 R(4) 4 732 scalar COM 361,539
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 372,545
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:23 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 372,545
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 374,547
H21M12 R(4) 4 460 scalar COM
H21M31 R(4) 4 556 scalar COM
H235M3 R(4) 4 724 scalar COM 361,539
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:23 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 374,547
H327M8 R(4) 4 276 scalar COM
H32M11 R(4) 4 348 scalar COM
H35E1 R(4) 4 140 scalar COM 279
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 374,547
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:23 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 373,546
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 300
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 362,540
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 286,291
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 381,382,425,427,584,585,780,784
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:23 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 372,545
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 305,511
HP228 R(4) 4 672 scalar COM
HP26 R(4) 4 728 scalar COM 361,539
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 305,511
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 828 1 21 COM
IE_GLB_TABLE I(4) 4 84 1 21 COM
IE_LOC_TABLE I(4) 4 660 1 21 COM
ILCOL I(4) 4 396 scalar COM
ILPAD1 I(4) 4 412 scalar COM
ILPAD2 I(4) 4 416 scalar COM
ILPAD3 I(4) 4 420 scalar COM
ILPAD4 I(4) 4 424 scalar COM
ILPAD5 I(4) 4 428 scalar COM
INUMQ I(4) 4 1324 1 100 COM
IQUILT_GROUP I(4) 4 1320 scalar COM
IRCOL I(4) 4 400 scalar COM
IRPAD1 I(4) 4 432 scalar COM
IRPAD2 I(4) 4 436 scalar COM
IRPAD3 I(4) 4 440 scalar COM
IRPAD4 I(4) 4 444 scalar COM
IRPAD5 I(4) 4 448 scalar COM
IS_GLB_TABLE I(4) 4 0 1 21 COM
IS_LOC_TABLE I(4) 4 492 1 21 COM
ITEMP I(4) 4 5870016 2 91719 COM
ITEMP2 I(4) 4 6236892 2 91719 COM
ITROW I(4) 4 408 scalar COM
JBPAD1 I(4) 4 452 scalar COM
JBPAD2 I(4) 4 456 scalar COM
JBPAD3 I(4) 4 460 scalar COM
JBPAD4 I(4) 4 464 scalar COM
JBPAD5 I(4) 4 468 scalar COM
JE_GLB_TABLE I(4) 4 252 1 21 COM
JE_LOC_TABLE I(4) 4 744 1 21 COM
JS_GLB_TABLE I(4) 4 168 1 21 COM
Page 28 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Type Bytes Offset Dimen Elements Attributes References
JS_LOC_TABLE I(4) 4 576 1 21 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 912 scalar COM
MPI_COMM_INTER I(4) 4 916 scalar COM
MPI_COMM_INTER_ARRAY I(4) 4 920 1 100 COM
MYIE I(4) 4 52 scalar COM 278,285,289,294,299,304,309,318,33
0,336,341,347,360,370,380,386,394,
398,415,424,433,440,444,449,453,46
2,466,472,479,482,498,503,510,515,
521,526,532,538,543,553,560,572,58
3,589,603,610,620,633,646,653,669,
677,686,695,705,730,738,753,765,77
6,791,807,816,825,833,845,852,868,
874,903,911,920,942,950,956,964,97
1
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 278,285,289,294,299,304,309,318,33
0,336,341,347,360,370,380,386,394,
398,415,424,433,440,444,449,453,46
1,466,472,479,482,498,503,510,515,
521,526,532,538,543,553,560,572,58
3,589,603,610,620,633,646,653,669,
677,686,695,705,730,738,753,765,77
6,791,807,816,825,833,845,852,868,
874,903,911,920,942,950,956,964,97
1
MYIS1 I(4) 4 44 scalar COM 462
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:23 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:23 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 350,501,535
ONE R(4) 4 40 scalar COM 279,282,305,306,307,310,311,312,36
1,371,372,373,374,387,400,401,434,
486,488,490,511,512,513,539,544,54
5,546,547,565,578,590,607,624,625,
658,741,758,795,813,829,830,856
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 300
QUARTR R(4) 4 48 scalar COM
RADCON R(4) 4 40 scalar COM 454,975
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 306,307,512,513
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 361,371,539,544
Page 31 Source Listing SWR93
2025-03-12 18:23 Symbol Table SWR93.F
Name Type Bytes Offset Dimen Elements Attributes References
ZERO R(4) 4 52 scalar COM 280,331,332,333,402,403,484,491,51
6,517,518,650,849
Page 32 Source Listing SWR93
2025-03-12 18:23 Subprograms/Common Blocks SWR93.F
SUBPROGRAMS/COMMON BLOCKS
Name Object Declared Type Bytes Dimen Elements Attributes References
GLB_TABLE Common 176 336 SAVE
HCON Common 57 872 SAVE
MAPPINGS Common 195 5024 SAVE
MPPCOM Common 140 1724 SAVE
PHYCON Common 51 80 SAVE
SWR93 Subr 2
SWRSAV Common 264 112 SAVE
TEMPCOM Common 181 6603768 SAVE
TOPO Common 189 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:23 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:23 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