1. C**** 2. C**** C090R.S Fortran source code for radiation routines 2003/11/19 3. C**** 4. C**** NASA/GISS Climate Model III Programmed by: 5. C**** Andrew A. Lacis 6. C**** NASA/Goddard Space Flight Center Gary L. Russell 7. C**** Institute for Space Studies 8. C**** 2880 Broadway, New York, NY 10025 9. C**** U.S.A. 10. C**** 11. SUBROUTINE RCOMP0 12. C**** 13. INCLUDE 'RA99M9.COM' 15. CHARACTER TITLE*80 16. C CHARACTER LABLAT*16,LABLON*16,LABLEV*16,LABTIM*16, NASAGISS*8 17. REAL*4 GHGC(6,1950:1990) 18. C REAL*4 VAOT(JM,12,1880:1991), RLON,RLAT(JM),RLEV,RTIM(12) 19. REAL*4 SAP1(IM,JM),SAP2(IM,JM), SAOTGA(2,1990:2100), SAOT(IM,JM) 20. C REAL*8 TKOFPF,WAVNA,WAVNB,PFWI 21. C**** 22. C**** RCOMP0 is called once at the beginning of a run 23. C**** 24. C**** Read in distributed PLANCK function and gas TAU table from disk 25. READ (12) PLANCK 26. CLOSE (12) 27. WRITE (6,*) 'Distributed PLANCK function read in on unit 12.' 28. READ (13) TAUTBL 29. CLOSE (13) 30. WRITE (6,*) 'Gas TAU table read in on unit 13.' 31. C**** Set permanent greenhouse gas concentrations (ppm) 32. PPMV(2) = 311.87d0 ! 1950 CO2 33. PPMV(6) = 285.27d-3 ! 1950 N2O = 283.54d-3 34. PPMV(7) = 1181.71d-3 ! 1950 CH4 = 1150.00d-3 35. PPMV(8) = .72d-6 ! 1950 F11 + other freons 36. PPMV(9) = 9.30d-6 ! 1950 F12 37. CALL GASES0 38. C CALL GASESY ! to be called for constant greenhouse gases 39. C**** Read in annual greenhouse gas concentrations for all years 40. READ (14,900) TITLE 41. WRITE (6,*) 'Greenhouse gas concentrations read in: ',TITLE 42. READ (14,900) 43. READ (14,900) 44. READ (14,900) 45. 41 READ (14,*,END=841) KYEAR,(GHGC(K,1950),K=1,6) 46. IF(KYEAR.lt.1950) GO TO 41 47. DO 42 IYEAR=1951,1990 48. READ (14,*,END=841) KYEAR,(GHGC(K,IYEAR),K=1,6) 48.1 42 IF(KYEAR.ne.IYEAR) GO TO 842 49. CLOSE (14) 50. C**** Set permanent aerosol optical thicknesses 51. FGOLDH(1) = .012d0 52. CALL AEROS0 53. C**** Read in monthly volcanic aerosol optical thicknesses 53.1 C 50 READ (15,END=850) TITLE 53.2 C READ (TITLE(69:72),905) KYEAR 53.3 C IF(KYEAR.lt.1879) GO TO 850 54. C DO 51 IYEAR=1880,1991 55. C READ (15,END=850) TITLE, IMX,JMX,LMX,NMX, 56. C * (VAOT(J,1,IYEAR),J=1,JM*12), RLON,RLAT,RLEV,RTIM, 57. C * LABLON,LABLAT,LABLEV,LABTIM, NASAGISS 58. C IF(NASAGISS.ne.'NASAGISS') GO TO 851 59. C READ (TITLE(69:72),903) KYEAR 60. C 51 IF(KYEAR.ne.IYEAR) GO TO 851 62. C CLOSE (15) 63. C**** Read in 2 tropospheric sulfate aerosol loading patterns 64. READ (16) TITLE,SAP1 65. WRITE (6,*) 'Sulfate aerosol pattern 1 read in: ',TITLE 66. C READ (16) TITLE,SAP2 67. C WRITE (6,*) 'Sulfate aerosol pattern 2 read in: ',TITLE 68. CLOSE (16) 68.1 DO 61 I=1,IM*JM 68.2 61 SAOT(I,1) = .7876 * .01806 * SAP1(I,1) 68.3 FGOLDH(4) = .7876 * .01806 ! 1950 tropospheric sulfates 69. C**** Read in global annual tropo sulfate aerosol optical thicknesses 70. C READ (17,900) TITLE 71. C WRITE (6,*) 'Sulfate aerosols optical thickness read in: ',TITLE 72. C READ (17,900) 73. C READ (17,900) 74. C READ (17,900) 75. C 71 READ (17,*,END=871) KYEAR,(SAOTGA(K,1990),K=1,2) 76. C IF(KYEAR.lt.1990) GO TO 71 77. C DO 72 IYEAR=1991,2100 78. C READ (17,*,END=871) KYEAR,(SAOTGA(K,IYEAR),K=1,2) 78.1 C 72 IF(KYEAR.ne.IYEAR) GO TO 872 79. C CLOSE (17) 80. C**** Define window flux to brightness temperature conversion factors 81. C WAVNA = 850.0 82. C WAVNB = 900.0 83. C DO 91 I=1,100 84. C PFWI = .001*I 85. C 91 TKPFW(I) = TKOFPF (WAVNA,WAVNB,PFWI) 86. C DO 92 I=101,190 87. C PFWI = .01*(I-100)+.1 88. C 92 TKPFW(I) = TKOFPF (WAVNA,WAVNB,PFWI) 89. C DO 93 I=191,630 90. C PFWI = .1*(I-190)+1. 91. C 93 TKPFW(I) = TKOFPF (WAVNA,WAVNB,PFWI) 92. RETURN 93. C**** 94. C**** 95. C**** 96. ENTRY RCOMPY 97. C**** 98. C**** RCOMPY is called once each year, trace gases may be changed 99. C**** 100. IF(JYEAR.lt.1950 .or. JYEAR.gt.2100) STOP 'RCOMPY: bad JYEAR' 101. C**** Change annual greenhouse gas concentrations (ppm) 102. IF(JYEAR.le.1990) then 103. PPMV(2) = GHGC(1,JYEAR) 104. PPMV(6) = GHGC(2,JYEAR)*1.e-3 105. PPMV(7) = GHGC(3,JYEAR)*1.e-3 106. PPMV(8) = GHGC(4,JYEAR)*1.e-6 ! what about GHGC(6) 107. PPMV(9) = GHGC(5,JYEAR)*1.e-6 107.1 else 107.2 PPMV(2) = GHGC(1,1990) * 1.005d0**(JYEAR-1990) 107.3 PPMV(6) = GHGC(2,1990)*1.e-3 107.4 PPMV(7) = GHGC(3,1990)*1.e-3 107.5 PPMV(8) = GHGC(4,1990)*1.e-6 ! what about GHGC(6) 107.6 PPMV(9) = GHGC(5,1990)*1.e-6 107.7 endif 108. CALL GASESY ! to be called for varying greenhouse gases 109. WRITE (6,910) JYEAR, PPMV(2),PPMV(6)*1.e3,PPMV(7)*1.e3, 109.1 * PPMV(8)*1.e6,PPMV(9)*1.e6 110. C**** Change annual tropospheric sulphate aerosol optical thicknesses 111. C DO 110 I=1,IM*JM 112. C 110 SAOT(I,1) = ((SAOTGA(1,JYEAR)-SAOTGA(1,1990))*SAP1(I,1) + 113. C + SAOTGA(2,JYEAR)*SAP2(I,1)) * .01806d0 114. RETURN 115. C**** 116. C**** 117. C**** 118. ENTRY RCOMPD 119. C**** 120. C**** RCOMPD is called once each day 121. C**** 122. C**** Change daily surface albedo parameters 123. CALL ALBEDD 124. C**** Change daily ozone concentrations 125. CALL OZONED 126. RETURN 127. C**** 128. C**** 129. C**** 130. ENTRY RCOMPJ 131. C**** 132. C**** RCOMPJ is called once each latitude 133. C**** 134. C**** Change volcanic aerosol optical thicknesses 135. C FGOLDH(1) = ??? 136. C**** Change latitudinal distribution of ozone concentration 137. CALL OZONEJ 138. RETURN 139. C**** 140. C**** 141. C**** 142. ENTRY RCOMP 143. C**** 144. C**** RCOMP is called once for each GCM grid box 145. C**** 146. CALL ALBED 147. CALL GASES 147.9 FGOLDH(4) = SAOT(ILON,JLAT) ! tropospheric sulphate aerosol 148. CALL AEROS 149. CALL THERML 150. CALL SOLAR 151. RETURN 152. C**** 153. 841 WRITE (6,*) 'Unexpected End-of-File on unit 14: IYEAR=',IYEAR 155. STOP 'RCOMP0 841' 155.1 842 WRITE (6,*) 'Inconsistant years on unit 14: IYEAR=',IYEAR,KYEAR 155.2 STOP 'RCOMP0 842' 156. 850 WRITE (6,*) 'Unexpected End-of-File on unit 15: IYEAR=',IYEAR 158. STOP 'RCOMP0 850' 159. 851 WRITE (6,*) 'Inconsistant years on unit 15: IYEAR=',IYEAR,KYEAR 161. STOP 'RCOMP0 851' 162. 871 WRITE (6,*) 'Unexpected End-of-File on unit 17: IYEAR=',IYEAR 164. STOP 'RCOMP0 871' 164.1 872 WRITE (6,*) 'Inconsistant years on unit 17: IYEAR=',IYEAR,KYEAR 164.3 STOP 'RCOMP0 872' 165. C**** 166. 900 FORMAT (A80) 167. 910 FORMAT ('0Greenhouse gases: JYEAR CO2*E6 N2O*E9 CH4*E9 ', 168. * 'F11*E12 F12*E12' / I24,5F9.2) 169. END 200. 201. BLOCK DATA RADNBD 202. INCLUDE 'RA99M9.COM' 203. C**** 204. C**** Control and input parameters 205. C**** 206. DATA VADATA / 207. C**** DESRT TNDRA GRASS SHRUB TREES DECID EVERG RAINF ROCKS CROPS ALGAE 208. C**** Visual albedoes for 4 seasons 209. W .350, .067, .089, .089, .078, .100, .067, .061, .100, .070, .001, 210. S .350, .063, .100, .100, .073, .055, .067, .061, .100, .070, .001, 211. S .350, .085, .091, .139, .085, .058, .083, .061, .100, .070, .001, 212. A .350, .080, .090, .111, .064, .055, .061, .061, .100, .070, .001, 213. C**** Near-IR albedoes for 4 seasons 214. W .350, .200, .267, .267, .233, .300, .200, .183, .100, .070, .001, 215. S .350, .206, .350, .300, .241, .218, .200, .183, .100, .070, .001, 216. S .350, .298, .364, .417, .298, .288, .250, .183, .100, .070, .001, 217. A .350, .255, .315, .333, .204, .218, .183, .183, .100, .070, .001, 218. C**** Water field capacities (kg/m**2) for 3 layers 219. 1 10., 30., 30., 30., 30., 30., 30., 200., 10., 30., 999., 220. 2 10., 200., 200., 300., 300., 450., 450., 450., 10., 200., 999., 221. 3 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 10., 222. C**** Masking depths (kg/m**2) 223. 4 10., 20., 20., 50., 200., 500.,1000.,2500., 10., 30., .001/ 224. C**** 225. DATA FRACSL,TKCICE/ 226. * 1.E-2, 258./ 227. C**** GLOBAL OCEAN LAND DESERT HAZE TR1 TR2 TR3 TR4 228. DATA FGOLDH/.012, .68, .32, 0., 0., 0., 0., 0., 0. 229. + ,.012, .68, .32, 0., 0., 0., 0., 0., 0./ 230. DATA ID5, ITR/ 231. * 5*0, 4*0/ 232. DATA KFORCE/ 233. * -123456789/ 234. C**** 235. C**** Basic RADCOM input data 236. C**** 237. DATA PLE /1013.2500, 961.7485, 879.3460, 741.3219, 566.2166, 238. + 401.4117, 262.3575, 154.2043, 71.8018, 10.0000, 239. + 5.0000, 2.0000, 1.E-05, 27*0.00/ 240. DATA HLE /1.E-10, 1., 2., 3., 6., 11., 20., 32., 47., 51., 241. * 71., 84.852, 99.99, 27*99.999/ 242. DATA TLB, TLT, TLM, U0GAS/ 243. * 40*250.,40*250.,40*250.,360*0./ 244. DATA ULGAS,TRACER,CLDTAU, QL/ 245. * 360*0.,160*0., 40*0., 40*0./ 246. DATA POCEAN,PEARTH, POICE, PLICE, AGESN, SNOWE,SNOWOI,SNOWLI/ 247. * .7, .1, .1, .1, 1., .3, .1, .2/ 248. DATA TGO, TGE, TGOI, TGLI, TS, WS,WEARTH/ 249. * 288.15,288.15,288.15,288.15,288.15, 2., 0./ 250. DATA S0, COSZ, PVT , BXA,SRBXAL/ 251. * 1367., .5,8*.125, 3*0.,153*0., 30*0./ 252. DATA JYEAR, JDAY, JLAT, ILON/ 253. * 1958, 0, 18, 18/ 254. C**** 255. C**** BLOCKD initialized default data 256. C**** 257. C AEROSOL RADIATIVE PROPERTIES,COMPOSITION,TYPE & VERTICAL DISTRIBUTION 258. C**** 259. DIMENSION QACID1(25),QACID2(25),QSLFT1(25),QSLFT2(25) 260. T ,QBSLT1(25),QBSLT2(25),QSSALT(25),QDUST1(25) 261. T ,QDUST2(25),QCARB1(25),QCARB2(25) 262. T ,SACID1(25),SACID2(25),SSLFT1(25),SSLFT2(25) 263. T ,SBSLT1(25),SBSLT2(25),SSSALT(25),SDUST1(25) 264. T ,SDUST2(25),SCARB1(25),SCARB2(25) 265. T ,CACID1(25),CACID2(25),CSLFT1(25),CSLFT2(25) 266. T ,CBSLT1(25),CBSLT2(25),CSSALT(25),CDUST1(25) 267. T ,CDUST2(25),CCARB1(25),CCARB2(25) 268. T ,QWATER(25),QICE25(25),SWATER(25),SICE25(25) 269. T ,CWATER(25),CICE25(25) 270. S ,XACID1(6),XACID2(6),XSLFT1(6),XSLFT2(6),XBSLT1(6),XBSLT2(6) 271. S ,XSSALT(6),XDUST1(6),XDUST2(6),XCARB1(6),XCARB2(6) 272. S ,YACID1(6),YACID2(6),YSLFT1(6),YSLFT2(6),YBSLT1(6),YBSLT2(6) 273. S ,YSSALT(6),YDUST1(6),YDUST2(6),YCARB1(6),YCARB2(6) 274. S ,ZACID1(6),ZACID2(6),ZSLFT1(6),ZSLFT2(6),ZBSLT1(6),ZBSLT2(6) 275. S ,ZSSALT(6),ZDUST1(6),ZDUST2(6),ZCARB1(6),ZCARB2(6) 276. S ,XWATER(6),XICE25(6),YWATER(6),YICE25(6),ZWATER(6),ZICE25(6) 277. C**** 278. EQUIVALENCE (TRAQEX(1, 1),QACID1),(TRAQEX(1, 2),QACID2) 279. 1 ,(TRAQEX(1, 3),QSLFT1),(TRAQEX(1, 4),QSLFT2) 280. 2 ,(TRAQEX(1, 5),QBSLT1),(TRAQEX(1, 6),QBSLT2) 281. 3 ,(TRAQEX(1, 7),QSSALT),(TRAQEX(1, 8),QDUST1) 282. 4 ,(TRAQEX(1, 9),QDUST2),(TRAQEX(1,10),QCARB1) 283. 5 ,(TRAQEX(1,11),QCARB2) 284. EQUIVALENCE (TRAQSC(1, 1),SACID1),(TRAQSC(1, 2),SACID2) 285. 1 ,(TRAQSC(1, 3),SSLFT1),(TRAQSC(1, 4),SSLFT2) 286. 2 ,(TRAQSC(1, 5),SBSLT1),(TRAQSC(1, 6),SBSLT2) 287. 3 ,(TRAQSC(1, 7),SSSALT),(TRAQSC(1, 8),SDUST1) 288. 4 ,(TRAQSC(1, 9),SDUST2),(TRAQSC(1,10),SCARB1) 289. 5 ,(TRAQSC(1,11),SCARB2) 290. EQUIVALENCE (TRACOS(1, 1),CACID1),(TRACOS(1, 2),CACID2) 291. 1 ,(TRACOS(1, 3),CSLFT1),(TRACOS(1, 4),CSLFT2) 292. 2 ,(TRACOS(1, 5),CBSLT1),(TRACOS(1, 6),CBSLT2) 293. 3 ,(TRACOS(1, 7),CSSALT),(TRACOS(1, 8),CDUST1) 294. 4 ,(TRACOS(1, 9),CDUST2),(TRACOS(1,10),CCARB1) 295. 5 ,(TRACOS(1,11),CCARB2) 296. EQUIVALENCE (TRCQEX(1, 1),QWATER),(TRCQEX(1, 2),QICE25) 297. EQUIVALENCE (TRCQSC(1, 1),SWATER),(TRCQSC(1, 2),SICE25) 298. EQUIVALENCE (TRCCOS(1, 1),CWATER),(TRCCOS(1, 2),CICE25) 299. C**** 300. EQUIVALENCE (SRAQEX(1, 1),XACID1),(SRAQEX(1, 2),XACID2) 301. 1 ,(SRAQEX(1, 3),XSLFT1),(SRAQEX(1, 4),XSLFT2) 302. 2 ,(SRAQEX(1, 5),XBSLT1),(SRAQEX(1, 6),XBSLT2) 303. 3 ,(SRAQEX(1, 7),XSSALT),(SRAQEX(1, 8),XDUST1) 304. 4 ,(SRAQEX(1, 9),XDUST2),(SRAQEX(1,10),XCARB1) 305. 5 ,(SRAQEX(1,11),XCARB2) 306. EQUIVALENCE (SRAQSC(1, 1),YACID1),(SRAQSC(1, 2),YACID2) 307. 1 ,(SRAQSC(1, 3),YSLFT1),(SRAQSC(1, 4),YSLFT2) 308. 2 ,(SRAQSC(1, 5),YBSLT1),(SRAQSC(1, 6),YBSLT2) 309. 3 ,(SRAQSC(1, 7),YSSALT),(SRAQSC(1, 8),YDUST1) 310. 4 ,(SRAQSC(1, 9),YDUST2),(SRAQSC(1,10),YCARB1) 311. 5 ,(SRAQSC(1,11),YCARB2) 312. EQUIVALENCE (SRACOS(1, 1),ZACID1),(SRACOS(1, 2),ZACID2) 313. 1 ,(SRACOS(1, 3),ZSLFT1),(SRACOS(1, 4),ZSLFT2) 314. 2 ,(SRACOS(1, 5),ZBSLT1),(SRACOS(1, 6),ZBSLT2) 315. 3 ,(SRACOS(1, 7),ZSSALT),(SRACOS(1, 8),ZDUST1) 316. 4 ,(SRACOS(1, 9),ZDUST2),(SRACOS(1,10),ZCARB1) 317. 5 ,(SRACOS(1,11),ZCARB2) 318. EQUIVALENCE (SRCQEX(1, 1),XWATER),(SRCQEX(1, 2),XICE25) 319. EQUIVALENCE (SRCQSC(1, 1),YWATER),(SRCQSC(1, 2),YICE25) 320. EQUIVALENCE (SRCCOS(1, 1),ZWATER),(SRCCOS(1, 2),ZICE25) 321. C**** 322. DATA NGOLDH/5/,NAERO/11/ 323. C 324. C----------------------------------------------------------------------- 325. C COMPOSITION & VERTICAL DISTRIBUTION FOR 5 SPECIFIED AEROSOL TYPES 326. C----------------------------------------------------------------------- 327. C TYPE 328. C 1 STRATOSPHERIC GLOBAL AEROSOL A,B,C ARE GLOBAL AVERAGE VALUES 329. C 2 TROPOSPHERIC OCEAN AEROSOL A,B,C ARE GLOBAL AVERAGE VALUES 330. C 3 TROPOSPHERIC LAND AEROSOL A,B,C ARE GLOBAL AVERAGE VALUES 331. C 4 Tropospheric sulphate aerosols 332. C 5 Sea salt aerosols 333. C 334. C 1 2 3 4 5 6 7 8 9 10 11 335. C ACID1 OCT82 SLFT1 SLFT2 BSLT1 BSLT2 SSALT DUST1 DUST2 MAY82 CARB2 336. DATA AGOLDH/ 337. 11.000, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 338. 2 0., 0., .018, .010, .012, .023, 0., 0., 0., 0., 0., 339. 3 0., 0., .031, .010, .021, .042, 0., 0., 0., 0., .018, 340. 4 0., 0., 0.,1.000, 0., 0., 0., 0., 0., 0., 0., 341. 5 0., 0., 0., 0., 0., 0.,1.000, 0., 0., 0., 0./ 342. DATA BGOLDH/ 343. 1 20., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 344. 2 0., 0., 4., 0., 4., 1., 0., 0., 0., 0., 0., 345. 3 0., 0., 4., 0., 4., 0., 0., 0., 0., 0., 0., 346. 4 0., 0., 0., 0., 0., 0., 0., 3.5, 0., 0., 0., 347. 5 0., 0.00, 0., 0., 0., 0., 0., 3.5, 0., 0., 0./ 348. DATA CGOLDH/ 349. 1 3., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 350. 2 0., 0., 3., 1., 3., .5, 1., 0., 0., 0., 0., 351. 3 0., 0., 3., 1., 3., 1., 0., 0., 0., 0., 1., 352. 4 0., 0., 0., 0., 0., 0., 0., 1., 1., 0., 0., 353. 5 0., 1., 0., 0., 0., 0., 0., 1., 0., 0., 0./ 354. C 355. C----------------------------------------------------------------------- 356. C THERMAL RADIATION 25 K-INTERVAL MERGED AEROSOL DATA FOR QEXT,QSCA,COSB 357. C----------------------------------------------------------------------- 358. DATA QACID1/ 359. + .04052, .05895, .08506, .06673, .05160, .04437, .03864, 360. + .02719, .01668, .01146, .00705, .03286, .02449, .03017, 361. + .03198, .02891, .02634, .02366, .02300, .02271, .02159, 362. + .08516, .08825, .08982, .09284/ 363. DATA SACID1/ 364. + .00095, .00361, .00273, .00226, .00150, .00141, .00131, 365. + .00090, .00049, .00029, .00014, .00072, .00049, .00031, 366. + .00023, .00023, .00022, .00020, .00019, .00018, .00018, 367. + .00183, .00201, .00205, .00207/ 368. DATA CACID1/ 369. + .11030, .17256, .17138, .19696, .19510, .18945, .18874, 370. + .18795, .18313, .17814, .17075, .10583, .09756, .08388, 371. + .07246, .07266, .07099, .06873, .06754, .06661, .06674, 372. + .11197, .11068, .10998, .10852/ 373. C 374. DATA QACID2/ 375. + .05764, .15189, .06264, .04527, .03973, .03646, .03375, 376. + .02163, .01337, .00979, .00724, .04076, .03631, .04273, 377. + .04072, .03752, .03290, .03012, .02968, .02914, .02763, 378. + .10731, .12510, .12901, .13232/ 379. DATA SACID2/ 380. + .00367, .00752, .00264, .00172, .00188, .00221, .00225, 381. + .00134, .00066, .00034, .00012, .00237, .00121, .00084, 382. + .00080, .00081, .00074, .00069, .00067, .00065, .00064, 383. + .00674, .00807, .00825, .00837/ 384. DATA CACID2/ 385. + .05720, .11171, .11850, .11443, .12325, .13171, .13500, 386. + .13575, .13419, .12666, .10961, .05186, .04026, .03219, 387. + .03060, .03105, .03041, .02959, .02911, .02884, .02901, 388. + .07145, .07168, .07134, .07096/ 389. C 390. DATA QSLFT1/ 391. + .15555, .16333, .16406, .16396, .16070, .14074, .11920, 392. + .09140, .07341, .06645, .05871, .15301, .13456, .15809, 393. + .16264, .14805, .12798, .10588, .09960, .09604, .08844, 394. + .35895, .27430, .26964, .27183/ 395. DATA SSLFT1/ 396. + .13162, .13152, .11642, .12932, .10550, .08323, .07081, 397. + .05079, .03287, .02458, .01871, .12787, .11183, .09490, 398. + .08739, .08716, .08022, .07182, .06899, .06700, .06496, 399. + .13067, .12933, .12878, .12808/ 400. DATA CSLFT1/ 401. + .52508, .48102, .59654, .66259, .66566, .70224, .71546, 402. + .69308, .62819, .55963, .45811, .52840, .54500, .51620, 403. + .50685, .52475, .54985, .58351, .59484, .60203, .61652, 404. + .45926, .47060, .47243, .47178/ 405. C 406. DATA QSLFT2/ 407. + .44109, .37065, .38095, .40554, .37738, .32564, .27970, 408. + .21687, .17752, .16154, .14952, .43239, .38517, .39512, 409. + .39098, .36978, .32960, .28406, .27042, .26204, .24771, 410. + .63665, .59084, .58844, .59078/ 411. DATA SSLFT2/ 412. + .37818, .31549, .29505, .33810, .28074, .22692, .19562, 413. + .14289, .09653, .07449, .06008, .36685, .33089, .28296, 414. + .26185, .26286, .24369, .22019, .21220, .20647, .20093, 415. + .31870, .30963, .30762, .30507/ 416. DATA CSLFT2/ 417. + .54586, .50074, .62826, .69007, .69596, .73443, .74600, 418. + .71846, .64430, .57291, .47311, .54977, .56612, .53939, 419. + .53105, .54799, .57221, .60426, .61497, .62179, .63518, 420. + .51454, .52095, .52268, .52316/ 421. C 422. DATA QBSLT1/ 423. + .19787, .15206, .14808, .15505, .14132, .12508, .10931, 424. + .07946, .05659, .04675, .03801, .20081, .15823, .15732, 425. + .15377, .14273, .13163, .12005, .11684, .11523, .11121, 426. + .36601, .39099, .39240, .39274/ 427. DATA SBSLT1/ 428. + .09892, .12369, .09780, .11017, .08914, .08577, .07794, 429. + .05688, .03912, .03069, .02440, .09492, .08277, .05817, 430. + .04773, .04970, .04568, .04058, .03865, .03717, .03641, 431. + .07710, .08232, .08235, .08163/ 432. DATA CBSLT1/ 433. + .54090, .49369, .59375, .67539, .69444, .71623, .71674, 434. + .69425, .63125, .57379, .48766, .54072, .57272, .57215, 435. + .57655, .59243, .60616, .62323, .62911, .63253, .63934, 436. + .51632, .50380, .50414, .50666/ 437. C 438. DATA QBSLT2/ 439. + .49004, .35700, .34009, .38146, .35476, .32874, .29258, 440. + .21726, .16067, .13571, .11451, .48169, .40550, .37263, 441. + .35312, .33842, .31466, .28850, .28051, .27574, .26813, 442. + .59495, .63654, .63850, .63742/ 443. DATA SBSLT2/ 444. + .26833, .30862, .25309, .29334, .24644, .24238, .22164, 445. + .16459, .11742, .09480, .07809, .26006, .23936, .17265, 446. + .14418, .15103, .13960, .12488, .11925, .11488, .11275, 447. + .19766, .20963, .20969, .20807/ 448. DATA CBSLT2/ 449. + .57850, .51330, .62334, .70306, .72063, .74166, .74111, 450. + .71466, .64442, .58410, .49911, .58174, .60690, .60535, 451. + .60954, .62353, .63716, .65423, .66019, .66381, .67030, 452. + .58670, .57707, .57759, .58014/ 453. C 454. DATA QSSALT/ 455. + .27651, .36950, .40122, .39669, .34286, .33458, .29978, 456. + .26075, .26470, .26660, .28507, .27114, .23752, .18761, 457. + .16890, .17532, .17705, .17827, .17801, .17743, .17914, 458. + .34241, .33620, .33607, .33681/ 459. DATA SSSALT/ 460. + .27651, .36950, .40121, .39659, .34226, .33245, .29555, 461. + .22360, .16290, .13425, .11177, .27114, .23751, .18755, 462. + .16883, .17526, .17700, .17823, .17797, .17739, .17911, 463. + .34241, .33620, .33607, .33681/ 464. DATA CSSALT/ 465. + .66858, .50298, .60372, .65282, .66694, .67041, .66666, 466. + .62258, .52248, .44732, .32878, .66866, .66680, .66404, 467. + .66252, .66281, .66265, .66244, .66232, .66223, .66226, 468. + .67338, .67406, .67410, .67408/ 469. C 470. DATA QDUST1/ 471. + .60958, .65996, .59890, .73030, .64827, .55835, .48157, 472. + .34847, .23144, .18097, .13460, .59012, .47533, .39938, 473. + .36575, .35808, .33834, .31587, .30849, .30369, .29821, 474. + .91360,1.14613,1.16193,1.16619/ 475. DATA SDUST1/ 476. + .32015, .60541, .49800, .59591, .46651, .39745, .34242, 477. + .23468, .13039, .08473, .04350, .29084, .23940, .16410, 478. + .13070, .13267, .12095, .10691, .10167, .09788, .09578, 479. + .39128, .54469, .55555, .55942/ 480. DATA CDUST1/ 481. + .50425, .49645, .57736, .63615, .63373, .66224, .67205, 482. + .67034, .65137, .61767, .53600, .49640, .47921, .43825, 483. + .40760, .41364, .41120, .40706, .40418, .40149, .40315, 484. + .47280, .39308, .38801, .38670/ 485. C 486. DATA QDUST2/ 487. + .95483, .71515, .77676, .91847, .93699, .89565, .82979, 488. + .74871, .70959, .69272, .68748, .94632, .90846, .85600, 489. + .83350, .83544, .82317, .80807, .80270, .79879, .79577, 490. + 1.02427,1.12417,1.13054,1.13169/ 491. DATA SDUST2/ 492. + .49885, .58157, .55165, .64038, .59140, .55222, .50136, 493. + .42019, .36087, .33502, .31667, .49026, .47989, .42207, 494. + .39751, .40487, .39774, .38819, .38426, .38107, .38027, 495. + .49780, .59147, .59817, .60013/ 496. DATA CDUST2/ 497. + .74352, .54594, .68229, .72513, .73598, .75710, .75041, 498. + .70723, .65024, .61702, .58021, .74556, .74741, .75647, 499. + .76384, .76647, .77599, .78746, .79136, .79400, .79700, 500. + .71874, .62817, .62224, .62062/ 501. C 502. DATA QCARB1/ 503. + .44718, .51882, .26055, .20526, .19295, .18655, .17520, 504. + .11120, .06749, .04893, .03537, .32912, .25261, .24973, 505. + .23947, .22883, .20424, .18781, .18400, .18032, .17370, 506. + .57200, .64430, .65267, .65790/ 507. DATA SCARB1/ 508. + .17857, .12659, .06506, .05088, .05317, .05712, .05562, 509. + .03310, .01705, .01009, .00493, .13908, .08683, .06332, 510. + .06114, .06260, .05755, .05319, .05155, .05032, .04981, 511. + .19594, .21003, .20967, .20853/ 512. DATA CCARB1/ 513. + .40490, .48729, .43960, .40824, .46236, .51422, .53366, 514. + .53211, .51283, .46211, .32882, .40923, .35984, .30817, 515. + .30468, .31306, .31215, .30857, .30555, .30388, .30644, 516. + .43102, .40748, .40436, .40208/ 517. C 518. DATA QCARB2/ 519. + .09591, .22971, .21603, .21745, .17928, .17061, .15202, 520. + .10846, .06721, .04817, .03076, .09456, .08428, .07093, 521. + .06589, .06737, .06766, .06782, .06771, .06754, .06792, 522. + .12455, .12130, .12121, .12155/ 523. DATA SCARB2/ 524. + .00748, .06133, .05031, .04978, .03714, .03448, .03065, 525. + .02099, .01137, .00688, .00291, .00728, .00544, .00350, 526. + .00276, .00291, .00290, .00288, .00285, .00282, .00286, 527. + .01420, .01327, .01324, .01332/ 528. DATA CCARB2/ 529. + .14117, .25269, .27090, .30506, .29845, .28974, .28880, 530. + .28843, .28603, .28395, .29112, .14128, .12741, .11121, 531. + .09892, .09935, .09786, .09604, .09517, .09448, .09466, 532. + .18297, .17686, .17658, .17696/ 533. C 534. DATA QWATER/ 535. + .82334, .89509,1.13254,1.20762,1.24075,1.18580,1.07585, 536. + .95283, .89542, .86914, .85864, .87834, .94021,1.03878, 537. + 1.07876,1.06927,1.06987,1.07153,1.07327,1.07505,1.07280, 538. + 1.20709,1.20194,1.20383,1.20978/ 539. DATA SWATER/ 540. + .34695, .68566, .86748, .89010, .83121, .75556, .65338, 541. + .51441, .40925, .36469, .31873, .39396, .39368, .43707, 542. + .45625, .44997, .45039, .45146, .45251, .45357, .45227, 543. + .85537, .85478, .85718, .86370/ 544. DATA CWATER/ 545. + .91848, .65450, .79206, .82335, .83709, .84869, .84338, 546. + .77907, .68419, .62521, .54076, .91355, .89224, .85667, 547. + .84557, .85029, .85229, .85399, .85411, .85389, .85524, 548. + .91095, .91472, .91488, .91467/ 549. C 550. DATA QICE25/ 551. + 1.15210, .81551, .98885,1.10325,1.17652,1.14217,1.07777, 552. + 1.08252,1.14496,1.16939,1.22006,1.16194,1.16781,1.19342, 553. + 1.20279,1.19736,1.19435,1.19146,1.19097,1.19095,1.18924, 554. + 1.19321,1.21794,1.21959,1.21942/ 555. DATA SICE25/ 556. + .57392, .45452, .57278, .68806, .74580, .69171, .64662, 557. + .62884, .64120, .64892, .66105, .59403, .60241, .67853, 558. + .70399, .68299, .66547, .64731, .64301, .64122, .63321, 559. + .71867, .77122, .77524, .77622/ 560. DATA CICE25/ 561. + .93634, .72920, .86084, .88431, .87489, .88472, .86613, 562. + .82078, .79850, .79041, .78539, .93377, .91036, .85751, 563. + .84228, .85220, .86089, .87036, .87263, .87355, .87810, 564. + .94697, .94840, .94812, .94714/ 565. C 566. C----------------------------------------------------------------------- 567. C SOLAR RADIATION 6 K-INTERVAL MERGED AEROSOL DATA FOR QEXT,QSCA,COSB 568. C----------------------------------------------------------------------- 569. C 570. DATA XACID1 /.05776, .10033, .19099, .36614, .55931,1.04703/ 571. DATA YACID1 /.01880, .09956, .19090, .36613, .55931,1.04703/ 572. DATA ZACID1 /.36054, .51871, .57276, .62068, .65273, .68988/ 573. C 574. DATA XACID2 /.13360, .33875, .51498, .68359, .79939, .94494/ 575. DATA YACID2 /.07420, .33691, .51483, .68358, .79939, .94494/ 576. C$ DATA ZACID2 /.40248, .62259, .68524, .71328, .71195, .72894/ 577. DATA ZACID2 /.39821, .54835, .60846, .63637, .63503, .65221/ 578. C 579. DATA XSLFT1 /.31035, .44757, .54238, .66756, .78260,1.04454/ 580. DATA YSLFT1 /.24589, .44490, .54224, .66755, .78260,1.04454/ 581. DATA ZSLFT1 /.70591, .67557, .66832, .66438, .66199, .66008/ 582. C 583. DATA XSLFT2 /.60959, .74888, .81124, .87560, .92632,1.00936/ 584. DATA YSLFT2 /.50477, .74262, .81090, .87556, .92631,1.00935/ 585. DATA ZSLFT2 /.74067, .70281, .69748, .69922, .70070, .70754/ 586. C 587. DATA XBSLT1 /.30419, .46195, .54908, .66403, .77732,1.02644/ 588. DATA YBSLT1 /.28732, .44765, .53358, .64786, .76063,1.00769/ 589. DATA ZBSLT1 /.67768, .66588, .66785, .66932, .66671, .66818/ 590. C 591. DATA XBSLT2 /.62145, .76377, .81783, .87743, .92782,1.00765/ 592. DATA YBSLT2 /.58466, .73120, .78367, .84258, .89259, .96944/ 593. DATA ZBSLT2 /.70368, .69767, .70313, .70847, .70983, .71935/ 594. C 595. DATA XSSALT /.64091, .78294, .83066, .87490, .92554,1.00414/ 596. DATA YSSALT /.64091, .78294, .83066, .87490, .92554,1.00413/ 597. DATA ZSSALT /.67233, .68272, .68718, .69084, .69334, .69627/ 598. C 599. DATA XDUST1/1.17571,1.20282,1.13894,1.08190,1.04572, .99864/ 600. DATA YDUST1/1.04642,1.12320,1.04442, .97057, .93288, .78720/ 601. DATA ZDUST1 /.72235, .68164, .69516, .72361, .74315, .80409/ 602. C 603. DATA XDUST2/1.09335,1.12888,1.09512,1.05217,1.02411,1.00081/ 604. DATA YDUST2 /.83740, .93590, .88162, .81721, .78602, .68767/ 605. DATA ZDUST2 /.78776, .76447, .77511, .79364, .80840, .85594/ 606. C 607. DATA XCARB1 /.74444,1.11851,1.14599,1.09902,1.05179,1.00292/ 608. DATA YCARB1 /.53412,1.11290,1.14544,1.09899,1.05179,1.00292/ 609. C$ DATA ZCARB1 /.75767, .74553, .72950, .71977, .71968, .74073/ 610. DATA ZCARB1 /.71248, .66984, .65284, .64292, .64282, .66426/ 611. C 612. DATA XCARB2 /.54418, .82500, .91922, .97919,1.00345, .99476/ 613. DATA YCARB2 /.19636, .34820, .40558, .44719, .46860, .48132/ 614. DATA ZCARB2 /.45878, .59691, .65112, .70444, .74341, .79820/ 615. C 616. DATA XWATER/1.10372,1.05381,1.03792,1.02265,1.01285, .99989/ 617. DATA YWATER /.84758,1.03190,1.02896,1.02226,1.01282, .99988/ 618. DATA ZWATER /.87621, .84587, .84884, .85323, .85888, .86321/ 619. C 620. DATA XICE25/1.05394,1.02884,1.02030,1.01257,1.00706, .99981/ 621. DATA YICE25 /.75677, .96035,1.00797,1.01184,1.00702, .99981/ 622. DATA ZICE25 /.92708, .88645, .87975, .87906, .87391, .87623/ 623. C 624. C----------------------------------------------------------------------- 625. C THERMAL RADIATION 25 K-INTERVAL MERGED CLOUD & SURFACE ALBEDO DATA 626. C----------------------------------------------------------------------- 627. DATA AGSIDV / 628. S .01407, .01653, .03230, .08764, .10055, .09095, .08892, 629. S .07985, .06411, .05926, .05398, .01576, .02449, .05091, 630. S .05680, .05325, .04652, .03809, .03574, .03451, .03082, 631. S .01757, .02022, .02059, .02082, 632. I .01407, .01653, .03230, .08764, .10055, .09095, .08892, 633. I .07985, .06411, .05926, .05398, .01576, .02449, .05091, 634. I .05680, .05325, .04652, .03809, .03574, .03451, .03082, 635. I .01757, .02022, .02059, .02082, 636. D .04500, .10414, .06739, .08448, .08516, .06283, .05230, 637. D .03382, .01901, .01542, .01178, .05142, .04835, .05505, 638. D .05600, .05310, .04603, .03731, .03472, .03328, .03000, 639. D .16159, .17592, .17812, .17927, 640. V 25*0.0/ 641. DATA AOCEAN / 642. + .04000, .05965, .06124, .08339, .09235, .09510, .09908, 643. + .11117, .12263, .12577, .12931, .04700, .06894, .08970, 644. + .09574, .09565, .09619, .09672, .09703, .09723, .09700, 645. + .04645, .04487, .04482, .04493/ 646. C 647. DATA CLDALB / 648. + .01332, .08190, .07036, .05082, .04486, .04673, .04770, 649. + .05130, .05240, .05251, .05259, .01558, .01763, .02410, 650. + .02571, .02514, .02448, .02366, .02347, .02340, .02294, 651. + .04566, .04499, .04518, .04544, 652. + .01407, .01653, .03230, .08764, .10055, .09095, .08892, 653. + .07985, .06411, .05926, .05398, .01576, .02449, .05091, 654. + .05680, .05325, .04652, .03809, .03574, .03451, .03082, 655. + .01757, .02022, .02059, .02082/ 656. C 657. DATA ASNALB /.60, .35, 13*0.0/ 658. DATA AOIALB /.55, .30, 13*0.0/ 659. DATA ALIALB /.60, .35, 13*0.0/ 660. C 661. C----------------------------------------------------------------------- 662. C TRACE GAS VERTICAL DISTRIBUTION & 1958 MEAN CONCENTRATION 663. C----------------------------------------------------------------------- 664. C 665. DATA CMANO2/ 666. 1 8.66E-06,5.15E-06,2.85E-06,1.50E-06,9.89E-07,6.91E-07,7.17E-07, 667. 2 8.96E-07,3.67E-06,4.85E-06,5.82E-06,6.72E-06,7.77E-06,8.63E-06, 668. 3 8.77E-06,8.14E-06,6.91E-06,5.45E-06,4.00E-06,2.67E-06,1.60E-06, 669. 4 8.36E-07,3.81E-07,1.58E-07,6.35E-08,2.57E-08,1.03E-08,4.18E-09, 670. 5 1.66E-09,6.57E-10,2.58E-10,1.02E-10,4.11E-11,1.71E-11,7.73E-12, 671. 6 9.07E-12,4.63E-12,2.66E-12,1.73E-12,1.28E-12,1.02E-12,1.00E-30/ 672. C 673. C 674. C GAS 1 2 3 4 5 6 7 8 9 675. C H2O CO2 O3 O2 NO2 N2O CH4 CCL3F1 CCL2F2 676. DATA PPMV/0.,315., 0.,210000., 0., .295, 1.4, 8.E-6, 25.E-6/ 677. C 678. C----------------------------------------------------------------------- 679. C TRACE GAS ABSORPTION COEFFICIENTS FOR F11 & F12 680. C----------------------------------------------------------------------- 681. C 682. DIMENSION F11PCM(25),F12PCM(25) 683. EQUIVALENCE (TRACEG(1,1),F11PCM),(TRACEG(1,2),F12PCM) 684. DATA F11PCM/ 685. + 13.6000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 686. + 0.0000, 0.0000, 0.0000, 0.0000,11.9504, 2.5138, 0.5054, 687. + 0.1086, 0.0308, 0.0178, 0.0054, 0.0000, 0.0000, 0.0000, 688. + 2.5220, 1.1731, 0.8627, 0.7445/ 689. DATA F12PCM/ 690. + 5.4900, 1.3339, 0.7739, 0.1304, 0.0286, 0.0051, 0.0000, 691. + 0.0000, 0.0000, 0.0000, 0.0000, 9.0745, 2.3577, 0.4135, 692. + 0.0575, 0.0000, 0.2507, 0.6215, 0.7262, 0.7972, 0.9150, 693. + 13.1663, 1.1564, 0.0388, 0.0082/ 694. END 1500. 1501. SUBROUTINE PHDATM (P,H,D,T,O,Q,S,OCM,WCM,NPHD,NATM) 1502. C 1503. C ------------------------------------------------------------------ 1504. C ------------- MCCLATCHY (1972) ATMOSPHERE DATA ----------- 1505. C ------------------------------------------------------------------ 1506. C 1507. C INPUT DATA 1508. C------------------ 1509. C NATM=0 GIVES ABREVIATED DATA FOR STANDARD ATMOSPHER 1510. C (INPUT: P OR H) (RETURNS: H OR P & D,T) 1511. C 1512. C NATM=1 GIVES ATMOSPHERE DATA FOR TROPICAL LATITUDES 1513. C NATM=2 GIVES ATMOSPHERE DATA FOR MIDLATITUDE SUMMER 1514. C NATM=3 GIVES ATMOSPHERE DATA FOR MIDLATITUDE WINTER 1515. C NATM=4 GIVES ATMOSPHERE DATA FOR SUBARCTIC SUMMER 1516. C NATM=5 GIVES ATMOSPHERE DATA FOR SUBARCTIC WINTER 1517. C NATM=6 GIVES ATMOSPHERE DATA FOR STANDARD ATMOSPHER 1518. C 1519. C NPHD=1 RETURNS H,D,T,O,Q,S DATA FOR GIVEN PRESSURE P 1520. C NPHD=2 RETURNS P,D,T,O,Q,S DATA FOR GIVEN HEIGHT H 1521. C NPHD=3 RETURNS P,H,T,O,Q,S DATA FOR GIVEN DENSITY D 1522. C 1523. C OUTPUT DATA 1524. C------------------ 1525. C P = PRESSURE IN MILLIBARS 1526. C H = HEIGHT IN KILOMETERS 1527. C D = DENSITY IN GRAMS/METER**3 1528. C T = TEMPERATURE (ABSOLUTE) 1529. C O = OZONE MIXING RATIO (GRAMS OZONE)/(GRAMS AIR) 1530. C Q = SPECIFIC HUMIDITY (GRAMS WATER VAPOR)/(GRAMS AIR) 1531. C S = SATURATION RATIO (GRAMS WATER VAPOR)/(GRAMS AIR) 1532. C OCM = OZONE (CM-STP) ABOVE GIVEN HEIGHT 1533. C WCM = WATER VAPOR (CM-STP) ABOVE GIVEN HEIGHT 1534. C 1535. C REMARKS 1536. C------------------ 1537. C INPUT P,H,D PARAMETERS ARE NOT ALTERED 1538. C P,D INTERPOLATION IS EXPONENTIAL WITH HEIGHT 1539. C NO EXTRAPOLATION IS MADE OUTSIDE 0-100 KM INTERVAL 1540. C S IS NOT COMPUTED ABOVE 40 KM (FORMULA NOT ACCURATE) 1541. C 1542. C R = Q/S GIVES RELATIVE HUMIDITY 1543. C W = Q/(1-Q) GIVES WATER VAPOR MIXING RATIO 1544. C N = D*2.079E 16 GIVES NUMBER DENSITY PER CM**3 1545. C 1546. C 1547. C 1548. C 1549. C 1550. DIMENSION PRS1(33),PRS2(33),PRS3(33),PRS4(33),PRS5(33),PRS6(33) 1551. 1 ,DNS1(33),DNS2(33),DNS3(33),DNS4(33),DNS5(33),DNS6(33) 1552. 2 ,TMP1(33),TMP2(33),TMP3(33),TMP4(33),TMP5(33),TMP6(33) 1553. 3 ,WVP1(33),WVP2(33),WVP3(33),WVP4(33),WVP5(33),WVP6(33) 1554. 4 ,OZO1(33),OZO2(33),OZO3(33),OZO4(33),OZO5(33),OZO6(33) 1555. DIMENSION PRES(33,6),DENS(33,6),TEMP(33,6),WVAP(33,6),OZON(33,6) 1556. C 1557. EQUIVALENCE 1558. + (PRES(1,1),PRS1(1)),(DENS(1,1),DNS1(1)),(TEMP(1,1),TMP1(1)) 1559. + ,(PRES(1,2),PRS2(1)),(DENS(1,2),DNS2(1)),(TEMP(1,2),TMP2(1)) 1560. + ,(PRES(1,3),PRS3(1)),(DENS(1,3),DNS3(1)),(TEMP(1,3),TMP3(1)) 1561. + ,(PRES(1,4),PRS4(1)),(DENS(1,4),DNS4(1)),(TEMP(1,4),TMP4(1)) 1562. + ,(PRES(1,5),PRS5(1)),(DENS(1,5),DNS5(1)),(TEMP(1,5),TMP5(1)) 1563. + ,(PRES(1,6),PRS6(1)),(DENS(1,6),DNS6(1)),(TEMP(1,6),TMP6(1)) 1564. EQUIVALENCE (WVAP(1,1),WVP1(1)),(OZON(1,1),OZO1(1)) 1565. EQUIVALENCE (WVAP(1,2),WVP2(1)),(OZON(1,2),OZO2(1)) 1566. EQUIVALENCE (WVAP(1,3),WVP3(1)),(OZON(1,3),OZO3(1)) 1567. EQUIVALENCE (WVAP(1,4),WVP4(1)),(OZON(1,4),OZO4(1)) 1568. EQUIVALENCE (WVAP(1,5),WVP5(1)),(OZON(1,5),OZO5(1)) 1569. EQUIVALENCE (WVAP(1,6),WVP6(1)),(OZON(1,6),OZO6(1)) 1570. C 1571. C 1572. DIMENSION HTKM(33) 1573. DATA HTKM/1.0E-09, 1., 2., 3., 4., 5., 6., 7., 8., 9.,10.,11. 1574. 1 ,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24. 1575. 2 ,25.,30.,35.,40.,45.,50.,70.,99.9/ 1576. C 1577. C 1578. C---------------------------------------------------------------------- 1579. C0000 GLOBAL U.S. (1976) STANDARD ATMOSPHERE P, T, GEO H PARAMETERS 1580. C---------------------------------------------------------------------- 1581. C 1582. DIMENSION SPLB(8),STLB(8),SHLB(8),SDLB(8) 1583. DATA SPLB/1013.25,226.32,54.748,8.6801,1.109,.66938,.039564 1584. + ,3.7338E-03/ 1585. DATA STLB/288.15,216.65,216.65,228.65,270.65,270.65,214.65,186.87/ 1586. DATA SHLB/0.0,11.0,20.0,32.0,47.0,51.0,71.0,84.852/ 1587. DATA SDLB/-6.5,0.0,1.0,2.8,0.0,-2.8,-2.0,0.0/ 1588. DATA HPCON/34.16319/ 1589. C 1590. C 1591. C----------------------------------------------------------------------- 1592. C1111 TROPICAL LATITUDES MCCLATCHY (1972) ATMOSPHERE DATA VS HEIGHT 1593. C----------------------------------------------------------------------- 1594. C 1595. DATA PRS1/ 1.013E 03,9.040E 02,8.050E 02,7.150E 02,6.330E 02, 1596. 1 5.590E 02,4.920E 02,4.320E 02,3.780E 02,3.290E 02,2.860E 02, 1597. 2 2.470E 02,2.130E 02,1.820E 02,1.560E 02,1.320E 02,1.110E 02, 1598. 3 9.370E 01,7.890E 01,6.660E 01,5.650E 01,4.800E 01,4.090E 01, 1599. 4 3.500E 01,3.000E 01,2.570E 01,1.220E 01,6.000E 00,3.050E 00, 1600. 5 1.590E 00,8.540E-01,5.790E-02,3.000E-04/ 1601. DATA DNS1/ 1.167E 03,1.064E 03,9.689E 02,8.756E 02,7.951E 02, 1602. 1 7.199E 02,6.501E 02,5.855E 02,5.258E 02,4.708E 02,4.202E 02, 1603. 2 3.740E 02,3.316E 02,2.929E 02,2.578E 02,2.260E 02,1.972E 02, 1604. 3 1.676E 02,1.382E 02,1.145E 02,9.515E 01,7.938E 01,6.645E 01, 1605. 4 5.618E 01,4.763E 01,4.045E 01,1.831E 01,8.600E 00,4.181E 00, 1606. 5 2.097E 00,1.101E 00,9.210E-02,5.000E-04/ 1607. DATA TMP1/ 300.0,294.0,288.0,284.0,277.0,270.0,264.0,257.0,250.0, 1608. 1244.0,237.0,230.0,224.0,217.0,210.0,204.0,197.0,195.0,199.0,203.0, 1609. 2207.0,211.0,215.0,217.0,219.0,221.0,232.0,243.0,254.0,265.0,270.0, 1610. 3 219.0,210.0/ 1611. DATA WVP1/1.9E 01,1.3E 01,9.3E 00,4.7E 00,2.2E 00,1.5E 00,8.5E-01, 1612. 1 4.7E-01,2.5E-01,1.2E-01,5.0E-02,1.7E-02,6.0E-03,1.8E-03,1.0E-03, 1613. 2 7.6E-04,6.4E-04,5.6E-04,5.0E-04,4.9E-04,4.5E-04,5.1E-04,5.1E-04, 1614. 3 5.4E-04,6.0E-04,6.7E-04,3.6E-04,1.1E-04,4.3E-05,1.9E-05,6.3E-06, 1615. 4 1.4E-07,1.0E-09/ 1616. DATA OZO1/5.6E-05,5.6E-05,5.4E-05,5.1E-05,4.7E-05,4.5E-05,4.3E-05, 1617. 1 4.1E-05,3.9E-05,3.9E-05,3.9E-05,4.1E-05,4.3E-05,4.5E-05,4.5E-05, 1618. 2 4.7E-05,4.7E-05,6.9E-05,9.0E-05,1.4E-04,1.9E-04,2.4E-04,2.8E-04, 1619. 3 3.2E-04,3.4E-04,3.4E-04,2.4E-04,9.2E-05,4.1E-05,1.3E-05,4.3E-06, 1620. 4 8.6E-08,4.3E-11/ 1621. C 1622. C----------------------------------------------------------------------- 1623. C2222 MIDLATITUDE SUMMER MCCLATCHY (1972) ATMOSPHERE DATA VS HEIGHT 1624. C----------------------------------------------------------------------- 1625. C 1626. DATA PRS2/ 1.013E 03,9.020E 02,8.020E 02,7.100E 02,6.280E 02, 1627. 1 5.540E 02,4.870E 02,4.260E 02,3.720E 02,3.240E 02,2.810E 02, 1628. 2 2.430E 02,2.090E 02,1.790E 02,1.530E 02,1.300E 02,1.110E 02, 1629. 3 9.500E 01,8.120E 01,6.950E 01,5.950E 01,5.100E 01,4.370E 01, 1630. 4 3.760E 01,3.220E 01,2.770E 01,1.320E 01,6.520E 00,3.330E 00, 1631. 5 1.760E 00,9.510E-01,6.710E-02,3.000E-04/ 1632. DATA DNS2/ 1.191E 03,1.080E 03,9.757E 02,8.846E 02,7.998E 02, 1633. 1 7.211E 02,6.487E 02,5.830E 02,5.225E 02,4.669E 02,4.159E 02, 1634. 2 3.693E 02,3.269E 02,2.882E 02,2.464E 02,2.104E 02,1.797E 02, 1635. 3 1.535E 02,1.305E 02,1.110E 02,9.453E 01,8.056E 01,6.872E 01, 1636. 4 5.867E 01,5.014E 01,4.288E 01,1.322E 01,6.519E 00,3.330E 00, 1637. 5 1.757E 00,9.512E-01,6.706E-02,5.000E-04/ 1638. DATA TMP2/ 294.0,290.0,285.0,279.0,273.0,267.0,261.0,255.0,248.0, 1639. 1242.0,235.0,229.0,222.0,216.0,216.0,216.0,216.0,216.0,216.0,217.0, 1640. 2218.0,219.0,220.0,222.0,223.0,224.0,234.0,245.0,258.0,270.0,276.0, 1641. 3 218.0,210.0/ 1642. DATA WVP2/1.4E 01,9.3E 00,5.9E 00,3.3E 00,1.9E 00,1.0E 00,6.1E-01, 1643. 1 3.7E-01,2.1E-01,1.2E-01,6.4E-02,2.2E-02,6.0E-03,1.8E-03,1.0E-03, 1644. 2 7.6E-04,6.4E-04,5.6E-04,5.0E-04,4.9E-04,4.5E-04,5.1E-04,5.1E-04, 1645. 3 5.4E-04,6.0E-04,6.7E-04,3.6E-04,1.1E-04,4.3E-05,1.9E-05,6.3E-06, 1646. 4 1.4E-07,1.0E-09/ 1647. DATA OZO2/6.0E-05,6.0E-05,6.0E-05,6.2E-05,6.4E-05,6.6E-05,6.9E-05, 1648. 1 7.5E-05,7.9E-05,8.6E-05,9.0E-05,1.1E-04,1.2E-04,1.5E-04,1.8E-04, 1649. 2 1.9E-04,2.1E-04,2.4E-04,2.8E-04,3.2E-04,3.4E-04,3.6E-04,3.6E-04, 1650. 3 3.4E-04,3.2E-04,3.0E-04,2.0E-04,9.2E-05,4.1E-05,1.3E-05,4.3E-06, 1651. 4 8.6E-08,4.3E-11/ 1652. C 1653. C----------------------------------------------------------------------- 1654. C3333 MIDLATITUDE WINTER MCCLATCHY (1972) ATMOSPHERE DATA VS HEIGHT 1655. C----------------------------------------------------------------------- 1656. C 1657. DATA PRS3/ 1.018E 03,8.973E 02,7.897E 02,6.938E 02,6.081E 02, 1658. 1 5.313E 02,4.627E 02,4.016E 02,3.473E 02,2.992E 02,2.568E 02, 1659. 2 2.199E 02,1.882E 02,1.610E 02,1.378E 02,1.178E 02,1.007E 02, 1660. 3 8.610E 01,7.350E 01,6.280E 01,5.370E 01,4.580E 01,3.910E 01, 1661. 4 3.340E 01,2.860E 01,2.430E 01,1.110E 01,5.180E 00,2.530E 00, 1662. 5 1.290E 00,6.820E-01,4.670E-02,3.000E-04/ 1663. DATA DNS3/ 1.301E 03,1.162E 03,1.037E 03,9.230E 02,8.282E 02, 1664. 1 7.411E 02,6.614E 02,5.886E 02,5.222E 02,4.619E 02,4.072E 02, 1665. 2 3.496E 02,2.999E 02,2.572E 02,2.206E 02,1.890E 02,1.620E 02, 1666. 3 1.388E 02,1.188E 02,1.017E 02,8.690E 01,7.421E 01,6.338E 01, 1667. 4 5.415E 01,4.624E 01,3.950E 01,1.783E 01,7.924E 00,3.625E 00, 1668. 5 1.741E 00,8.954E-01,7.051E-02,5.000E-04/ 1669. DATA TMP3/ 272.2,268.7,265.2,261.7,255.7,249.7,243.7,237.7,231.7, 1670. 1225.7,219.7,219.2,218.7,218.2,217.7,217.2,216.7,216.2,215.7,215.2, 1671. 2215.2,215.2,215.2,215.2,215.2,215.2,217.4,227.8,243.2,258.5,265.7, 1672. 3 230.7,210.2/ 1673. DATA WVP3/3.5E 00,2.5E 00,1.8E 00,1.2E 00,6.6E-01,3.8E-01,2.1E-01, 1674. 1 8.5E-02,3.5E-02,1.6E-02,7.5E-03,6.9E-03,6.0E-03,1.8E-03,1.0E-03, 1675. 2 7.6E-04,6.4E-04,5.6E-04,5.0E-04,4.9E-04,4.5E-04,5.1E-04,5.1E-04, 1676. 3 5.4E-04,6.0E-04,6.7E-04,3.6E-04,1.1E-04,4.3E-05,1.9E-05,6.3E-06, 1677. 4 1.4E-07,1.0E-09/ 1678. DATA OZO3/6.0E-05,5.4E-05,4.9E-05,4.9E-05,4.9E-05,5.8E-05,6.4E-05, 1679. 1 7.7E-05,9.0E-05,1.2E-04,1.6E-04,2.1E-04,2.6E-04,3.0E-04,3.2E-04, 1680. 2 3.4E-04,3.6E-04,3.9E-04,4.1E-04,4.3E-04,4.5E-04,4.3E-04,4.3E-04, 1681. 3 3.9E-04,3.6E-04,3.4E-04,1.9E-04,9.2E-05,4.1E-05,1.3E-05,4.3E-06, 1682. 4 8.6E-08,4.3E-11/ 1683. C 1684. C----------------------------------------------------------------------- 1685. C4444 SUBARCTIC SUMMER MCCLATCHY (1972) ATMOSPHERE DATA VS HEIGHT 1686. C----------------------------------------------------------------------- 1687. C 1688. DATA PRS4/ 1.010E 03,8.960E 02,7.929E 02,7.000E 02,6.160E 02, 1689. 1 5.410E 02,4.730E 02,4.130E 02,3.590E 02,3.107E 02,2.677E 02, 1690. 2 2.300E 02,1.977E 02,1.700E 02,1.460E 02,1.250E 02,1.080E 02, 1691. 3 9.280E 01,7.980E 01,6.860E 01,5.890E 01,5.070E 01,4.360E 01, 1692. 4 3.750E 01,3.227E 01,2.780E 01,1.340E 01,6.610E 00,3.400E 00, 1693. 5 1.810E 00,9.870E-01,7.070E-02,3.000E-04/ 1694. DATA DNS4/ 1.220E 03,1.110E 03,9.971E 02,8.985E 02,8.077E 02, 1695. 1 7.244E 02,6.519E 02,5.849E 02,5.231E 02,4.663E 02,4.142E 02, 1696. 2 3.559E 02,3.059E 02,2.630E 02,2.260E 02,1.943E 02,1.671E 02, 1697. 3 1.436E 02,1.235E 02,1.062E 02,9.128E 01,7.849E 01,6.750E 01, 1698. 4 5.805E 01,4.963E 01,4.247E 01,1.338E 01,6.614E 00,3.404E 00, 1699. 5 1.817E 00,9.868E-01,7.071E-02,5.000E-04/ 1700. DATA TMP4/ 287.0,282.0,276.0,271.0,266.0,260.0,253.0,246.0,239.0, 1701. 1232.0,225.0,225.0,225.0,225.0,225.0,225.0,225.0,225.0,225.0,225.0, 1702. 2225.0,225.0,225.0,225.0,226.0,228.0,235.0,247.0,262.0,274.0,277.0, 1703. 3 216.0,210.0/ 1704. DATA WVP4/9.1E 00,6.0E 00,4.2E 00,2.7E 00,1.7E 00,1.0E 00,5.4E-01, 1705. 1 2.9E-01,1.3E-02,4.2E-02,1.5E-02,9.4E-03,6.0E-03,1.8E-03,1.0E-03, 1706. 2 7.6E-04,6.4E-04,5.6E-04,5.0E-04,4.9E-04,4.5E-04,5.1E-04,5.1E-04, 1707. 3 5.4E-04,6.0E-04,6.7E-04,3.6E-04,1.1E-04,4.3E-05,1.9E-05,6.3E-06, 1708. 4 1.4E-07,1.0E-09/ 1709. DATA OZO4/4.9E-05,5.4E-05,5.6E-05,5.8E-05,6.0E-05,6.4E-05,7.1E-05, 1710. 1 7.5E-05,7.9E-05,1.1E-04,1.3E-04,1.8E-04,2.1E-04,2.6E-04,2.8E-04, 1711. 2 3.2E-04,3.4E-04,3.9E-04,4.1E-04,4.1E-04,3.9E-04,3.6E-04,3.2E-04, 1712. 3 3.0E-04,2.8E-04,2.6E-04,1.4E-04,9.2E-05,4.1E-05,1.3E-05,4.3E-06, 1713. 4 8.6E-08,4.3E-11/ 1714. C 1715. C----------------------------------------------------------------------- 1716. C5555 SUBARCTIC WINTER MCCLATCHY (1972) ATMOSPHERE DATA VS HEIGHT 1717. C----------------------------------------------------------------------- 1718. C 1719. DATA PRS5/ 1.013E 03,8.878E 02,7.775E 02,6.798E 02,5.932E 02, 1720. 1 5.158E 02,4.467E 02,3.853E 02,3.308E 02,2.829E 02,2.418E 02, 1721. 2 2.067E 02,1.766E 02,1.510E 02,1.291E 02,1.103E 02,9.431E 01, 1722. 3 8.058E 01,6.882E 01,5.875E 01,5.014E 01,4.277E 01,3.647E 01, 1723. 4 3.109E 01,2.649E 01,2.256E 01,1.020E 01,4.701E 00,2.243E 00, 1724. 5 1.113E 00,5.719E-01,4.016E-02,3.000E-04/ 1725. DATA DNS5/ 1.372E 03,1.193E 03,1.058E 03,9.366E 02,8.339E 02, 1726. 1 7.457E 02,6.646E 02,5.904E 02,5.226E 02,4.538E 02,3.879E 02, 1727. 2 3.315E 02,2.834E 02,2.422E 02,2.071E 02,1.770E 02,1.517E 02, 1728. 3 1.300E 02,1.113E 02,9.529E 01,8.155E 01,6.976E 01,5.966E 01, 1729. 4 5.100E 01,4.358E 01,3.722E 01,1.645E 01,7.368E 00,3.330E 00, 1730. 5 1.569E 00,7.682E-01,5.695E-02,5.000E-04/ 1731. DATA TMP5/ 257.1,259.1,255.9,252.7,247.7,240.9,234.1,227.3,220.6, 1732. 1217.2,217.2,217.2,217.2,217.2,217.2,217.2,216.6,216.0,215.4,214.8, 1733. 2214.1,213.6,213.0,212.4,211.8,211.2,216.0,222.2,234.7,247.0,259.3, 1734. 3 245.7,210.0/ 1735. DATA WVP5/1.2E 00,1.2E 00,9.4E-01,6.8E-01,4.1E-01,2.0E-01,9.8E-02, 1736. 1 5.4E-02,1.1E-02,8.4E-03,5.5E-03,3.8E-03,2.6E-03,1.8E-03,1.0E-03, 1737. 2 7.6E-04,6.4E-04,5.6E-04,5.0E-04,4.9E-04,4.5E-04,5.1E-04,5.1E-04, 1738. 3 5.4E-04,6.0E-04,6.7E-04,3.6E-04,1.1E-04,4.3E-05,1.9E-05,6.3E-06, 1739. 4 1.4E-07,1.0E-09/ 1740. DATA OZO5/4.1E-05,4.1E-05,4.1E-05,4.3E-05,4.5E-05,4.7E-05,4.9E-05, 1741. 1 7.1E-05,9.0E-05,1.6E-04,2.4E-04,3.2E-04,4.3E-04,4.7E-04,4.9E-04, 1742. 2 5.6E-04,6.2E-04,6.2E-04,6.2E-04,6.0E-04,5.6E-04,5.1E-04,4.7E-04, 1743. 3 4.3E-04,3.6E-04,3.2E-04,1.5E-04,9.2E-05,4.1E-05,1.3E-05,4.3E-06, 1744. 4 8.6E-08,4.3E-11/ 1745. C 1746. C---------------------------------------------------------------------- 1747. C6666 GLOBAL U.S. (1976) STANDARD ATMOSPHERE P, T, GEO H PARAMETERS 1748. C---------------------------------------------------------------------- 1749. C 1750. DATA PRS6/ 1.01325E+03,8.987E+02,7.950E+02,7.011E+02,6.164E+02, 1751. 1 5.402E+02,4.718E+02,4.106E+02,3.560E+02,3.074E+02,2.644E+02, 1752. 2 2.263E+02,1.933E+02,1.651E+02,1.410E+02,1.204E+02,1.029E+02, 1753. 3 8.787E+01,7.505E+01,6.410E+01,5.475E+01,4.678E+01,4.000E+01, 1754. 4 3.422E+01,2.931E+01,2.511E+01,1.172E+01,5.589E+00,2.775E+00, 1755. 5 1.431E+00,7.594E-01,4.634E-02,2.384E-04/ 1756. DATA DNS6/ 1.225E+03,1.112E+03,1.006E+03,9.091E+02,8.191E+02, 1757. 1 7.361E+02,6.597E+02,5.895E+02,5.252E+02,4.663E+02,4.127E+02, 1758. 2 3.639E+02,3.108E+02,2.655E+02,2.268E+02,1.937E+02,1.654E+02, 1759. 3 1.413E+02,1.207E+02,1.031E+02,8.803E+01,7.487E+01,6.373E+01, 1760. 4 5.428E+01,4.627E+01,3.947E+01,1.801E+01,8.214E+00,3.851E+00, 1761. 5 1.881E+00,9.775E-01,7.424E-02,4.445E-04/ 1762. DATA TMP6/ 1763. 1 288.150,281.650,275.150,268.650,262.150,255.650,249.150, 1764. 2 242.650,236.150,229.650,223.150,216.650,216.650,216.650, 1765. 3 216.650,216.650,216.650,216.650,216.650,216.650,216.650, 1766. 4 217.650,218.650,219.650,220.650,221.650,226.650,237.050, 1767. 5 251.050,265.050,270.650,217.450,186.870/ 1768. DATA WVP6/ 1.083E+01,6.323E+00,3.612E+00,2.015E+00,1.095E+00, 1769. 1 5.786E-01,2.965E-01,1.469E-01,7.021E-02,3.226E-02,1.419E-02, 1770. 2 5.956E-03,5.002E-03,4.186E-03,3.490E-03,2.896E-03,2.388E-03, 1771. 3 1.954E-03,1.583E-03,1.267E-03,9.967E-04,8.557E-04,7.104E-04, 1772. 4 5.600E-04,4.037E-04,2.406E-04,5.404E-05,2.464E-05,1.155E-05, 1773. 5 5.644E-06,2.932E-06,2.227E-07,1.334E-09/ 1774. DATA OZO6/ 7.526E-05,3.781E-05,6.203E-05,3.417E-05,5.694E-05, 1775. 1 3.759E-05,5.970E-05,4.841E-05,7.102E-05,6.784E-05,9.237E-05, 1776. 2 9.768E-05,1.251E-04,1.399E-04,1.715E-04,1.946E-04,2.300E-04, 1777. 3 2.585E-04,2.943E-04,3.224E-04,3.519E-04,3.714E-04,3.868E-04, 1778. 4 3.904E-04,3.872E-04,3.728E-04,2.344E-04,9.932E-05,3.677E-05, 1779. 5 1.227E-05,4.324E-06,5.294E-08,1.262E-10/ 1780. C 1781. C 1782. IF(NATM.GT.0) GO TO 200 1783. O=1.E-10 1784. Q=1.E-10 1785. S=1.E-10 1786. OCM=1.E-10 1787. WCM=1.E-10 1788. IF(NPHD.LT.2) GO TO 150 1789. DO 110 N=2,8 1790. IF(H.LT.SHLB(N)) GO TO 120 1791. 110 CONTINUE 1792. N=9 1793. 120 N=N-1 1794. IF(ABS(SDLB(N)).LT.1.E-04) GO TO 130 1795. P=SPLB(N)*(1.E0+SDLB(N)/STLB(N)*(H-SHLB(N)))**(-HPCON/SDLB(N)) 1796. GO TO 140 1797. 130 P=SPLB(N)*EXP(-HPCON/STLB(N)*(H-SHLB(N))) 1798. 140 T=STLB(N)+SDLB(N)*(H-SHLB(N)) 1799. D=P/T*28.9644E 05/8.31432E 03 1800. RETURN 1801. C 1802. 150 CONTINUE 1803. DO 160 N=2,8 1804. IF(P.GT.SPLB(N)) GO TO 170 1804.2 160 CONTINUE 1805. N=9 1806. 170 N=N-1 1807. IF(ABS(SDLB(N)).LT.1.E-04) GO TO 180 1808. H=SHLB(N)+STLB(N)/SDLB(N)*((SPLB(N)/P)**(SDLB(N)/HPCON)-1.E0) 1809. GO TO 190 1810. 180 H=SHLB(N)+STLB(N)/HPCON*ALOG(SPLB(N)/P) 1811. 190 T=STLB(N)+SDLB(N)*(H-SHLB(N)) 1812. D=P/T*28.9644E 05/8.31432E 03 1813. RETURN 1814. C 1815. 200 CONTINUE 1816. IF(NPHD.EQ.1) GO TO 240 1817. IF(NPHD.EQ.2) GO TO 220 1818. XX=D 1819. XI=DENS(1,NATM) 1820. IF(D.GT.XI) XX=XI 1821. IF(D.LT.5.0E-04) GO TO 280 1822. DO 210 J=2,33 1823. XJ=DENS(J,NATM) 1824. IF(XX.GT.XJ) GO TO 260 1825. 210 XI=XJ 1826. 220 XX=H 1827. XI=HTKM(1) 1828. IF(H.LT.XI) XX=XI 1829. IF(H.GT.99.9E0) GO TO 280 1830. DO 230 J=2,33 1831. XJ=HTKM(J) 1832. IF(XX.LT.XJ) GO TO 260 1833. 230 XI=XJ 1834. 240 XX=P 1835. XI=PRES(1,NATM) 1836. IF(P.GT.XI) XX=XI 1837. IF(P.LT.3.0E-04) GO TO 280 1838. DO 250 J=2,33 1839. XJ=PRES(J,NATM) 1840. IF(XX.GT.XJ) GO TO 260 1841. 250 XI=XJ 1842. 260 DELTA=(XX-XI)/(XJ-XI) 1843. I=J-1 1844. IF(NPHD.NE.2) H=HTKM(I)+(HTKM(J)-HTKM(I))*ALOG(XX/XI)/ALOG(XJ/XI) 1845. PI=PRES(I,NATM) 1846. PJ=PRES(J,NATM) 1847. DI=DENS(I,NATM) 1848. DJ=DENS(J,NATM) 1849. IF(NPHD.NE.1) P=PI+DELTA*(PJ-PI) 1850. IF(NPHD.NE.3) D=DI+DELTA*(DJ-DI) 1851. T=TEMP(I,NATM)+DELTA*(TEMP(J,NATM)-TEMP(I,NATM)) 1852. O=OZON(I,NATM)/DI+DELTA*(OZON(J,NATM)/DJ-OZON(I,NATM)/DI) 1853. Q=WVAP(I,NATM)/DI+DELTA*(WVAP(J,NATM)/DJ-WVAP(I,NATM)/DI) 1854. ES=10.**(9.4051-2353./T) 1855. IF(P.LT.PI) PI=P 1856. S=1.E+06 1857. RS=(PI-ES+0.622E0*ES)/(0.622E0*ES) 1858. IF(RS.GT.1.E-06) S=1.E0/RS 1859. OI=O 1860. QI=Q 1861. OCM=0.E0 1862. WCM=0.E0 1863. DO 270 K=J,33 1864. PJ=PRES(K,NATM) 1865. DJ=DENS(K,NATM) 1866. OJ=OZON(K,NATM)/DJ 1867. QJ=WVAP(K,NATM)/DJ 1868. DP=PI-PJ 1869. OCM=OCM+0.5E0*(OI+OJ)*DP 1870. WCM=WCM+0.5E0*(QI+QJ)*DP 1871. OI=OJ 1872. QI=QJ 1873. 270 PI=PJ 1874. WCM=WCM/0.980E0*22420.7E0/18.0E0 1875. OCM=OCM/0.980E0*22420.7E0/48.0E0 1876. RETURN 1877. 280 T=210.0E0 1878. IF(NATM.EQ.6) T=186.87 1879. O=1.E-10 1880. Q=1.E-10 1881. S=1.E-10 1882. OCM=1.E-10 1883. WCM=1.E-10 1884. IF(NPHD.NE.1) P=1.E-05 1885. IF(NPHD.NE.2) H=99.99E0 1886. IF(NPHD.NE.3) D=2.E-05 1887. RETURN 1888. END 2000. 2001. SUBROUTINE ALBEDD 2003. C**** 2004. C**** ALBEDD calculates the surface albedo for differnt K-intervals 2005. C**** 2006. INCLUDE 'RA99M9.COM' 2007. C**** Working common block 2008. COMMON /WORK04/ PL(40),DPL(40), 2009. T TRCALB(40),TRGALB(40),BGFEMT(40),BGFEMD(40), 2010. T TRAEXT(40,25),TAUN(LMR,NKTR),TAUSL(25),FTAUSL(25), 2011. T DFLB(40,25),UFLB(40,25),WFLB(40,25),DFSL(25),UFSL(25),WFSL(25), 2012. S EXTAER(40,6),SCTAER(40,6),COSAER(40,6) 2013. PARAMETER (NV=10,WETTRA=1.,WETSRA=1.,DMOICE=10.,DMLICE=10., 2014. * ECLTRA=1.,EOCTRA=1.,ESNTRA=1.,EICTRA=1.,EDSTRA=1.,EVGTRA=1., 2015. * ZCLSRA=1.,ZOCSRA=1.,ZSNSRA=1.,ZICSRA=1.,ZDSSRA=1.,ZVGSRA=1.) 2018. REAL*4 ALVISK(11,4),ALNIRK(11,4),FIELDC(11,3),VTMASK(11), 2019. * ALVISD(11,2),ALNIRD(11,2), SRBALB(6),SRXALB(6) 2020. REAL*8 FVNOS(NV) 2021. INTEGER*4 ISEAS(0:5) 2022. EQUIVALENCE (VADATA(1,1,1),ALVISK),(VADATA(1,1,2),ALNIRK), 2023. * (VADATA(1,1,3),FIELDC),(VADATA(1,4,3),VTMASK), 2024. * (ASNALB(1),ASNVIS),(ASNALB(2),ASNNIR), 2025. * (AOIALB(1),AOIVIS),(AOIALB(2),AOINIR), 2026. * (ALIALB(1),ALIVIS),(ALIALB(2),ALINIR), 2027. * (SRBXAL(1,1),SRBALB),(SRBXAL(1,2),SRXALB) 2027.1 EQUIVALENCE 2027.2 + (BXA(1),BOCVIS),(BXA(5),BEAVIS),(BXA( 9),BOIVIS),(BXA(13),BLIVIS) 2027.3 +,(BXA(2),BOCNIR),(BXA(6),BEANIR),(BXA(10),BOINIR),(BXA(14),BLINIR) 2027.4 +,(BXA(3),XOCVIS),(BXA(7),XEAVIS),(BXA(11),XOIVIS),(BXA(15),XLIVIS) 2027.5 +,(BXA(4),XOCNIR),(BXA(8),XEANIR),(BXA(12),XOINIR),(BXA(16),XLINIR) 2028. C**** WIN SPR SUM ATM 2029. DATA ISEAS / -77, 15, 105, 196, 288, 380 / 2038. JEQ = JM/2+1 2039. C**** 2040. C**** Define VISual and Near IR ALbedoes, ALVISD and ALNIRD, 2041. C**** for different vegetation types 2042. C**** 2045. DO 10 KS=1,4 2048. 10 IF(ISEAS(KS).ge.JDAY) GO TO 20 2049. C KS=5 2051. 20 WT1 = DFLOAT(ISEAS(KS)-JDAY) / (ISEAS(KS)-ISEAS(KS-1)) 2053. KS1 = 1+MOD(KS,4) 2054. KS2 = 1+MOD(KS+1,4) 2055. KN1 = 1+MOD(KS+2,4) 2056. KN2 = 1+MOD(KS+3,4) 2057. DO 30 K=1,NV 2058. C**** Southern hemisphere 2059. ALVISD(K,1) = ALVISK(K,KS1)*WT1 + ALVISK(K,KS2)*(1.-WT1) 2060. ALNIRD(K,1) = ALNIRK(K,KS1)*WT1 + ALNIRK(K,KS2)*(1.-WT1) 2061. C**** Northern hemisphere 2062. ALVISD(K,2) = ALVISK(K,KN1)*WT1 + ALVISK(K,KN2)*(1.-WT1) 2063. 30 ALNIRD(K,2) = ALNIRK(K,KN1)*WT1 + ALNIRK(K,KN2)*(1.-WT1) 2064. RETURN 2065. C**** 2066. C**** 2067. ENTRY ALBED 2068. C**** 2069. C**** ALBED calculates the solar albedoes, the thermal albedoes, and 2070. C**** derivative of Planck function with respect to ground temperature 2071. C**** 2072. IHEM = 2 2073. IF(JLAT.LT.JEQ) IHEM = 1 2074. C**** 2075. C**** Thermal albedoes and Planck function derivative 2076. C**** 2077. C**** Snow 2078. ASNAGE = 0.35*EXP(-0.2*AGESN) 2079. DO 110 K=1,NKTR 2080. TRGALB(K) = 0. 2081. BGFEMD(K) = 0. 2082. 110 BGFEMT(K) = 0. 2083. BOCSUM = 0. 2084. BEASUM = 0. 2085. BOISUM = 0. 2086. BLISUM = 0. 2087. DO 120 ITYPE=1,4 2088. 120 DTRUFG(ITYPE) = 0. 2089. C**** Ocean albedo specification 2090. IF(POCEAN.le.0.) GO TO 220 2091. X = 1./(1.+WS) 2092. AV = (-.0147087*X*X + .0292266*X - .0081079)*EOCTRA 2093. BV = (1.01673 - .0083652*WS)*EOCTRA 2094. ITOC = TGO 2095. WTOC = TGO-ITOC 2096. BOCM = 0. 2097. BOCP = 0. 2098. BOCSUM = 0. 2099. DO 210 K=1,NKTR 2100. TRAPOC = AV + BV*AOCEAN(K) 2101. BOCM1 = (PLANCK(ITOC-1,K) 2102. - -(PLANCK(ITOC-1,K)-PLANCK(ITOC ,K))*WTOC)*(1.-TRAPOC) 2103. BOCM = BOCM+BOCM1 2104. BOCP1 = (PLANCK(ITOC+1,K) 2105. - -(PLANCK(ITOC+1,K)-PLANCK(ITOC+2,K))*WTOC)*(1.-TRAPOC) 2106. BOCP = BOCP+BOCP1 2107. BOC = (PLANCK(ITOC ,K) 2108. - -(PLANCK(ITOC ,K)-PLANCK(ITOC+1,K))*WTOC)*(1.-TRAPOC) 2109. BOCSUM = BOCSUM+BOC 2110. TRGALB(K) = TRGALB(K) + POCEAN*TRAPOC 2111. BGFEMD(K) = BGFEMD(K) + POCEAN*(BOCP1-BOCM1) 2112. 210 BGFEMT(K) = BGFEMT(K) + POCEAN*BOC 2113. DTRUFG(1) = .5*(BOCP-BOCM) 2114. C**** Ocean ice 2115. 220 IF(POICE.le.0.) GO TO 240 2116. EXPSNO = EXP(-SNOWOI/DMOICE) 2117. ITOI = TGOI 2118. WTOI = TGOI-ITOI 2119. BOIM = 0. 2120. BOIP = 0. 2121. BOISUM = 0. 2122. DO 230 K=1,NKTR 2123. TRAPOI = AGSIDV(K,1)*ESNTRA*(1.-EXPSNO) 2124. + +AGSIDV(K,2)*EICTRA*EXPSNO 2125. BOIM1 = (PLANCK(ITOI-1,K) 2126. - -(PLANCK(ITOI-1,K)-PLANCK(ITOI ,K))*WTOI)*(1.-TRAPOI) 2127. BOIM = BOIM+BOIM1 2128. BOIP1 = (PLANCK(ITOI+1,K) 2129. - -(PLANCK(ITOI+1,K)-PLANCK(ITOI+2,K))*WTOI)*(1.-TRAPOI) 2130. BOIP = BOIP+BOIP1 2131. BOI = (PLANCK(ITOI ,K) 2132. - -(PLANCK(ITOI ,K)-PLANCK(ITOI+1,K))*WTOI)*(1.-TRAPOI) 2133. BOISUM = BOISUM+BOI 2134. TRGALB(K) = TRGALB(K) + POICE*TRAPOI 2135. BGFEMD(K) = BGFEMD(K) + POICE*(BOIP1-BOIM1) 2136. 230 BGFEMT(K) = BGFEMT(K) + POICE*BOI 2137. DTRUFG(3) = .5*(BOIP-BOIM) 2138. C**** Land ice 2139. 240 IF(PLICE.le.0.) GO TO 300 2140. EXPSNL = EXP(-SNOWLI/DMLICE) 2141. ITLI = TGLI 2142. WTLI = TGLI-ITLI 2143. BLIM = 0. 2144. BLIP = 0. 2145. BGF = 0. 2146. BLISUM = 0. 2147. DO 250 K=1,NKTR 2148. TRAPLI = AGSIDV(K,1)*ESNTRA*(1.-EXPSNL) 2149. + +AGSIDV(K,2)*EICTRA*EXPSNL 2150. BLIM1 = (PLANCK(ITLI-1,K) 2151. - -(PLANCK(ITLI-1,K)-PLANCK(ITLI ,K))*WTLI)*(1.-TRAPLI) 2152. BLIM = BLIM+BLIM1 2153. BLIP1 = (PLANCK(ITLI+1,K) 2154. - -(PLANCK(ITLI+1,K)-PLANCK(ITLI+2,K))*WTLI)*(1.-TRAPLI) 2155. BLIP = BLIP+BLIP1 2156. BLI = (PLANCK(ITLI ,K) 2157. - -(PLANCK(ITLI ,K)-PLANCK(ITLI+1,K))*WTLI)*(1.-TRAPLI) 2158. BLISUM = BLISUM+BLI 2159. TRGALB(K) = TRGALB(K) + PLICE*TRAPLI 2160. BGFEMD(K) = BGFEMD(K) + PLICE*(BLIP1-BLIM1) 2161. 250 BGFEMT(K) = BGFEMT(K) + PLICE*BLI 2162. DTRUFG(4) = .5*(BLIP-BLIM) 2163. C**** 2164. C**** Ground: FSNOW + FBARE + FVEGE = 1 2164.1 C**** If there is snow, then FVNOS is the no snow vegetation fraction 2165. C**** 2166. 300 IF(PEARTH.le.0.) GO TO 400 2167. FBARE = PVT(1) + PVT(10) 2168. FSNOW = 0. 2169. IF(SNOWE .le.1.E-4) GO TO 320 2170. FVNOS( 1) = PVT( 1)*EXP(-SNOWE/VTMASK( 1)) 2170.1 FVNOS(10) = PVT(10)*EXP(-SNOWE/VTMASK(10)) 2170.2 FBARE = FVNOS(1) + FVNOS(10) 2171. FSNOW = 1. - FBARE 2172. DO 310 K=2,NV-1 2173. FVNOS(K) = PVT(K)*EXP(-SNOWE/VTMASK(K)) 2174. 310 FSNOW = FSNOW - FVNOS(K) 2175. 320 FVEGE = 1. - FBARE - FSNOW 2176. ITEA = TGE 2177. WTEA = TGE-ITEA 2178. BEAM = 0. 2179. BEAP = 0. 2180. BEASUM = 0. 2181. DO 330 K=1,NKTR 2182. TRAPEA = AGSIDV(K,1)*FSNOW 2183. + +AGSIDV(K,3)*FBARE*(1.-WETTRA*WEARTH) 2184. + +AGSIDV(K,4)*FVEGE 2185. BEAM1 = (PLANCK(ITEA-1,K) 2186. - -(PLANCK(ITEA-1,K)-PLANCK(ITEA ,K))*WTEA)*(1.-TRAPEA) 2187. BEAM = BEAM+BEAM1 2188. BEAP1 = (PLANCK(ITEA+1,K) 2189. - -(PLANCK(ITEA+1,K)-PLANCK(ITEA+2,K))*WTEA)*(1.-TRAPEA) 2190. BEAP = BEAP+BEAP1 2191. BEA = (PLANCK(ITEA ,K) 2192. - -(PLANCK(ITEA ,K)-PLANCK(ITEA+1,K))*WTEA)*(1.-TRAPEA) 2193. BEASUM = BEASUM+BEA 2194. TRGALB(K) = TRGALB(K) + PEARTH*TRAPEA 2195. BGFEMD(K) = BGFEMD(K) + PEARTH*(BEAP1-BEAM1) 2196. 330 BGFEMT(K) = BGFEMT(K) + PEARTH*BEA 2197. DTRUFG(2) = .5*(BEAP-BEAM) 2198. C**** 2199. C**** Calculate composite thermal albedoes at the surface 2200. C**** 2201. 400 BGF = 0. 2202. DO 410 K=1,NKTR 2203. BGFEMD(K) = BGFEMD(K)*0.5 2204. 410 BGF = BGF + BGFEMT(K) 2205. BGM = BOCM*POCEAN + BEAM*PEARTH + BOIM*POICE + BLIM*PLICE 2206. BGP = BOCP*POCEAN + BEAP*PEARTH + BOIP*POICE + BLIP*PLICE 2207. TTRUFG = .5*(BGP-BGM) 2208. FTRUFG(1) = BOCSUM/BGF 2209. FTRUFG(2) = BEASUM/BGF 2210. FTRUFG(3) = BOISUM/BGF 2211. FTRUFG(4) = BLISUM/BGF 2212. C**** 2213. C**** Solar albedoes 2214. C**** 2215. IF(COSZ.lt..01) RETURN 2216. C**** Snow albedo 2217. BSNVIS = ASNVIS + ASNAGE 2218. BSNNIR = ASNNIR + ASNAGE 2219. C XSNVIS = BSNVIS 2220. C XSNNIR = BSNNIR 2221. DO 510 I=1,16 2222. 510 BXA(I) = 0. 2223. C**** Open water 2224. IF(POCEAN.le.0.) GO TO 620 2228. BOCNIR = .08 2229. XOCNIR = .05 / (1.1*COSZ**1.4 + .15) 2230. BOCVIS = .08 2231. XOCVIS = XOCNIR 2240. C**** Ocean ice 2241. 620 IF(POICE.le.0.) GO TO 640 2242. BOIVIS = AOIVIS*EXPSNO + BSNVIS*(1.-EXPSNO) 2243. BOINIR = AOINIR*EXPSNO + BSNNIR*(1.-EXPSNO) 2244. XOIVIS = BOIVIS 2245. XOINIR = BOINIR 2246. C**** Glacial ice 2247. 640 IF(PLICE.le.0.) GO TO 700 2248. BLIVIS = ALIVIS*EXPSNL + BSNVIS*(1.-EXPSNL) 2249. BLINIR = ALINIR*EXPSNL + BSNNIR*(1.-EXPSNL) 2250. XLIVIS = BLIVIS 2251. XLINIR = BLINIR 2252. C**** 2253. C**** Ground: sum PVT = 1; or FSNOW + sum FVNOS = 1 2254. C**** 2255. 700 IF(PEARTH.le.0.) GO TO 800 2256. IF(SNOWE.gt.1.E-4) GO TO 720 2257. C**** Ground without snow 2258. BEAVIS = (PVT( 1)*ALVISD( 1,IHEM) 2258.1 * +PVT(10)*ALVISD(10,IHEM))*(1.-.5*WEARTH*WETSRA) 2259. BEANIR = (PVT( 1)*ALNIRD( 1,IHEM) 2259.1 * +PVT(10)*ALNIRD(10,IHEM))*(1.-.5*WEARTH*WETSRA) 2260. DO 710 K=2,NV-1 2261. BEAVIS = BEAVIS + PVT(K)*ALVISD(K,IHEM) 2262. 710 BEANIR = BEANIR + PVT(K)*ALNIRD(K,IHEM) 2263. GO TO 740 2264. C**** Snow on ground 2265. 720 BEAVIS = (FVNOS( 1)*ALVISD( 1,IHEM) 2265.1 * +FVNOS(10)*ALVISD(10,IHEM))*(1.-.5*WEARTH*WETSRA) 2266. BEANIR = (FVNOS( 1)*ALNIRD( 1,IHEM) 2266.1 * +FVNOS(10)*ALNIRD(10,IHEM))*(1.-.5*WEARTH*WETSRA) 2267. DO 730 K=2,NV-1 2268. BEAVIS = BEAVIS + FVNOS(K)*ALVISD(K,IHEM) 2269. 730 BEANIR = BEANIR + FVNOS(K)*ALNIRD(K,IHEM) 2270. BEAVIS = BEAVIS + FSNOW*BSNVIS 2271. BEANIR = BEANIR + FSNOW*BSNNIR 2272. 740 XEAVIS = BEAVIS 2273. XEANIR = BEANIR 2274. C**** 2275. C**** Calculate composite solar albedoes at the surface 2276. C**** 2277. 800 SRBALB(6) = POCEAN*BOCVIS+PEARTH*BEAVIS+POICE*BOIVIS+PLICE*BLIVIS 2278. SRXALB(6) = POCEAN*XOCVIS+PEARTH*XEAVIS+POICE*XOIVIS+PLICE*XLIVIS 2279. SRBALB(1) = POCEAN*BOCNIR+PEARTH*BEANIR+POICE*BOINIR+PLICE*BLINIR 2280. SRXALB(1) = POCEAN*XOCNIR+PEARTH*XEANIR+POICE*XOINIR+PLICE*XLINIR 2281. DO 810 K=2,5 2282. SRBALB(K) = SRBALB(1) 2283. 810 SRXALB(K) = SRXALB(1) 2284. RETURN 2285. END 2500. 2501. SUBROUTINE AEROS0 2502. C**** 2503. INCLUDE 'RA99M9.COM' 2504. C**** Working common block 2505. COMMON /WORK04/ PL(40),DPL(40), 2506. T TRCALB(40),TRGALB(40),BGFEMT(40),BGFEMD(40), 2507. T TRAEXT(40,25),TAUN(LMR,NKTR),TAUSL(25),FTAUSL(25), 2508. T DFLB(40,25),UFLB(40,25),WFLB(40,25),DFSL(25),UFSL(25),WFSL(25), 2509. S EXTAER(40,6),SCTAER(40,6),COSAER(40,6) 2510. C**** 2511. PARAMETER (FOGTSL=0., FCLDTR=1., ECLTRA=1.) 2512. REAL*4 SRAX(LMR,NKSR,5),SRAS(LMR,NKSR,5),SRAC(LMR,NKSR,5), 2513. * OTN(LMR,11,5) 2514. C**** 2515. IF(NGOLDH.gt.5 .or. NAERO.gt.11) GO TO 800 2516. C**** 2517. C**** Calculate normalized optical thicknesses for each layer 2518. C**** 2519. WRITE (6,900) (NINT(PLE(L)),L=1,LMR+1) 2520. WRITE (6,901) (HLE(L),L=1,LMR+1) 2521. WRITE (6,902) ( L ,L=1,LMR) 2522. WRITE (6,903) 2523. DO 100 NG=1,NGOLDH 2524. DO 100 NA=1,NAERO 2525. DO 100 L=1,LMR 2526. 100 OTN(L,NA,NG) = 0. 2527. C**** Stratospheric volcanic aerosol (NA=1,NG=1), layers 8 and 9 2528. OTN(8,1,1) = AGOLDH(1,1)*.375 2529. OTN(9,1,1) = AGOLDH(1,1)*.625 2530. WRITE (6,910) 1,1, (OTN(L,1,1),L=1,9) 2531. C**** Global ocean aerosol (NG=2), global land aerosol (NG=3) 2532. DO 130 NG=2,3 2533. DO 130 NA=1,NAERO 2534. IF(AGOLDH(NA,NG).lt.1.e-6) GO TO 130 2535. C = CGOLDH(NA,NG) 2536. BC = EXP(-BGOLDH(NA,NG)/C) 2537. ABC = AGOLDH(NA,NG)*(1.+BC) 2538. DO 120 L=1,LMR 2539. 120 OTN(L,NA,NG) = ABC/(1.+BC*EXP(AMIN1(HLE(L )/C,80.))) - 2540. - ABC/(1.+BC*EXP(AMIN1(HLE(L+1)/C,80.))) 2541. WRITE (6,910) NG,NA, (OTN(L,NA,NG),L=1,LMR) 2542. 130 continue 2543. C**** Tropospheric sulphate aerosol (NA=4,NG=4), layers 1 to 3 2544. IF(AGOLDH(4,4).le.0. .or. NGOLDH.lt.4) GO TO 200 2545. DO 140 L=1,3 2546. 140 OTN(L,4,4) = AGOLDH(4,4)/3. 2547. WRITE (6,910) 4,4, (OTN(L,4,4),L=1,3) 2547.1 C**** Sea salt aerosol (NA=7,NG=5), layers 1 to 3 2547.2 IF(AGOLDH(7,5).le.0. .or. NGOLDH.lt.5) GO TO 200 2547.3 OTN(1,7,5) = AGOLDH(7,5)*.500 2547.4 OTN(2,7,5) = AGOLDH(7,5)*.375 2547.5 OTN(3,7,5) = AGOLDH(7,5)*.125 2547.6 WRITE (6,910) 5,7, (OTN(L,7,5),L=1,3) 2548. C**** 2549. C**** THERMAL: SET (5) AEROSOL TYPE COMPOSITIONS & VERTICAL DISTRIBUTION 2550. C**** 2551. 200 DO 210 NG=1,NGOLDH 2552. DO 210 K=1,NKTR 2553. DO 210 L=1,LMR 2554. 210 TRAX(L,K,NG) = 0. 2555. C**** Dry aerosols 2556. DO 230 NG=1,NGOLDH 2557. DO 230 NA=1,NAERO 2558. IF(AGOLDH(NA,NG).lt.1.e-6) GO TO 230 2559. DO 220 K=1,NKTR 2560. DO 220 L=1,LMR 2561. 220 TRAX(L,K,NG) = TRAX(L,K,NG) + 2562. + OTN(L,NA,NG)*(TRAQEX(K,NA)-TRACOS(K,NA)*TRAQSC(K,NA)) 2563. 230 continue 2564. C**** Cloud aerosols 2565. DO 250 J=1,2 2566. DO 250 K=1,NKTR 2567. 250 TRCX(K,J) = TRCQEX(K,J)-TRCCOS(K,J)*TRCQSC(K,J) 2568. C**** 2569. C**** SOLAR: SET (5) AEROSOL TYPE COMPOSITIONS & VERTICAL DISTRIBUTION 2570. C**** 2571. DO 310 NG=1,NGOLDH 2572. DO 310 K=1,NKSR 2573. DO 310 L=1,LMR 2574. SRAX(L,K,NG) = 1.e-30 2575. SRAS(L,K,NG) = 1.e-30 2576. 310 SRAC(L,K,NG) = 0. 2577. C**** Dry aerosols 2578. DO 330 NG=1,NGOLDH 2579. DO 330 NA=1,NAERO 2580. IF(AGOLDH(NA,NG).lt.1.e-6) GO TO 330 2581. DO 320 K=1,NKSR 2582. DO 320 L=1,LMR 2583. SRAX(L,K,NG) = SRAX(L,K,NG) + OTN(L,NA,NG)*SRAQEX(K,NA) 2584. SRAS(L,K,NG) = SRAS(L,K,NG) + OTN(L,NA,NG)*SRAQSC(K,NA) 2585. 320 SRAC(L,K,NG) = SRAC(L,K,NG) + 2586. + OTN(L,NA,NG)*SRACOS(K,NA)*SRAQSC(K,NA) 2587. 330 continue 2588. C**** 2589. DO 350 NG=1,NGOLDH 2590. DO 350 K=1,NKSR 2591. DO 350 L=1,LMR 2592. 350 SRAC(L,K,NG) = SRAC(L,K,NG) / SRAS(L,K,NG) 2593. RETURN 2594. C**** 2595. C**** 2596. C**** 2597. ENTRY AEROS 2598. C**** 2599. C**** CLOUD AEROSOL AMOUNTS & DISTRIBUTIONS 2600. C**** 2602. LTOPCL=0 2603. DO 430 L=1,LMR 2604. KCLD=1 2605. IF(TLM(L).lt.TKCICE) KCLD=2 2609. IF(CLDTAU(L).gt..1) LTOPCL=L 2610. C**** THERMAL 2611. DO 420 K=1,NKTR 2612. SUMEXT = 1.e-30 2613. DO 410 NG=1,NGOLDH 2614. 410 SUMEXT = SUMEXT + FGOLDH(NG)*TRAX(L,K,NG) 2615. TRAEXT(L,K) = SUMEXT + CLDTAU(L)*TRCX(K,KCLD)*FCLDTR 2616. 420 TAUN(L,K) = TAUN(L,K) + TRAEXT(L,K) 2617. 430 continue 2618. C**** 2619. C**** CLOUD ALBEDO & SURFACE LAYER FOG SPECIFICATION 2620. C**** 2621. C DO 510 K=1,NKTR 2622. C 510 FTAUSL(K) = FOGTSL*TRCX(K,1)*FCLDTR 2623. IF(LTOPCL.gt.0) GO TO 530 2624. C**** No clouds 2625. DO 520 K=1,NKTR 2626. 520 TRCALB(K) = 0. 2627. GO TO 600 2628. C**** Clouds: set albedoes 2629. 530 KCLD=1 2630. IF(TLM(LTOPCL).lt.TKCICE) KCLD=2 2631. DO 540 K=1,NKTR 2632. 540 TRCALB(K) = (1.-EXP(-CLDTAU(LTOPCL)*TRCX(K,KCLD)))*CLDALB(K,KCLD) 2633. * *ECLTRA*FCLDTR 2634. C**** 2635. C**** SOLAR 2636. C**** 2637. 600 DO 620 K=1,NKSR 2638. DO 620 L=1,LMR 2639. EXTSUM = 1.e-30 2640. SCTSUM = 1.e-30 2641. COSSUM = 0. 2642. DO 610 NG=1,NGOLDH 2643. EXTSUM = EXTSUM + FGOLDH(NG)*SRAX(L,K,NG) 2644. SCTSUM = SCTSUM + FGOLDH(NG)*SRAS(L,K,NG) 2645. 610 COSSUM = COSSUM + FGOLDH(NG)*SRAS(L,K,NG)*SRAC(L,K,NG) 2646. EXTAER(L,K) = EXTSUM 2647. SCTAER(L,K) = SCTSUM 2648. 620 COSAER(L,K) = COSSUM/SCTSUM 2649. RETURN 2650. C**** 2651. 800 WRITE (6,*) ' Either NGOLDH or NAERO exceeds internal', 2652. * ' dimentions of OTN. NGOLDH,NAERO=',NGOLDH,NAERO 2653. STOP 'AEROS0 800:' 2654. C**** 2655. 900 FORMAT ( 2656. *'0AEROS0 NGOLD NAERO' / 2657. *' ----- -----' / 2658. *' 1. Volcano 1. Stratospheric sulphate' / 2659. *' 2. Ocean 3. Natural tropospheric sulphate' / 2660. *' 3. Land 4. Anthropogenic tropospheric sulphate' / 2661. *' 4. Sulphate 5. Basalt 1' / 2662. *' 5. Sea salt 6. Basalt 2' / 2663. *' 7. Sea salt' / 2664. *' 8. Dust 1' / 2665. *' 11. Carbon' // 2666. * ' Edge P: ',20I6 / (10X,20I6)) 2667. 901 FORMAT (' Edge H: ',20F6.1 / (10X,20F6.1)) 2668. 902 FORMAT (' Layers: ',20I6 / (12X,20I6)) 2669. 903 FORMAT ( 2670. *' NGOLD NAERO Normalized optical thickness for each layer' / 2671. *' ----- ----- -------------------------------------------') 2672. 910 FORMAT (I4,I7,2X,20(1X,F5.4) / (13X,20(1X,F5.4))) 2673. END 3000. 3001. SUBROUTINE GASES0 3002. C**** 3003. C**** GASES0 calculates gas amounts and calls TAUGAS. 3004. C**** Output: ULGAS = gas amounts for 9 different gasses 3005. C**** 3006. INCLUDE 'RA99M9.COM' 3006.4 C**** Working common block 3006.5 COMMON /WORK04/ PL(40),DPL(40), 3006.6 T TRCALB(40),TRGALB(40),BGFEMT(40),BGFEMD(40), 3006.7 T TRAEXT(40,25),TAUN(LMR,NKTR),TAUSL(25),FTAUSL(25), 3006.8 T DFLB(40,25),UFLB(40,25),WFLB(40,25),DFSL(25),UFSL(25),WFSL(25), 3006.9 S EXTAER(40,6),SCTAER(40,6),COSAER(40,6) 3007. REAL*4 SPLB(8),STLB(8),SHLB(8),SDLB(8), Z0(9),ZH(9), 3008. * PLE0(LMR+1),TLM0(LMR), TAUN0(LMR,NKTR) 3009. LOGICAL*4 Q2F(15),Q2V(15) 3010. DATA HPCON /34.16319/, P0 /1013.25/ 3011. DATA Q2F /.FALSE.,.FALSE.,.FALSE.,.FALSE.,.FALSE., 3011.1 * .TRUE. ,.TRUE. ,.TRUE. ,.TRUE. ,.TRUE. , 3011.2 * .TRUE. ,.TRUE. ,.TRUE. ,.FALSE.,.FALSE./ 3012. DATA Q2V /.TRUE. ,.TRUE. ,.TRUE. ,.TRUE. ,.TRUE. , 3012.1 * .FALSE.,.FALSE.,.FALSE.,.FALSE.,.FALSE., 3012.2 * .FALSE.,.FALSE.,.FALSE.,.TRUE. ,.TRUE. / 3013. C**** GLOBAL U.S. (1976) STANDARD ATMOSPHERE P, T, GEO H PARAMETERS 3014. DATA SPLB,STLB,SHLB,SDLB / 3015. * 1013.25,226.32,54.748,8.6801, 1.109,.66938,.039564,3.7338E-03, 3016. * 288.15,216.65,216.65,228.65,270.65,270.65, 214.65,186.87 , 3017. * 0. , 11. , 20. , 32. , 47. , 51. , 71. , 84.852, 3018. * -6.5 , 0. , 1. , 2.8 , 0. , -2.8 , -2. , 0./ 3019. C**** Gasses: H2O CO2 O3 O2 NO2 N2O CH4 F11 F12 3020. C$ DATA Z0 / 0., 0., 0., 0., 0., 15., 10., 12., 12./ 3021. DATA Z0 / 0., 0., 0., 0., 0.,915.,910., 12., 12./ 3022. DATA ZH / 8., 8., 8., 8., 8., 10., 30., 3., 3./ 3023. C**** 3024. C**** Use initial PLE to fix standard heights for gas distributions 3025. C**** 3026. byPDYN = 1. / (PLE(1)-PLE(LMA+1)) 3027. PTOP = PLE(LMA+1) 3028. DO 110 L=1,LMR 3029. DPL(L) = PLE(L)-PLE(L+1) 3030. 110 PL(L) = (PLE(L)+PLE(L+1))*.5 3031. C**** 3032. DO 160 L=1,LMR 3033. P = PLE(L) 3034. DO 120 N=2,8 3035. IF(P.GT.SPLB(N)) GO TO 130 3036. 120 CONTINUE 3037. N=9 3038. 130 N=N-1 3039. IF(ABS(SDLB(N)).LT.1.E-04) GO TO 140 3040. H = SHLB(N)+STLB(N)/SDLB(N)*((SPLB(N)/P)**(SDLB(N)/HPCON)-1.) 3041. GO TO 150 3042. 140 H = SHLB(N)+STLB(N)/HPCON*ALOG(SPLB(N)/P) 3043. 150 T = STLB(N)+SDLB(N)*(H-SHLB(N)) 3044. TLB(L) = T 3045. 160 HLE(L) = H 3046. HLE(1) = 1.E-10 3047. HLE(LMR+1) = 99.99 3048. TLB(LMR+1) = STLB(8) 3049. DO 170 L=1,LMR 3050. TLT(L) = TLB(L+1) 3051. 170 TLM(L) = .5*(TLB(L)+TLT(L)) 3052. C**** 3053. C**** DEFINE GLOBAL MEAN GAS AMOUNTS FOR TRACEGAS 3054. C**** 3055. C**** Water vapor 3056. RHP = .77 3057. EST = 10.**(9.4051-2353./TLB(1)) 3058. FWB = .662*RHP*EST/(PLE(1)-RHP*EST) 3059. DO 211 L=1,LMR 3060. PLT = PLE(L+1) 3061. DP = PLE(L)-PLT 3062. RHP = .77*(PLT/P0-.02)/.98 3063. EST = 10.**(9.4051-2353./TLT(L)) 3064. FWT = .662*RHP*EST/(PLT-RHP*EST) 3065. IF(FWT.GT.3.E-06) GO TO 210 3066. FWT = 3.E-06 3067. RHP = FWT*PLT/(EST*(FWT+0.662)) 3068. 210 ULGASL = .5*(FWB+FWT)*DP*1270. 3069. C$210 ULGASL = .5*(FWB+FWT)*DP*1268.75 3070. U0GAS(L,1) = ULGASL 3071. QL(L) = ULGASL/(ULGASL+1268.75*DP) 3072. 211 FWB = FWT 3073. C**** O2 is uniformly mixed 3075. DO 230 L=1,LMR 3076. 230 U0GAS(L,4) = PPMV(4)*.8*DPL(L)/P0 3077. CALL SETAO2 (U0GAS(1,4),LMR) 3078. C**** Ozone 3079. CALL OZONE0 3080. XX = SETAO3 (OCM) 3081. C**** NO2 3082. HDN = 0. 3083. FDN = CMANO2(0) 3084. M=1 3085. DO 254 L=1,LMR 3086. FSM = 0. 3087. HUP = HLE(L+1) 3088. 251 IF(2.*M.GE.HUP) GO TO 253 3089. FSM = FSM + (2.*M-HDN)*(CMANO2(M)+FDN)*.5 3090. HDN = 2.*M 3091. FDN = CMANO2(M) 3092. M=M+1 3093. IF(M.LE.41) GO TO 251 3094. U0GAS(L,5) = FSM 3095. DO 252 LS=L+1,LMR 3096. 252 U0GAS(LS,5) = 0. 3097. GO TO 260 3098. 253 FUP = CMANO2(M) + (CMANO2(M-1)-CMANO2(M))*(2.*M-HUP)/2. 3099. C$ U0GAS(L,5) = FSM + (HUP-HDN)*(FUP+FDN)*.5 correct 3100. U0GAS(L,5) = FSM + (HUP-HDN)*(CMANO2(M)+FDN)*.5 3101. HDN = HUP 3102. 254 FDN = FUP 3102.1 C**** Save PLE and TLM for future calls to GASESY 3102.2 260 DO 290 L=1,LMR 3102.3 PLE0(L) = PLE(L) 3102.4 290 TLM0(L) = TLM(L) 3102.5 PLE0(LMR+1) = PLE(LMR+1) 3102.6 RETURN 3102.7 C**** 3102.8 C**** 3102.9 C**** 3103. ENTRY GASESY 3103.1 C**** 3103.2 C**** Restore initial values of PLE, DPL, PL and TLM from GASES0 3103.3 DO 310 L=1,LMR 3103.4 PLE(L) = PLE0(L) 3103.5 DPL(L) = PLE0(L) - PLE0(L+1) 3103.6 PL(L) = (PLE0(L) + PLE0(L+1))*.5 3103.7 310 TLM(L) = TLM0(L) 3103.8 C**** CO2 is uniformly mixed 3103.9 DO 320 L=1,LMR 3104. 320 U0GAS(L,2) = PPMV(2)*.8*DPL(L)/P0 3104.1 C**** Restore initial ozone distribution 3104.2 CALL OZONE0 3104.3 C**** Specify vertical gas distribution of N20, CH4, F11 and F12 3104.4 DO 330 N=6,9 3105. DO 330 L=1,LMR 3106. U0GAS(L,N) = PPMV(N)*.8*DPL(L)/P0 3107. ZT = (HLE(L+1)-Z0(N))/ZH(N) 3108. IF(ZT.le.0.) GO TO 330 3109. ZB = (HLE(L)-Z0(N))/ZH(N) 3110. EXPZT = EXP(-ZT) 3111. EXPZB = EXP(-ZB) 3112. IF(ZB.lt.0.) EXPZB = 1.-ZB 3113. U0GAS(L,N) = U0GAS(L,N)*(EXPZB-EXPZT)/(ZT-ZB) 3114. 330 continue 3115. C**** 3116. DO 340 N=1,9 3117. DO 340 L=1,LMR 3118. 340 ULGAS(L,N) = U0GAS(L,N) 3119. C**** 3120. C**** Calculate optical depths for fixed gases that overlap 3121. C**** SEC: H2O CO2 O3 H2O H2O CO2 CO2 O3 N2O N2O CH4 F11 F12 H2O H2O 3122. C**** PRI: H2O CO2 O3 CO2 O3 H2O O3 CO2 H2O CO2 H2O WIN 3123. C**** Q2F: T T T T T T T T 3124. C**** 3125. DO 350 L=1,LMR*NKTR 3126. 350 TAUN(L,1) = 0. 3127. CALL TAUGAS (Q2F) 3128. DO 360 L=1,LMR*NKTR 3129. 360 TAUN0(L,1) = TAUN(L,1) 3130. RETURN 3131. C**** 3132. C**** 3133. ENTRY GASES 3134. C**** 3135. C**** Calculate time varying gas amounts 3136. C**** 3137. DO 410 L=1,LMR 3138. DPL(L) = PLE(L)-PLE(L+1) 3139. 410 PL(L) = (PLE(L)+PLE(L+1))*.5 3140. C**** Water vapor 3141. DO 420 L=1,LMR 3142. 420 ULGAS(L,1) = 1268.75*DPL(L)*QL(L) 3143. C$420 ULGAS(L,1) = 1268.75*DPL(L)*QL(L)/(1.-QL(L)) ********CORRECT 3144. C**** Ozone 3145. C CALL OZONE 3146. DO 430 L=1,LMR 3147. 430 ULGAS(L,3) = U0GAS(L,3) 3148. C**** CO2, O2, N2O, CH4, F11, F12 vary in time below radiation layers 3149. RMASS = (PLE(1)-PTOP) * byPDYN 3150. DO 440 L=1,LMA 3151. ULGAS(L,2) = U0GAS(L,2)*RMASS 3152. ULGAS(L,4) = U0GAS(L,4)*RMASS 3153. ULGAS(L,6) = U0GAS(L,6)*RMASS 3154. ULGAS(L,7) = U0GAS(L,7)*RMASS 3155. ULGAS(L,8) = U0GAS(L,8)*RMASS 3156. 440 ULGAS(L,9) = U0GAS(L,9)*RMASS 3157. C**** 3158. C**** Calculate optical depths for time varying gases that overlap 3159. C**** SEC: H2O CO2 O3 H2O H2O CO2 CO2 O3 N2O N2O CH4 F11 F12 H2O H2O 3160. C**** PRI: H2O CO2 O3 CO2 O3 H2O O3 CO2 H2O CO2 H2O WIN CON 3161. C**** Q2V: T T T T T T T 3162. C**** 3163. DO 510 L=1,LMR*NKTR 3164. 510 TAUN(L,1) = TAUN0(L,1) 3165. CALL TAUGAS (Q2V) 3166. C**** 3167. RETURN 3168. END 3500. 3501. SUBROUTINE TAUGAS (Q2) 3502. C**** 3503. C**** TAUGAS calulates minus the logarithm of the fraction of the 3504. C**** radiative energy transmitted through each layer for each 3505. C**** K-interval. 3506. C**** 3507. C**** Input: LMR,TLM,ULGAS,TRACEG,PL,DPL,TAUTBL 3508. C**** Output: TAUN = -ln F(trans)/F(input) 3509. INCLUDE 'RA99M9.COM' 3509.4 C**** Working common block 3509.5 COMMON /WORK04/ PL(40),DPL(40), 3509.6 T TRCALB(40),TRGALB(40),BGFEMT(40),BGFEMD(40), 3509.7 T TRAEXT(40,25),TAUN(LMR,NKTR),TAUSL(25),FTAUSL(25), 3509.8 T DFLB(40,25),UFLB(40,25),WFLB(40,25),DFSL(25),UFSL(25),WFSL(25), 3509.9 S EXTAER(40,6),SCTAER(40,6),COSAER(40,6) 3510. PARAMETER (NPX=15,NTX=8,NGUX=652, TLOX=181.,DTX=23.) 3511. INTEGER*4 IGASS(11),K0(11),NK(11),ITBL0(11),NTBL(11) 3512. REAL*4 ULOX(NPX,11),DUX(NPX,11),PX(NPX),H2OCON(25) 3513. LOGICAL*4 Q2(15) 3514. C**** 3515. C**** Secondary: H2O CO2 O3 H2O H2O CO2 CO2 O3 N2O N2O CH4 3516. C**** Principal: H2O CO2 O3 CO2 O3 H2O O3 CO2 H2O CO2 H2O 3517. DATA IGASS / 1, 2, 3, 1, 1, 2, 2, 3, 6, 6, 7/, 3518. C * IGASP / 1, 2, 3, 2, 3, 1, 3, 2, 1, 2, 1/, 3519. * K0 / 1, 11, 21, 11, 21, 1, 21, 11, 1, 11, 1/, 3520. * NK /10, 10, 4, 10, 4, 10, 4, 10, 10, 10, 10/, 3521. * ITBL0 / 0,250,340,376,466,502,552,572,622,632,642/, 3522. * NTBL /25, 9, 9, 9, 9, 5, 5, 5, 1, 1, 1/ 3523. C**** 3524. DATA PX /1000., 975., 910., 800., 645., 480., 330., 3525. * 205., 110., 40., 7.5, 3.5, 1.0, 0.1, .001/ 3526. DATA ULOX/ .25E+2,.25E+2,.50E+2,.50E+2,.25E+2,.50E+1,.10E+1, 3527. *.25E+0,.10E+0,.50E-1,.10E-1,.10E-1,.10E-3,.10E-5,.10E-5, 3528. *.50E+1,.50E+1,.80E+1,.10E+2,.20E+2,.20E+2,.10E+1,.10E+2,.80E+1, 3529. *.10E+1,.10E+1,.25E+0,.50E-1,.50E-2,.50E-3, .10E-3,.10E-3, 3530. *.40E-3,.60E-3,.10E-2,.24E-2,.48E-2,.48E-2,.64E-2,.64E-2,.64E-2, 3531. *.40E-2,.10E-4,.80E-7,.40E-7, .25E+2,.25E+2,.50E+2,.50E+2, 3532. *.25E+2,.50E+1,.10E+1,.25E+0,.10E+0,.50E-1,.10E-1,.10E-1,.10E-3, 3533. *.10E-5,.10E-5, .25E+2,.25E+2,.50E+2,.50E+2,.25E+2,.50E+1, 3534. *.10E+1,.25E+0,.10E+0,.50E-1,.10E-1,.10E-1,.10E-3,.10E-5,.10E-5, 3535. * .50E+1,.50E+1,.80E+1,.10E+2,.20E+2,.20E+2,.10E+2,.10E+2, 3536. *.80E+1,.10E+1,.10E+1,.25E+0,.50E-1,.50E-2,.50E-3, .50E+1, 3537. *.50E+1,.80E+1,.10E+2,.20E+2,.20E+2,.10E+2,.10E+2,.80E+1,.10E+1, 3538. *.10E+1,.25E+0,.50E-1,.50E-2,.50E-3, .10E-3,.10E-3,.40E-3, 3539. *.60E-3,.10E-2,.24E-2,.48E-2,.48E-2,.64E-2,.64E-2,.64E-2,.40E-2, 3540. *.10E-4,.80E-7,.40E-7, .11E-1,.11E-1,.18E-1,.31E-1,.37E-1, 3541. *.35E-1,.31E-1,.24E-1,.18E-1,.13E-1,.11E-2,.66E-3,.44E-3,.44E-4, 3542. *.44E-6, .11E-1,.11E-1,.18E-1,.31E-1,.37E-1,.35E-1,.31E-1, 3543. *.24E-1,.18E-1,.13E-1,.11E-2,.66E-3,.44E-3,.44E-4,.44E-6, 3544. *.64E-1,.64E-1,.10E+0,.18E+0,.22E+0,.20E+0,.18E+0,.14E+0,.10E+0, 3545. *.77E-1,.64E-2,.38E-2,.26E-2,.26E-3,.26E-5/ 3546. C**** 3547. DATA DUX/ .75E+2,.75E+2,.10E+3,.10E+3,.75E+2,.50E+2,.10E+2, 3548. *.20E+1,.20E+0,.10E+0,.50E-1,.10E-1,.40E-2,.40E-3,.40E-4, 3549. *.50E+1,.50E+1,.80E+1,.10E+2,.10E+2,.10E+2,.10E+2,.10E+2,.80E+1, 3550. *.50E+1,.35E+1,.25E+0,.25E+0,.10E+0,.10E-1, .30E-3,.30E-3, 3551. *.50E-3,.80E-3,.10E-2,.16E-2,.64E-2,.16E-2,.25E-1,.25E-1,.25E-1, 3552. *.45E-2,.25E-2,.10E-2,.25E-4, .24E+3,.24E+3,.30E+3,.30E+3, 3553. *.24E+3,.15E+3,.30E+2,.60E+1,.60E+0,.30E+0,.15E+0,.30E-1,.12E-1, 3554. *.12E-2,.12E-3, .24E+3,.24E+3,.30E+3,.30E+3,.24E+3,.15E+3, 3555. *.30E+2,.60E+1,.60E+0,.30E+0,.15E+0,.30E-1,.12E-1,.12E-2,.12E-3, 3556. * .10E+2,.10E+2,.16E+2,.20E+2,.20E+2,.20E+2,.20E+2,.20E+2, 3557. *.16E+2,.10E+2,.70E+1,.50E+0,.50E+0,.20E+0,.20E-1, .10E+2, 3558. *.10E+2,.16E+2,.20E+2,.20E+2,.20E+2,.20E+2,.20E+2,.16E+2,.10E+2, 3559. *.70E+1,.50E+0,.50E+0,.20E+0,.20E-1, .60E-3,.60E-3,.10E-2, 3560. *.16E-2,.20E-2,.32E-2,.13E-1,.32E-1,.50E-1,.50E-1,.50E-1,.90E-2, 3561. *.50E-2,.20E-2,.50E-4, 45*0./ 3562. C**** 3563. DATA H2OCON/ .767116, .322401, .572299,.58537, .48869, 3564. * .43539, .44322, .64072, .89293, 1.12733,1.65550, .865210, 3565. * 1.38403,1.80159,1.99196, 2.03403, 2.20561,2.42859,2.56883, 3566. * 2.67157,2.71888, .45534, .44735, .44534, .44365/ 3567. C**** 3568. C**** ABSORPTION (TAU) INTERPOLATION FOR GAS AMOUNTS IN ULGAS(N,K) 3569. C**** 3570. IPX=2 3571. DO 300 L=1,LMR 3572. C**** Locate model layer pressure between IPX and IPX-1 3573. 20 WPB = (PL(L)-PX(IPX))/(PX(IPX-1)-PX(IPX)) 3574. IF(WPB.GE.0. .OR. IPX.GE.NPX) GO TO 30 3575. IPX = IPX+1 3576. GO TO 20 3577. C**** Locate model layer temperature between ITX and ITX+1 3578. 30 WTB = (TLM(L)-TLOX)/DTX + 1. 3579. ITX = WTB 3580. IF(ITX.LE.0) ITX = 1 3581. IF(ITX.GE.NTX) ITX = NTX-1 3582. WTB = WTB-ITX 3583. C**** 3584. WBB = WPB*WTB 3585. WBA = WPB-WBB 3586. WAB = WTB-WBB 3587. WAA = 1.-(WBB+WBA+WAB) 3588. C**** 3589. C**** N2=1,8: absorption is non-linear function of gas amount 3590. C**** 3591. DO 180 N2=1,8 3592. IF(.NOT.Q2(N2)) GO TO 180 3593. UGAS = ULGAS(L,IGASS(N2)) 3594. IF(UGAS.LT.1.E-10) GO TO 180 3595. C**** 3596. NU = NTBL(N2) 3597. IF(UGAS.LE.ULOX(IPX,N2)) GO TO 110 3598. IF(UGAS.GE.ULOX(IPX,N2)+DUX(IPX,N2)*(NU-1)) GO TO 120 3599. C**** UGAS(IPX) is located within the interpolation limits 3600. XUA = (UGAS-ULOX(IPX,N2))/DUX(IPX,N2) 3601. IUA = XUA 3602. WAAA = WAA*(IUA+1-XUA) 3603. WABA = WAB*(IUA+1-XUA) 3604. WAAB = WAA*(XUA-IUA) 3605. WABB = WAB*(XUA-IUA) 3606. GO TO 130 3607. C**** UGAS(IPX) is less than the minimum interpolation value 3608. 110 IUA = 0 3609. WAAA = WAA*(UGAS/ULOX(IPX,N2)) 3610. WABA = WAB*(UGAS/ULOX(IPX,N2)) 3611. WAAB = 0. 3612. WABB = 0. 3613. GO TO 130 3614. C**** UGAS(IPX) exceeds the maximum interpolation value 3615. 120 IUA = NU-2 3616. WAAA = 0. 3617. WABA = 0. 3618. WAAB = WAA*(UGAS/(ULOX(IPX,N2)+DUX(IPX,N2)*(NU-1))) 3619. WABB = WAB*(UGAS/(ULOX(IPX,N2)+DUX(IPX,N2)*(NU-1))) 3620. C**** 3621. 130 IF(UGAS.LE.ULOX(IPX-1,N2)) GO TO 140 3622. IF(UGAS.GE.ULOX(IPX-1,N2)+(NU-1)*DUX(IPX-1,N2)) GO TO 150 3623. C**** UGAS(IPX-1) is located within the interpolation limits 3624. XUB = (UGAS-ULOX(IPX-1,N2))/DUX(IPX-1,N2) 3625. IUB = XUB 3626. WBAA = WBA*(IUB+1-XUB) 3627. WBBA = WBB*(IUB+1-XUB) 3628. WBAB = WBA*(XUB-IUB) 3629. WBBB = WBB*(XUB-IUB) 3630. GO TO 160 3631. C**** UGAS(IPX-1) is less than the minimum interpolation value 3632. 140 IUB = 0 3633. WBAA = WBA*(UGAS/ULOX(IPX-1,N2)) 3634. WBBA = WBB*(UGAS/ULOX(IPX-1,N2)) 3635. WBAB = 0. 3636. WBBB = 0. 3637. GO TO 160 3638. C**** UGAS(IPX-1) exceeds the maximum interpolation value 3639. 150 IUB = NU-2 3640. WBAA = 0. 3641. WBBA = 0. 3642. WBAB = WBA*(UGAS/(ULOX(IPX-1,N2)+DUX(IPX-1,N2)*(NU-1))) 3643. WBBB = WBB*(UGAS/(ULOX(IPX-1,N2)+DUX(IPX-1,N2)*(NU-1))) 3644. C**** 3645. 160 NKJ = NK(N2) 3646. IAAA = ITBL0(N2) + NKJ*IUA 3647. IAAB = IAAA + NKJ 3648. IBAA = ITBL0(N2) + NKJ*IUB 3649. IBAB = IBAA + NKJ 3650. C**** 3651. DO 170 K=1,NKJ 3652. 170 TAUN(L,K+K0(N2)) = TAUN(L,K+K0(N2)) 3653. * + WAAA*TAUTBL(K+IAAA,ITX ,IPX ) 3654. * + WAAB*TAUTBL(K+IAAB,ITX ,IPX ) 3655. * + WABA*TAUTBL(K+IAAA,ITX+1,IPX ) 3656. * + WABB*TAUTBL(K+IAAB,ITX+1,IPX ) 3657. * + WBAA*TAUTBL(K+IBAA,ITX ,IPX-1) 3658. * + WBAB*TAUTBL(K+IBAB,ITX ,IPX-1) 3659. * + WBBA*TAUTBL(K+IBAA,ITX+1,IPX-1) 3660. * + WBBB*TAUTBL(K+IBAB,ITX+1,IPX-1) 3661. 180 CONTINUE 3662. C**** 3663. C**** N2=9,11: absorption is a linear function of gas amount 3664. C**** 3665. DO 280 N2=9,11 3666. IF(.NOT.Q2(N2)) GO TO 280 3667. UGAS = ULGAS(L,IGASS(N2)) 3668. IF(UGAS.LT.1.E-10) GO TO 280 3669. C**** 3670. WAAB = WAA*(UGAS/ULOX(IPX ,N2)) 3671. WABB = WAB*(UGAS/ULOX(IPX ,N2)) 3672. WBAB = WBA*(UGAS/ULOX(IPX-1,N2)) 3673. WBBB = WBB*(UGAS/ULOX(IPX-1,N2)) 3674. C**** 3675. DO 270 K=1,NK(N2) 3676. 270 TAUN(L,K+K0(N2)) = TAUN(L,K+K0(N2)) 3677. * + WAAB*TAUTBL(K+ITBL0(N2),ITX ,IPX ) 3678. * + WABB*TAUTBL(K+ITBL0(N2),ITX+1,IPX ) 3679. * + WBAB*TAUTBL(K+ITBL0(N2),ITX ,IPX-1) 3680. * + WBBB*TAUTBL(K+ITBL0(N2),ITX+1,IPX-1) 3681. 280 CONTINUE 3682. 300 CONTINUE 3683. C**** 3684. C**** PICK UP CCL3F1 (F11) ABSORPTION 3685. C**** 3686. IF(.NOT.Q2(12)) GO TO 410 3687. DO 402 K=1,25 3688. IF(TRACEG(K,1).LT.1.E-10) GO TO 402 3689. DO 401 L=1,LMR 3690. 401 TAUN(L,K) = TAUN(L,K) + ULGAS(L,8)*TRACEG(K,1) 3691. 402 CONTINUE 3692. C**** 3693. C**** PICK UP CCL2F2 (F12) ABSORPTION 3694. C**** 3695. 410 IF(.NOT.Q2(13)) GO TO 420 3696. DO 412 K=1,25 3697. IF(TRACEG(K,2).LT.1.E-10) GO TO 412 3698. DO 411 L=1,LMR 3699. 411 TAUN(L,K) = TAUN(L,K) + ULGAS(L,9)*TRACEG(K,2) 3700. 412 CONTINUE 3701. C**** 3702. C**** PICK UP WINDOW H2O GASEOUS ABSORPTION 3703. C**** 3704. C 420 IF(.NOT.Q2(14)) GO TO 430 3705. C DO 421 L=1,LMR 3706. C 421 TAUN(L,1) = TAUN(L,1) 3707. C**** 3708. C**** PICK UP H2O CONTINUUM ABSORPTION 3709. C**** 3710. 420 IF(.NOT.Q2(15)) GO TO 440 3711. DO 431 L=1,LMR 3712. 431 TAUN(L,1) = TAUN(L,1) + 2.21866E-11* 3713. * PL(L)*ULGAS(L,1)*EXP(1800./TLM(L))* 3714. * (ULGAS(L,1)/DPL(L)+.808563) 3715. C$*** REMOVE FOLLOWING STATEMENT TO CORRECT 3716. C IF(LMR.GT.0) RETURN 3717. C DO 433 L=1,LMR 3718. C PH2O = 12.38E-4*ULGAS(L,1)*PL(L)/DPL(L) 3719. C TH2O = EXP(1800./TLM(L)-6.081081) 3720. C COEC = PH2O*TH2O + .0015*(PL(L)-PH2O) 3721. C DO 432 K=2,25 3722. C COEF = H2OCON(K)*1.E-5 3723. C 432 TAUN(L,K) = TAUN(L,K) + ULGAS(L,1)*COEC*COEF 3724. C 433 CONTINUE 3725. 440 RETURN 3726. END 4000. 4001. SUBROUTINE OZONE0 4002. INCLUDE 'RA99M9.COM' 4003. C LOGICAL SKIPI 4004. C 4005. C 4006. C----------------------------------------------------------------------- 4007. C 4008. C LONDON ET AL (1976) JUL,1957-DEC,1970 NCAR ATLAS OF TOTAL OZONE 4009. C 4010. C AVERAGE GLOBAL COLUMN AMOUNT -- O3AVE(MONTH,LATITUDE,LONGITUDE) 4011. C 4012. C MONTH=1-12 JAN,FEB,...,DEC 4013. C LAT =1-18 -85,-75,..., 85 4014. C 4015. C----------------------------------------------------------------------- 4016. PARAMETER (FLONO3=0.) 4017. REAL O3AVEA(216),O3AVEB(216),O3AVEC(216),O3AVED(216),O3AVEE(216) 4018. REAL O3AVEF(216),O3AVEG(216),O3AVEH(216),O3AVEI(216),O3AVEJ(216) 4019. REAL O3AVEK(216),O3AVEL(216),O3AVEM(216),O3AVEN(216),O3AVEO(216) 4020. REAL O3AVEP(216),O3AVEQ(216),O3AVER(216),O3AVE(12,18,18) 4021. EQUIVALENCE (O3AVE(1,1,10),O3AVEA(1)),(O3AVE(1,1,11),O3AVEB(1)) 4022. 1 ,(O3AVE(1,1,12),O3AVEC(1)),(O3AVE(1,1,13),O3AVED(1)) 4023. 2 ,(O3AVE(1,1,14),O3AVEE(1)),(O3AVE(1,1,15),O3AVEF(1)) 4024. 3 ,(O3AVE(1,1,16),O3AVEG(1)),(O3AVE(1,1,17),O3AVEH(1)) 4025. 4 ,(O3AVE(1,1,18),O3AVEI(1)),(O3AVE(1,1,01),O3AVEJ(1)) 4026. 5 ,(O3AVE(1,1,02),O3AVEK(1)),(O3AVE(1,1,03),O3AVEL(1)) 4027. 6 ,(O3AVE(1,1,04),O3AVEM(1)),(O3AVE(1,1,05),O3AVEN(1)) 4028. 7 ,(O3AVE(1,1,06),O3AVEO(1)),(O3AVE(1,1,07),O3AVEP(1)) 4029. 8 ,(O3AVE(1,1,08),O3AVEQ(1)),(O3AVE(1,1,09),O3AVER(1)) 4030. DATA O3AVEA/ 4031. A .317,.295,.291,.292,.293,.298,.300,.305,.313,.324,.369,.355, 4032. B .319,.300,.296,.292,.291,.300,.301,.304,.314,.322,.358,.350, 4033. C .312,.301,.295,.287,.286,.298,.302,.305,.316,.322,.343,.335, 4034. D .299,.291,.285,.280,.279,.290,.295,.300,.307,.319,.327,.316, 4035. E .281,.275,.279,.268,.266,.278,.282,.290,.295,.306,.306,.296, 4036. F .266,.261,.259,.256,.252,.261,.267,.277,.280,.289,.285,.277, 4037. G .252,.249,.248,.246,.240,.249,.252,.262,.264,.273,.265,.258, 4038. H .240,.238,.240,.242,.237,.242,.240,.249,.252,.258,.251,.245, 4039. I .232,.230,.238,.241,.240,.238,.234,.241,.241,.245,.239,.236, 4040. J .235,.235,.244,.252,.253,.244,.236,.237,.232,.230,.230,.232, 4041. K .249,.256,.264,.269,.267,.261,.245,.245,.238,.234,.233,.237, 4042. L .278,.289,.294,.300,.294,.284,.265,.265,.256,.249,.248,.261, 4043. M .318,.338,.343,.351,.342,.324,.300,.296,.287,.275,.279,.299, 4044. N .347,.368,.383,.383,.370,.351,.335,.319,.304,.288,.296,.321, 4045. O .364,.394,.418,.410,.402,.371,.358,.340,.312,.298,.302,.325, 4046. P .356,.388,.421,.414,.394,.360,.337,.319,.299,.285,.292,.313, 4047. Q .364,.403,.431,.426,.398,.358,.328,.303,.292,.287,.297,.324, 4048. R .373,.421,.447,.440,.408,.355,.323,.295,.289,.291,.305,.329/ 4049. DATA O3AVEB/ 4050. A .318,.295,.291,.293,.293,.299,.301,.305,.314,.326,.372,.358, 4051. B .321,.300,.295,.293,.291,.301,.301,.306,.314,.326,.361,.353, 4052. C .315,.302,.296,.291,.288,.300,.303,.306,.318,.328,.348,.340, 4053. D .307,.296,.291,.284,.278,.298,.299,.305,.314,.326,.335,.324, 4054. E .294,.285,.286,.272,.270,.286,.288,.296,.302,.315,.315,.304, 4055. F .278,.271,.265,.260,.258,.270,.273,.283,.287,.298,.293,.284, 4056. G .262,.259,.254,.250,.247,.255,.259,.268,.270,.282,.274,.266, 4057. H .247,.246,.244,.245,.239,.245,.247,.255,.255,.266,.257,.250, 4058. I .235,.235,.239,.244,.240,.238,.236,.244,.244,.249,.244,.239, 4059. J .233,.234,.243,.251,.249,.240,.234,.235,.232,.231,.231,.231, 4060. K .247,.254,.263,.267,.262,.253,.242,.240,.237,.232,.232,.237, 4061. L .279,.287,.296,.282,.286,.275,.260,.257,.253,.246,.246,.258, 4062. M .320,.336,.345,.348,.325,.309,.293,.282,.279,.267,.272,.294, 4063. N .346,.369,.379,.377,.348,.330,.317,.299,.286,.280,.288,.312, 4064. O .368,.406,.412,.401,.373,.345,.332,.312,.293,.284,.293,.316, 4065. P .366,.409,.423,.418,.386,.349,.326,.307,.290,.278,.295,.312, 4066. Q .366,.407,.428,.429,.396,.352,.323,.296,.287,.282,.298,.318, 4067. R .372,.420,.446,.441,.407,.352,.320,.292,.286,.290,.305,.327/ 4068. DATA O3AVEC/ 4069. A .319,.296,.292,.294,.294,.299,.302,.306,.316,.328,.372,.359, 4070. B .321,.300,.295,.297,.293,.303,.305,.309,.319,.332,.367,.359, 4071. C .322,.309,.302,.297,.293,.309,.309,.314,.326,.338,.362,.353, 4072. D .324,.313,.303,.294,.295,.314,.311,.318,.330,.342,.353,.343, 4073. E .315,.308,.296,.286,.287,.305,.306,.314,.326,.335,.338,.326, 4074. F .294,.290,.281,.271,.273,.287,.290,.299,.307,.319,.312,.303, 4075. G .274,.272,.264,.258,.258,.268,.272,.281,.286,.297,.290,.281, 4076. H .254,.254,.251,.248,.248,.254,.257,.263,.267,.276,.271,.262, 4077. I .240,.239,.241,.245,.241,.243,.244,.250,.251,.256,.250,.246, 4078. J .230,.231,.238,.249,.246,.237,.234,.233,.234,.233,.230,.228, 4079. K .238,.244,.251,.258,.253,.244,.236,.235,.233,.228,.228,.230, 4080. L .259,.269,.276,.279,.268,.254,.246,.241,.238,.235,.237,.246, 4081. M .289,.305,.312,.306,.289,.270,.261,.255,.249,.246,.252,.268, 4082. N .321,.347,.354,.343,.315,.291,.281,.273,.262,.259,.268,.285, 4083. O .351,.394,.396,.384,.353,.315,.300,.288,.275,.271,.282,.296, 4084. P .363,.414,.422,.415,.382,.333,.313,.292,.281,.276,.292,.306, 4085. Q .366,.415,.430,.433,.398,.346,.313,.288,.282,.280,.299,.317, 4086. R .372,.421,.445,.441,.406,.348,.316,.289,.285,.289,.306,.327/ 4087. DATA O3AVED/ 4088. A .320,.296,.293,.294,.295,.300,.303,.308,.317,.330,.374,.361, 4089. B .322,.300,.297,.299,.296,.307,.310,.314,.323,.339,.373,.366, 4090. C .329,.313,.310,.304,.302,.320,.318,.326,.338,.352,.373,.367, 4091. D .343,.330,.318,.306,.315,.333,.329,.337,.354,.366,.370,.366, 4092. E .334,.324,.311,.299,.312,.326,.329,.333,.352,.357,.354,.342, 4093. F .304,.300,.291,.279,.285,.302,.308,.315,.324,.328,.325,.312, 4094. G .277,.276,.268,.262,.266,.279,.283,.289,.296,.303,.299,.283, 4095. H .256,.257,.253,.249,.252,.259,.266,.269,.274,.278,.273,.263, 4096. I .242,.243,.243,.248,.247,.251,.255,.256,.258,.260,.253,.249, 4097. J .231,.234,.238,.250,.255,.251,.250,.246,.248,.244,.237,.229, 4098. K .235,.241,.248,.257,.259,.257,.248,.246,.245,.244,.233,.230, 4099. L .256,.261,.267,.270,.269,.262,.251,.247,.247,.248,.239,.248, 4100. M .293,.304,.306,.302,.288,.272,.259,.256,.256,.256,.254,.269, 4101. N .327,.344,.356,.346,.319,.291,.272,.270,.264,.267,.270,.285, 4102. O .356,.392,.402,.388,.359,.312,.289,.281,.276,.281,.285,.297, 4103. P .368,.416,.424,.415,.388,.328,.304,.285,.279,.284,.295,.309, 4104. Q .370,.418,.436,.436,.402,.338,.306,.283,.278,.284,.301,.320, 4105. R .373,.422,.446,.441,.407,.345,.312,.286,.275,.291,.307,.328/ 4106. DATA O3AVEE/ 4107. A .319,.295,.293,.295,.296,.300,.304,.309,.318,.332,.375,.362, 4108. B .325,.301,.300,.302,.300,.309,.313,.319,.328,.345,.378,.370, 4109. C .332,.314,.312,.310,.310,.327,.329,.335,.347,.362,.381,.375, 4110. D .348,.334,.324,.312,.328,.346,.366,.352,.372,.381,.377,.373, 4111. E .337,.327,.318,.303,.322,.335,.342,.347,.363,.366,.358,.344, 4112. F .301,.297,.292,.282,.291,.307,.314,.321,.331,.332,.324,.309, 4113. G .275,.271,.269,.264,.270,.279,.286,.292,.299,.301,.293,.281, 4114. H .255,.253,.252,.251,.253,.258,.265,.269,.275,.277,.268,.262, 4115. I .245,.244,.246,.250,.249,.253,.254,.257,.259,.260,.252,.249, 4116. J .240,.239,.245,.255,.256,.260,.256,.253,.253,.251,.243,.237, 4117. K .247,.248,.252,.263,.270,.268,.258,.256,.256,.252,.244,.238, 4118. L .263,.263,.268,.277,.282,.276,.261,.259,.259,.258,.251,.251, 4119. M .299,.304,.309,.309,.302,.291,.269,.266,.268,.269,.269,.275, 4120. N .346,.358,.365,.353,.335,.307,.276,.272,.276,.283,.289,.300, 4121. O .379,.400,.414,.401,.373,.319,.286,.280,.283,.293,.303,.314, 4122. P .382,.421,.437,.427,.398,.323,.293,.280,.280,.293,.308,.321, 4123. Q .375,.424,.444,.440,.405,.334,.298,.278,.276,.290,.306,.326, 4124. R .374,.424,.448,.443,.406,.345,.310,.284,.281,.292,.309,.328/ 4125. DATA O3AVEF/ 4126. A .318,.294,.294,.295,.298,.301,.304,.311,.320,.333,.377,.361, 4127. B .324,.298,.300,.304,.305,.310,.315,.323,.331,.348,.383,.371, 4128. C .337,.311,.314,.313,.317,.330,.333,.344,.354,.369,.386,.377, 4129. D .350,.330,.324,.317,.332,.349,.351,.362,.378,.390,.380,.372, 4130. E .333,.322,.314,.307,.323,.339,.345,.358,.369,.372,.357,.340, 4131. F .300,.292,.286,.284,.294,.307,.316,.327,.335,.334,.323,.307, 4132. G .275,.269,.264,.263,.269,.277,.285,.292,.300,.303,.290,.279, 4133. H .254,.251,.250,.251,.254,.256,.261,.267,.271,.276,.266,.261, 4134. I .243,.242,.242,.247,.248,.250,.247,.251,.252,.258,.253,.247, 4135. J .237,.239,.243,.253,.255,.255,.246,.243,.244,.245,.239,.236, 4136. K .246,.247,.253,.263,.265,.265,.253,.245,.247,.247,.239,.238, 4137. L .265,.265,.276,.283,.284,.280,.261,.254,.253,.258,.250,.250, 4138. M .306,.309,.321,.316,.318,.292,.273,.259,.265,.271,.273,.277, 4139. N .365,.369,.381,.363,.347,.313,.278,.264,.275,.290,.302,.307, 4140. O .396,.416,.431,.415,.405,.322,.282,.271,.288,.303,.321,.328, 4141. P .397,.433,.455,.436,.404,.322,.287,.273,.276,.302,.320,.333, 4142. Q .382,.429,.451,.442,.408,.331,.297,.274,.273,.295,.311,.333, 4143. R .375,.427,.450,.445,.407,.343,.309,.283,.280,.295,.311,.330/ 4144. DATA O3AVEG/ 4145. A .317,.293,.293,.295,.299,.299,.305,.311,.320,.335,.378,.360, 4146. B .323,.296,.300,.304,.306,.310,.317,.325,.334,.353,.385,.367, 4147. C .335,.307,.310,.312,.318,.328,.335,.347,.357,.376,.390,.372, 4148. D .346,.324,.320,.317,.332,.349,.354,.367,.384,.393,.384,.368, 4149. E .331,.318,.311,.305,.324,.339,.349,.365,.378,.377,.360,.339, 4150. F .301,.293,.286,.285,.296,.309,.321,.334,.344,.339,.325,.309, 4151. G .276,.270,.266,.267,.271,.280,.287,.295,.303,.308,.294,.282, 4152. H .257,.253,.250,.252,.254,.257,.261,.266,.271,.279,.268,.261, 4153. I .240,.241,.241,.246,.246,.250,.246,.249,.253,.259,.254,.248, 4154. J .234,.238,.245,.256,.258,.259,.244,.243,.241,.243,.237,.235, 4155. K .244,.249,.259,.271,.274,.274,.257,.251,.248,.248,.238,.237, 4156. L .270,.272,.289,.297,.298,.294,.277,.267,.260,.262,.251,.254, 4157. M .329,.338,.353,.338,.333,.313,.296,.275,.273,.282,.281,.296, 4158. N .401,.414,.424,.392,.369,.329,.298,.272,.282,.303,.321,.341, 4159. O .420,.451,.461,.432,.389,.331,.291,.272,.279,.313,.343,.358, 4160. P .411,.451,.468,.447,.403,.320,.289,.271,.277,.308,.334,.349, 4161. Q .386,.434,.456,.443,.404,.332,.297,.273,.273,.300,.317,.339, 4162. R .378,.430,.453,.446,.407,.342,.310,.282,.279,.296,.314,.332/ 4163. DATA O3AVEH/ 4164. A .315,.292,.293,.295,.299,.297,.303,.311,.320,.334,.378,.358, 4165. B .320,.294,.298,.303,.306,.308,.316,.325,.337,.355,.387,.362, 4166. C .330,.304,.307,.311,.315,.323,.334,.345,.360,.381,.389,.366, 4167. D .339,.318,.312,.314,.328,.344,.355,.368,.388,.401,.384,.360, 4168. E .325,.313,.302,.300,.318,.339,.354,.369,.381,.380,.360,.337, 4169. F .299,.291,.285,.284,.296,.313,.326,.340,.350,.343,.328,.312, 4170. G .277,.271,.269,.269,.272,.281,.288,.296,.308,.311,.298,.289, 4171. H .257,.253,.252,.254,.253,.257,.262,.267,.272,.281,.272,.265, 4172. I .241,.241,.241,.246,.245,.248,.246,.248,.253,.260,.255,.250, 4173. J .234,.236,.242,.256,.260,.260,.246,.244,.240,.241,.237,.237, 4174. K .243,.246,.257,.273,.279,.276,.261,.258,.251,.246,.238,.238, 4175. L .270,.269,.288,.299,.308,.299,.283,.276,.269,.263,.252,.257, 4176. M .327,.339,.358,.349,.351,.337,.313,.292,.288,.280,.284,.302, 4177. N .407,.419,.432,.407,.390,.356,.324,.298,.300,.304,.327,.368, 4178. O .421,.455,.459,.439,.393,.333,.306,.287,.289,.311,.345,.377, 4179. P .408,.452,.465,.443,.399,.323,.296,.276,.279,.309,.338,.362, 4180. Q .387,.437,.459,.444,.404,.334,.301,.276,.277,.302,.320,.345, 4181. R .379,.433,.455,.447,.408,.343,.313,.282,.279,.298,.315,.336/ 4182. DATA O3AVEI/ 4183. A .313,.291,.291,.293,.299,.296,.302,.310,.319,.333,.379,.354, 4184. B .316,.292,.295,.300,.307,.306,.315,.322,.333,.354,.384,.354, 4185. C .322,.302,.301,.307,.309,.319,.331,.340,.357,.379,.385,.356, 4186. D .328,.310,.301,.306,.316,.332,.347,.359,.380,.397,.379,.348, 4187. E .315,.304,.293,.296,.308,.328,.345,.360,.374,.376,.356,.329, 4188. F .292,.285,.277,.278,.288,.304,.318,.330,.340,.340,.324,.306, 4189. G .271,.266,.262,.263,.266,.277,.283,.291,.301,.307,.293,.284, 4190. H .253,.249,.249,.252,.250,.256,.261,.267,.271,.278,.267,.263, 4191. I .240,.238,.240,.247,.244,.248,.247,.250,.254,.258,.251,.249, 4192. J .233,.236,.243,.254,.259,.258,.248,.246,.241,.243,.238,.238, 4193. K .242,.246,.256,.268,.273,.271,.260,.255,.250,.244,.240,.239, 4194. L .258,.266,.278,.290,.295,.288,.277,.269,.265,.257,.253,.256, 4195. M .294,.308,.325,.326,.322,.308,.297,.284,.278,.271,.277,.287, 4196. N .338,.368,.383,.371,.357,.329,.316,.294,.287,.288,.303,.324, 4197. O .375,.420,.429,.411,.382,.328,.312,.293,.287,.299,.322,.354, 4198. P .388,.440,.454,.437,.396,.328,.307,.285,.282,.305,.330,.359, 4199. Q .386,.439,.457,.444,.404,.338,.309,.283,.280,.304,.321,.349, 4200. R .379,.435,.456,.448,.408,.345,.316,.286,.281,.300,.317,.337/ 4201. DATA O3AVEJ/ 4202. A .313,.290,.290,.291,.298,.294,.301,.309,.318,.331,.378,.353, 4203. B .313,.291,.291,.296,.304,.302,.311,.318,.330,.348,.382,.350, 4204. C .315,.297,.294,.300,.306,.310,.325,.334,.348,.364,.378,.346, 4205. D .316,.301,.292,.297,.305,.317,.334,.346,.360,.371,.366,.335, 4206. E .304,.293,.283,.286,.295,.313,.330,.344,.356,.359,.346,.316, 4207. F .284,.276,.268,.271,.279,.297,.309,.320,.325,.330,.317,.296, 4208. G .265,.258,.254,.257,.261,.273,.280,.288,.289,.296,.287,.274, 4209. H .250,.245,.244,.249,.247,.255,.260,.265,.268,.273,.263,.257, 4210. I .237,.235,.238,.246,.246,.249,.247,.249,.251,.257,.249,.247, 4211. J .234,.236,.245,.256,.259,.255,.248,.249,.244,.245,.242,.238, 4212. K .244,.249,.259,.271,.273,.270,.258,.256,.253,.247,.243,.242, 4213. L .261,.273,.283,.291,.292,.284,.271,.269,.263,.257,.254,.257, 4214. M .289,.305,.319,.321,.315,.301,.287,.281,.273,.268,.272,.282, 4215. N .321,.347,.364,.358,.344,.319,.305,.293,.282,.281,.291,.313, 4216. O .357,.400,.409,.397,.373,.332,.314,.295,.286,.293,.309,.333, 4217. P .377,.429,.442,.429,.396,.338,.317,.294,.287,.302,.321,.351, 4218. Q .385,.439,.458,.443,.407,.345,.318,.292,.284,.304,.322,.349, 4219. R .380,.437,.458,.449,.408,.348,.319,.289,.283,.301,.319,.340/ 4220. DATA O3AVEK/ 4221. A .311,.289,.289,.290,.298,.293,.300,.308,.317,.329,.377,.352, 4222. B .308,.290,.288,.291,.301,.296,.307,.315,.326,.340,.377,.344, 4223. C .305,.291,.287,.293,.297,.302,.315,.325,.335,.346,.369,.333, 4224. D .299,.289,.281,.287,.293,.302,.317,.327,.335,.344,.353,.318, 4225. E .287,.279,.272,.277,.281,.295,.309,.320,.325,.332,.331,.301, 4226. F .272,.264,.259,.262,.268,.281,.292,.300,.300,.309,.305,.282, 4227. G .257,.249,.246,.250,.254,.264,.271,.278,.279,.285,.278,.263, 4228. H .246,.239,.239,.245,.245,.252,.255,.261,.262,.267,.259,.250, 4229. I .234,.231,.239,.245,.245,.248,.245,.249,.248,.254,.246,.243, 4230. J .235,.237,.247,.258,.260,.257,.250,.250,.245,.246,.241,.240, 4231. K .248,.254,.264,.276,.276,.272,.262,.258,.255,.250,.248,.246, 4232. L .267,.278,.289,.300,.296,.286,.272,.270,.263,.258,.258,.262, 4233. M .292,.310,.325,.329,.319,.302,.288,.280,.273,.268,.274,.281, 4234. N .323,.346,.365,.365,.347,.320,.305,.291,.282,.281,.292,.305, 4235. O .352,.390,.405,.398,.378,.338,.316,.300,.290,.294,.309,.330, 4236. P .376,.424,.440,.431,.404,.350,.323,.303,.293,.303,.321,.349, 4237. Q .386,.442,.462,.448,.411,.354,.324,.298,.289,.306,.325,.349, 4238. R .381,.441,.459,.452,.410,.352,.322,.293,.286,.301,.320,.342/ 4239. DATA O3AVEL/ 4240. A .309,.290,.288,.288,.295,.292,.299,.307,.315,.327,.375,.350, 4241. B .306,.289,.287,.288,.298,.293,.304,.311,.320,.333,.372,.340, 4242. C .298,.286,.282,.288,.290,.294,.308,.316,.322,.332,.362,.325, 4243. D .289,.280,.274,.281,.282,.290,.304,.312,.317,.325,.342,.309, 4244. E .276,.269,.264,.268,.271,.281,.293,.300,.304,.313,.318,.290, 4245. F .262,.256,.253,.255,.258,.267,.278,.283,.283,.293,.294,.272, 4246. G .250,.245,.241,.245,.246,.255,.261,.267,.265,.282,.272,.256, 4247. H .240,.235,.236,.243,.240,.245,.249,.254,.253,.260,.254,.247, 4248. I .232,.229,.239,.245,.244,.247,.241,.245,.241,.246,.243,.241, 4249. J .235,.236,.247,.258,.258,.254,.246,.246,.239,.240,.238,.240, 4250. K .248,.253,.263,.273,.271,.267,.256,.253,.245,.243,.243,.244, 4251. L .265,.274,.287,.293,.290,.281,.267,.262,.256,.251,.253,.258, 4252. M .293,.307,.324,.323,.315,.298,.284,.275,.268,.263,.271,.278, 4253. N .326,.348,.370,.363,.347,.320,.304,.290,.281,.278,.291,.306, 4254. O .357,.391,.412,.404,.380,.347,.322,.303,.296,.296,.313,.334, 4255. P .381,.431,.447,.439,.412,.363,.331,.311,.301,.308,.331,.353, 4256. Q .389,.449,.470,.456,.417,.363,.329,.306,.296,.308,.331,.354, 4257. R .382,.441,.462,.454,.413,.354,.325,.296,.289,.301,.319,.343/ 4258. DATA O3AVEM/ 4259. A .309,.290,.288,.289,.293,.292,.299,.306,.313,.325,.374,.350, 4260. B .306,.289,.286,.285,.296,.291,.300,.308,.316,.326,.369,.339, 4261. C .297,.284,.281,.285,.288,.290,.302,.308,.315,.324,.355,.323, 4262. D .287,.278,.272,.275,.277,.284,.295,.300,.306,.316,.333,.304, 4263. E .273,.266,.261,.263,.267,.274,.284,.288,.292,.302,.311,.286, 4264. F .260,.253,.250,.252,.253,.261,.268,.273,.275,.284,.288,.269, 4265. G .247,.244,.241,.245,.243,.250,.254,.260,.260,.270,.268,.254, 4266. H .238,.234,.235,.242,.239,.243,.244,.250,.249,.255,.253,.245, 4267. I .231,.231,.238,.244,.242,.246,.238,.242,.239,.243,.242,.239, 4268. J .236,.238,.247,.257,.254,.253,.245,.244,.237,.235,.235,.236, 4269. K .250,.254,.263,.270,.266,.264,.254,.250,.244,.239,.237,.243, 4270. L .270,.279,.289,.290,.285,.279,.267,.261,.256,.250,.251,.258, 4271. M .301,.317,.329,.322,.314,.298,.285,.277,.270,.263,.270,.282, 4272. N .342,.367,.380,.369,.351,.326,.309,.294,.286,.284,.295,.314, 4273. O .380,.412,.424,.411,.388,.357,.331,.311,.303,.302,.325,.347, 4274. P .398,.448,.457,.449,.419,.373,.343,.318,.309,.314,.341,.366, 4275. Q .396,.456,.480,.466,.424,.370,.338,.311,.303,.311,.336,.363, 4276. R .384,.442,.464,.456,.414,.358,.327,.297,.290,.302,.322,.344/ 4277. DATA O3AVEN/ 4278. A .311,.291,.287,.288,.293,.292,.297,.305,.312,.325,.373,.350, 4279. B .307,.290,.286,.285,.293,.292,.300,.305,.315,.326,.366,.341, 4280. C .300,.287,.283,.282,.288,.292,.300,.306,.313,.324,.351,.323, 4281. D .290,.281,.274,.276,.279,.285,.293,.298,.303,.315,.330,.308, 4282. E .276,.272,.265,.264,.267,.274,.281,.287,.288,.302,.309,.289, 4283. F .263,.259,.254,.253,.257,.262,.267,.272,.274,.285,.287,.273, 4284. G .252,.247,.244,.248,.247,.252,.254,.260,.262,.270,.268,.259, 4285. H .243,.238,.239,.244,.241,.245,.245,.251,.251,.257,.253,.249, 4286. I .236,.233,.238,.244,.244,.246,.238,.243,.242,.245,.243,.242, 4287. J .237,.241,.247,.256,.255,.254,.245,.245,.242,.234,.234,.236, 4288. K .252,.259,.266,.271,.269,.269,.257,.256,.251,.242,.240,.245, 4289. L .277,.286,.296,.298,.292,.290,.276,.275,.267,.259,.259,.267, 4290. M .323,.342,.352,.339,.333,.319,.303,.298,.288,.280,.285,.296, 4291. N .374,.403,.413,.392,.376,.351,.332,.319,.306,.303,.317,.340, 4292. O .408,.448,.448,.433,.410,.375,.351,.330,.317,.318,.343,.368, 4293. P .418,.467,.473,.464,.426,.383,.347,.328,.316,.319,.347,.376, 4294. Q .402,.459,.482,.474,.426,.374,.343,.313,.306,.313,.338,.368, 4295. R .384,.440,.463,.458,.415,.360,.328,.299,.291,.301,.319,.344/ 4296. DATA O3AVEO/ 4297. A .313,.291,.288,.288,.292,.292,.298,.305,.312,.324,.364,.351, 4298. B .311,.294,.289,.286,.294,.293,.302,.306,.316,.326,.358,.345, 4299. C .308,.296,.291,.286,.294,.297,.303,.310,.316,.330,.354,.331, 4300. D .301,.292,.284,.282,.286,.295,.301,.307,.310,.326,.334,.318, 4301. E .290,.283,.274,.273,.276,.286,.291,.297,.299,.314,.314,.302, 4302. F .280,.272,.266,.263,.264,.272,.277,.283,.286,.297,.295,.286, 4303. G .267,.261,.256,.254,.255,.260,.263,.268,.272,.280,.276,.271, 4304. H .254,.250,.249,.249,.247,.251,.251,.256,.259,.264,.261,.258, 4305. I .242,.242,.243,.245,.244,.248,.242,.247,.248,.252,.248,.248, 4306. J .237,.242,.249,.256,.255,.255,.245,.244,.243,.237,.236,.236, 4307. K .253,.256,.267,.271,.270,.270,.259,.258,.252,.245,.242,.248, 4308. L .279,.283,.296,.296,.294,.292,.280,.279,.269,.260,.260,.268, 4309. M .327,.339,.357,.345,.338,.328,.319,.309,.293,.284,.285,.302, 4310. N .386,.409,.421,.405,.388,.363,.346,.332,.314,.311,.319,.348, 4311. O .419,.450,.459,.445,.418,.384,.361,.338,.322,.320,.340,.373, 4312. P .419,.461,.473,.468,.423,.358,.358,.331,.316,.319,.343,.376, 4313. Q .401,.453,.477,.469,.423,.375,.345,.314,.307,.312,.333,.361, 4314. R .382,.437,.461,.455,.415,.361,.329,.299,.291,.301,.316,.341/ 4315. DATA O3AVEP/ 4316. A .314,.293,.289,.290,.292,.294,.299,.305,.312,.323,.363,.352, 4317. B .315,.298,.293,.290,.294,.299,.303,.307,.316,.324,.365,.350, 4318. C .315,.303,.296,.291,.300,.306,.311,.316,.323,.336,.360,.341, 4319. D .308,.301,.293,.291,.297,.308,.312,.318,.324,.337,.345,.329, 4320. E .299,.292,.284,.283,.285,.299,.306,.311,.317,.326,.327,.314, 4321. F .285,.280,.272,.272,.274,.284,.293,.296,.301,.308,.306,.297, 4322. G .272,.266,.262,.261,.262,.269,.275,.280,.283,.289,.284,.280, 4323. H .256,.253,.251,.251,.251,.255,.256,.264,.266,.271,.267,.263, 4324. I .241,.242,.244,.245,.245,.248,.245,.251,.251,.255,.252,.251, 4325. J .236,.239,.247,.253,.253,.251,.242,.244,.239,.237,.235,.236, 4326. K .248,.250,.262,.267,.264,.262,.254,.250,.244,.240,.235,.239, 4327. L .268,.270,.286,.287,.284,.278,.267,.264,.256,.250,.245,.256, 4328. M .301,.308,.329,.322,.317,.300,.297,.281,.272,.264,.263,.279, 4329. N .351,.362,.380,.372,.360,.337,.320,.305,.295,.285,.287,.316, 4330. O .383,.406,.427,.415,.391,.365,.345,.324,.310,.304,.310,.342, 4331. P .393,.428,.450,.441,.404,.373,.353,.324,.310,.310,.321,.356, 4332. Q .387,.435,.461,.456,.412,.370,.341,.313,.303,.306,.321,.353, 4333. R .381,.432,.457,.452,.413,.361,.328,.299,.291,.298,.314,.338/ 4334. DATA O3AVEQ/ 4335. A .315,.293,.289,.291,.293,.295,.298,.305,.312,.323,.362,.354, 4336. B .316,.301,.295,.291,.294,.300,.303,.307,.316,.322,.361,.350, 4337. C .318,.305,.297,.292,.298,.306,.311,.314,.324,.334,.354,.340, 4338. D .309,.301,.292,.289,.295,.305,.312,.317,.326,.335,.343,.326, 4339. E .295,.288,.279,.279,.284,.297,.305,.305,.316,.321,.324,.310, 4340. F .279,.272,.266,.269,.272,.281,.289,.291,.299,.303,.305,.293, 4341. G .263,.259,.254,.257,.259,.266,.273,.276,.281,.285,.284,.277, 4342. H .247,.246,.244,.248,.247,.252,.253,.261,.265,.269,.267,.259, 4343. I .235,.236,.239,.244,.243,.246,.243,.247,.251,.253,.249,.246, 4344. J .231,.234,.243,.250,.251,.247,.240,.238,.233,.234,.232,.233, 4345. K .242,.244,.257,.262,.260,.255,.247,.243,.235,.235,.228,.233, 4346. L .257,.263,.278,.280,.275,.269,.258,.252,.242,.239,.235,.243, 4347. M .280,.288,.308,.307,.299,.287,.274,.267,.255,.250,.246,.259, 4348. N .309,.319,.348,.340,.332,.309,.293,.286,.273,.264,.261,.282, 4349. O .339,.357,.388,.376,.360,.334,.320,.305,.289,.282,.279,.306, 4350. P .365,.393,.424,.411,.386,.355,.340,.316,.300,.303,.297,.329, 4351. Q .375,.415,.445,.439,.404,.365,.336,.310,.298,.299,.306,.338, 4352. R .379,.428,.453,.447,.412,.360,.326,.298,.291,.296,.310,.335/ 4353. DATA O3AVER/ 4354. A .316,.295,.291,.292,.292,.296,.299,.305,.313,.323,.361,.355, 4355. B .317,.301,.296,.292,.292,.300,.302,.305,.314,.319,.358,.348, 4356. C .316,.303,.295,.289,.291,.301,.306,.307,.317,.324,.348,.336, 4357. D .303,.294,.286,.283,.285,.296,.304,.304,.313,.322,.333,.318, 4358. E .283,.277,.272,.272,.273,.284,.290,.296,.302,.309,.314,.299, 4359. F .265,.262,.259,.259,.259,.268,.274,.282,.286,.293,.293,.279, 4360. G .252,.249,.248,.249,.247,.253,.258,.265,.272,.277,.273,.265, 4361. H .241,.238,.240,.242,.241,.244,.246,.252,.257,.260,.256,.249, 4362. I .231,.229,.238,.241,.241,.242,.237,.242,.244,.247,.242,.239, 4363. J .231,.233,.242,.249,.251,.246,.237,.235,.230,.230,.229,.230, 4364. K .241,.250,.257,.265,.262,.257,.245,.243,.234,.230,.229,.231, 4365. L .260,.273,.281,.285,.280,.272,.257,.256,.245,.238,.237,.245, 4366. M .285,.302,.312,.314,.305,.294,.278,.277,.262,.252,.251,.262, 4367. N .310,.331,.347,.346,.336,.320,.303,.298,.281,.267,.267,.283, 4368. O .331,.354,.383,.378,.364,.342,.324,.315,.293,.278,.279,.297, 4369. P .350,.379,.414,.398,.381,.343,.335,.317,.299,.287,.285,.311, 4370. Q .367,.404,.436,.428,.399,.361,.332,.307,.295,.293,.298,.327, 4371. R .376,.424,.450,.442,.409,.358,.326,.296,.290,.294,.306,.332/ 4372. C 4373. DIMENSION AO3AVE(18,12),SO3JF(11,19),SO3SO(11,19) 4374. DATA AO3AVE/ .3148,.3160,.3171,.3159,.3027,.2824,.2645, 4375. A.2493,.2376,.2344,.2455,.2667,.3038,.3467,.3753,.3842,.3817,.3780, 4376. B.2926,.2959,.3008,.3035,.2943,.2763,.2600,.2463,.2366,.2366,.2500, 4377. C.2735,.3166,.3661,.4076,.4270,.4310,.4309,.2904,.2937,.2974,.2959, 4378. D.2869,.2704,.2561,.2454,.2403,.2443,.2590,.2844,.3293,.3803,.4210, 4379. E.4439,.4534,.4539,.2918,.2943,.2965,.2940,.2834,.2687,.2561,.2476, 4380. F.2450,.2538,.2676,.2888,.3259,.3692,.4077,.4325,.4454,.4476,.2951, 4381. G.2979,.2994,.3001,.2904,.2731,.2575,.2467,.2441,.2548,.2675,.2873, 4382. H.3181,.3517,.3828,.4002,.4080,.4096,.2960,.3012,.3084,.3132,.3044, 4383. I.2852,.2660,.2515,.2465,.2521,.2641,.2802,.3023,.3257,.3417,.3457, 4384. J.3521,.3517,.3008,.3070,.3153,.3211,.3127,.2934,.2714,.2545,.2437, 4385. K.2440,.2528,.2665,.2875,.3064,.3191,.3222,.3210,.3201,.3074,.3126, 4386. L.3221,.3276,.3211,.3015,.2783,.2603,.2478,.2431,.2499,.2624,.2784, 4387. M.2928,.3024,.3017,.2954,.2914,.3156,.3224,.3326,.3391,.3300,.3071, 4388. N.2827,.2632,.2489,.2399,.2455,.2566,.2720,.2854,.2939,.2931,.2889, 4389. O.2854,.3282,.3354,.3456,.3504,.3368,.3124,.2899,.2692,.2532,.2389, 4390. P.2415,.2521,.2672,.2844,.2967,.3003,.2986,.2966,.3723,.3713,.3661, 4391. Q.3538,.3332,.3072,.2826,.2626,.2481,.2359,.2373,.2489,.2700,.2936, 4392. R.3113,.3172,.3154,.3130,.3554,.3533,.3467,.3353,.3146,.2925,.2723, 4393. S.2562,.2450,.2350,.2387,.2554,.2828,.3140,.3331,.3406,.3408,.3351/ 4394. C 4395. DATA SO3JF/ 4396. A 13.0,12.3,11.7,10.5,8.90,6.20,4.50,3.30,2.20,1.80,1.00, 4397. B 13.6,12.9,11.9,10.3,8.30,6.10,4.45,3.40,2.50,1.85,1.00, 4398. C 14.8,13.9,12.8,10.3,8.00,6.00,4.55,3.60,2.70,1.90,1.00, 4399. D 16.6,15.1,14.0,11.0,7.95,6.00,4.65,3.70,2.95,1.95,1.00, 4400. E 18.1,16.0,14.6,12.0,8.00,6.00,4.80,3.75,3.00,1.98,1.00, 4401. F 18.3,16.3,14.8,12.6,8.20,6.15,4.80,3.80,3.05,2.00,1.00, 4402. G 17.3,16.1,14.7,12.7,9.10,6.10,4.70,3.75,3.00,2.00,1.00, 4403. H 16.3,15.5,14.5,12.6,9.00,6.00,4.55,3.65,2.95,1.98,1.00, 4404. I 15.7,14.9,14.1,12.4,8.70,5.90,4.40,3.45,2.80,1.96,1.00, 4405. J 15.3,14.1,13.5,12.2,8.30,5.85,4.25,3.40,2.75,1.95,1.00, 4406. K 15.6,14.9,14.0,12.4,9.00,6.10,4.55,3.50,2.85,1.96,1.00, 4407. L 17.4,16.6,16.0,14.0,10.0,7.30,5.10,3.90,3.00,1.97,1.00, 4408. M 17.6,18.3,17.8,15.8,12.3,9.00,6.05,4.40,3.20,1.97,1.00, 4409. N 16.0,16.9,17.8,16.8,15.2,12.0,7.90,5.10,3.65,1.97,1.00, 4410. O 12.3,13.8,15.7,16.2,16.2,14.8,10.0,6.00,4.00,1.96,1.00, 4411. P 12.0,11.9,12.0,13.8,14.3,14.3,12.0,6.80,4.30,1.95,1.00, 4412. Q 11.9,11.8,11.7,11.6,11.8,12.0,10.3,7.20,4.50,1.90,1.00, 4413. R 11.6,11.5,11.4,11.2,11.0,10.4,9.00,7.20,4.15,1.85,1.00, 4414. S 11.2,10.9,10.7,10.5,10.0,9.75,8.60,7.00,3.80,1.80,1.00/ 4415. DATA SO3SO/ 4416. A 10.5,10.5,10.5,10.6,10.5,10.3,8.20,4.80,3.10,1.90,1.00, 4417. B 11.5,11.5,11.6,12.1,12.1,10.8,8.05,4.95,3.40,1.92,1.00, 4418. C 12.7,13.8,14.0,14.1,12.9,10.9,7.95,5.10,3.70,1.96,1.00, 4419. D 15.4,15.9,16.0,15.4,13.2,10.7,7.40,5.15,3.85,1.98,1.00, 4420. E 17.9,18.0,17.4,16.1,13.0,10.0,6.70,4.90,3.80,1.99,1.00, 4421. F 18.3,18.6,17.8,16.1,12.1,9.10,5.95,4.80,3.70,2.00,1.00, 4422. G 18.6,18.5,17.8,15.9,11.1,8.00,5.55,4.40,3.45,2.00,1.00, 4423. H 18.2,18.1,17.2,15.1,10.3,7.40,5.10,4.00,3.10,1.99,1.00, 4424. I 17.5,16.8,16.2,14.0,9.90,7.00,4.90,3.85,2.95,1.98,1.00, 4425. J 16.5,15.8,15.0,12.9,9.40,6.65,4.80,3.70,2.90,1.96,1.00, 4426. K 16.3,15.8,15.0,12.9,9.20,6.80,5.00,3.85,2.95,1.96,1.00, 4427. L 16.4,16.2,15.8,14.0,9.80,7.10,5.10,3.95,3.00,1.96,1.00, 4428. M 16.6,16.5,16.2,14.8,10.8,7.75,5.50,4.05,3.05,1.97,1.00, 4429. N 16.5,16.6,16.5,16.0,12.1,9.00,6.00,4.40,3.10,1.97,1.00, 4430. O 15.8,16.2,16.4,16.1,14.2,10.9,6.60,4.50,3.20,1.97,1.00, 4431. P 12.2,14.2,15.5,15.3,14.7,12.4,7.40,4.70,3.10,1.96,1.00, 4432. Q 11.6,11.9,12.1,14.0,13.9,12.3,8.00,4.40,2.95,1.90,1.00, 4433. R 11.2,11.2,11.4,11.6,11.8,10.9,8.00,3.95,2.60,1.87,1.00, 4434. S 11.0,10.8,10.5,10.3,10.1,9.70,7.00,3.65,2.20,1.80,1.00/ 4435. C 4436. DIMENSION XJDMO(14),HKMSPR(14),HKMAUT(14) 4437. DIMENSION CNCAUT(14),CNCSPR(14),DEGLAT(14) 4438. DATA DEGLAT/-85.0,-71.0,-59.0,-47.0,-35.0,-22.0,-9.0, 4439. + 9.0,22.0,35.0,47.0,59.0,71.0,85.0/ 4440. DATA XJDMO/-15.0,16.0,45.0,75.0,105.0,136.0,166.0,197.0,228.0 4441. + ,258.0,289.0,319.0,350.0,381.0/ 4442. DATA HKMSPR/18.5,18.5,19.0,23.5,24.0,24.5,26.5, 4443. + 26.5,25.0,22.5,21.0,20.0,18.5,16.5/ 4444. DATA HKMAUT/16.5,18.5,20.0,21.0,22.5,25.0,26.5, 4445. + 26.5,24.5,24.0,23.5,19.0,18.5,18.5/ 4446. DATA CNCSPR/0.0181,0.0212,0.0187,0.0167,0.0162,0.0183,0.0175, 4447. + 0.0187,0.0200,0.0196,0.0225,0.0291,0.0287,0.0300/ 4448. DATA CNCAUT/0.0300,0.0287,0.0291,0.0225,0.0196,0.0200,0.0187, 4449. + 0.0175,0.0183,0.0162,0.0167,0.0187,0.0212,0.0181/ 4450. C 4451. DIMENSION PLBSO3(11),SOJDAY(6),PMLAT(6) 4452. DATA PLBSO3/10.0,7.0,5.0,3.0,2.0,1.5,1.0,0.7,0.5,0.3,0.1/ 4453. DATA SOJDAY/-91.,31.,92.,213.,274.,396./ 4454. DATA PMLAT/1.,1.,-1.,-1.,1.,1./ 4455. DIMENSION AO3JIM(99),O3LB(40),PLE0(40) 4456. DIMENSION CONCS(99),CONCA(99),BHKMS(99),BHKMA(99) 4457. DIMENSION WTJLAT(99),WTJLON(99),ILATIJ(99),ILONIJ(99) 4458. DIMENSION WTLSEP(99),WTLJAN(99),LSEPJ(99),LJANJ(99) 4459. DATA ACMMGG/2.37251E-4/,ACMPKM/7.1509E-4/,H10MB/31.05467/ 4460. DATA A,B,C,D/0.331,23.0,4.553,5.23/ 4461. C 4462. C----------------------------------------------------------------------- 4463. C----SET O3 VERTICAL PROFILE PARAMETERS FOR LATITUDE GCM GRID POINTS 4464. C----------------------------------------------------------------------- 4465. C SKIPI =.FALSE. 4466. C IF(ABS(FLONO3).LT.1.E-04) SKIPI =.TRUE. 4467. DO 100 L=1,LMR 4468. 100 PLE0(L)=PLE(L) 4469. DO 103 J=1,JM 4470. DGLATJ=DGLAT(J) 4471. ILATI=(DGLATJ+95.001)/10. 4472. IF(ILATI.LT. 1) ILATI= 1 4473. IF(ILATI.GT.17