Page 1 Source Listing LWR88
2025-03-12 18:21 /tmp/ifortNM9FTa.i
1 # 1 "LWR88.F"
2 SUBROUTINE LWR88(HEATRA,GRNFLX,TOPFLX,
3 1 PRESS,TEMP,RH2O,QO3,CLDFAC,
4 2 CAMT,NCLDS,KTOP,KBTM)
5 C---------------------------------------------------------------------
6 C SUBROUTINE LWR88 COMPUTES TEMPERATURE-CORRECTED CO2 TRANSMISSION
7 C FUNCTIONS AND ALSO COMPUTES THE PRESSURE GRID AND LAYER OPTICAL
8 C PATHS.
9 C INPUTS: (COMMON BLOCKS)
10 C CLDFAC CLDCOM
11 C PRESS,TEMP,RH2O,QO3 RADISW
12 C CAMT,NCLDS,KTOP,KBTM RADISW
13 C CO251,CO258,CDT51,CDT58 CO2BD3
14 C C2D51,C2D58,CO2M51,CO2M58 CO2BD3
15 C CDTM51,CDTM58,C2DM51,C2DM58 CO2BD3
16 C STEMP,GTEMP CO2BD3
17 C CO231,CO238,CDT31,CDT38 CO2BD2
18 C C2D31,C2D38 CO2BD2
19 C CO271,CO278,CDT71,CDT78 CO2BD4
20 C C2D71,C2D78 CO2BD4
21 C BETINW BDWIDE
22 C OUTPUTS:
23 C HEATRA,GRNFLX,TOPFLX LWOUT
24 C CALLED BY:
25 C RADMN OR INPUT ROUTINE OF MODEL
26 C CALLS:
27 C FST88
28 C
29 COMMON/PHYCON/AMOLWT,CSUBP,DIFFCTR,G,GRAVDR,O3DIFCTR,P0,
30 * P0XZP2,P0XZP8,P0X2,RADCON,RGAS,RGASSP,SECPDA
31 COMMON/PHYCON/RATCO2MW,RATH2OMW
32 COMMON/PHYCON/RADCON1
33 COMMON/PHYCON/GINV,P0INV,GP0INV
34 COMMON/HCON/HUNDRED,HNINETY,SIXTY,FIFTY,TEN,EIGHT,FIVE,
35 * FOUR,THREE,TWO,ONE,HAF,QUARTR,ZERO
36 COMMON/HCON/H83E26,H71E26,H1E15,H1E13,H1E11,H1E8,H4E5,
37 * H165E5,H5725E4,H488E4,H1E4,H24E3,H20788E3,
38 * H2075E3,H1224E3,H5E2,H3082E2,H3E2,H2945E2,
39 * H23E2,H15E2,H35E1,H3P6,H181E1,H18E1,H2P9,H2P8,
40 * H2P5,H1P8,H1P4387,H1P4,H1P25892,HP8,HP518,
41 * HP369,HP1
42 COMMON/HCON/H44871M2,H559M3,H1M3,H987M4,H285M4,H1M4,
43 * H6938M5,H394M5,H37412M5,H1439M5,H128M5,H1M5,
44 * H7M6,H4999M6,H25452M6,H1M6,H391M7,H1174M7,
45 * H8725M8,H327M8,H257M8,H1M8,H23M10,H14M10,
46 * H11M10,H1M10,H83M11,H82M11,H8M11,H77M11,
47 * H72M11,H53M11,H48M11,H44M11,H42M11,H37M11,
48 * H35M11,H32M11,H3M11,H28M11,H24M11,H23M11,
49 * H2M11,H18M11,H15M11,H14M11,H114M11,H11M11,
50 * H1M11,H96M12,H93M12,H77M12,H74M12,H65M12,
51 * H62M12,H6M12,H45M12,H44M12,H4M12,H38M12,
52 * H37M12,H3M12,H29M12,H28M12,H24M12,H21M12,
53 * H16M12,H14M12,H12M12,H8M13,H46M13,H36M13,
54 * H135M13,H12M13,H1M13,H3M14,H15M14,H14M14,
55 * H1M17,H1M18,H1M19,H1M20,H1M21,H1M22,H1M23,
56 * H1M24,H26M30,H14M30,H25M31,H21M31,H12M31,
57 * H9M32,H55M32,H45M32,H4M33,H62M34,H1M60
Page 2 Source Listing LWR88
2025-03-12 18:21 LWR88.F
58 COMMON/HCON/HMP575,HM13EZ,HM19EZ,HM1E1,HM181E1,HM1E2
59 COMMON/HCON/H1E6,H2E6,H1M2,HMP66667,HM6666M2,HP166666,
60 * H41666M2,HMP5,HM2M2,H29316E2,H1226E1,H3116E1,
61 * H9P94,HP6,H625M2,HP228,HP60241,HM1797E1,
62 * H8121E1,H2E2,HM1EZ,H26E2,H44194M2,H1P41819
63 COMMON/HCON/HP219,HP144,HP816,H69766E5,H235M3,HP26,
64 * H129M2,H75826M4,H1P082,HP805,H1386E2,
65 * H658M2,H1036E2,H2118M2,H42M2,H323M4,
66 * H67390E2,HP3795,HP5048,H102M5,H451M6
67 COMMON/HCON/H16E1,HM161E1,H161E1,H3M3,H101M16,
68 * HM1597E1,H25E2,HP118666,H15M5,H3P5,H18E3,
69 * H6P08108,HMP805,HP602409,HP526315,
70 * H28571M2,H1M16
71 COMMON/HCON/H3M4
72 COMMON/HCON/HM8E1
73 COMMON/HCON/H28E1
74 C-----------------------------------------------------------------------
75 INCLUDE "parmeta"
98 INCLUDE "mpp.h"
99 # 1 "./sp.h" 1
100 # 4
101
102 # 77 "LWR88.F" 2
103 C-----------------------------------------------------------------------
104 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE:
105 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL
106 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE
107 C BANDTA FOR DEFINITION
108 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS
109 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE
110 C BDCOMB FOR DEFINITION
111 C INLTE = NO. LEVELS USED FOR NLTE CALCS.
112 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS.
113 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED
114 C FROM THE ABOVE PARAMETERS.
173 PARAMETER (L=LM)
174 PARAMETER (NBLW=163,NBLX=47,NBLY=15)
175 PARAMETER (NBLM=NBLY-1)
176 PARAMETER (LP1=L+1,LP2=L+2,LP3=L+3)
177 PARAMETER (LM1=L-1,LM2=L-2,LM3=L-3)
178 PARAMETER (LL=2*L,LLP1=LL+1,LLP2=LL+2,LLP3=LL+3)
179 PARAMETER (LLM1=LL-1,LLM2=LL-2,LLM3=LL-3)
180 PARAMETER (LP1M=LP1*LP1,LP1M1=LP1M-1)
181 PARAMETER (LP1V=LP1*(1+2*L/2))
182 PARAMETER (LP121=LP1*NBLY)
183 PARAMETER (LL3P=3*L+2)
184 PARAMETER (NB=12)
185 PARAMETER (INLTE=3,INLTEP=INLTE+1,NNLTE=56)
186 PARAMETER (NB1=NB-1)
187 PARAMETER (KO2=12)
188 PARAMETER (KO21=KO2+1,KO2M=KO2-1)
189 C PARAMETER SETTINGS FOR THE LONGWAVE AND SHORTWAVE RADIATION CODE:
190 C L = NO. VERTICAL LEVELS (ALSO LAYERS) IN MODEL
191 C NBLW = NO. FREQ. BANDS FOR APPROX COMPUTATIONS. SEE
192 C BANDTA FOR DEFINITION
193 C NBLX = NO. FREQ BANDS FOR APPROX CTS COMPUTATIONS
194 C NBLY = NO. FREQ. BANDS FOR EXACT CTS COMPUTATIONS. SEE
Page 3 Source Listing LWR88
2025-03-12 18:21 LWR88.F
195 C BDCOMB FOR DEFINITION
196 C INLTE = NO. LEVELS USED FOR NLTE CALCS.
197 C NNLTE = INDEX NO. OF FREQ. BAND IN NLTE CALCS.
198 C NB,KO2 ARE SHORTWAVE PARAMETERS; OTHER QUANTITIES ARE DERIVED
199 C FROM THE ABOVE PARAMETERS.
200 C THE FOLLOWING COMMON BLOCKS CONTAIN PRETABULATED CO2 TRANSMISSION
201 C FUNCTIONS, EVALUATED USING THE METHODS OF FELS AND
202 C SCHWARZKOPF (1981) AND SCHWARZKOPF AND FELS (1985),
203 C***COMMON CO2BD3 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE
204 C AND PRESSURE DERIVATIVES FOR THE 560-800 CM-1 BAND. ALSO INCLUDED
205 C ARE THE STANDARD TEMPERATURES AND THE WEIGHTING FUNCTION. THESE
206 C DATA ARE IN BLOCK DATA BD3:
207 C CO251 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE)
208 C WITH P(SFC)=1013.25 MB
209 C CO258 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE)
210 C WITH P(SFC)= ^810 MB
211 C CDT51 = FIRST TEMPERATURE DERIVATIVE OF CO251
212 C CDT58 = FIRST TEMPERATURE DERIVATIVE OF CO258
213 C C2D51 = SECOND TEMPERATURE DERIVATIVE OF CO251
214 C C2D58 = SECOND TEMPERATURE DERIVATIVE OF CO251
215 C CO2M51 = TRANSMISSION FCTNS FOR T0 FOR ADJACENT PRESSURE
216 C LEVELS, WITH NO PRESSURE QUADRATURE. USED FOR
217 C NEARBY LAYER COMPUTATIONS. P(SFC)=1013.25 MB
218 C CO2M58 = SAME AS CO2M51,WITH P(SFC)= ^810 MB
219 C CDTM51 = FIRST TEMPERATURE DERIVATIVE OF CO2M51
220 C CDTM58 = FIRST TEMPERATURE DERIVATIVE OF CO2M58
221 C C2DM51 = SECOND TEMPERATURE DERIVATIVE OF CO2M51
222 C C2DM58 = SECOND TEMPERATURE DERIVATIVE OF CO2M58
223 C STEMP = STANDARD TEMPERATURES FOR MODEL PRESSURE LEVEL
224 C STRUCTURE WITH P(SFC)=1013.25 MB
225 C GTEMP = WEIGHTING FUNCTION FOR MODEL PRESSURE LEVEL
226 C STRUCTURE WITH P(SFC)=1013.25 MB.
227 C B0 = TEMP. COEFFICIENT USED FOR CO2 TRANS. FCTN.
228 C CORRECTION FOR T(K). (SEE REF. 4 AND BD3)
229 C B1 = TEMP. COEFFICIENT, USED ALONG WITH B0
230 C B2 = TEMP. COEFFICIENT, USED ALONG WITH B0
231 C B3 = TEMP. COEFFICIENT, USED ALONG WITH B0
232 C
233 COMMON/CO2BD3/CO251(LP1,LP1),CO258(LP1,LP1),CDT51(LP1,LP1),
234 1 CDT58(LP1,LP1),C2D51(LP1,LP1),C2D58(LP1,LP1),CO2M51(L),
235 2 CO2M58(L),CDTM51(L),CDTM58(L),C2DM51(L),C2DM58(L),
236 3 STEMP(LP1),GTEMP(LP1),B0,B1,B2,B3
237 C
238 C***COMMON CO2BD2 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE
239 C AND PRESSURE DERIVATIVES FOR THE 560-670 CM-1 PART OF THE 15 UM
240 C CO2 BAND. THESE DATA ARE IN BLOCK DATA BD2.
241 C CO231 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE)
242 C WITH P(SFC)=1013.25 MB
243 C CO238 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE)
244 C WITH P(SFC)= ^810 MB
245 C CDT31 = FIRST TEMPERATURE DERIVATIVE OF CO231
246 C CDT38 = FIRST TEMPERATURE DERIVATIVE OF CO238
247 C C2D31 = SECOND TEMPERATURE DERIVATIVE OF CO231
248 C C2D38 = SECOND TEMPERATURE DERIVATIVE OF CO231
249 C
250 COMMON / CO2BD2 / CO231(LP1),CO238(LP1),CDT31(LP1),
251 1 CDT38(LP1),C2D31(LP1),C2D38(LP1)
Page 4 Source Listing LWR88
2025-03-12 18:21 LWR88.F
252 C
253 C***COMMON CO2BD4 CONTAINS CO2 TRANSMISSION FUNCTIONS AND TEMPERATURE
254 C AND PRESSURE DERIVATIVES FOR THE 670-800 CM-1 PART OF THE 15 UM
255 C CO2 BAND. THESE DATA ARE IN BLOCK DATA BD4.
256 C CO271 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE)
257 C WITH P(SFC)=1013.25 MB
258 C CO278 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE)
259 C WITH P(SFC)= ^810 MB
260 C CDT71 = FIRST TEMPERATURE DERIVATIVE OF CO271
261 C CDT78 = FIRST TEMPERATURE DERIVATIVE OF CO278
262 C C2D71 = SECOND TEMPERATURE DERIVATIVE OF CO271
263 C C2D78 = SECOND TEMPERATURE DERIVATIVE OF CO271
264 C
265 COMMON / CO2BD4 / CO271(LP1),CO278(LP1),CDT71(LP1),
266 1 CDT78(LP1),C2D71(LP1),C2D78(LP1)
267 C
268 C***COMMON CO2BD5 CONTAINS CO2 TRANSMISSION FUNCTIONS FOR THE 2270-
269 C 2380 PART OF THE 4.3 UM CO2 BAND. THESE DATA ARE IN BLOCK DATA BD5.
270 C CO211 = TRANSMISSION FCTNS FOR T0 (STD. PROFILE)
271 C WITH P(SFC)=1013.25 MB
272 C CO218 = TRANSMISSION FCTNS. FOR T0 (STD. PROFILE)
273 C WITH P(SFC)= ^810 MB
274 C
275 COMMON / CO2BD5 / CO211(LP1),CO218(LP1)
276 C
277 C COMMON BLOCK BANDTA CONTAINS RANDOM BAND PARAMETERS FOR THE LW
278 C CALCULATIONS USING 10 CM-1 WIDE BANDS.THE 15 UM CO2 COMPLEX
279 C IS 2 BANDS,560-670 AND 670-800 CM-1. OZONE COEFFICIENTS ARE
280 C IN 3 BANDS,670-800 (14.1 UM),990-1070 AND 1070-1200 (9.6 UM).
281 C THE (NBLW) BANDS NOW INCLUDE:
282 C 56 BANDS, 10 CM-1 WIDE 0 - 560 CM-1
283 C 2 BANDS, 15 UM COMPLEX 560 - 670 CM-1
284 C 670 - 800 CM-1
285 C 3 "CONTINUUM" BANDS 800 - 900 CM-1
286 C 900 - 990 CM-1
287 C 1070 - 1200 CM-1
288 C 1 BAND FOR 9.6 UM BAND 990 - 1070 CM-1
289 C 100 BANDS, 10 CM-1 WIDE 1200 - 2200 CM-1
290 C 1 BAND FOR 4.3 UM SRC 2270 - 2380 CM-1
291 C THUS NBLW PRESENTLY EQUALS 163
292 C ALL BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER
293 C
294 C ARNDM = RANDOM "A" PARAMETER FOR (NBLW) BANDS
295 C BRNDM = RANDOM "B" PARAMETER FOR (NBLW) BANDS
296 C BETAD = CONTINUUM COEFFICIENTS FOR (NBLW) BANDS
297 C AP,BP = CAPPHI COEFFICIENTS FOR (NBLW) BANDS
298 C ATP,BTP = CAPPSI COEFFICIENTS FOR (NBLW) BANDS
299 C BANDLO = LOWEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS
300 C BANDHI = HIGHEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS
301 C AO3RND = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE
302 C BANDS
303 C BO3RND = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE
304 C BANDS
305 C AB15 = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS
306 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2
307 C DATA FOR ARNDM,BRNDM,AP,BP,ATP,BTP,AO3RND,BO3RND ARE OBTAINED BY
308 C USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS ARE FROM
Page 5 Source Listing LWR88
2025-03-12 18:21 LWR88.F
309 C ROBERTS (1976).
310 COMMON / BANDTA / ARNDM(NBLW),BRNDM(NBLW),BETAD(NBLW),AP(NBLW),
311 1 BP(NBLW),ATP(NBLW),BTP(NBLW),BANDLO(NBLW),
312 2 BANDHI(NBLW),AO3RND(3),BO3RND(3),AB15(2)
313 C
314 C COMMON BLOCK BDWIDE CONTAINS RANDOM BAND PARAMETERS FOR SPECIFIC
315 C WIDE BANDS. AT PRESENT,THE INFORMATION CONSISTS OF 1) RANDOM
316 C MODEL PARAMETERS FOR THE 15 UM BAND,560-800 CM-1; 2) THE
317 C CONTINUUM COEFFICIENT FOR THE 800-990,1070-1200 CM-1 BAND
318 C SPECIFICALLY:
319 C AWIDE = RANDOM "A" PARAMETER FOR BAND
320 C BWIDE = RANDOM "B" PARAMETER FOR BAND
321 C BETAWD = CONTINUUM COEFFICIENTS FOR BAND
322 C APWD,BPWD = CAPPHI COEFFICIENTS FOR BAND
323 C ATPWD,BTPWD = CAPPSI COEFFICIENTS FOR BAND
324 C BDLOWD = LOWEST FREQUENCY IN EACH FREQ BAND
325 C BDHIWD = HIGHEST FREQUENCY IN EACH FREQ BAND
326 C AB15WD = THE PRODUCT ARNDM*BRNDM FOR THE ONE BAND
327 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2
328 C BETINW = CONT.COEFFICIENT FOR A SPECIFIED WIDE
329 C FREQ.BAND (800-990 AND 1070-1200 CM-1).
330 C SKO2D = 1./BETINW, USED IN SPA88 FOR CONT. COEFFS
331 C SKC1R = BETAWD/BETINW, USED FOR CONT. COEFF. FOR
332 C 15 UM BAND IN FST88
333 C SKO3R = RATIO OF CONT. COEFF. FOR 9.9 UM BAND TO
334 C BETINW, USED FOR 9.6 UM CONT COEFF IN FST88
335 C DATA FOR AWIDE,BWIDE,APWD,BPWD,ATPWD,BTPWD,AO3WD,BO3WD ARE
336 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS
337 C ARE FROM ROBERTS (1976).
338 COMMON / BDWIDE / AWIDE,BWIDE,BETAWD,
339 1 APWD,BPWD,ATPWD,BTPWD,
340 2 BDLOWD,BDHIWD,BETINW,
341 3 AB15WD,SKO2D,SKC1R,SKO3R
342 C
343 C COMMON BLOCK BDCOMB CONTAINS RANDOM BAND PARAMETERS FOR THE LW
344 C CALCULATIONS USING COMBINED WIDE FREQUENCY BANDS BETWEEN 160 AND
345 C 1200 CM-1,AS WELL AS THE 2270-2380 BAND FOR SOURCE CALC.
346 C BANDS 1-8: COMBINED WIDE FREQUENCY BANDS FOR 160-560 CM-1
347 C BANDS 9-14: FREQUENCY BANDS,AS IN BANDTA (NARROW BANDS)
348 C FOR 560-1200 CM-1
349 C BAND 15: FREQUENCY BAND 2270-2380 CM-1,USED FOR SOURCE
350 C CALCULATION ONLY
351 C THUS NBLY PRESENTLY EQUALS 15
352 C
353 C BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER
354 C ACOMB = RANDOM "A" PARAMETER FOR (NBLY) BANDS
355 C BCOMB = RANDOM "B" PARAMETER FOR (NBLY) BANDS
356 C BETACM = CONTINUUM COEFFICIENTS FOR (NBLY) BANDS
357 C APCM,BPCM = CAPPHI COEFFICIENTS FOR (NBLY) BANDS
358 C ATPCM,BTPCM = CAPPSI COEFFICIENTS FOR (NBLY) BANDS
359 C BDLOCM = LOWEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS
360 C BDHICM = HIGHEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS
361 C AO3CM = RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE
362 C BANDS
363 C BO3CM = RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE
364 C BANDS
365 C AB15CM = THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS
Page 6 Source Listing LWR88
2025-03-12 18:21 LWR88.F
366 C REPRESENTING THE 15 UM BAND COMPLEX OF CO2
367 C BETINC = CONT.COEFFICIENT FOR A SPECIFIED WIDE
368 C FREQ.BAND (800-990 AND 1070-1200 CM-1).
369 C IBAND = INDEX NO OF THE 40 WIDE BANDS USED IN
370 C COMBINED WIDE BAND CALCULATIONS. IN OTHER
371 C WORDS,INDEX TELLING WHICH OF THE 40 WIDE
372 C BANDS BETWEEN 160-560 CM-1 ARE INCLUDED IN
373 C EACH OF THE FIRST 8 COMBINED WIDE BANDS
374 C DATA FOR ACOMB,BCOMB,APCM,BPCM,ATPCM,BTPCM,AO3CM,BO3CM ARE
375 C OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS
376 C ARE FROM ROBERTS (1976). IBAND INDEX VALUES ARE OBTAINED BY
377 C EXPERIMENTATION.
378 COMMON / BDCOMB / IBAND(40),ACOMB(NBLY),BCOMB(NBLY),
379 1 BETACM(NBLY),APCM(NBLY),BPCM(NBLY),ATPCM(NBLY),
380 2 BTPCM(NBLY),BDLOCM(NBLY),BDHICM(NBLY),BETINC,
381 3 AO3CM(3),BO3CM(3),AB15CM(2)
382 C
383 C
384 DIMENSION PRESS(IDIM1:IDIM2,LP1),TEMP(IDIM1:IDIM2,LP1),
385 1 RH2O(IDIM1:IDIM2,L),QO3(IDIM1:IDIM2,L)
386 DIMENSION CLDFAC(IDIM1:IDIM2,LP1,LP1),CAMT(IDIM1:IDIM2,LP1)
387 DIMENSION NCLDS(IDIM1:IDIM2),KTOP(IDIM1:IDIM2,LP1),
388 1 KBTM(IDIM1:IDIM2,LP1)
389 DIMENSION HEATRA(IDIM1:IDIM2,L),GRNFLX(IDIM1:IDIM2),
390 1 TOPFLX(IDIM1:IDIM2)
391 DIMENSION DELP2(IDIM1:IDIM2,L)
392 C
393 DIMENSION QH2O(IDIM1:IDIM2,L),T(IDIM1:IDIM2,LP1)
394 DIMENSION P(IDIM1:IDIM2,LP1),DELP(IDIM1:IDIM2,L)
395 DIMENSION CO21(IDIM1:IDIM2,LP1,LP1),CO2NBL(IDIM1:IDIM2,L)
396 DIMENSION CO2SP1(IDIM1:IDIM2,LP1),CO2SP2(IDIM1:IDIM2,LP1)
397 DIMENSION VAR1(IDIM1:IDIM2,L),VAR2(IDIM1:IDIM2,L),
398 1 VAR3(IDIM1:IDIM2,L),VAR4(IDIM1:IDIM2,L)
399 DIMENSION CNTVAL(IDIM1:IDIM2,LP1)
400 DIMENSION TOTO3(IDIM1:IDIM2,LP1),TPHIO3(IDIM1:IDIM2,LP1),
401 1 TOTPHI(IDIM1:IDIM2,LP1)
402 DIMENSION TOTVO2(IDIM1:IDIM2,LP1),EMX1(IDIM1:IDIM2),
403 1 EMX2(IDIM1:IDIM2),EMPL(IDIM1:IDIM2,LLP1)
404 C
405 DIMENSION CO2R(IDIM1:IDIM2,LP1),DIFT(IDIM1:IDIM2,LP1)
406 DIMENSION CO2R1(IDIM1:IDIM2,LP1),DCO2D1(IDIM1:IDIM2,LP1)
407 DIMENSION D2CD21(IDIM1:IDIM2,LP1),D2CD22(IDIM1:IDIM2,LP1)
408 DIMENSION CO2R2(IDIM1:IDIM2,LP1),DCO2D2(IDIM1:IDIM2,LP1)
409 DIMENSION CO2MR(IDIM1:IDIM2,L),CO2MD(IDIM1:IDIM2,L),
410 1 CO2M2D(IDIM1:IDIM2,L)
411 DIMENSION TDAV(IDIM1:IDIM2,LP1),TSTDAV(IDIM1:IDIM2,LP1),
412 1 VV(IDIM1:IDIM2,L),VSUM3(IDIM1:IDIM2,LP1),VSUM1(IDIM1:IDIM2),
413 2 VSUM2(IDIM1:IDIM2)
414 DIMENSION A1(IDIM1:IDIM2),A2(IDIM1:IDIM2)
415 DIMENSION DCO2DT(IDIM1:IDIM2,LP1),D2CDT2(IDIM1:IDIM2,LP1)
416 C
417 DIMENSION TEXPSL(IDIM1:IDIM2,LP1),TLSQU(IDIM1:IDIM2,LP1)
418 DIMENSION VSUM4(IDIM1:IDIM2,L)
419 EQUIVALENCE (VSUM3,TLSQU,TEXPSL)
420 EQUIVALENCE (VV,VSUM4)
421 c
422 C
Page 7 Source Listing LWR88
2025-03-12 18:21 LWR88.F
423 C****COMPUTE FLUX PRESSURES (P) AND DIFFERENCES (DELP2,DELP)
424 C****COMPUTE FLUX LEVEL TEMPERATURES (T) AND CONTINUUM TEMPERATURE
425 C CORRECTIONS (TEXPSL)
426 DO 103 K=2,L
427 DO 103 I=MYIS,MYIE
428 P(I,K)=HAF*(PRESS(I,K-1)+PRESS(I,K))
429 T(I,K)=HAF*(TEMP(I,K-1)+TEMP(I,K))
430 103 CONTINUE
431 DO 105 I=MYIS,MYIE
432 P(I,1)=ZERO
433 P(I,LP1)=PRESS(I,LP1)
434 T(I,1)=TEMP(I,1)
435 T(I,LP1)=TEMP(I,LP1)
436 105 CONTINUE
437 DO 107 K=1,L
438 DO 107 I=MYIS,MYIE
439 DELP2(I,K)=P(I,K+1)-P(I,K)
440 DELP(I,K)=ONE/DELP2(I,K)
441 107 CONTINUE
442 C****COMPUTE ARGUMENT FOR CONT.TEMP.COEFF.
443 C (THIS IS 1800.(1./TEMP-1./296.))
444 DO 125 K=1,LP1
445 DO 125 I=MYIS,MYIE
446 TEXPSL(I,K)=H18E3/TEMP(I,K)-H6P08108
447 C...THEN TAKE EXPONENTIAL
448 TEXPSL(I,K)=EXP(TEXPSL(I,K))
449 125 CONTINUE
450 C***COMPUTE OPTICAL PATHS FOR H2O AND O3, USING THE DIFFUSIVITY
451 C APPROXIMATION FOR THE ANGULAR INTEGRATION (1.66). OBTAIN THE
452 C UNWEIGHTED VALUES(VAR1,VAR3) AND THE WEIGHTED VALUES(VAR2,VAR4).
453 C THE QUANTITIES H3M4(.0003) AND H3M3(.003) APPEARING IN THE VAR2 AND
454 C VAR4 EXPRESSIONS ARE THE APPROXIMATE VOIGT CORRECTIONS FOR H2O AND
455 C O3,RESPECTIVELY.
456 C
457 DO 131 K=1,L
458 DO 131 I=MYIS,MYIE
459 QH2O(I,K)=RH2O(I,K)*DIFFCTR
460 C---VV IS THE LAYER-MEAN PRESSURE (IN ATM),WHICH IS NOT THE SAME AS
461 C THE LEVEL PRESSURE (PRESS)
462 VV(I,K)=HAF*(P(I,K+1)+P(I,K))*P0INV
463 VAR1(I,K)=DELP2(I,K)*QH2O(I,K)*GINV
464 VAR3(I,K)=DELP2(I,K)*QO3(I,K)*DIFFCTR*GINV
465 VAR2(I,K)=VAR1(I,K)*(VV(I,K)+H3M4)
466 VAR4(I,K)=VAR3(I,K)*(VV(I,K)+H3M3)
467 C COMPUTE OPTICAL PATH FOR THE H2O CONTINUUM, USING ROBERTS COEFFS.
468 C (BETINW),AND TEMP. CORRECTION (TEXPSL). THE DIFFUSIVITY FACTOR
469 C (WHICH CANCELS OUT IN THIS EXPRESSION) IS ASSUMED TO BE 1.66. THE
470 C USE OF THE DIFFUSIVITY FACTOR HAS BEEN SHOWN TO BE A SIGNIFICANT
471 C SOURCE OF ERROR IN THE CONTINUUM CALCS.,BUT THE TIME PENALTY OF
472 C AN ANGULAR INTEGRATION IS SEVERE.
473 C
474 CNTVAL(I,K)=TEXPSL(I,K)*RH2O(I,K)*VAR2(I,K)*BETINW/
475 1 (RH2O(I,K)+RATH2OMW)
476 131 CONTINUE
477 C COMPUTE SUMMED OPTICAL PATHS FOR H2O,O3 AND CONTINUUM
478 DO 201 I=MYIS,MYIE
479 TOTPHI(I,1)=ZERO
Page 8 Source Listing LWR88
2025-03-12 18:21 LWR88.F
480 TOTO3(I,1)=ZERO
481 TPHIO3(I,1)=ZERO
482 TOTVO2(I,1)=ZERO
483 201 CONTINUE
484 DO 203 K=2,LP1
485 DO 203 I=MYIS,MYIE
486 TOTPHI(I,K)=TOTPHI(I,K-1)+VAR2(I,K-1)
487 TOTO3(I,K)=TOTO3(I,K-1)+VAR3(I,K-1)
488 TPHIO3(I,K)=TPHIO3(I,K-1)+VAR4(I,K-1)
489 TOTVO2(I,K)=TOTVO2(I,K-1)+CNTVAL(I,K-1)
490 203 CONTINUE
491 C---EMX1 IS THE ADDITIONAL PRESSURE-SCALED MASS FROM PRESS(L) TO
492 C P(L). IT IS USED IN NEARBY LAYER AND EMISS CALCULATIONS.
493 C---EMX2 IS THE ADDITIONAL PRESSURE-SCALED MASS FROM PRESS(L) TO
494 C P(LP1). IT IS USED IN CALCULATIONS BETWEEN FLUX LEVELS L AND LP1.
495 C
496 DO 801 I=MYIS,MYIE
497 EMX1(I)=QH2O(I,L)*PRESS(I,L)*(PRESS(I,L)-P(I,L))*GP0INV
498 EMX2(I)=QH2O(I,L)*PRESS(I,L)*(P(I,LP1)-PRESS(I,L))*GP0INV
499 801 CONTINUE
500 C---EMPL IS THE PRESSURE SCALED MASS FROM P(K) TO PRESS(K) (INDEX 2-LP1)
501 C OR TO PRESS(K+1) (INDEX LP2-LL)
502 DO 811 K=1,L
503 DO 811 I=MYIS,MYIE
504 EMPL(I,K+1)=QH2O(I,K)*P(I,K+1)*(P(I,K+1)-PRESS(I,K))*GP0INV
505 811 CONTINUE
506 DO 812 K=1,LM1
507 DO 812 I=MYIS,MYIE
508 EMPL(I,LP2+K-1)=QH2O(I,K+1)*P(I,K+1)*(PRESS(I,K+1)-P(I,K+1))
509 1 *GP0INV
510 812 CONTINUE
511 DO 821 I=MYIS,MYIE
512 EMPL(I,1)=VAR2(I,L)
513 EMPL(I,LLP1)=EMPL(I,LL)
514 821 CONTINUE
515 C***COMPUTE WEIGHTED TEMPERATURE (TDAV) AND PRESSURE (TSTDAV) INTEGRALS
516 C FOR USE IN OBTAINING TEMP. DIFFERENCE BET. SOUNDING AND STD.
517 C TEMP. SOUNDING (DIFT)
518 DO 161 I=MYIS,MYIE
519 TSTDAV(I,1)=ZERO
520 TDAV(I,1)=ZERO
521 161 CONTINUE
522 DO 162 K=1,LP1
523 DO 162 I=MYIS,MYIE
524 VSUM3(I,K)=TEMP(I,K)-STEMP(K)
525 162 CONTINUE
526 DO 163 K=1,L
527 DO 165 I=MYIS,MYIE
528 VSUM2(I)=GTEMP(K)*DELP2(I,K)
529 VSUM1(I)=VSUM2(I)*VSUM3(I,K)
530 TSTDAV(I,K+1)=TSTDAV(I,K)+VSUM2(I)
531 TDAV(I,K+1)=TDAV(I,K)+VSUM1(I)
532 165 CONTINUE
533 163 CONTINUE
534 C
535 C****EVALUATE COEFFICIENTS FOR CO2 PRESSURE INTERPOLATION (A1,A2)
536 DO 171 I=MYIS,MYIE
Page 9 Source Listing LWR88
2025-03-12 18:21 LWR88.F
537 A1(I)=(PRESS(I,LP1)-P0XZP8)/P0XZP2
538 A2(I)=(P0-PRESS(I,LP1))/P0XZP2
539 171 CONTINUE
540 C***PERFORM CO2 PRESSURE INTERPOLATION ON ALL INPUTTED TRANSMISSION
541 C FUNCTIONS AND TEMP. DERIVATIVES
542 C---SUCCESSIVELY COMPUTING CO2R,DCO2DT AND D2CDT2 IS DONE TO SAVE
543 C STORAGE (AT A SLIGHT LOSS IN COMPUTATION TIME)
544 DO 184 K=1,LP1
545 DO 184 I=MYIS,MYIE
546 CO2R1(I,K)=A1(I)*CO231(K)+A2(I)*CO238(K)
547 D2CD21(I,K)=H1M3*(A1(I)*C2D31(K)+A2(I)*C2D38(K))
548 DCO2D1(I,K)=H1M2*(A1(I)*CDT31(K)+A2(I)*CDT38(K))
549 CO2R2(I,K)=A1(I)*CO271(K)+A2(I)*CO278(K)
550 D2CD22(I,K)=H1M3*(A1(I)*C2D71(K)+A2(I)*C2D78(K))
551 DCO2D2(I,K)=H1M2*(A1(I)*CDT71(K)+A2(I)*CDT78(K))
552 184 CONTINUE
553 DO 190 K=1,L
554 DO 190 I=MYIS,MYIE
555 CO2MR(I,K)=A1(I)*CO2M51(K)+A2(I)*CO2M58(K)
556 CO2MD(I,K)=H1M2*(A1(I)*CDTM51(K)+A2(I)*CDTM58(K))
557 CO2M2D(I,K)=H1M3*(A1(I)*C2DM51(K)+A2(I)*C2DM58(K))
558 190 CONTINUE
559 C***COMPUTE CO2 TEMPERATURE INTERPOLATIONS FOR ALL BANDS,USING DIFT
560 C
561 C THE CASE WHERE K=1 IS HANDLED FIRST. WE ARE NOW REPLACING
562 C 3-DIMENSIONAL ARRAYS BY 2-D ARRAYS, TO SAVE SPACE. THUS THIS
563 C CALCULATION IS FOR (I,KP,1)
564 DO 211 KP=2,LP1
565 DO 211 I=MYIS,MYIE
566 DIFT(I,KP)=TDAV(I,KP)/TSTDAV(I,KP)
567 211 CONTINUE
568 DO 212 I=MYIS,MYIE
569 CO21(I,1,1)=1.0
570 CO2SP1(I,1)=1.0
571 CO2SP2(I,1)=1.0
572 212 CONTINUE
573 DO 215 KP=2,LP1
574 DO 215 I=MYIS,MYIE
575 C---CALCULATIONS FOR KP>1 FOR K=1
576 CO2R(I,KP)=A1(I)*CO251(KP,1)+A2(I)*CO258(KP,1)
577 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(KP,1)+A2(I)*CDT58(KP,1))
578 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(KP,1)+A2(I)*C2D58(KP,1))
579 CO21(I,KP,1)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+
580 1 HAF*DIFT(I,KP)*D2CDT2(I,KP))
581 C---CALCULATIONS FOR (EFFECTIVELY) KP=1,K>KP. THESE USE THE
582 C SAME VALUE OF DIFT DUE TO SYMMETRY
583 CO2R(I,KP)=A1(I)*CO251(1,KP)+A2(I)*CO258(1,KP)
584 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(1,KP)+A2(I)*CDT58(1,KP))
585 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(1,KP)+A2(I)*C2D58(1,KP))
586 CO21(I,1,KP)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+
587 1 HAF*DIFT(I,KP)*D2CDT2(I,KP))
588 215 CONTINUE
589 C THE TRANSMISSION FUNCTIONS USED IN SPA88 MAY BE COMPUTED NOW.
590 C---(IN THE 250 LOOP,DIFT REALLY SHOULD BE (I,1,K), BUT DIFT IS
591 C INVARIANT WITH RESPECT TO K,KP,AND SO (I,1,K)=(I,K,1))
592 DO 250 K=2,LP1
593 DO 250 I=MYIS,MYIE
Page 10 Source Listing LWR88
2025-03-12 18:21 LWR88.F
594 CO2SP1(I,K)=CO2R1(I,K)+DIFT(I,K)*(DCO2D1(I,K)+HAF*DIFT(I,K)*
595 1 D2CD21(I,K))
596 CO2SP2(I,K)=CO2R2(I,K)+DIFT(I,K)*(DCO2D2(I,K)+HAF*DIFT(I,K)*
597 1 D2CD22(I,K))
598 250 CONTINUE
599 C
600 C NEXT THE CASE WHEN K=2...L
601 DO 220 K=2,L
602 DO 222 KP=K+1,LP1
603 DO 222 I=MYIS,MYIE
604 DIFT(I,KP)=(TDAV(I,KP)-TDAV(I,K))/
605 1 (TSTDAV(I,KP)-TSTDAV(I,K))
606 CO2R(I,KP)=A1(I)*CO251(KP,K)+A2(I)*CO258(KP,K)
607 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(KP,K)+A2(I)*CDT58(KP,K))
608 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(KP,K)+A2(I)*C2D58(KP,K))
609 CO21(I,KP,K)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+
610 1 HAF*DIFT(I,KP)*D2CDT2(I,KP))
611 CO2R(I,KP)=A1(I)*CO251(K,KP)+A2(I)*CO258(K,KP)
612 DCO2DT(I,KP)=H1M2*(A1(I)*CDT51(K,KP)+A2(I)*CDT58(K,KP))
613 D2CDT2(I,KP)=H1M3*(A1(I)*C2D51(K,KP)+A2(I)*C2D58(K,KP))
614 CO21(I,K,KP)=CO2R(I,KP)+DIFT(I,KP)*(DCO2DT(I,KP)+
615 1 HAF*DIFT(I,KP)*D2CDT2(I,KP))
616 222 CONTINUE
617 220 CONTINUE
618 C FINALLY THE CASE WHEN K=KP,K=2..LP1
619 DO 206 K=2,LP1
620 DO 206 I=MYIS,MYIE
621 DIFT(I,K)=HAF*(VSUM3(I,K)+VSUM3(I,K-1))
622 CO2R(I,K)=A1(I)*CO251(K,K)+A2(I)*CO258(K,K)
623 DCO2DT(I,K)=H1M2*(A1(I)*CDT51(K,K)+A2(I)*CDT58(K,K))
624 D2CDT2(I,K)=H1M3*(A1(I)*C2D51(K,K)+A2(I)*C2D58(K,K))
625 CO21(I,K,K)=CO2R(I,K)+DIFT(I,K)*(DCO2DT(I,K)+
626 1 HAF*DIFT(I,K)*D2CDT2(I,K))
627 206 CONTINUE
628 C--- WE AREN'T DOING NBL TFS ON THE 100 CM-1 BANDS .
629 DO 260 K=1,L
630 DO 260 I=MYIS,MYIE
631 CO2NBL(I,K)=CO2MR(I,K)+VSUM3(I,K)*(CO2MD(I,K)+HAF*
632 1 VSUM3(I,K)*CO2M2D(I,K))
633 260 CONTINUE
634 C***COMPUTE TEMP. COEFFICIENT BASED ON T(K) (SEE REF.2)
635 DO 264 K=1,LP1
636 DO 264 I=MYIS,MYIE
637 IF (T(I,K).LE.H25E2) THEN
638 TLSQU(I,K)=B0+(T(I,K)-H25E2)*
639 1 (B1+(T(I,K)-H25E2)*
640 2 (B2+B3*(T(I,K)-H25E2)))
641 ELSE
642 TLSQU(I,K)=B0
643 ENDIF
644 264 CONTINUE
645 C***APPLY TO ALL CO2 TFS
646 DO 280 K=1,LP1
647 DO 282 KP=1,LP1
648 DO 282 I=MYIS,MYIE
649 CO21(I,KP,K)=CO21(I,KP,K)*(ONE-TLSQU(I,KP))+TLSQU(I,KP)
650 282 CONTINUE
Page 11 Source Listing LWR88
2025-03-12 18:21 LWR88.F
651 280 CONTINUE
652 DO 284 K=1,LP1
653 DO 286 I=MYIS,MYIE
654 CO2SP1(I,K)=CO2SP1(I,K)*(ONE-TLSQU(I,1))+TLSQU(I,1)
655 CO2SP2(I,K)=CO2SP2(I,K)*(ONE-TLSQU(I,1))+TLSQU(I,1)
656 286 CONTINUE
657 284 CONTINUE
658 DO 288 K=1,L
659 DO 290 I=MYIS,MYIE
660 CO2NBL(I,K)=CO2NBL(I,K)*(ONE-TLSQU(I,K))+TLSQU(I,K)
661 290 CONTINUE
662 288 CONTINUE
663 CALL FST88(HEATRA,GRNFLX,TOPFLX,
664 1 QH2O,PRESS,P,DELP,DELP2,TEMP,T,
665 2 CLDFAC,NCLDS,KTOP,KBTM,CAMT,
666 3 CO21,CO2NBL,CO2SP1,CO2SP2,
667 4 VAR1,VAR2,VAR3,VAR4,CNTVAL,
668 5 TOTO3,TPHIO3,TOTPHI,TOTVO2,
669 6 EMX1,EMX2,EMPL)
670 RETURN
671 END
ENTRY POINTS
Name
lwr88_
Page 12 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
SYMBOL CROSS REFERENCE
Name Object Declared Type Bytes Dimen Elements Attributes References
103 Label 430 426,427
105 Label 436 431
107 Label 441 437,438
125 Label 449 444,445
131 Label 476 457,458
161 Label 521 518
162 Label 525 522,523
163 Label 533 526
165 Label 532 527
171 Label 539 536
184 Label 552 544,545
190 Label 558 553,554
201 Label 483 478
203 Label 490 484,485
206 Label 627 619,620
211 Label 567 564,565
212 Label 572 568
215 Label 588 573,574
220 Label 617 601
222 Label 616 602,603
250 Label 598 592,593
260 Label 633 629,630
264 Label 644 635,636
280 Label 651 646
282 Label 650 647,648
284 Label 657 652
286 Label 656 653
288 Label 662 658
290 Label 661 659
801 Label 499 496
811 Label 505 502,503
812 Label 510 506,507
821 Label 514 511
A1 Local 414 R(4) 4 1 45 537,546,547,548,549,550,551,555,55
6,557,576,577,578,583,584,585,606,
607,608,611,612,613,622,623,624
A2 Local 414 R(4) 4 1 45 538,546,547,548,549,550,551,555,55
6,557,576,577,578,583,584,585,606,
607,608,611,612,613,622,623,624
BANDTA Common 310 5900
BDCOMB Common 378 736
BDWIDE Common 338 56
CAMT Dummy 4 R(4) 4 2 2070 ARG,INOUT 665
CLDFAC Dummy 3 R(4) 4 3 95220 ARG,INOUT 665
CNTVAL Local 399 R(4) 4 2 2070 474,489,667
CO21 Local 395 R(4) 4 3 95220 569,579,586,609,614,625,649,666
CO2BD2 Common 250 1104
CO2BD3 Common 233 52248
CO2BD4 Common 265 1104
CO2BD5 Common 275 368
CO2M2D Local 410 R(4) 4 2 2025 557,632
Page 13 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Object Declared Type Bytes Dimen Elements Attributes References
CO2MD Local 409 R(4) 4 2 2025 556,631
CO2MR Local 409 R(4) 4 2 2025 555,631
CO2NBL Local 395 R(4) 4 2 2025 631,660,666
CO2R Local 405 R(4) 4 2 2070 576,579,583,586,606,609,611,614,62
2,625
CO2R1 Local 406 R(4) 4 2 2070 546,594
CO2R2 Local 408 R(4) 4 2 2070 549,596
CO2SP1 Local 396 R(4) 4 2 2070 570,594,654,666
CO2SP2 Local 396 R(4) 4 2 2070 571,596,655,666
D2CD21 Local 407 R(4) 4 2 2070 547,595
D2CD22 Local 407 R(4) 4 2 2070 550,597
D2CDT2 Local 415 R(4) 4 2 2070 578,580,585,587,608,610,613,615,62
4,626
DCO2D1 Local 406 R(4) 4 2 2070 548,594
DCO2D2 Local 408 R(4) 4 2 2070 551,596
DCO2DT Local 415 R(4) 4 2 2070 577,579,584,586,607,609,612,614,62
3,625
DELP Local 394 R(4) 4 2 2025 440,664
DELP2 Local 391 R(4) 4 2 2025 439,440,463,464,528,664
DIFT Local 405 R(4) 4 2 2070 566,579,580,586,587,594,596,604,60
9,610,614,615,621,625,626
EMPL Local 403 R(4) 4 2 4095 504,508,512,513,669
EMX1 Local 402 R(4) 4 1 45 497,669
EMX2 Local 403 R(4) 4 1 45 498,669
EXP Func 448 scalar 448
FST88 Subr 663 663
GLB_TABLE Common 151 128 SAVE
GRNFLX Dummy 2 R(4) 4 1 45 ARG,INOUT 663
HCON Common 34 872
HEATRA Dummy 2 R(4) 4 2 2025 ARG,INOUT 663
I Local 427 I(4) 4 scalar 427,428,429,431,432,433,434,435,43
8,439,440,445,446,448,458,459,462,
463,464,465,466,474,475,478,479,48
0,481,482,485,486,487,488,489,496,
497,498,503,504,507,508,511,512,51
3,518,519,520,523,524,527,528,529,
530,531,536,537,538,545,546,547,54
8,549,550,551,554,555,556,557,565,
566,568,569,570,571,574,576,577,57
8,579,580,583,584,585,586,587,593,
594,595,596,597,603,604,605,606,60
7,608,609,610,611,612,613,614,615,
620,621,622,623,624,625,626,630,63
1,632,636,637,638,639,640,642,648,
649,653,654,655,659,660
IDIM1 Param 94 I(4) 4 scalar 384,385,386,387,388,389,390,391,39
3,394,395,396,397,398,399,400,401,
402,403,405,406,407,408,409,410,41
1,412,413,414,415,417,418
IDIM2 Param 94 I(4) 4 scalar 384,385,386,387,388,389,390,391,39
3,394,395,396,397,398,399,400,401,
402,403,405,406,407,408,409,410,41
1,412,413,414,415,417,418
IGSTL Param 89 I(4) 4 scalar 94
IGSTR Param 89 I(4) 4 scalar 94
Page 14 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Object Declared Type Bytes Dimen Elements Attributes References
IM Param 78 I(4) 4 scalar 91,94,157,158,159,160,161,165,166,
171
INLTE Param 185 I(4) 4 scalar 185
INLTEP Param 185 I(4) 4 scalar
INPES Param 86 I(4) 4 scalar 91,94,141,142,143,152,153
ITAIL Param 91 I(4) 4 scalar
JDIM1 Param 95 I(4) 4 scalar
JDIM2 Param 95 I(4) 4 scalar
JGSTL Param 90 I(4) 4 scalar 95
JGSTR Param 90 I(4) 4 scalar 95
JM Param 78 I(4) 4 scalar 92,95,157,158,159,160,161,165,166,
171
JNPES Param 86 I(4) 4 scalar 92,95,141,142,143,152,153
JTAIL Param 92 I(4) 4 scalar
K Local 426 I(4) 4 scalar 426,428,429,437,439,440,444,446,44
8,457,459,462,463,464,465,466,474,
475,484,486,487,488,489,502,504,50
6,508,522,524,526,528,529,530,531,
544,546,547,548,549,550,551,553,55
5,556,557,592,594,595,596,597,601,
602,604,605,606,607,608,609,611,61
2,613,614,619,621,622,623,624,625,
626,629,631,632,635,637,638,639,64
0,642,646,649,652,654,655,658,660
KBTM Dummy 4 I(4) 4 2 2070 ARG,INOUT 665
KO2 Param 187 I(4) 4 scalar 188
KO21 Param 188 I(4) 4 scalar
KO2M Param 188 I(4) 4 scalar
KP Local 564 I(4) 4 scalar 564,566,573,576,577,578,579,580,58
3,584,585,586,587,602,604,605,606,
607,608,609,610,611,612,613,614,61
5,647,649
KTOP Dummy 4 I(4) 4 2 2070 ARG,INOUT 665
L Param 173 I(4) 4 scalar 176,177,178,181,183,234,235,385,38
9,391,393,394,395,397,398,409,410,
412,418,426,437,457,497,498,502,51
2,526,553,601,629,658
LL Param 178 I(4) 4 scalar 178,179,513
LL3P Param 183 I(4) 4 scalar
LLM1 Param 179 I(4) 4 scalar
LLM2 Param 179 I(4) 4 scalar
LLM3 Param 179 I(4) 4 scalar
LLP1 Param 178 I(4) 4 scalar 403,513
LLP2 Param 178 I(4) 4 scalar
LLP3 Param 178 I(4) 4 scalar
LM Param 78 I(4) 4 scalar 166,173
LM1 Param 177 I(4) 4 scalar 506
LM2 Param 177 I(4) 4 scalar
LM3 Param 177 I(4) 4 scalar
LP1 Param 176 I(4) 4 scalar 180,181,182,233,234,236,250,251,26
5,266,275,384,386,387,388,393,394,
395,396,399,400,401,402,405,406,40
7,408,411,412,415,417,433,435,444,
484,498,522,537,538,544,564,573,59
2,602,619,635,646,647,652
Page 15 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Object Declared Type Bytes Dimen Elements Attributes References
LP121 Param 182 I(4) 4 scalar
LP1M Param 180 I(4) 4 scalar 180
LP1M1 Param 180 I(4) 4 scalar
LP1V Param 181 I(4) 4 scalar
LP2 Param 176 I(4) 4 scalar 508
LP3 Param 176 I(4) 4 scalar
LSM Param 78 I(4) 4 scalar
LWR88 Subr 2
MAPPINGS Common 170 5024 SAVE
MPPCOM Common 115 1464 SAVE
NB Param 184 I(4) 4 scalar 186
NB1 Param 186 I(4) 4 scalar
NBLM Param 175 I(4) 4 scalar
NBLW Param 174 I(4) 4 scalar 310,311,312
NBLX Param 174 I(4) 4 scalar
NBLY Param 174 I(4) 4 scalar 175,182,378,379,380
NCLDS Dummy 4 I(4) 4 1 45 ARG,INOUT 665
NNLTE Param 185 I(4) 4 scalar
P Local 394 R(4) 4 2 2070 428,432,433,439,462,497,498,504,50
8,664
PHYCON Common 29 80
PRESS Dummy 3 R(4) 4 2 2070 ARG,INOUT 428,433,497,498,504,508,537,538,66
4
QH2O Local 393 R(4) 4 2 2025 459,463,497,498,504,508,664
QO3 Dummy 3 R(4) 4 2 2025 ARG,INOUT 464
RH2O Dummy 3 R(4) 4 2 2025 ARG,INOUT 459,474,475
T Local 393 R(4) 4 2 2070 429,434,435,637,638,639,640,664
TDAV Local 411 R(4) 4 2 2070 520,531,566,604
TEMP Dummy 3 R(4) 4 2 2070 ARG,INOUT 429,434,435,446,524,664
TEMPCOM Common 156 6603768 SAVE
TEXPSL Local 417 R(4) 4 2 2070 446,448,474
TLSQU Local 417 R(4) 4 2 2070 638,642,649,654,655,660
TOPFLX Dummy 2 R(4) 4 1 45 ARG,INOUT 663
TOPO Common 164 17478548 SAVE
TOTO3 Local 400 R(4) 4 2 2070 480,487,668
TOTPHI Local 401 R(4) 4 2 2070 479,486,668
TOTVO2 Local 402 R(4) 4 2 2070 482,489,668
TPHIO3 Local 400 R(4) 4 2 2070 481,488,668
TSTDAV Local 411 R(4) 4 2 2070 519,530,566,605
VAR1 Local 397 R(4) 4 2 2025 463,465,667
VAR2 Local 397 R(4) 4 2 2025 465,474,486,512,667
VAR3 Local 398 R(4) 4 2 2025 464,466,487,667
VAR4 Local 398 R(4) 4 2 2025 466,488,667
VSUM1 Local 412 R(4) 4 1 45 529,531
VSUM2 Local 413 R(4) 4 1 45 528,529,530
VSUM3 Local 412 R(4) 4 2 2070 524,529,621,631,632
VSUM4 Local 418 R(4) 4 2 2025
VV Local 412 R(4) 4 2 2025 462,465,466
Page 16 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
TYPE COMPONENTS/COMMON VARIABLES
Name Type Bytes Offset Dimen Elements Attributes References
AB15 R(4) 4 5892 1 2 COM
AB15CM R(4) 4 728 1 2 COM
AB15WD R(4) 4 40 scalar COM
ACOMB R(4) 4 160 1 15 COM
AMOLWT R(4) 4 0 scalar COM
AO3CM R(4) 4 704 1 3 COM
AO3RND R(4) 4 5868 1 3 COM
AP R(4) 4 1956 1 163 COM
APCM R(4) 4 340 1 15 COM
APWD R(4) 4 12 scalar COM
ARNDM R(4) 4 0 1 163 COM
ATP R(4) 4 3260 1 163 COM
ATPCM R(4) 4 460 1 15 COM
ATPWD R(4) 4 20 scalar COM
AWIDE R(4) 4 0 scalar COM
B0 R(4) 4 52232 scalar COM 638,642
B1 R(4) 4 52236 scalar COM 639
B2 R(4) 4 52240 scalar COM 640
B3 R(4) 4 52244 scalar COM 640
BANDHI R(4) 4 5216 1 163 COM
BANDLO R(4) 4 4564 1 163 COM
BCOMB R(4) 4 220 1 15 COM
BDHICM R(4) 4 640 1 15 COM
BDHIWD R(4) 4 32 scalar COM
BDLOCM R(4) 4 580 1 15 COM
BDLOWD R(4) 4 28 scalar COM
BETACM R(4) 4 280 1 15 COM
BETAD R(4) 4 1304 1 163 COM
BETAWD R(4) 4 8 scalar COM
BETINC R(4) 4 700 scalar COM
BETINW R(4) 4 36 scalar COM 474
BO3CM R(4) 4 716 1 3 COM
BO3RND R(4) 4 5880 1 3 COM
BP R(4) 4 2608 1 163 COM
BPCM R(4) 4 400 1 15 COM
BPWD R(4) 4 16 scalar COM
BRNDM R(4) 4 652 1 163 COM
BTP R(4) 4 3912 1 163 COM
BTPCM R(4) 4 520 1 15 COM
BTPWD R(4) 4 24 scalar COM
BWIDE R(4) 4 4 scalar COM
C2D31 R(4) 4 736 1 46 COM 547
C2D38 R(4) 4 920 1 46 COM 547
C2D51 R(4) 4 33856 2 2116 COM 578,585,608,613,624
C2D58 R(4) 4 42320 2 2116 COM 578,585,608,613,624
C2D71 R(4) 4 736 1 46 COM 550
C2D78 R(4) 4 920 1 46 COM 550
C2DM51 R(4) 4 51504 1 45 COM 557
C2DM58 R(4) 4 51684 1 45 COM 557
CDT31 R(4) 4 368 1 46 COM 548
CDT38 R(4) 4 552 1 46 COM 548
CDT51 R(4) 4 16928 2 2116 COM 577,584,607,612,623
CDT58 R(4) 4 25392 2 2116 COM 577,584,607,612,623
Page 17 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
CDT71 R(4) 4 368 1 46 COM 551
CDT78 R(4) 4 552 1 46 COM 551
CDTM51 R(4) 4 51144 1 45 COM 556
CDTM58 R(4) 4 51324 1 45 COM 556
CO211 R(4) 4 0 1 46 COM
CO218 R(4) 4 184 1 46 COM
CO231 R(4) 4 0 1 46 COM 546
CO238 R(4) 4 184 1 46 COM 546
CO251 R(4) 4 0 2 2116 COM 576,583,606,611,622
CO258 R(4) 4 8464 2 2116 COM 576,583,606,611,622
CO271 R(4) 4 0 1 46 COM 549
CO278 R(4) 4 184 1 46 COM 549
CO2M51 R(4) 4 50784 1 45 COM 555
CO2M58 R(4) 4 50964 1 45 COM 555
CSUBP R(4) 4 4 scalar COM
DIFFCTR R(4) 4 8 scalar COM 459,464
EIGHT R(4) 4 20 scalar COM
FIFTY R(4) 4 12 scalar COM
FIVE R(4) 4 24 scalar COM
FOUR R(4) 4 28 scalar COM
G R(4) 4 12 scalar COM
G2LI I(4) 4 0 1 239 COM
G2LJ I(4) 4 1912 1 389 COM
GINV R(4) 4 68 scalar COM 463,464
GP0INV R(4) 4 76 scalar COM 497,498,504,509
GRAVDR R(4) 4 16 scalar COM
GTEMP R(4) 4 52048 1 46 COM 528
H101M16 R(4) 4 808 scalar COM
H102M5 R(4) 4 784 scalar COM
H1036E2 R(4) 4 756 scalar COM
H114M11 R(4) 4 384 scalar COM
H1174M7 R(4) 4 268 scalar COM
H11M10 R(4) 4 296 scalar COM
H11M11 R(4) 4 388 scalar COM
H1224E3 R(4) 4 112 scalar COM
H1226E1 R(4) 4 652 scalar COM
H128M5 R(4) 4 240 scalar COM
H129M2 R(4) 4 732 scalar COM
H12M12 R(4) 4 472 scalar COM
H12M13 R(4) 4 492 scalar COM
H12M31 R(4) 4 560 scalar COM
H135M13 R(4) 4 488 scalar COM
H1386E2 R(4) 4 748 scalar COM
H1439M5 R(4) 4 236 scalar COM
H14M10 R(4) 4 292 scalar COM
H14M11 R(4) 4 380 scalar COM
H14M12 R(4) 4 468 scalar COM
H14M14 R(4) 4 508 scalar COM
H14M30 R(4) 4 548 scalar COM
H15E2 R(4) 4 136 scalar COM
H15M11 R(4) 4 376 scalar COM
H15M14 R(4) 4 504 scalar COM
H15M5 R(4) 4 824 scalar COM
H161E1 R(4) 4 800 scalar COM
H165E5 R(4) 4 84 scalar COM
Page 18 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
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 446
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 548,551,556,577,584,607,612,623
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 547,550,557,578,585,608,613,624
H1M4 R(4) 4 220 scalar COM
H1M5 R(4) 4 244 scalar COM
H1M6 R(4) 4 260 scalar COM
H1M60 R(4) 4 584 scalar COM
H1M8 R(4) 4 284 scalar COM
H1P082 R(4) 4 740 scalar COM
H1P25892 R(4) 4 180 scalar COM
H1P4 R(4) 4 176 scalar COM
H1P41819 R(4) 4 704 scalar COM
H1P4387 R(4) 4 172 scalar COM
H1P8 R(4) 4 168 scalar COM
H2075E3 R(4) 4 108 scalar COM
H20788E3 R(4) 4 104 scalar COM
H2118M2 R(4) 4 760 scalar COM
H21M12 R(4) 4 460 scalar COM
H21M31 R(4) 4 556 scalar COM
H235M3 R(4) 4 724 scalar COM
H23E2 R(4) 4 132 scalar COM
H23M10 R(4) 4 288 scalar COM
H23M11 R(4) 4 364 scalar COM
H24E3 R(4) 4 100 scalar COM
H24M11 R(4) 4 360 scalar COM
H24M12 R(4) 4 456 scalar COM
H25452M6 R(4) 4 256 scalar COM
H257M8 R(4) 4 280 scalar COM
H25E2 R(4) 4 816 scalar COM 637,638,639,640
H25M31 R(4) 4 552 scalar COM
H26E2 R(4) 4 696 scalar COM
H26M30 R(4) 4 544 scalar COM
Page 19 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
H28571M2 R(4) 4 852 scalar COM
H285M4 R(4) 4 216 scalar COM
H28E1 R(4) 4 868 scalar COM
H28M11 R(4) 4 356 scalar COM
H28M12 R(4) 4 452 scalar COM
H29316E2 R(4) 4 648 scalar COM
H2945E2 R(4) 4 128 scalar COM
H29M12 R(4) 4 448 scalar COM
H2E2 R(4) 4 688 scalar COM
H2E6 R(4) 4 616 scalar COM
H2M11 R(4) 4 368 scalar COM
H2P5 R(4) 4 164 scalar COM
H2P8 R(4) 4 160 scalar COM
H2P9 R(4) 4 156 scalar COM
H3082E2 R(4) 4 120 scalar COM
H3116E1 R(4) 4 656 scalar COM
H323M4 R(4) 4 768 scalar COM
H327M8 R(4) 4 276 scalar COM
H32M11 R(4) 4 348 scalar COM
H35E1 R(4) 4 140 scalar COM
H35M11 R(4) 4 344 scalar COM
H36M13 R(4) 4 484 scalar COM
H37412M5 R(4) 4 232 scalar COM
H37M11 R(4) 4 340 scalar COM
H37M12 R(4) 4 440 scalar COM
H38M12 R(4) 4 436 scalar COM
H391M7 R(4) 4 264 scalar COM
H394M5 R(4) 4 228 scalar COM
H3E2 R(4) 4 124 scalar COM
H3M11 R(4) 4 352 scalar COM
H3M12 R(4) 4 444 scalar COM
H3M14 R(4) 4 500 scalar COM
H3M3 R(4) 4 804 scalar COM 466
H3M4 R(4) 4 860 scalar COM 465
H3P5 R(4) 4 828 scalar COM
H3P6 R(4) 4 144 scalar COM
H41666M2 R(4) 4 636 scalar COM
H42M11 R(4) 4 336 scalar COM
H42M2 R(4) 4 764 scalar COM
H44194M2 R(4) 4 700 scalar COM
H44871M2 R(4) 4 200 scalar COM
H44M11 R(4) 4 332 scalar COM
H44M12 R(4) 4 428 scalar COM
H451M6 R(4) 4 788 scalar COM
H45M12 R(4) 4 424 scalar COM
H45M32 R(4) 4 572 scalar COM
H46M13 R(4) 4 480 scalar COM
H488E4 R(4) 4 92 scalar COM
H48M11 R(4) 4 328 scalar COM
H4999M6 R(4) 4 252 scalar COM
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
Page 20 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
H55M32 R(4) 4 568 scalar COM
H5725E4 R(4) 4 88 scalar COM
H5E2 R(4) 4 116 scalar COM
H625M2 R(4) 4 668 scalar COM
H62M12 R(4) 4 416 scalar COM
H62M34 R(4) 4 580 scalar COM
H658M2 R(4) 4 752 scalar COM
H65M12 R(4) 4 412 scalar COM
H67390E2 R(4) 4 772 scalar COM
H6938M5 R(4) 4 224 scalar COM
H69766E5 R(4) 4 720 scalar COM
H6M12 R(4) 4 420 scalar COM
H6P08108 R(4) 4 836 scalar COM 446
H71E26 R(4) 4 60 scalar COM
H72M11 R(4) 4 320 scalar COM
H74M12 R(4) 4 408 scalar COM
H75826M4 R(4) 4 736 scalar COM
H77M11 R(4) 4 316 scalar COM
H77M12 R(4) 4 404 scalar COM
H7M6 R(4) 4 248 scalar COM
H8121E1 R(4) 4 684 scalar COM
H82M11 R(4) 4 308 scalar COM
H83E26 R(4) 4 56 scalar COM
H83M11 R(4) 4 304 scalar COM
H8725M8 R(4) 4 272 scalar COM
H8M11 R(4) 4 312 scalar COM
H8M13 R(4) 4 476 scalar COM
H93M12 R(4) 4 400 scalar COM
H96M12 R(4) 4 396 scalar COM
H987M4 R(4) 4 212 scalar COM
H9M32 R(4) 4 564 scalar COM
H9P94 R(4) 4 660 scalar COM
HAF R(4) 4 44 scalar COM 428,429,462,580,587,594,596,610,61
5,621,626,631
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
HM2M2 R(4) 4 644 scalar COM
HM6666M2 R(4) 4 628 scalar COM
HM8E1 R(4) 4 864 scalar COM
HMP5 R(4) 4 640 scalar COM
HMP575 R(4) 4 588 scalar COM
HMP66667 R(4) 4 624 scalar COM
HMP805 R(4) 4 840 scalar COM
HNINETY R(4) 4 4 scalar COM
HP1 R(4) 4 196 scalar COM
HP118666 R(4) 4 820 scalar COM
HP144 R(4) 4 712 scalar COM
HP166666 R(4) 4 632 scalar COM
Page 21 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
HP219 R(4) 4 708 scalar COM
HP228 R(4) 4 672 scalar COM
HP26 R(4) 4 728 scalar COM
HP369 R(4) 4 192 scalar COM
HP3795 R(4) 4 776 scalar COM
HP5048 R(4) 4 780 scalar COM
HP518 R(4) 4 188 scalar COM
HP526315 R(4) 4 848 scalar COM
HP6 R(4) 4 664 scalar COM
HP602409 R(4) 4 844 scalar COM
HP60241 R(4) 4 676 scalar COM
HP8 R(4) 4 184 scalar COM
HP805 R(4) 4 744 scalar COM
HP816 R(4) 4 716 scalar COM
HTMG R(4) 4 743768 3 4183695 COM
HUNDRED R(4) 4 0 scalar COM
IBAND I(4) 4 0 1 40 COM
IBROW I(4) 4 404 scalar COM
ICHUNKTAB I(4) 4 620 1 8 COM
IE_GLB_TABLE I(4) 4 32 1 8 COM
IE_LOC_TABLE I(4) 4 556 1 8 COM
ILCOL I(4) 4 396 scalar COM
ILPAD1 I(4) 4 412 scalar COM
ILPAD2 I(4) 4 416 scalar COM
ILPAD3 I(4) 4 420 scalar COM
ILPAD4 I(4) 4 424 scalar COM
ILPAD5 I(4) 4 428 scalar COM
INUMQ I(4) 4 1064 1 100 COM
IQUILT_GROUP I(4) 4 1060 scalar COM
IRCOL I(4) 4 400 scalar COM
IRPAD1 I(4) 4 432 scalar COM
IRPAD2 I(4) 4 436 scalar COM
IRPAD3 I(4) 4 440 scalar COM
IRPAD4 I(4) 4 444 scalar COM
IRPAD5 I(4) 4 448 scalar COM
IS_GLB_TABLE I(4) 4 0 1 8 COM
IS_LOC_TABLE I(4) 4 492 1 8 COM
ITEMP I(4) 4 5870016 2 91719 COM
ITEMP2 I(4) 4 6236892 2 91719 COM
ITROW I(4) 4 408 scalar COM
JBPAD1 I(4) 4 452 scalar COM
JBPAD2 I(4) 4 456 scalar COM
JBPAD3 I(4) 4 460 scalar COM
JBPAD4 I(4) 4 464 scalar COM
JBPAD5 I(4) 4 468 scalar COM
JE_GLB_TABLE I(4) 4 96 1 8 COM
JE_LOC_TABLE I(4) 4 588 1 8 COM
JS_GLB_TABLE I(4) 4 64 1 8 COM
JS_LOC_TABLE I(4) 4 524 1 8 COM
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
Page 22 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
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 427,431,438,445,458,478,485,496,50
3,507,511,518,523,527,536,545,554,
565,568,574,593,603,620,630,636,64
8,653,659
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 427,431,438,445,458,478,485,496,50
3,507,511,518,523,527,536,545,554,
565,568,574,593,603,620,630,636,64
8,653,659
MYIS1 I(4) 4 44 scalar COM
MYIS1_P1 I(4) 4 84 scalar COM
MYIS1_P2 I(4) 4 88 scalar COM
MYIS1_P3 I(4) 4 92 scalar COM
MYIS1_P4 I(4) 4 96 scalar COM
MYIS2 I(4) 4 48 scalar COM
MYIS_P1 I(4) 4 64 scalar COM
MYIS_P2 I(4) 4 68 scalar COM
MYIS_P3 I(4) 4 72 scalar COM
MYIS_P4 I(4) 4 76 scalar COM
MYIS_P5 I(4) 4 80 scalar COM
MYJE I(4) 4 236 scalar COM
MYJE1 I(4) 4 240 scalar COM
MYJE1_P1 I(4) 4 280 scalar COM
MYJE1_P2 I(4) 4 284 scalar COM
MYJE1_P3 I(4) 4 288 scalar COM
MYJE1_P4 I(4) 4 292 scalar COM
MYJE2 I(4) 4 244 scalar COM
MYJE2_P1 I(4) 4 296 scalar COM
MYJE2_P2 I(4) 4 300 scalar COM
MYJE2_P3 I(4) 4 304 scalar COM
MYJE2_P4 I(4) 4 308 scalar COM
MYJE3 I(4) 4 248 scalar COM
MYJE3_P4 I(4) 4 312 scalar COM
MYJE4 I(4) 4 252 scalar COM
MYJE4_P1 I(4) 4 316 scalar COM
MYJE4_P4 I(4) 4 320 scalar COM
MYJE5 I(4) 4 256 scalar COM
MYJE5_P1 I(4) 4 324 scalar COM
MYJE5_P2 I(4) 4 328 scalar COM
MYJE_P1 I(4) 4 260 scalar COM
Page 23 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
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
MY_JS_GLB I(4) 4 16 scalar COM
MY_JS_LOC I(4) 4 32 scalar COM
MY_N I(4) 4 332 scalar COM
MY_NE I(4) 4 348 scalar COM
MY_NEB I(4) 4 364 1 8 COM
MY_NW I(4) 4 360 scalar COM
MY_S I(4) 4 340 scalar COM
MY_SE I(4) 4 352 scalar COM
MY_SW I(4) 4 356 scalar COM
MY_W I(4) 4 344 scalar COM
NPES I(4) 4 4 scalar COM
O3DIFCTR R(4) 4 20 scalar COM
ONE R(4) 4 40 scalar COM 440,649,654,655,660
P0 R(4) 4 24 scalar COM 538
P0INV R(4) 4 72 scalar COM 462
P0X2 R(4) 4 36 scalar COM
P0XZP2 R(4) 4 28 scalar COM 537,538
P0XZP8 R(4) 4 32 scalar COM 537
QUARTR R(4) 4 48 scalar COM
Page 24 Source Listing LWR88
2025-03-12 18:21 Symbol Table LWR88.F
Name Type Bytes Offset Dimen Elements Attributes References
RADCON R(4) 4 40 scalar COM
RADCON1 R(4) 4 64 scalar COM
RATCO2MW R(4) 4 56 scalar COM
RATH2OMW R(4) 4 60 scalar COM 475
RGAS R(4) 4 44 scalar COM
RGASSP R(4) 4 48 scalar COM
SECPDA R(4) 4 52 scalar COM
SIXTY R(4) 4 8 scalar COM
SKC1R R(4) 4 48 scalar COM
SKO2D R(4) 4 44 scalar COM
SKO3R R(4) 4 52 scalar COM
STEMP R(4) 4 51864 1 46 COM 524
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
ZERO R(4) 4 52 scalar COM 432,479,480,481,482,519,520
Page 25 Source Listing LWR88
2025-03-12 18:21 Subprograms/Common Blocks LWR88.F
SUBPROGRAMS/COMMON BLOCKS
Name Object Declared Type Bytes Dimen Elements Attributes References
BANDTA Common 310 5900
BDCOMB Common 378 736
BDWIDE Common 338 56
CO2BD2 Common 250 1104
CO2BD3 Common 233 52248
CO2BD4 Common 265 1104
CO2BD5 Common 275 368
GLB_TABLE Common 151 128 SAVE
HCON Common 34 872
LWR88 Subr 2
MAPPINGS Common 170 5024 SAVE
MPPCOM Common 115 1464 SAVE
PHYCON Common 29 80
TEMPCOM Common 156 6603768 SAVE
TOPO Common 164 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
Page 26 Source Listing LWR88
2025-03-12 18:21 LWR88.F
-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
-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
Page 27 Source Listing LWR88
2025-03-12 18:21 LWR88.F
-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
-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 : LWR88.lst
no -o
COMPILER: Intel(R) Fortran 19.1-1655