Context Navigation

suswn.F90 @ 5408

Last change on this file since 5408 was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM

Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence

Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

File size: 43.6 KB

Line
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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format