LMDZ
suswn.F90
Go to the documentation of this file.
1 SUBROUTINE suswn (KTSW, KSW)
2 
3 !**** *SUSW* - INITIALIZE COMMON YOESW
4 
5 ! PURPOSE.
6 ! --------
7 ! INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS
8 ! NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES
9 
10 !** INTERFACE.
11 ! ----------
12 ! *CALL* *SUSW
13 
14 ! EXPLICIT ARGUMENTS :
15 ! --------------------
16 ! NONE
17 
18 ! IMPLICIT ARGUMENTS :
19 ! --------------------
20 ! COMMON YOESW
21 
22 ! METHOD.
23 ! -------
24 ! SEE DOCUMENTATION
25 
26 ! EXTERNALS.
27 ! ----------
28 
29 ! REFERENCE.
30 ! ----------
31 ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
32 
33 ! AUTHOR.
34 ! -------
35 ! JEAN-JACQUES MORCRETTE *ECMWF*
36 
37 ! MODIFICATIONS.
38 ! --------------
39 ! ORIGINAL : 88-12-15
40 ! 97-04-16 JJ Morcrette 2 and 4 interval spectral resolution
41 ! 00-10-24 JJ Morcrette sea-ice albedo revisited
42 ! 00-12-14 JJ Morcrette
43 ! and Ph.Dubuisson B.Bonnel 6 spectral interval resolution
44 ! 01-04-17 Ph.Dubuisson, B.Bonnel, JJ.Morcrette 6 sp.int.resolu'n
45 ! 01-06-28 B.Bonnel, JJ.Morcrette, Ph.Dubuisson Rayleigh (2/4/6)
46 ! 01-11-05 Ph.Dubuisson, JJMorcrette (new 2 intervals for TL/AD)
47 ! 03-04-24 JJ Morcrette SRTM with 14 spectral intervals
48 ! M.Hamrud 01-Oct-2003 CY28 Cleaning
49 
50 ! ------------------------------------------------------------------
51 
52 USE parkind1 ,ONLY : jpim ,jprb
53 USE yomhook ,ONLY : lhook, dr_hook
54 
55 USE yomcst , ONLY : rg
56 USE yoerad , ONLY : ntsw ,lonewsw ,lrayl
57 USE yoesw , ONLY : apad ,bpad ,rray ,rsun ,&
58  & rpdh1 ,rpdu1 ,rpnh ,rpnu ,rswce ,&
59  & rswcp ,rtdh2o ,rtdumg ,rth2o ,rtumg ,&
60  & d ,rroma ,rromb ,rrasy ,rhsra ,&
61  & rhsrb ,rhsrc ,rhsrd ,rhsre ,rhsrf ,&
62  & rhsrta ,rhsrtb ,rexpo3 ,&
63  & radjust ,nexpo3
64 
65 IMPLICIT NONE
66 
67 INTEGER(KIND=JPIM),INTENT(IN) :: KTSW
68 INTEGER(KIND=JPIM),INTENT(IN) :: KSW
69 INTEGER(KIND=JPIM) :: IMPSRTM(14)
70 
71 ! ----------------------------------------------------------------
72 REAL(KIND=JPRB) :: ZAPAD2(2,3,7) , ZBPAD2(2,3,7) , ZD2(2,3)&
73  & , ZRAY2(2,6) , ZSUN2(2) , ZSWCE2(2) , ZSWCP2(2)&
74  & , ZROMA2(2) , ZROMB2(2) , ZRASY2(2)&
75  & , ZRA2(2) , ZRB2(2) , ZRC2(2)&
76  & , ZRD2(2) , ZRE2(2) , ZRF2(2)
77 
78 REAL(KIND=JPRB) :: ZAPAD4(4,3,7) , ZBPAD4(4,3,7) , ZD4(4,3)&
79  & , ZRAY4(4,6) , ZSUN4(4) , ZSWCE4(4) , ZSWCP4(4)&
80  & , ZROMA4(4) , ZROMB4(4) , ZRASY4(4)&
81  & , ZRA4(4) , ZRB4(4) , ZRC4(4)&
82  & , ZRD4(4) , ZRE4(4) , ZRF4(4)
83 
84 REAL(KIND=JPRB) :: ZAPAD6(6,3,7) , ZBPAD6(6,3,7) , ZD6(6,3)&
85  & , ZRAY6(6,6) , ZSUN6(6) , ZSWCE6(6) , ZSWCP6(6)&
86  & , ZROMA6(6) , ZROMB6(6) , ZRASY6(6)&
87  & , ZRA6(6) , ZRB6(6) , ZRC6(6)&
88  & , ZRD6(6) , ZRE6(6) , ZRF6(6)&
89  & , ZRBY6(6,6)
90 
91 
92 INTEGER(KIND=JPIM) :: JC3, JC6, JI, JJ, JW, I_K, ISW6
93 
94 REAL(KIND=JPRB) :: ZH2O, ZPDHNIS, ZPDH2O, ZPDUNIS,&
95  & ZPDUMG, ZPRHNIS, ZPRH2O, &
96  & ZPRUNIS, ZPRUMG, ZRTO1, ZRTO2, ZTDHNIS, &
97  & ZTDUNIS, ZTHNIS, &
98  & ZTUNIS, ZUMG, ZADJUS2, ZADJUS4, ZADJUS6
99 
100 REAL(KIND=JPRB) :: ZTDH2IS, ZTDU2IS, ZTH2IS, ZTU2IS, ZPDH2IS, ZPDU2IS, ZPRH2IS, ZPRU2IS
101 REAL(KIND=JPRB) :: ZTDH4IS, ZTDU4IS, ZTH4IS, ZTU4IS, ZPDH4IS, ZPDU4IS, ZPRH4IS, ZPRU4IS
102 REAL(KIND=JPRB) :: ZTDH6IS, ZTDU6IS, ZTH6IS, ZTU6IS, ZPDH6IS, ZPDU6IS, ZPRH6IS, ZPRU6IS
103 REAL(KIND=JPRB) :: ZHOOK_HANDLE
104 
105 #include "abor1.intfb.h"
106 
107 ! ----------------------------------------------------------------
108 
109 !* 1. CLEAR-SKY ABSORPTION COEFFICIENTS FOR N SPECTRAL INTERVALS
110 ! --------------------------------------------------------
111 
112 !* spectrally independent
113 
114 !-- New coefficients
115 
116 IF (lhook) CALL dr_hook('SUSWN',0,zhook_handle)
117 ztdhnis = 0.450_jprb
118 ztdunis = 0.375_jprb
119 zthnis = 273._jprb
120 ztunis = 273._jprb
121 zpdhnis = 0.90_jprb
122 zpdunis = 0.75_jprb
123 zprhnis = 101300._jprb
124 zprunis = 101300._jprb
125 isw6=6
126 
127 ! ----------------------------------------------------------------
128 
129 !* 1.1 COEFFICIENTS FOR TWO SPECTRAL INTERVALS
130 ! ---------------------------------------
131 
132 !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
133 ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
134 ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
135 ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
136 
137 zadjus2 = 1.0_jprb
138 ztdh2is = 0.450_jprb
139 ztdu2is = 0.375_jprb
140 zth2is = 273._jprb
141 ztu2is = 273._jprb
142 zpdh2is = 0.90_jprb
143 zpdu2is = 0.75_jprb
144 zprh2is = 101300._jprb
145 zpru2is = 101300._jprb
146 
147 !* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
148 
149 zsun2(1) = 0.459760_jprb
150 
151 !-- New coefficients
152 zd2(1,:)= (/ 0.800000000_jprb, 0.900000000_jprb, 0.000000000_jprb /)
153 
154 zapad2(1, 1, :) = (/&
155  & 0.1095121e+03_jprb,&
156  & 0.8906647e+03_jprb,&
157  & 0.1030601e+03_jprb,&
158  & 0.000000000e-00_jprb,&
159  & 0.000000000e-00_jprb,&
160  & 0.000000000e-00_jprb,&
161  & 0.000000000e-00_jprb/)
162 zapad2(1, 2, :) = (/&
163  & 0.5124541e+03_jprb,&
164  & 0.3478800e+04_jprb,&
165  & 0.1244272e+04_jprb,&
166  & 0.000000000e-00_jprb,&
167  & 0.000000000e-00_jprb,&
168  & 0.000000000e-00_jprb,&
169  & 0.000000000e-00_jprb/)
170 zapad2(1, 3, :) = (/&
171  & 0.4766857e-04_jprb,&
172  & 0.4019806e-02_jprb,&
173  & 0.1936247e+00_jprb,&
174  & 0.1886522e+02_jprb,&
175  & 0.6498215e+02_jprb,&
176  & 0.1799387e+02_jprb,&
177  & 0.000000000e+00_jprb/)
178 
179 zbpad2(1, 1, :) = (/&
180  & 0.1095121e+03_jprb,&
181  & 0.8919103e+03_jprb,&
182  & 0.1142218e+03_jprb,&
183  & 0.100000000e+01_jprb,&
184  & 0.000000000e-00_jprb,&
185  & 0.000000000e-00_jprb,&
186  & 0.000000000e-00_jprb/)
187 zbpad2(1, 2, :) = (/&
188  & 0.5124541e+03_jprb,&
189  & 0.3499516e+04_jprb,&
190  & 0.1278186e+04_jprb,&
191  & 0.100000000e+01_jprb,&
192  & 0.000000000e-00_jprb,&
193  & 0.000000000e-00_jprb,&
194  & 0.000000000e-00_jprb/)
195 zbpad2(1, 3, :) = (/&
196  & 0.4766857e-04_jprb,&
197  & 0.4153343e-02_jprb,&
198  & 0.1925810e+00_jprb,&
199  & 0.1938938e+02_jprb,&
200  & 0.6918191e+02_jprb,&
201  & 0.2278795e+02_jprb,&
202  & 0.100000000e+01_jprb/)
203 
204 zray2(1,:)= (/&
205  & .389107e-01_jprb, .101608e+01_jprb,-.354921e+01_jprb,&
206  & .659510e+01_jprb,-.592423e+01_jprb, .202348e+01_jprb/)
207 
208 !-- Original coefficients of ECMWF model (before 27 June 2000)
209 
210 !ZTDH2IS = 0.450_JPRB
211 !ZTDU2IS = 0.375_JPRB
212 !ZTH2IS = 296._JPRB
213 !ZTU2IS = 296._JPRB
214 !ZPDH2IS = 0.90_JPRB
215 !ZPDU2IS = 0.75_JPRB
216 !ZPRH2IS = 100000._JPRB
217 !ZPRU2IS = 100000._JPRB
218 
219 !ZD2(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
220 
221 !ZAPAD2(1, 1, :) = (/&
222 ! &0.184678379E+06_JPRB,&
223 ! &0.553080884E+05_JPRB,&
224 ! &0.248143712E+04_JPRB,&
225 ! &0.000000000E-00_JPRB,&
226 ! &0.000000000E-00_JPRB,&
227 ! &0.000000000E-00_JPRB,&
228 ! &0.000000000E-00_JPRB/)
229 !ZAPAD2(1, 2, :) = (/&
230 ! &0.715303869E+01_JPRB,&
231 ! &0.219386847E+03_JPRB,&
232 ! &0.830001089E+03_JPRB,&
233 ! &0.000000000E-00_JPRB,&
234 ! &0.000000000E-00_JPRB,&
235 ! &0.000000000E-00_JPRB,&
236 ! &0.000000000E-00_JPRB/)
237 !ZAPAD2(1, 3, :) = (/&
238 ! &0.925887084E-04_JPRB,&
239 ! &0.129353723E-01_JPRB,&
240 ! &0.800821928E+00_JPRB,&
241 ! &0.242715973E+02_JPRB,&
242 ! &0.878331486E+02_JPRB,&
243 ! &0.191559725E+02_JPRB,&
244 ! &0.000000000E+00_JPRB/)
245 
246 !ZBPAD2(1, 1, :) = (/&
247 ! &0.184678379E+06_JPRB,&
248 ! &0.555188347E+05_JPRB,&
249 ! &0.253257443E+04_JPRB,&
250 ! &0.100000000E+01_JPRB,&
251 ! &0.000000000E-00_JPRB,&
252 ! &0.000000000E-00_JPRB,&
253 ! &0.000000000E-00_JPRB/)
254 !ZBPAD2(1, 2, :) = (/&
255 ! &0.715303869E+01_JPRB,&
256 ! &0.219441875E+03_JPRB,&
257 ! &0.831119997E+03_JPRB,&
258 ! &0.100000000E+01_JPRB,&
259 ! &0.000000000E-00_JPRB,&
260 ! &0.000000000E-00_JPRB,&
261 ! &0.000000000E-00_JPRB/)
262 !ZBPAD2(1, 3, :) = (/&
263 ! &0.925887084E-04_JPRB,&
264 ! &0.131812683E-01_JPRB,&
265 ! &0.812706117E+00_JPRB,&
266 ! &0.249863591E+02_JPRB,&
267 ! &0.931071925E+02_JPRB,&
268 ! &0.252233437E+02_JPRB,&
269 ! &0.100000000E+01_JPRB/)
270 
271 zray2(1,:)= (/&
272  & .428937e-01_jprb, .890743e+00_jprb,-.288555e+01_jprb,&
273  & .522744e+01_jprb,-.469173e+01_jprb, .161645e+01_jprb/)
274 
275 zswce2(1) = 0.0_jprb
276 zswcp2(1) = 0.0_jprb
277 
278 ! ----------------------------------------------------------------
279 
280 !* Near-Infrared (0.69 - 4.0 Microns) is not sub-divided
281 
282 ! ----------------------------------------------------------------
283 
284 !* 0.68 - 4.00 Micron
285 
286 zsun2(2) = 0.540240_jprb
287 
288 !-- New coefficients
289 
290 zd2(2,:)= (/ 0.000000000_jprb, 0.900000000_jprb, 0.000000000_jprb/)
291 
292 zapad2(2, 1, :) = (/&
293  & 0.1345066e-03_jprb,&
294  & 0.9276888e+00_jprb,&
295  & 0.3449105e+03_jprb,&
296  & 0.7143025e+04_jprb,&
297  & 0.6531248e+04_jprb,&
298  & 0.2210815e+03_jprb,&
299  & 0.0000000e+00_jprb/)
300 zapad2(2, 2, :) = (/&
301  & 0.7776703e-05_jprb,&
302  & 0.6851592e-02_jprb,&
303  & 0.5533947e+00_jprb,&
304  & 0.2538585e+01_jprb,&
305  & 0.000000000e+00_jprb,&
306  & 0.000000000e+00_jprb,&
307  & 0.000000000e+00_jprb/)
308 zapad2(2, 3, :) = (/&
309  & 0.932012720e+02_jprb,&
310  & 0.437866983e+03_jprb,&
311  & 0.000000000e-00_jprb,&
312  & 0.000000000e-00_jprb,&
313  & 0.000000000e-00_jprb,&
314  & 0.000000000e-00_jprb,&
315  & 0.000000000e+00_jprb/)
316 
317 zbpad2(2, 1, :) = (/&
318  & 0.1345066e-03_jprb,&
319  & 0.9297728e+00_jprb,&
320  & 0.3523821e+03_jprb,&
321  & 0.7931590e+04_jprb,&
322  & 0.8979086e+04_jprb,&
323  & 0.5050864e+03_jprb,&
324  & 0.1000000e+01_jprb/)
325 zbpad2(2, 2, :) = (/&
326  & 0.7776703e-05_jprb,&
327  & 0.6976301e-02_jprb,&
328  & 0.6074521e+00_jprb,&
329  & 0.3517521e+01_jprb,&
330  & 0.100000000e+01_jprb,&
331  & 0.000000000e+00_jprb,&
332  & 0.000000000e+00_jprb/)
333 zbpad2(2, 3, :) = (/&
334  & 0.932012720e+02_jprb,&
335  & 0.439078599e+03_jprb,&
336  & 0.100000000e+01_jprb,&
337  & 0.000000000e+00_jprb,&
338  & 0.000000000e+00_jprb,&
339  & 0.000000000e+00_jprb,&
340  & 0.000000000e+00_jprb/)
341 
342 zray2(2,:)= (/&
343  & .104635e-01_jprb, .102535e+00_jprb,-.536771e+00_jprb,&
344  & .123640e+01_jprb,-.128930e+01_jprb, .502426e+00_jprb/)
345 
346 !-- Original coefficients of ECMWF model (before 27 June 2000)
347 
348 !ZD2(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB/)
349 
350 !ZAPAD2(2, 1, :) = (/&
351 ! &0.577160808E-06_JPRB,&
352 ! &0.681760991E-02_JPRB,&
353 ! &0.502416370E+01_JPRB,&
354 ! &0.209771270E+03_JPRB,&
355 ! &0.644493139E+03_JPRB,&
356 ! &0.120654163E+03_JPRB,&
357 ! &0.000000000E+00_JPRB/)
358 !ZAPAD2(2, 2, :) = (/&
359 ! &0.422874202E-09_JPRB,&
360 ! &0.257499477E-05_JPRB,&
361 ! &0.349570458E-02_JPRB,&
362 ! &0.917758876E+00_JPRB,&
363 ! &0.327669989E+02_JPRB,&
364 ! &0.874155115E+02_JPRB,&
365 ! &0.000000000E+00_JPRB/)
366 !ZAPAD2(2, 3, :) = (/&
367 ! &0.932012720E+02_JPRB,&
368 ! &0.437866983E+03_JPRB,&
369 ! &0.000000000E-00_JPRB,&
370 ! &0.000000000E-00_JPRB,&
371 ! &0.000000000E-00_JPRB,&
372 ! &0.000000000E-00_JPRB,&
373 ! &0.000000000E+00_JPRB/)
374 
375 !ZBPAD2(2, 1, :) = (/&
376 ! &0.577160808E-06_JPRB,&
377 ! &0.682595589E-02_JPRB,&
378 ! &0.508797283E+01_JPRB,&
379 ! &0.223901000E+03_JPRB,&
380 ! &0.770640908E+03_JPRB,&
381 ! &0.186812480E+03_JPRB,&
382 ! &0.100000000E+01_JPRB/)
383 !ZBPAD2(2, 2, :) = (/&
384 ! &0.422874202E-09_JPRB,&
385 ! &0.257567242E-05_JPRB,&
386 ! &0.349912408E-02_JPRB,&
387 ! &0.920770481E+00_JPRB,&
388 ! &0.330990960E+02_JPRB,&
389 ! &0.903249595E+02_JPRB,&
390 ! &0.100000000E+01_JPRB/)
391 !ZBPAD2(2, 3, :) = (/&
392 ! &0.932012720E+02_JPRB,&
393 ! &0.439078599E+03_JPRB,&
394 ! &0.100000000E+01_JPRB,&
395 ! &0.000000000E+00_JPRB,&
396 ! &0.000000000E+00_JPRB,&
397 ! &0.000000000E+00_JPRB,&
398 ! &0.000000000E+00_JPRB/)
399 
400 !ZRAY2(2,:)= (/&
401 ! &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,&
402 ! &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/)
403 
404 zswce2(2) = 0.0_jprb
405 zswcp2(2) = 0.0_jprb
406 
407 !=====================================================================
408 
409 !* 1.2 COEFFICIENTS FOR FOUR SPECTRAL INTERVALS
410 ! ----------------------------------------
411 
412 !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
413 ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
414 ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
415 ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
416 
417 zadjus4 = 1.0_jprb
418 
419 !* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
420 
421 zsun4(1) = 0.459760_jprb
422 
423 !-- New coefficients
424 
425 !ZD4(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
426 
427 !ZAPAD4(1, 1, :) = (/&
428 ! &0.1095121E+03_JPRB,&
429 ! &0.8906647E+03_JPRB,&
430 ! &0.1030601E+03_JPRB,&
431 ! &0.0000000E-00_JPRB,&
432 ! &0.0000000E-00_JPRB,&
433 ! &0.0000000E-00_JPRB,&
434 ! &0.0000000E-00_JPRB/)
435 !ZAPAD4(1, 2, :) = (/&
436 ! &0.5124541E+03_JPRB,&
437 ! &0.3478800E+04_JPRB,&
438 ! &0.1244272E+04_JPRB,&
439 ! &0.0000000E-00_JPRB,&
440 ! &0.0000000E-00_JPRB,&
441 ! &0.0000000E-00_JPRB,&
442 ! &0.0000000E-00_JPRB/)
443 !ZAPAD4(1, 3, :) = (/&
444 ! &0.4766857E-04_JPRB,&
445 ! &0.4019806E-02_JPRB,&
446 ! &0.1936247E+00_JPRB,&
447 ! &0.1886522E+02_JPRB,&
448 ! &0.6498215E+02_JPRB,&
449 ! &0.1799387E+02_JPRB,&
450 ! &0.0000000E+00_JPRB/)
451 
452 !ZBPAD4(1, 1, :) = (/&
453 ! &0.1095121E+03_JPRB,&
454 ! &0.8919103E+03_JPRB,&
455 ! &0.1142218E+03_JPRB,&
456 ! &0.1000000E+01_JPRB,&
457 ! &0.0000000E-00_JPRB,&
458 ! &0.0000000E-00_JPRB,&
459 ! &0.0000000E-00_JPRB/)
460 !ZBPAD4(1, 2, :) = (/&
461 ! &0.5124541E+03_JPRB,&
462 ! &0.3499516E+04_JPRB,&
463 ! &0.1278186E+04_JPRB,&
464 ! &0.1000000E+01_JPRB,&
465 ! &0.0000000E-00_JPRB,&
466 ! &0.0000000E-00_JPRB,&
467 ! &0.0000000E-00_JPRB/)
468 !ZBPAD4(1, 3, :) = (/&
469 ! &0.4766857E-04_JPRB,&
470 ! &0.4153343E-02_JPRB,&
471 ! &0.1925810E+00_JPRB,&
472 ! &0.1938938E+02_JPRB,&
473 ! &0.6918191E+02_JPRB,&
474 ! &0.2278795E+02_JPRB,&
475 ! &0.1000000E+01_JPRB/)
476 
477 !ZRAY4(1,:)= (/&
478 ! &.389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,&
479 ! &.659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/)
480 
481 !-- Original coefficients in ECMWF scheme
482 
483 ztdh4is = 0.450_jprb
484 ztdu4is = 0.375_jprb
485 zth4is = 296._jprb
486 ztu4is = 296._jprb
487 zpdh4is = 0.90_jprb
488 zpdu4is = 0.75_jprb
489 zprh4is = 100000._jprb
490 zpru4is = 100000._jprb
491 
492 zd4(1,:)= (/ 0.000000000_jprb, 0.000000000_jprb, 0.000000000_jprb /)
493 
494 zapad4(1, 1, :) = (/&
495  & 0.184678379e+06_jprb,&
496  & 0.553080884e+05_jprb,&
497  & 0.248143712e+04_jprb,&
498  & 0.000000000e-00_jprb,&
499  & 0.000000000e-00_jprb,&
500  & 0.000000000e-00_jprb,&
501  & 0.000000000e-00_jprb/)
502 zapad4(1, 2, :) = (/&
503  & 0.715303869e+01_jprb,&
504  & 0.219386847e+03_jprb,&
505  & 0.830001089e+03_jprb,&
506  & 0.000000000e-00_jprb,&
507  & 0.000000000e-00_jprb,&
508  & 0.000000000e-00_jprb,&
509  & 0.000000000e-00_jprb/)
510 zapad4(1, 3, :) = (/&
511  & 0.925887084e-04_jprb,&
512  & 0.129353723e-01_jprb,&
513  & 0.800821928e+00_jprb,&
514  & 0.242715973e+02_jprb,&
515  & 0.878331486e+02_jprb,&
516  & 0.191559725e+02_jprb,&
517  & 0.000000000e+00_jprb/)
518 
519 zbpad4(1, 1, :) = (/&
520  & 0.184678379e+06_jprb,&
521  & 0.555188347e+05_jprb,&
522  & 0.253257443e+04_jprb,&
523  & 0.100000000e+01_jprb,&
524  & 0.000000000e-00_jprb,&
525  & 0.000000000e-00_jprb,&
526  & 0.000000000e-00_jprb/)
527 zbpad4(1, 2, :) = (/&
528  & 0.715303869e+01_jprb,&
529  & 0.219441875e+03_jprb,&
530  & 0.831119997e+03_jprb,&
531  & 0.100000000e+01_jprb,&
532  & 0.000000000e-00_jprb,&
533  & 0.000000000e-00_jprb,&
534  & 0.000000000e-00_jprb/)
535 zbpad4(1, 3, :) = (/&
536  & 0.925887084e-04_jprb,&
537  & 0.131812683e-01_jprb,&
538  & 0.812706117e+00_jprb,&
539  & 0.249863591e+02_jprb,&
540  & 0.931071925e+02_jprb,&
541  & 0.252233437e+02_jprb,&
542  & 0.100000000e+01_jprb/)
543 
544 zray4(1,:)= (/&
545  & .428937e-01_jprb, .890743e+00_jprb,-.288555e+01_jprb,&
546  & .522744e+01_jprb,-.469173e+01_jprb, .161645e+01_jprb/)
547 
548 zswce4(1) = 0.0_jprb
549 zswcp4(1) = 0.0_jprb
550 
551 ! ----------------------------------------------------------------
552 
553 !* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
554 
555 ! ----------------------------------------------------------------
556 
557 !* 0.69 - 1.19 Micron
558 
559 zsun4(2) = 0.326158_jprb
560 
561 !-- New coefficients
562 
563 !ZD4(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /)
564 
565 !ZAPAD4(2, 1, :) = (/&
566 ! &0.1335726E+02_JPRB,&
567 ! &0.2939136E+04_JPRB,&
568 ! &0.4010585E+05_JPRB,&
569 ! &0.7195030E+05_JPRB,&
570 ! &0.1648338E+05_JPRB,&
571 ! &0.3373738E+03_JPRB,&
572 ! &0.0000000E+00_JPRB/)
573 !ZAPAD4(2, 2, :) = (/&
574 ! &0.2001271E-01_JPRB,&
575 ! &0.2480831E+01_JPRB,&
576 ! &0.3444162E+02_JPRB,&
577 ! &0.4788946E+02_JPRB,&
578 ! &0.0000000E+00_JPRB,&
579 ! &0.0000000E+00_JPRB,&
580 ! &0.0000000E+00_JPRB/)
581 !ZAPAD4(2, 3, :) = (/&
582 ! &0.000000000E+00_JPRB,&
583 ! &0.000000000E+00_JPRB,&
584 ! &0.000000000E+00_JPRB,&
585 ! &0.000000000E+00_JPRB,&
586 ! &0.000000000E+00_JPRB,&
587 ! &0.000000000E+00_JPRB,&
588 ! &0.000000000E+00_JPRB/)
589 
590 !ZBPAD4(2, 1, :) = (/&
591 ! &0.1335726E+02_JPRB,&
592 ! &0.2942327E+04_JPRB,&
593 ! &0.4077237E+05_JPRB,&
594 ! &0.7749017E+05_JPRB,&
595 ! &0.2123132E+05_JPRB,&
596 ! &0.6659687E+03_JPRB,&
597 ! &0.1000000E+01_JPRB/)
598 !ZBPAD4(2, 2, :) = (/&
599 ! &0.2001271E-01_JPRB,&
600 ! &0.2549067E+01_JPRB,&
601 ! &0.3752433E+02_JPRB,&
602 ! &0.6276637E+02_JPRB,&
603 ! &0.1000000E+01_JPRB,&
604 ! &0.0000000E+00_JPRB,&
605 ! &0.0000000E+00_JPRB/)
606 !ZBPAD4(2, 3, :) = (/&
607 ! &1.000000000E+00_JPRB,&
608 ! &0.000000000E+00_JPRB,&
609 ! &0.000000000E+00_JPRB,&
610 ! &0.000000000E+00_JPRB,&
611 ! &0.000000000E+00_JPRB,&
612 ! &0.000000000E+00_JPRB,&
613 ! &0.000000000E+00_JPRB/)
614 
615 !DATA (ZRAY4(2,K),K=1,6) /&
616 ! & .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,&
617 ! &-.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/
618 
619 !-- Original coefficients in the ECMWF model
620 
621 zd4(2,:)= (/ 0.000000000_jprb, 0.000000000_jprb, 1.000000000_jprb /)
622 
623 zapad4(2, 1, :) = (/&
624  & 0.690730834e-02_jprb,&
625  & 0.151704275e+01_jprb,&
626  & 0.751477543e+02_jprb,&
627  & 0.759770236e+03_jprb,&
628  & 0.109800326e+04_jprb,&
629  & 0.148407574e+03_jprb,&
630  & 0.000000000e+00_jprb/)
631 zapad4(2, 2, :) = (/&
632  & 0.863790752e-03_jprb,&
633  & 0.448762291e+00_jprb,&
634  & 0.332530367e+02_jprb,&
635  & 0.190914146e+03_jprb,&
636  & 0.000000000e+00_jprb,&
637  & 0.000000000e+00_jprb,&
638  & 0.000000000e+00_jprb/)
639 zapad4(2, 3, :) = (/&
640  & 0.000000000e+00_jprb,&
641  & 0.000000000e+00_jprb,&
642  & 0.000000000e+00_jprb,&
643  & 0.000000000e+00_jprb,&
644  & 0.000000000e+00_jprb,&
645  & 0.000000000e+00_jprb,&
646  & 0.000000000e+00_jprb/)
647 
648 zbpad4(2, 1, :) = (/&
649  & 0.690730834e-02_jprb,&
650  & 0.151954406e+01_jprb,&
651  & 0.756512527e+02_jprb,&
652  & 0.779384997e+03_jprb,&
653  & 0.121113108e+04_jprb,&
654  & 0.207678436e+03_jprb,&
655  & 0.100000000e+01_jprb/)
656 zbpad4(2, 2, :) = (/&
657  & 0.863790752e-03_jprb,&
658  & 0.448948107e+00_jprb,&
659  & 0.333186750e+02_jprb,&
660  & 0.192727216e+03_jprb,&
661  & 0.100000000e+01_jprb,&
662  & 0.000000000e+00_jprb,&
663  & 0.000000000e+00_jprb/)
664 zbpad4(2, 3, :) = (/&
665  & 1.000000000e+00_jprb,&
666  & 0.000000000e+00_jprb,&
667  & 0.000000000e+00_jprb,&
668  & 0.000000000e+00_jprb,&
669  & 0.000000000e+00_jprb,&
670  & 0.000000000e+00_jprb,&
671  & 0.000000000e+00_jprb/)
672 
673 DATA (zray4(2,i_k),i_k=1,6) /&
674  & .164261e-01_jprb, .000000e+00_jprb, .000000e+00_jprb,&
675  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/
676 
677 zswce4(2) = 0.0_jprb
678 zswcp4(2) = 0.0_jprb
679 
680 ! ----------------------------------------------------------------
681 
682 !* 1.19 - 2.38 Microns
683 
684 zsun4(3) = 0.180608_jprb
685 
686 !-- New coefficients
687 
688 !ZD4(3,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
689 
690 !ZAPAD4(3, 1, :) = (/&
691 ! &0.3325841E-03_JPRB,&
692 ! &0.6194496E+00_JPRB,&
693 ! &0.1497138E+03_JPRB,&
694 ! &0.2314864E+04_JPRB,&
695 ! &0.2380109E+04_JPRB,&
696 ! &0.9553823E+02_JPRB,&
697 ! &0.0000000E+00_JPRB/)
698 !ZAPAD4(3, 2, :) = (/&
699 ! &0.4552471E-03_JPRB,&
700 ! &0.4084154E+00_JPRB,&
701 ! &0.6114905E+01_JPRB,&
702 ! &0.7102540E+01_JPRB,&
703 ! &0.0000000E+00_JPRB,&
704 ! &0.0000000E+00_JPRB,&
705 ! &0.0000000E+00_JPRB/)
706 !ZAPAD4(3, 3, :) = (/&
707 ! &0.000000000E+00_JPRB,&
708 ! &0.000000000E+00_JPRB,&
709 ! &0.000000000E+00_JPRB,&
710 ! &0.000000000E+00_JPRB,&
711 ! &0.000000000E+00_JPRB,&
712 ! &0.000000000E+00_JPRB,&
713 ! &0.000000000E+00_JPRB/)
714 
715 !ZBPAD4(3, 1, :) = (/&
716 ! &0.3325841E-03_JPRB,&
717 ! &0.6231947E+00_JPRB,&
718 ! &0.1553098E+03_JPRB,&
719 ! &0.2822458E+04_JPRB,&
720 ! &0.3885194E+04_JPRB,&
721 ! &0.2700235E+03_JPRB,&
722 ! &0.1000000E+01_JPRB/)
723 !ZBPAD4(3, 2, :) = (/&
724 ! &0.4552471E-03_JPRB,&
725 ! &0.4088242E+00_JPRB,&
726 ! &0.6411905E+01_JPRB,&
727 ! &0.9444439E+01_JPRB,&
728 ! &0.1000000E+01_JPRB,&
729 ! &0.0000000E+00_JPRB,&
730 ! &0.0000000E+00_JPRB/)
731 !ZBPAD4(3, 3, :) = (/&
732 ! &1.000000000E+00_JPRB,&
733 ! &0.000000000E+00_JPRB,&
734 ! &0.000000000E+00_JPRB,&
735 ! &0.000000000E+00_JPRB,&
736 ! &0.000000000E+00_JPRB,&
737 ! &0.000000000E+00_JPRB,&
738 ! &0.000000000E+00_JPRB/)
739 
740 !ZRAY4(3,:)= (/&
741 ! &.288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,&
742 ! &.540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/)
743 
744 zd4(3,:)= (/ 0.000000000_jprb, 0.000000000_jprb, 1.000000000_jprb /)
745 
746 zapad4(3, 1, :) = (/&
747  & 0.837531303e-05_jprb,&
748  & 0.173886341e-01_jprb,&
749  & 0.518852799e+01_jprb,&
750  & 0.159078416e+03_jprb,&
751  & 0.493273523e+03_jprb,&
752  & 0.102567293e+03_jprb,&
753  & 0.000000000e+00_jprb/)
754 zapad4(3, 2, :) = (/&
755  & 0.657978575e-02_jprb,&
756  & 0.752617872e+00_jprb,&
757  & 0.158209734e+02_jprb,&
758  & 0.410274915e+02_jprb,&
759  & 0.000000000e+00_jprb,&
760  & 0.000000000e+00_jprb,&
761  & 0.000000000e+00_jprb/)
762 zapad4(3, 3, :) = (/&
763  & 0.000000000e+00_jprb,&
764  & 0.000000000e+00_jprb,&
765  & 0.000000000e+00_jprb,&
766  & 0.000000000e+00_jprb,&
767  & 0.000000000e+00_jprb,&
768  & 0.000000000e+00_jprb,&
769  & 0.000000000e+00_jprb/)
770 
771 zbpad4(3, 1, :) = (/&
772  & 0.837531303e-05_jprb,&
773  & 0.174882536e-01_jprb,&
774  & 0.534536580e+01_jprb,&
775  & 0.180351767e+03_jprb,&
776  & 0.673126838e+03_jprb,&
777  & 0.182718543e+03_jprb,&
778  & 0.100000000e+01_jprb/)
779 zbpad4(3, 2, :) = (/&
780  & 0.657978575e-02_jprb,&
781  & 0.753752065e+00_jprb,&
782  & 0.159286262e+02_jprb,&
783  & 0.424278450e+02_jprb,&
784  & 0.100000000e+01_jprb,&
785  & 0.000000000e+00_jprb,&
786  & 0.000000000e+00_jprb/)
787 zbpad4(3, 3, :) = (/&
788  & 1.000000000e+00_jprb,&
789  & 0.000000000e+00_jprb,&
790  & 0.000000000e+00_jprb,&
791  & 0.000000000e+00_jprb,&
792  & 0.000000000e+00_jprb,&
793  & 0.000000000e+00_jprb,&
794  & 0.000000000e+00_jprb/)
795 
796 zray4(3,:)= (/&
797  & .180438e-02_jprb, .000000e+00_jprb, .000000e+00_jprb,&
798  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
799 
800 zswce4(3) = 0.0_jprb
801 zswcp4(3) = 0.0_jprb
802 
803 ! ----------------------------------------------------------------
804 
805 !* 2.38 - 4.00 Microns
806 
807 zsun4(4) = 0.033474_jprb
808 
809 !-- New coefficients
810 
811 !ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
812 
813 !ZAPAD4(4, 1, :) = (/&
814 ! &0.2122889E-06_JPRB,&
815 ! &0.9030576E-03_JPRB,&
816 ! &0.2431282E+00_JPRB,&
817 ! &0.4901345E+01_JPRB,&
818 ! &0.3996347E+01_JPRB,&
819 ! &0.3910227E+01_JPRB,&
820 ! &0.0000000E+00_JPRB/)
821 !ZAPAD4(4, 2, :) = (/&
822 ! &0.1215163E-03_JPRB,&
823 ! &0.1222574E+00_JPRB,&
824 ! &0.9382420E+01_JPRB,&
825 ! &0.6875727E+02_JPRB,&
826 ! &0.2746421E+02_JPRB,&
827 ! &0.0000000E+00_JPRB,&
828 ! &0.0000000E+00_JPRB/)
829 !ZAPAD4(4, 3, :) = (/&
830 ! &0.263068898E+02_JPRB,&
831 ! &0.146425875E+03_JPRB,&
832 ! &0.860137809E+02_JPRB,&
833 ! &0.000000000E+00_JPRB,&
834 ! &0.000000000E+00_JPRB,&
835 ! &0.000000000E+00_JPRB,&
836 ! &0.000000000E+00_JPRB/)
837 
838 !ZBPAD4(4, 1, :) = (/&
839 ! &0.2122889E-06_JPRB,&
840 ! &0.9379083E-03_JPRB,&
841 ! &0.2957335E+00_JPRB,&
842 ! &0.8747190E+01_JPRB,&
843 ! &0.1015794E+02_JPRB,&
844 ! &0.1361277E+02_JPRB,&
845 ! &0.1000000E+01_JPRB/)
846 !ZBPAD4(4, 2, :) = (/&
847 ! &0.1215163E-03_JPRB,&
848 ! &0.1255648E+00_JPRB,&
849 ! &0.1060119E+02_JPRB,&
850 ! &0.8414439E+02_JPRB,&
851 ! &0.4299438E+02_JPRB,&
852 ! &0.1000000E+01_JPRB,&
853 ! &0.0000000E+00_JPRB/)
854 !ZBPAD4(4, 3, :) = (/&
855 ! &0.263068898E+02_JPRB,&
856 ! &0.152569217E+03_JPRB,&
857 ! &0.976791971E+02_JPRB,&
858 ! &0.100000000E+01_JPRB,&
859 ! &0.000000000E+00_JPRB,&
860 ! &0.000000000E+00_JPRB,&
861 ! &0.000000000E+00_JPRB/)
862 
863 !ZRAY4(4,:)= (/&
864 ! &.493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,&
865 ! &.741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/)
866 
867 zd4(4,:)= (/ 0.000000000_jprb, 0.000000000_jprb, 0.000000000_jprb /)
868 
869 zapad4(4, 1, :) = (/&
870  & 0.122118185e-06_jprb,&
871  & 0.154042531e-02_jprb,&
872  & 0.141152193e+01_jprb,&
873  & 0.685368761e+02_jprb,&
874  & 0.216522281e+03_jprb,&
875  & 0.421228746e+02_jprb,&
876  & 0.000000000e+00_jprb/)
877 zapad4(4, 2, :) = (/&
878  & 0.364233560e-10_jprb,&
879  & 0.217340835e-06_jprb,&
880  & 0.292623386e-03_jprb,&
881  & 0.797100631e-01_jprb,&
882  & 0.319103672e+01_jprb,&
883  & 0.110530283e+02_jprb,&
884  & 0.000000000e+00_jprb/)
885 zapad4(4, 3, :) = (/&
886  & 0.263068898e+02_jprb,&
887  & 0.146425875e+03_jprb,&
888  & 0.860137809e+02_jprb,&
889  & 0.000000000e+00_jprb,&
890  & 0.000000000e+00_jprb,&
891  & 0.000000000e+00_jprb,&
892  & 0.000000000e+00_jprb/)
893 
894 zbpad4(4, 1, :) = (/&
895  & 0.122118185e-06_jprb,&
896  & 0.156287582e-02_jprb,&
897  & 0.156994562e+01_jprb,&
898  & 0.102304103e+03_jprb,&
899  & 0.475479878e+03_jprb,&
900  & 0.188714799e+03_jprb,&
901  & 0.100000000e+01_jprb/)
902 zbpad4(4, 2, :) = (/&
903  & 0.364233560e-10_jprb,&
904  & 0.218265806e-06_jprb,&
905  & 0.297085917e-03_jprb,&
906  & 0.834253403e-01_jprb,&
907  & 0.358290002e+01_jprb,&
908  & 0.139206164e+02_jprb,&
909  & 0.100000000e+01_jprb/)
910 zbpad4(4, 3, :) = (/&
911  & 0.263068898e+02_jprb,&
912  & 0.152569217e+03_jprb,&
913  & 0.976791971e+02_jprb,&
914  & 0.100000000e+01_jprb,&
915  & 0.000000000e+00_jprb,&
916  & 0.000000000e+00_jprb,&
917  & 0.000000000e+00_jprb/)
918 
919 zray4(4,:)= (/&
920  & .136020e-03_jprb, .000000e+00_jprb, .000000e+00_jprb,&
921  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
922 
923 zswce4(4) = 0.0_jprb
924 zswcp4(4) = 0.0_jprb
925 
926 !=====================================================================
927 
928 !* 1.3 COEFFICIENTS FOR SIX SPECTRAL INTERVALS
929 ! ---------------------------------------
930 
931 !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
932 ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
933 ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
934 ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
935 
936 zadjus6 = 0.9905709_jprb
937 zadjus6 = 1.0_jprb
938 
939 ztdh6is = 0.450_jprb
940 ztdu6is = 0.375_jprb
941 zth6is = 273._jprb
942 ztu6is = 273._jprb
943 zpdh6is = 0.90_jprb
944 zpdu6is = 0.75_jprb
945 zprh6is = 101300._jprb
946 zpru6is = 101300._jprb
947 
948 !* 1st spectral interval: U.V. (0.18 - 0.25 Micron)
949 
950 zsun6(1) = 0.001917_jprb
951 
952 zd6(1,:)= (/ 1.000000000_jprb, 1.000000000_jprb, 0.000000000_jprb /)
953 
954 zapad6(1, 1, :) = (/&
955  & 0.000000000e-00_jprb,&
956  & 0.000000000e-00_jprb,&
957  & 0.000000000e-00_jprb,&
958  & 0.000000000e-00_jprb,&
959  & 0.000000000e-00_jprb,&
960  & 0.000000000e-00_jprb,&
961  & 0.000000000e-00_jprb/)
962 zapad6(1, 2, :) = (/&
963  & 0.000000000e-00_jprb,&
964  & 0.000000000e-00_jprb,&
965  & 0.000000000e-00_jprb,&
966  & 0.000000000e-00_jprb,&
967  & 0.000000000e-00_jprb,&
968  & 0.000000000e-00_jprb,&
969  & 0.000000000e-00_jprb/)
970 zapad6(1, 3, :) = (/&
971  & 0.000000000e-00_jprb,&
972  & 0.000000000e-00_jprb,&
973  & 0.000000000e-00_jprb,&
974  & 0.000000000e-00_jprb,&
975  & 0.000000000e-00_jprb,&
976  & 0.000000000e-00_jprb,&
977  & 0.000000000e-00_jprb/)
978 
979 zbpad6(1, 1, :) = (/&
980  & 0.100000000e+01_jprb,&
981  & 0.000000000e-00_jprb,&
982  & 0.000000000e-00_jprb,&
983  & 0.000000000e-00_jprb,&
984  & 0.000000000e-00_jprb,&
985  & 0.000000000e-00_jprb,&
986  & 0.000000000e-00_jprb/)
987 zbpad6(1, 2, :) = (/&
988  & 0.100000000e+01_jprb,&
989  & 0.000000000e-00_jprb,&
990  & 0.000000000e-00_jprb,&
991  & 0.000000000e-00_jprb,&
992  & 0.000000000e-00_jprb,&
993  & 0.000000000e-00_jprb,&
994  & 0.000000000e-00_jprb/)
995 zbpad6(1, 3, :) = (/&
996  & 0.100000000e+01_jprb,&
997  & 0.000000000e-00_jprb,&
998  & 0.000000000e-00_jprb,&
999  & 0.000000000e-00_jprb,&
1000  & 0.000000000e-00_jprb,&
1001  & 0.000000000e-00_jprb,&
1002  & 0.000000000e-00_jprb/)
1003 
1004 zray6(1,:)= (/&
1005  & .334018e+00_jprb,-.298776e+01_jprb, .613911e+02_jprb,&
1006  & -.164331e+03_jprb, .197946e+03_jprb,-.652206e+02_jprb/)
1007 zrby6(1,:)= (/&
1008  & .43959002e+01_jprb, .000000e+00_jprb, .000000e+00_jprb,&
1009  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
1010 
1011 zswce6(1) = 0.0_jprb
1012 zswcp6(1) = 0.0_jprb
1013 
1014 nexpo3(1) = 7
1015 rexpo3(1, 1, :) =(/&
1016  & 0.051395e+00_jprb, 0.048250e+00_jprb, 0.112339e+00_jprb,&
1017  & 0.101426e+00_jprb, 0.007700e+00_jprb, 0.441320e+00_jprb,&
1018  & 0.237571e+00_jprb /)
1019 rexpo3(1, 2, :) =(/&
1020  & 0.100022e+02_jprb, 0.851159e+02_jprb, 0.346737e+03_jprb,&
1021  & 0.158501e+02_jprb, 0.724223e+01_jprb, 0.177828e+03_jprb,&
1022  & 0.467708e+02_jprb /)
1023 
1024 !* 2nd spectral interval: U.V. (0.25 - 0.44 Micron)
1025 
1026 zsun6(2) = 0.135708_jprb
1027 
1028 zd6(2,:)= (/ 1.000000000_jprb, 1.000000000_jprb, 0.000000000_jprb /)
1029 
1030 zapad6(2, 1, :) = (/&
1031  & 0.000000000e-00_jprb,&
1032  & 0.000000000e-00_jprb,&
1033  & 0.000000000e-00_jprb,&
1034  & 0.000000000e-00_jprb,&
1035  & 0.000000000e-00_jprb,&
1036  & 0.000000000e-00_jprb,&
1037  & 0.000000000e-00_jprb/)
1038 zapad6(2, 2, :) = (/&
1039  & 0.000000000e-00_jprb,&
1040  & 0.000000000e-00_jprb,&
1041  & 0.000000000e-00_jprb,&
1042  & 0.000000000e-00_jprb,&
1043  & 0.000000000e-00_jprb,&
1044  & 0.000000000e-00_jprb,&
1045  & 0.000000000e-00_jprb/)
1046 zapad6(2, 3, :) = (/&
1047  & 0.000000000e-00_jprb,&
1048  & 0.000000000e-00_jprb,&
1049  & 0.000000000e-00_jprb,&
1050  & 0.000000000e-00_jprb,&
1051  & 0.000000000e-00_jprb,&
1052  & 0.000000000e-00_jprb,&
1053  & 0.000000000e-00_jprb/)
1054 
1055 zbpad6(2, 1, :) = (/&
1056  & 0.100000000e+01_jprb,&
1057  & 0.000000000e-00_jprb,&
1058  & 0.000000000e-00_jprb,&
1059  & 0.000000000e-00_jprb,&
1060  & 0.000000000e-00_jprb,&
1061  & 0.000000000e-00_jprb,&
1062  & 0.000000000e-00_jprb/)
1063 zbpad6(2, 2, :) = (/&
1064  & 0.100000000e+01_jprb,&
1065  & 0.000000000e-00_jprb,&
1066  & 0.000000000e-00_jprb,&
1067  & 0.000000000e-00_jprb,&
1068  & 0.000000000e-00_jprb,&
1069  & 0.000000000e-00_jprb,&
1070  & 0.000000000e-00_jprb/)
1071 zbpad6(2, 3, :) = (/&
1072  & 0.100000000e+01_jprb,&
1073  & 0.000000000e-00_jprb,&
1074  & 0.000000000e-00_jprb,&
1075  & 0.000000000e-00_jprb,&
1076  & 0.000000000e-00_jprb,&
1077  & 0.000000000e-00_jprb,&
1078  & 0.000000000e-00_jprb/)
1079 
1080 zray6(2,:)= (/&
1081  & -.895408e-03_jprb, .283130e+01_jprb,-.867917e+01_jprb,&
1082  & .148672e+02_jprb,-.127282e+02_jprb, .422212e+01_jprb/)
1083 zrby6(2,:)= (/&
1084  & .55503070e+00_jprb, .000000e+00_jprb, .000000e+00_jprb,&
1085  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
1086 
1087 zswce6(2) = 0.0_jprb
1088 zswcp6(2) = 0.0_jprb
1089 
1090 nexpo3(2) = 7
1091 rexpo3(2, 1, :) =(/&
1092  & 0.043801e+00_jprb, 0.078893e+00_jprb, 0.036839e+00_jprb,&
1093  & 0.022503e+00_jprb, 0.042333e+00_jprb, 0.037870e+00_jprb,&
1094  & 0.737762e+00_jprb /)
1095 rexpo3(2, 2, :) =(/&
1096  & 0.234249e+01_jprb, 0.125170e+00_jprb, 0.549527e+02_jprb,&
1097  & 0.257041e+03_jprb, 0.476838e+00_jprb, 0.911993e+01_jprb,&
1098  & 0.000000e+00_jprb /)
1099 
1100 !* 3rd spectral interval: Visible (0.44 - 0.69 Micron)
1101 
1102 zsun6(3) = 0.322135_jprb
1103 
1104 zd6(3,:)= (/ 0.800000000_jprb, 0.900000000_jprb, 0.000000000_jprb /)
1105 
1106 zapad6(3, 1, :) = (/&
1107  & 0.1762097e+03_jprb,&
1108  & 0.1641762e+03_jprb,&
1109  & 0.8687919e+02_jprb,&
1110  & 0.0000000e-00_jprb,&
1111  & 0.0000000e-00_jprb,&
1112  & 0.0000000e-00_jprb,&
1113  & 0.0000000e-00_jprb/)
1114 zapad6(3, 2, :) = (/&
1115  & 0.5581224e+00_jprb,&
1116  & 0.1748430e+03_jprb,&
1117  & 0.1134123e+04_jprb,&
1118  & 0.3490429e+03_jprb,&
1119  & 0.0000000e-00_jprb,&
1120  & 0.0000000e-00_jprb,&
1121  & 0.0000000e-00_jprb/)
1122 zapad6(3, 3, :) = (/&
1123  & 0.000000000e-00_jprb,&
1124  & 0.000000000e-00_jprb,&
1125  & 0.000000000e-00_jprb,&
1126  & 0.000000000e-00_jprb,&
1127  & 0.000000000e-00_jprb,&
1128  & 0.000000000e-00_jprb,&
1129  & 0.000000000e-00_jprb/)
1130 
1131 zbpad6(3, 1, :) = (/&
1132  & 0.1762097e+03_jprb,&
1133  & 0.1663950e+03_jprb,&
1134  & 0.8939724e+02_jprb,&
1135  & 0.1000000e+01_jprb,&
1136  & 0.0000000e-00_jprb,&
1137  & 0.0000000e-00_jprb,&
1138  & 0.0000000e-00_jprb/)
1139 zbpad6(3, 2, :) = (/&
1140  & 0.5581224e+00_jprb,&
1141  & 0.1749251e+03_jprb,&
1142  & 0.1159910e+04_jprb,&
1143  & 0.3893268e+03_jprb,&
1144  & 0.1000000e+01_jprb,&
1145  & 0.0000000e-00_jprb,&
1146  & 0.0000000e-00_jprb/)
1147 zbpad6(3, 3, :) = (/&
1148  & 0.100000000e+01_jprb,&
1149  & 0.000000000e-00_jprb,&
1150  & 0.000000000e-00_jprb,&
1151  & 0.000000000e-00_jprb,&
1152  & 0.000000000e-00_jprb,&
1153  & 0.000000000e-00_jprb,&
1154  & 0.000000000e-00_jprb/)
1155 
1156 zray6(3,:)= (/&
1157  & .474852e-01_jprb, .605110e+00_jprb,-.248070e+01_jprb,&
1158  & .485580e+01_jprb,-.444043e+01_jprb, .152369e+01_jprb/)
1159 zrby6(3,:)= (/&
1160  & .10528199e+00_jprb, .000000e+00_jprb, .000000e+00_jprb,&
1161  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
1162 
1163 zswce6(3) = 0.0_jprb
1164 zswcp6(3) = 0.0_jprb
1165 
1166 nexpo3(3) = 6
1167 rexpo3(3, 1, :) =(/&
1168  & 0.063442e+00_jprb, 0.058550e+00_jprb, 0.237534e+00_jprb,&
1169  & 0.412292e+00_jprb, 0.126141e+00_jprb, 0.102041e+00_jprb,&
1170  & 0.000000e+00_jprb /)
1171 rexpo3(3, 2, :) =(/&
1172  & 0.125170e+00_jprb, 0.119209e-01_jprb, 0.119209e+00_jprb,&
1173  & 0.417233e-01_jprb, 0.894070e-01_jprb, 0.000000e+00_jprb,&
1174  & 0.000000e+00_jprb /)
1175 
1176 ! ----------------------------------------------------------------
1177 
1178 !* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
1179 
1180 ! ----------------------------------------------------------------
1181 
1182 !* 0.69 - 1.19 Micron
1183 
1184 !* UMG is O2 only
1185 
1186 zsun6(4) = 0.326158_jprb
1187 
1188 zd6(4,:)= (/ 0.000000000_jprb, 0.900000000_jprb, 1.000000000_jprb /)
1189 
1190 zapad6(4, 1, :) = (/&
1191  & 0.1335726e+02_jprb,&
1192  & 0.2939136e+04_jprb,&
1193  & 0.4010585e+05_jprb,&
1194  & 0.7195030e+05_jprb,&
1195  & 0.1648338e+05_jprb,&
1196  & 0.3373738e+03_jprb,&
1197  & 0.0000000e+00_jprb/)
1198 zapad6(4, 2, :) = (/&
1199  & 0.2001271e-01_jprb,&
1200  & 0.2480831e+01_jprb,&
1201  & 0.3444162e+02_jprb,&
1202  & 0.4788946e+02_jprb,&
1203  & 0.0000000e+00_jprb,&
1204  & 0.0000000e+00_jprb,&
1205  & 0.0000000e+00_jprb/)
1206 zapad6(4, 3, :) = (/&
1207  & 0.000000000e+00_jprb,&
1208  & 0.000000000e+00_jprb,&
1209  & 0.000000000e+00_jprb,&
1210  & 0.000000000e+00_jprb,&
1211  & 0.000000000e+00_jprb,&
1212  & 0.000000000e+00_jprb,&
1213  & 0.000000000e+00_jprb/)
1214 
1215 zbpad6(4, 1, :) = (/&
1216  & 0.1335726e+02_jprb,&
1217  & 0.2942327e+04_jprb,&
1218  & 0.4077237e+05_jprb,&
1219  & 0.7749017e+05_jprb,&
1220  & 0.2123132e+05_jprb,&
1221  & 0.6659687e+03_jprb,&
1222  & 0.1000000e+01_jprb/)
1223 zbpad6(4, 2, :) = (/&
1224  & 0.2001271e-01_jprb,&
1225  & 0.2549067e+01_jprb,&
1226  & 0.3752433e+02_jprb,&
1227  & 0.6276637e+02_jprb,&
1228  & 0.1000000e+01_jprb,&
1229  & 0.0000000e+00_jprb,&
1230  & 0.0000000e+00_jprb/)
1231 zbpad6(4, 3, :) = (/&
1232  & 1.000000000e+00_jprb,&
1233  & 0.000000000e+00_jprb,&
1234  & 0.000000000e+00_jprb,&
1235  & 0.000000000e+00_jprb,&
1236  & 0.000000000e+00_jprb,&
1237  & 0.000000000e+00_jprb,&
1238  & 0.000000000e+00_jprb/)
1239 
1240 zray6(4,:)= (/&
1241  & .519710e-01_jprb,-.150667e+00_jprb, .154170e+01_jprb,&
1242  & -.518901e+01_jprb, .668011e+01_jprb,-.294843e+01_jprb/)
1243 zrby6(4,:)= (/&
1244  & .16436996e-01_jprb, .000000e+00_jprb, .000000e+00_jprb,&
1245  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
1246 
1247 zswce6(4) = 0.0_jprb
1248 zswcp6(4) = 0.0_jprb
1249 
1250 nexpo3(4) = 4
1251 rexpo3(4, 1, :) =(/&
1252  & 0.000074e+00_jprb, 0.320194e+00_jprb, 0.082915e+00_jprb,&
1253  & 0.596816e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1254  & 0.000000e+00_jprb /)
1255 rexpo3(4, 2, :) =(/&
1256  & 0.232458e+00_jprb, 0.119209e-01_jprb, 0.178814e-01_jprb,&
1257  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1258  & 0.000000e+00_jprb /)
1259 
1260 ! ----------------------------------------------------------------
1261 
1262 !* 1.19 - 2.38 Microns
1263 
1264 !* UMG is CO2 only
1265 
1266 zsun6(5) = 0.180608_jprb
1267 
1268 zd6(5,:)= (/ 0.000000000_jprb, 0.800000000_jprb, 1.000000000_jprb /)
1269 
1270 zapad6(5, 1, :) = (/&
1271  & 0.3325841e-03_jprb,&
1272  & 0.6194496e+00_jprb,&
1273  & 0.1497138e+03_jprb,&
1274  & 0.2314864e+04_jprb,&
1275  & 0.2380109e+04_jprb,&
1276  & 0.9553823e+02_jprb,&
1277  & 0.0000000e+00_jprb/)
1278 zapad6(5, 2, :) = (/&
1279  & 0.4552471e-03_jprb,&
1280  & 0.4084154e+00_jprb,&
1281  & 0.6114905e+01_jprb,&
1282  & 0.7102540e+01_jprb,&
1283  & 0.0000000e+00_jprb,&
1284  & 0.0000000e+00_jprb,&
1285  & 0.0000000e+00_jprb/)
1286 zapad6(5, 3, :) = (/&
1287  & 0.000000000e+00_jprb,&
1288  & 0.000000000e+00_jprb,&
1289  & 0.000000000e+00_jprb,&
1290  & 0.000000000e+00_jprb,&
1291  & 0.000000000e+00_jprb,&
1292  & 0.000000000e+00_jprb,&
1293  & 0.000000000e+00_jprb/)
1294 
1295 zbpad6(5, 1, :) = (/&
1296  & 0.3325841e-03_jprb,&
1297  & 0.6231947e+00_jprb,&
1298  & 0.1553098e+03_jprb,&
1299  & 0.2822458e+04_jprb,&
1300  & 0.3885194e+04_jprb,&
1301  & 0.2700235e+03_jprb,&
1302  & 0.1000000e+01_jprb/)
1303 zbpad6(5, 2, :) = (/&
1304  & 0.4552471e-03_jprb,&
1305  & 0.4088242e+00_jprb,&
1306  & 0.6411905e+01_jprb,&
1307  & 0.9444439e+01_jprb,&
1308  & 0.1000000e+01_jprb,&
1309  & 0.0000000e+00_jprb,&
1310  & 0.0000000e+00_jprb/)
1311 zbpad6(5, 3, :) = (/&
1312  & 1.000000000e+00_jprb,&
1313  & 0.000000000e+00_jprb,&
1314  & 0.000000000e+00_jprb,&
1315  & 0.000000000e+00_jprb,&
1316  & 0.000000000e+00_jprb,&
1317  & 0.000000000e+00_jprb,&
1318  & 0.000000000e+00_jprb/)
1319 
1320 zray6(5,:)= (/&
1321  & .288474e-02_jprb, .133984e-02_jprb,-.152812e-01_jprb,&
1322  & .540845e-01_jprb,-.747444e-01_jprb, .364285e-01_jprb/)
1323 zrby6(5,:)= (/&
1324  & .18073079e-02_jprb, .000000e+00_jprb, .000000e+00_jprb,&
1325  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
1326 
1327 zswce6(5) = 0.0_jprb
1328 zswcp6(5) = 0.0_jprb
1329 
1330 nexpo3(5) = 0
1331 rexpo3(5, 1, :) =(/&
1332  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1333  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1334  & 0.000000e+00_jprb /)
1335 rexpo3(5, 2, :) =(/&
1336  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1337  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1338  & 0.000000e+00_jprb /)
1339 
1340 ! ----------------------------------------------------------------
1341 
1342 !* 2.38 - 4.00 Microns
1343 
1344 zsun6(6) = 0.033474_jprb
1345 
1346 zd6(6,:)= (/ 0.000000000_jprb, 0.000000000_jprb, 0.000000000_jprb /)
1347 
1348 zapad6(6, 1, :) = (/&
1349  & 0.2122889e-06_jprb,&
1350  & 0.9030576e-03_jprb,&
1351  & 0.2431282e+00_jprb,&
1352  & 0.4901345e+01_jprb,&
1353  & 0.3996347e+01_jprb,&
1354  & 0.3910227e+01_jprb,&
1355  & 0.0000000e+00_jprb/)
1356 zapad6(6, 2, :) = (/&
1357  & 0.1215163e-03_jprb,&
1358  & 0.1222574e+00_jprb,&
1359  & 0.9382420e+01_jprb,&
1360  & 0.6875727e+02_jprb,&
1361  & 0.2746421e+02_jprb,&
1362  & 0.0000000e+00_jprb,&
1363  & 0.0000000e+00_jprb/)
1364 zapad6(6, 3, :) = (/&
1365  & 0.263068898e+02_jprb,&
1366  & 0.146425875e+03_jprb,&
1367  & 0.860137809e+02_jprb,&
1368  & 0.000000000e+00_jprb,&
1369  & 0.000000000e+00_jprb,&
1370  & 0.000000000e+00_jprb,&
1371  & 0.000000000e+00_jprb/)
1372 
1373 zbpad6(6, 1, :) = (/&
1374  & 0.2122889e-06_jprb,&
1375  & 0.9379083e-03_jprb,&
1376  & 0.2957335e+00_jprb,&
1377  & 0.8747190e+01_jprb,&
1378  & 0.1015794e+02_jprb,&
1379  & 0.1361277e+02_jprb,&
1380  & 0.1000000e+01_jprb/)
1381 zbpad6(6, 2, :) = (/&
1382  & 0.1215163e-03_jprb,&
1383  & 0.1255648e+00_jprb,&
1384  & 0.1060119e+02_jprb,&
1385  & 0.8414439e+02_jprb,&
1386  & 0.4299438e+02_jprb,&
1387  & 0.1000000e+01_jprb,&
1388  & 0.0000000e+00_jprb/)
1389 zbpad6(6, 3, :) = (/&
1390  & 0.263068898e+02_jprb,&
1391  & 0.152569217e+03_jprb,&
1392  & 0.976791971e+02_jprb,&
1393  & 0.100000000e+01_jprb,&
1394  & 0.000000000e+00_jprb,&
1395  & 0.000000000e+00_jprb,&
1396  & 0.000000000e+00_jprb/)
1397 
1398 zray6(6,:)= (/&
1399  & .493545e-03_jprb, .355696e-03_jprb,-.284543e-02_jprb,&
1400  & .741248e-02_jprb,-.773691e-02_jprb, .370679e-02_jprb/)
1401 zrby6(6,:)= (/&
1402  & .13618247e-03_jprb, .000000e+00_jprb, .000000e+00_jprb,&
1403  & .000000e+00_jprb, .000000e+00_jprb, .000000e+00_jprb/)
1404 
1405 zswce6(6) = 0.0_jprb
1406 zswcp6(6) = 0.0_jprb
1407 
1408 nexpo3(6) = 0
1409 rexpo3(6, 1, :) =(/&
1410  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1411  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1412  & 0.000000e+00_jprb /)
1413 rexpo3(6, 2, :) =(/&
1414  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1415  & 0.000000e+00_jprb, 0.000000e+00_jprb, 0.000000e+00_jprb,&
1416  & 0.000000e+00_jprb /)
1417 
1418 !=====================================================================
1419 !* 2.4 OPTICAL PARAMETERS FOR RAIN DROPS
1420 ! Savijarvi et al. (1996)
1421 
1422 zrto1 = 0.003_jprb
1423 zrto2 = -0.22_jprb
1424 ! CAUTION JUST TEMPORARY PARAMETERS
1425 zroma2(1:2)= (/ 0.00008_jprb , 0.264_jprb /)
1426 zromb2(1:2)= (/ 0.23_jprb , 0.09_jprb /)
1427 zrasy2(1:2)= (/ 0.88_jprb , 0.94_jprb /)
1428 
1429 zra2(1:2)= (/ 1.5_jprb , 1.5_jprb /)
1430 zrb2(1:2)= (/ 0.50_jprb , 1.13_jprb /)
1431 zrc2(1:2)= (/ 5.58e-7_jprb , 8.55e-4_jprb/)
1432 zrd2(1:2)= (/ 1.25e-7_jprb , 1.28e-3_jprb/)
1433 zre2(1:2)= (/ 0.841_jprb , 0.786_jprb /)
1434 zrf2(1:2)= (/ 2.08e-3_jprb , 5.32e-3_jprb/)
1435 
1436 zroma4(1:4)= (/ 0.00008_jprb , 0.0105_jprb , 0.264_jprb , 0.465_jprb /)
1437 zromb4(1:4)= (/ 0.23_jprb , 0.22_jprb , 0.09_jprb , 0.001_jprb /)
1438 zrasy4(1:4)= (/ 0.88_jprb , 0.89_jprb , 0.94_jprb , 0.97_jprb /)
1439 
1440 zra4(1:4)= (/ 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb /)
1441 zrb4(1:4)= (/ 0.50_jprb , 0.78_jprb , 1.13_jprb , 2.00_jprb /)
1442 zrc4(1:4)= (/ 5.58e-7_jprb , 2.18e-5_jprb, 8.55e-4_jprb, 1.94e-1_jprb /)
1443 zrd4(1:4)= (/ 1.25e-7_jprb , 2.25e-5_jprb, 1.28e-3_jprb, 8.04e-3_jprb /)
1444 zre4(1:4)= (/ 0.841_jprb , 0.821_jprb , 0.786_jprb , 0.820_jprb /)
1445 zrf4(1:4)= (/ 2.08e-3_jprb , 3.06e-3_jprb, 5.32e-3_jprb, 5.59e-3_jprb /)
1446 
1447 zroma6(1:6)= (/ 0.00008_jprb , 0.00008_jprb , 0.00008_jprb , 0.0105_jprb , 0.264_jprb , 0.465_jprb /)
1448 zromb6(1:6)= (/ 0.23_jprb , 0.23_jprb , 0.23_jprb , 0.22_jprb , 0.09_jprb , 0.001_jprb /)
1449 zrasy6(1:6)= (/ 0.88_jprb , 0.88_jprb , 0.88_jprb , 0.89_jprb , 0.94_jprb , 0.97_jprb /)
1450 
1451 zra6(1:6)= (/ 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb /)
1452 zrb6(1:6)= (/ 0.50_jprb , 0.50_jprb , 0.50_jprb , 0.78_jprb , 1.13_jprb , 2.00_jprb /)
1453 zrc6(1:6)= (/ 5.58e-7_jprb , 5.58e-7_jprb , 5.58e-7_jprb , 2.18e-5_jprb, 8.55e-4_jprb, 1.94e-1_jprb /)
1454 zrd6(1:6)= (/ 1.25e-7_jprb , 1.25e-7_jprb , 1.25e-7_jprb , 2.25e-5_jprb, 1.28e-3_jprb, 8.04e-3_jprb /)
1455 zre6(1:6)= (/ 0.841_jprb , 0.841_jprb , 0.841_jprb , 0.821_jprb , 0.786_jprb , 0.820_jprb /)
1456 zrf6(1:6)= (/ 2.08e-3_jprb , 2.08e-3_jprb , 2.08e-3_jprb , 3.06e-3_jprb, 5.32e-3_jprb, 5.59e-3_jprb /)
1457 
1458 !=====================================================================
1459 ! 3. FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW
1460 
1461 ! SW : 14 spectral intervals : mapping to SW6 spectral intervals
1462 ! 0.185-0.25-0.44-0.69-1.19-2.38-4.00
1463 ! 3.846 - 3.077 6
1464 ! 3.077 - 2.500 6
1465 ! 2.500 - 2.150 5
1466 ! 2.150 - 1.942 5
1467 ! 1.942 - 1.626 5
1468 ! 1.626 - 1.299 5
1469 ! 1.299 - 1.242 5
1470 ! 1.242 - 0.7782 4
1471 ! 0.7782- 0.6250 4
1472 ! 0.6250- 0.4415 3
1473 ! 0.4415- 0.3448 2
1474 ! 0.3448- 0.2632 2
1475 ! 0.2632- 0.2000 1
1476 ! 12.195 - 3.846 6
1477 
1478 impsrtm(1:14)= (/ 6,6,5,5,5,5,5,4,4,3,2,2,1,6 /)
1479 !=====================================================================
1480 
1481 !* 4. SET VALUES.
1482 ! -----------
1483 
1484 IF (lonewsw) THEN
1485  zpdh2o = zpdhnis
1486  zpdumg = zpdunis
1487  zprh2o = zprhnis
1488  zprumg = zprunis
1489  rtdh2o = ztdhnis
1490  rtdumg = ztdunis
1491  rth2o = zthnis
1492  rtumg = ztunis
1493 ENDIF
1494 
1495 rpdh1=zpdh2o+1.0_jprb
1496 rpdu1=zpdumg+1.0_jprb
1497 zh2o=1.0_jprb/( 10._jprb* rg * rpdh1 )
1498 zumg=1.0_jprb/( 10._jprb* rg * rpdu1 )
1499 rpnu = zumg/(zprumg**zpdumg)
1500 rpnh = zh2o/(zprh2o**zpdh2o)
1501 
1502 IF(.NOT.ALLOCATED(rsun)) ALLOCATE (rsun(ntsw))
1503 
1504 DO jw=1,isw6
1505  rsun(jw)=0.0_jprb
1506  rswce(jw)=0.0_jprb
1507  rswcp(jw)=0.0_jprb
1508 
1509  DO jc3=1,3
1510  d(jw,jc3)=0.0_jprb
1511  ENDDO
1512  DO jc6=1,6
1513  rray(jw,jc6)=0.0_jprb
1514  ENDDO
1515  DO ji=1,3
1516  DO jj=1,7
1517  apad(jw,ji,jj)=0.0_jprb
1518  bpad(jw,ji,jj)=0.0_jprb
1519  ENDDO
1520  ENDDO
1521 ENDDO
1522 
1523 IF (ksw == 2) THEN
1524  zpdh2o = zpdh2is
1525  zpdumg = zpdu2is
1526  zprh2o = zprh2is
1527  zprumg = zpru2is
1528  rtdh2o = ztdh2is
1529  rtdumg = ztdu2is
1530  rth2o = zth2is
1531  rtumg = ztu2is
1532 
1533  radjust=zadjus2
1534  rhsrta=zrto1
1535  rhsrtb=zrto2
1536  DO jw=1,ksw
1537  rsun(jw)=zsun2(jw)
1538 
1539  rswce(jw)=zswce2(jw)
1540  rswcp(jw)=zswcp2(jw)
1541 
1542  rroma(jw)=zroma2(jw)
1543  rromb(jw)=zromb2(jw)
1544  rrasy(jw)=zrasy2(jw)
1545  rhsra(jw)=zra2(jw)
1546  rhsrb(jw)=zrb2(jw)
1547  rhsrc(jw)=zrc2(jw)
1548  rhsrd(jw)=zrd2(jw)
1549  rhsre(jw)=zre2(jw)
1550  rhsrf(jw)=zrf2(jw)
1551 
1552  DO jc3=1,3
1553  d(jw,jc3)=zd2(jw,jc3)
1554  ENDDO
1555  DO jc6=1,6
1556  rray(jw,jc6)=zray2(jw,jc6)
1557  ENDDO
1558  DO ji=1,3
1559  DO jj=1,7
1560  apad(jw,ji,jj)=zapad2(jw,ji,jj)
1561  bpad(jw,ji,jj)=zbpad2(jw,ji,jj)
1562  ENDDO
1563  ENDDO
1564  ENDDO
1565 
1566 ELSEIF (ksw == 4) THEN
1567  zpdh2o = zpdh4is
1568  zpdumg = zpdu4is
1569  zprh2o = zprh4is
1570  zprumg = zpru4is
1571  rtdh2o = ztdh4is
1572  rtdumg = ztdu4is
1573  rth2o = zth4is
1574  rtumg = ztu4is
1575 
1576  radjust=zadjus4
1577  rhsrta=zrto1
1578  rhsrtb=zrto2
1579  DO jw=1,ksw
1580  rsun(jw)=zsun4(jw)
1581 
1582  rswce(jw)=zswce4(jw)
1583  rswcp(jw)=zswcp4(jw)
1584 
1585  rroma(jw)=zroma4(jw)
1586  rromb(jw)=zromb4(jw)
1587  rrasy(jw)=zrasy4(jw)
1588  rhsra(jw)=zra4(jw)
1589  rhsrb(jw)=zrb4(jw)
1590  rhsrc(jw)=zrc4(jw)
1591  rhsrd(jw)=zrd4(jw)
1592  rhsre(jw)=zre4(jw)
1593  rhsrf(jw)=zrf4(jw)
1594 
1595  DO jc3=1,3
1596  d(jw,jc3)=zd4(jw,jc3)
1597  ENDDO
1598  DO jc6=1,6
1599  rray(jw,jc6)=zray4(jw,jc6)
1600  ENDDO
1601  DO ji=1,3
1602  DO jj=1,7
1603  apad(jw,ji,jj)=zapad4(jw,ji,jj)
1604  bpad(jw,ji,jj)=zbpad4(jw,ji,jj)
1605  ENDDO
1606  ENDDO
1607  ENDDO
1608 
1609 ELSEIF (ksw == 6) THEN
1610  zpdh2o = zpdh6is
1611  zpdumg = zpdu6is
1612  zprh2o = zprh6is
1613  zprumg = zpru6is
1614  rtdh2o = ztdh6is
1615  rtdumg = ztdu6is
1616  rth2o = zth6is
1617  rtumg = ztu6is
1618 
1619  radjust=zadjus6
1620  rhsrta=zrto1
1621  rhsrtb=zrto2
1622  DO jw=1,ksw
1623  rsun(jw)=zsun6(jw)
1624 
1625  rswce(jw)=zswce6(jw)
1626  rswcp(jw)=zswcp6(jw)
1627 
1628  rroma(jw)=zroma6(jw)
1629  rromb(jw)=zromb6(jw)
1630  rrasy(jw)=zrasy6(jw)
1631  rhsra(jw)=zra6(jw)
1632  rhsrb(jw)=zrb6(jw)
1633  rhsrc(jw)=zrc6(jw)
1634  rhsrd(jw)=zrd6(jw)
1635  rhsre(jw)=zre6(jw)
1636  rhsrf(jw)=zrf6(jw)
1637 
1638  DO jc3=1,3
1639  d(jw,jc3)=zd6(jw,jc3)
1640  ENDDO
1641  IF (lrayl) THEN
1642  DO jc6=1,6
1643  rray(jw,jc6)=zray6(jw,jc6)
1644  ENDDO
1645  ELSE
1646  DO jc6=1,6
1647  rray(jw,jc6)=zrby6(jw,jc6)
1648  ENDDO
1649  ENDIF
1650  DO ji=1,3
1651  DO jj=1,7
1652  apad(jw,ji,jj)=zapad6(jw,ji,jj)
1653  bpad(jw,ji,jj)=zbpad6(jw,ji,jj)
1654  ENDDO
1655  ENDDO
1656  ENDDO
1657 
1658 ELSEIF (ktsw == 14 .AND. ksw == 14) THEN
1659  print *,'SUSWN: 14-SPECTRAL INTERVALS --> RRTM_SW'
1660 
1661 ELSE
1662  print *,'SUSWN: NTSW, NSW: ',ktsw,ksw
1663  CALL abor1('SUSWN: WRONG NUMBER OF SW INTERVALS')
1664 ENDIF
1665 
1666 ! ----------------------------------------------------------------
1667 
1668 !* 5. PRINT VALUES.
1669 ! ------------
1670 
1671 ! PRINT 9700,RPDH1,RPDU1,RPNU,RPNH
1672 !9700 FORMAT(1X,'RPDH1,RPDU1,RPNU,RPNH ',4F10.5)
1673 ! DO JW=1,KSW
1674 ! PRINT *,' SW INTERVAL ',JW
1675 ! PRINT 9701,JW,RSUN(JW),RSWCE(JW),RSWCP(JW)
1676 !9701 FORMAT(1X,I3,F10.7,2F10.5)
1677 ! PRINT 9702,JW,(D(JW,JC3),JC3=1,3)
1678 !9702 FORMAT(1X,I3,3F10.5)
1679 ! PRINT 9703,JW,(RRAY(JW,JC6),JC6=1,6)
1680 !9703 FORMAT(1X,I3,6E15.8)
1681 ! PRINT 9704,JW,(JABS,(APAD(JW,JABS,JJ),JJ=1,7),JABS=1,3)
1682 !9704 FORMAT(1X,I3,(I3,7E15.8))
1683 ! PRINT 9706,JW,RROMA(JW),RROMB(JW),RRASY(JW),RHSRTA,RHSRTB
1684 !9706 FORMAT(1X,I3,' Rain Optics from Rain Rate: ',F8.5,F6.3,3F6.3)
1685 ! PRINT 9707,JW,RHSRA(JW),RHSRB(JW),RHSRC(JW),RHSRD(JW)
1686 ! & ,RHSRE(JW),RHSRF(JW)
1687 !9707 FORMAT(1X,I3,' Rain Optics as function of Reff: ',6E11.4)
1688 ! ENDDO
1689 
1690 ! ----------------------------------------------------------------
1691 
1692 IF (lhook) CALL dr_hook('SUSWN',1,zhook_handle)
1693 END SUBROUTINE suswn
yoesw::rhsre
real(kind=jprb), dimension(6) rhsre
Definition: yoesw.F90:110
yoesw
Definition: yoesw.F90:1
yoesw::d
real(kind=jprb), dimension(6, 3) d
Definition: yoesw.F90:28
yoesw::rexpo3
real(kind=jprb), dimension(6, 2, 7) rexpo3
Definition: yoesw.F90:29
yoesw::rth2o
real(kind=jprb) rth2o
Definition: yoesw.F90:26
yoerad::ntsw
integer(kind=jpim) ntsw
Definition: yoerad.F90:30
yoesw::nexpo3
integer(kind=jpim), dimension(6) nexpo3
Definition: yoesw.F90:30
yoesw::rpnu
real(kind=jprb) rpnu
Definition: yoesw.F90:21
yoesw::rhsrb
real(kind=jprb), dimension(6) rhsrb
Definition: yoesw.F90:107
yoesw::rsun
real(kind=jprb), dimension(:), allocatable rsun
Definition: yoesw.F90:16
yoesw::bpad
real(kind=jprb), dimension(6, 3, 7) bpad
Definition: yoesw.F90:14
yoesw::rtdumg
real(kind=jprb) rtdumg
Definition: yoesw.F90:25
abor1
subroutine abor1(CDTEXT)
Definition: abor1.F90:2
yoesw::rrasy
real(kind=jprb), dimension(6) rrasy
Definition: yoesw.F90:104
yomcst::rg
real(kind=jprb) rg
Definition: yomcst.F90:29
yoesw::rhsrc
real(kind=jprb), dimension(6) rhsrc
Definition: yoesw.F90:108
yoerad::lonewsw
logical lonewsw
Definition: yoerad.F90:50
yoesw::rromb
real(kind=jprb), dimension(6) rromb
Definition: yoesw.F90:103
yoesw::rpdh1
real(kind=jprb) rpdh1
Definition: yoesw.F90:18
yoesw::rhsrd
real(kind=jprb), dimension(6) rhsrd
Definition: yoesw.F90:109
yoesw::rtumg
real(kind=jprb) rtumg
Definition: yoesw.F90:27
yoesw::apad
real(kind=jprb), dimension(6, 3, 7) apad
Definition: yoesw.F90:13
yoesw::rswce
real(kind=jprb), dimension(6) rswce
Definition: yoesw.F90:22
parkind1::jprb
integer, parameter jprb
Definition: parkind1.F90:31
yoerad
Definition: yoerad.F90:1
yoesw::rhsra
real(kind=jprb), dimension(6) rhsra
Definition: yoesw.F90:106
yomhook::lhook
logical lhook
Definition: yomhook.F90:12
yoesw::rtdh2o
real(kind=jprb) rtdh2o
Definition: yoesw.F90:24
suswn
subroutine suswn(KTSW, KSW)
Definition: suswn.F90:2
yoesw::rroma
real(kind=jprb), dimension(6) rroma
Definition: yoesw.F90:102
yoesw::rhsrta
real(kind=jprb) rhsrta
Definition: yoesw.F90:112
yomhook::dr_hook
subroutine dr_hook(CDNAME, KSWITCH, PKEY)
Definition: yomhook.F90:17
yoesw::rpnh
real(kind=jprb) rpnh
Definition: yoesw.F90:20
yoesw::radjust
real(kind=jprb) radjust
Definition: yoesw.F90:127
parkind1::jpim
integer, parameter jpim
Definition: parkind1.F90:13
yoesw::rhsrf
real(kind=jprb), dimension(6) rhsrf
Definition: yoesw.F90:111
yoesw::rray
real(kind=jprb), dimension(6, 6) rray
Definition: yoesw.F90:15
yomcst
Definition: yomcst.F90:1
yoesw::rhsrtb
real(kind=jprb) rhsrtb
Definition: yoesw.F90:113
yoesw::rpdu1
real(kind=jprb) rpdu1
Definition: yoesw.F90:19
yoerad::lrayl
logical lrayl
Definition: yoerad.F90:58
yoesw::rswcp
real(kind=jprb), dimension(6) rswcp
Definition: yoesw.F90:23