Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

suclopn.F90 @ 5445

Last change on this file since 5445 was 2160, checked in by Laurent Fairhead, 10 years ago
Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand
File size: 31.7 KB

Line
1	SUBROUTINE SUCLOPN (KTSW,KSW , KLEV)
2
3	!**** SUCLOP - INITIALIZE COMMON YOECLOP
4
5	! PURPOSE.
6	! --------
7	! INITIALIZE YOMCLOP, WITH CLOUD OPTICAL PARAMETERS
8
9	!** INTERFACE.
10	! ----------
11	! CALL SUCLOPN
12	! FROM SUECRAD
13
14	! EXPLICIT ARGUMENTS :
15	! --------------------
16	! NONE
17
18	! IMPLICIT ARGUMENTS :
19	! --------------------
20	! COMMON YOECLOP
21
22	! METHOD.
23	! -------
24	! SEE DOCUMENTATION
25
26	! EXTERNALS.
27	! ----------
28	! NONE
29
30	! REFERENCE.
31	! ----------
32	! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
33	! "INTEGRATED FORECASTING SYSTEM"
34
35	! Fouquart, 1987: NATO ASI, 223-284
36	! A. Slingo, 1989: J. Atmos. Sci., 46, 1419-1427
37	! Ebert and Curry, 1992: J. Geophys. Res., 97D, 3831-3836
38	! Sun and Shine, 1994: Quart. J. Roy. Meteor. Soc., 120, 111-138
39	! Fu and Liou, 1993: J. Atmos. Sci., 50, 2008-2025
40	! Fu, 1996: J. Climate, 9, 2058-2082
41	! Fu et al., 1998: J. Climate, 11, 2223-2237
42
43	! AUTHOR.
44	! -------
45	! JEAN-JACQUES MORCRETTE ECMWF
46
47	! MODIFICATIONS.
48	! --------------
49	! ORIGINAL : 92-02-29
50	! 97-04-16 JJ Morcrette 2 and 4 spectral intervals
51	! 00-10-25 JJMorcrette 6 spectral intervals
52	! 01-01-16 JJ Morcrette ice cloud from Fu 96, Fu et al 98
53
54	! ------------------------------------------------------------------
55
56	#include "tsmbkind.h"
57
58	USE YOESW , ONLY : RYFWCA ,RYFWCB ,RYFWCC ,RYFWCD ,&
59	&RYFWCE ,RYFWCF ,REBCUA ,REBCUB ,REBCUC ,&
60	&REBCUD ,REBCUE ,REBCUF ,REBCUG ,REBCUH ,&
61	&REBCUI ,REBCUJ ,RASWCA ,RASWCB ,RASWCC ,&
62	&RASWCD ,RASWCE ,RASWCF ,RSUSHE ,RSUSHF ,&
63	&RSUSHH ,RSUSHK ,RSUSHA ,RSUSHG ,RSUSHFA ,&
64	&RSUSHC ,RSUSHD ,REFFIA ,REFFIB ,RHSAVI ,&
65	&RFUAA0 ,RFUAA1 ,RFUBB0 ,RFUBB1 ,RFUBB2 ,&
66	&RFUBB3 ,RFUCC0 ,RFUCC1 ,RFUCC2 ,RFUCC3 ,&
67	&RFUETA ,RFULIO ,RFLAA0 ,RFLAA1 ,RFLBB0 ,&
68	&RFLBB1 ,RFLBB2 ,RFLBB3 ,RFLCC0 ,RFLCC1 ,&
69	&RFLCC2 ,RFLCC3 ,RFLDD0 ,RFLDD1 ,RFLDD2 ,&
70	&RFLDD3 ,RLINLI ,RTIW ,RRIW
71	USE YOERAD , ONLY : RAOVLP ,RBOVLP
72
73	IMPLICIT NONE
74
75
76	! DUMMY INTEGER SCALARS
77	INTEGER_M :: KSW
78	INTEGER_M :: KTSW
79	INTEGER_M :: KLEV
80	INTEGER_M :: JC, JNU
81
82
83
84
85	! -----------------------------------------------------------------
86	REAL_B :: ZEBCUA2(2) ,ZEBCUB2(2) ,ZEBCUC2(2) ,ZEBCUD2(2)&
87	&, ZEBCUE2(2) ,ZEBCUF2(2) ,ZYFWCA2(2) ,ZYFWCB2(2)&
88	&, ZYFWCC2(2) ,ZYFWCD2(2) ,ZYFWCE2(2) ,ZYFWCF2(2)&
89	&, ZASWCA2(2) ,ZASWCB2(2) ,ZASWCC2(2) ,ZASWCD2(2)&
90	&, ZASWCE2(2) ,ZASWCF2(2) ,ZSUSHE2(2) ,ZSUSHF2(2)&
91	&, ZSUSHH2(2) ,ZSUSHK2(2) ,ZSUSHA2(2) ,ZSUSHG2(2)&
92	&, ZFLAA02(2) ,ZFLAA12(2) ,ZFLBB02(2) ,ZFLBB12(2)&
93	&, ZFLBB22(2) ,ZFLBB32(2) ,ZFLCC02(2) ,ZFLCC12(2)&
94	&, ZFLCC22(2) ,ZFLCC32(2) ,ZFLDD02(2) ,ZFLDD12(2)&
95	&, ZFLDD22(2) ,ZFLDD32(2)&
96	&, ZFUAA02(2) ,ZFUAA12(2) ,ZFUBB02(2) ,ZFUBB12(2)&
97	&, ZFUBB22(2) ,ZFUBB32(2) ,ZFUCC02(2) ,ZFUCC12(2)&
98	&, ZFUCC22(2) ,ZFUCC32(2) ,ZFUDD02(2) ,ZFUDD12(2)&
99	&, ZFUDD22(2) ,ZFUDD32(2)
100
101	REAL_B :: ZEBCUA4(4) ,ZEBCUB4(4) ,ZEBCUC4(4) ,ZEBCUD4(4)&
102	&, ZEBCUE4(4) ,ZEBCUF4(4) ,ZYFWCA4(4) ,ZYFWCB4(4)&
103	&, ZYFWCC4(4) ,ZYFWCD4(4) ,ZYFWCE4(4) ,ZYFWCF4(4)&
104	&, ZASWCA4(4) ,ZASWCB4(4) ,ZASWCC4(4) ,ZASWCD4(4)&
105	&, ZASWCE4(4) ,ZASWCF4(4) ,ZSUSHE4(4) ,ZSUSHF4(4)&
106	&, ZSUSHH4(4) ,ZSUSHK4(4) ,ZSUSHA4(4) ,ZSUSHG4(4)&
107	&, ZFLAA04(4) ,ZFLAA14(4) ,ZFLBB04(4) ,ZFLBB14(4)&
108	&, ZFLBB24(4) ,ZFLBB34(4) ,ZFLCC04(4) ,ZFLCC14(4)&
109	&, ZFLCC24(4) ,ZFLCC34(4) ,ZFLDD04(4) ,ZFLDD14(4)&
110	&, ZFLDD24(4) ,ZFLDD34(4)&
111	&, ZFUAA04(4) ,ZFUAA14(4) ,ZFUBB04(4) ,ZFUBB14(4)&
112	&, ZFUBB24(4) ,ZFUBB34(4) ,ZFUCC04(4) ,ZFUCC14(4)&
113	&, ZFUCC24(4) ,ZFUCC34(4) ,ZFUDD04(4) ,ZFUDD14(4)&
114	&, ZFUDD24(4) ,ZFUDD34(4)
115
116	REAL_B :: ZEBCUA6(6) ,ZEBCUB6(6) ,ZEBCUC6(6) ,ZEBCUD6(6)&
117	&, ZEBCUE6(6) ,ZEBCUF6(6) ,ZYFWCA6(6) ,ZYFWCB6(6)&
118	&, ZYFWCC6(6) ,ZYFWCD6(6) ,ZYFWCE6(6) ,ZYFWCF6(6)&
119	&, ZASWCA6(6) ,ZASWCB6(6) ,ZASWCC6(6) ,ZASWCD6(6)&
120	&, ZASWCE6(6) ,ZASWCF6(6) ,ZSUSHE6(6) ,ZSUSHF6(6)&
121	&, ZSUSHH6(6) ,ZSUSHK6(6) ,ZSUSHA6(6) ,ZSUSHG6(6)&
122	&, ZFLAA06(6) ,ZFLAA16(6) ,ZFLBB06(6) ,ZFLBB16(6)&
123	&, ZFLBB26(6) ,ZFLBB36(6) ,ZFLCC06(6) ,ZFLCC16(6)&
124	&, ZFLCC26(6) ,ZFLCC36(6) ,ZFLDD06(6) ,ZFLDD16(6)&
125	&, ZFLDD26(6) ,ZFLDD36(6)&
126	&, ZFUAA06(6) ,ZFUAA16(6) ,ZFUBB06(6) ,ZFUBB16(6)&
127	&, ZFUBB26(6) ,ZFUBB36(6) ,ZFUCC06(6) ,ZFUCC16(6)&
128	&, ZFUCC26(6) ,ZFUCC36(6) ,ZFUDD06(6) ,ZFUDD16(6)&
129	&, ZFUDD26(6) ,ZFUDD36(6)
130
131
132	REAL_B :: ZAOVLP(3), ZBOVLP(3)
133
134	! -----------------------------------------------------------------
135
136	!* 1. SHORTWAVE CLOUD OPTICAL PROPERTIES
137	! ----------------------------------
138
139	! ------------------------------------------------------------------
140
141	!* 1.1 TWO SPECTRAL INTERVALS
142	! ----------------------
143
144	! SW : 0.25 - 0.68 - 4.00 microns
145
146	!* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
147
148	! optical properties
149	ZEBCUA2 = (/ 3.448E-04_JPRB , 3.448E-04_JPRB /)
150	ZEBCUB2 = (/ 2.431_JPRB , 2.431_JPRB /)
151	ZEBCUC2 = (/ 0.00001_JPRB , 0.024366_JPRB /)
152	ZEBCUD2 = (/ _ZERO_ , 2.487E-04_JPRB /)
153	ZEBCUE2 = (/ 0.7661_JPRB , 0.7866_JPRB /)
154	ZEBCUF2 = (/ 5.851E-04_JPRB , 5.937E-04_JPRB /)
155
156	! optical properties
157	! ZEBCUA2 = (/ 3.448E-04 , 3.448E-04 /)
158	! ZEBCUB2 = (/ 2.431 , 2.431 /)
159	! ZEBCUC2 = (/ 0.00001 , 0.035589 /)
160	! ZEBCUD2 = (/ 0. , 2.757E-04 /)
161	! ZEBCUE2 = (/ 0.7661 , 0.7921 /)
162	! ZEBCUF2 = (/ 5.851E-04 , 5.893E-04 /)
163
164	!* Water cloud properties - from Fouquart (1987)
165
166	ZYFWCA2 = (/ _ZERO_ , _ZERO_ /)
167	ZYFWCB2 = (/ 1.5_JPRB , 1.5_JPRB /)
168	ZYFWCC2 = (/ 0.9999_JPRB , 0.9988_JPRB /)
169	ZYFWCD2 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB /)
170	ZYFWCE2 = (/ _HALF_ , 0.05_JPRB /)
171	ZYFWCF2 = (/ 0.865_JPRB , 0.910_JPRB /)
172
173	!* Water cloud properties - from Slingo (1989)
174
175	ZASWCA2 = (/ 2.817_JPRB , 2.4584_JPRB /)
176	ZASWCB2 = (/ 1.305_JPRB , 1.3994_JPRB /)
177	ZASWCC2 = (/-5.62E-08_JPRB , 1.111E-02_JPRB /)
178	ZASWCD2 = (/ 1.63E-07_JPRB , 8.613E-04_JPRB /)
179	ZASWCE2 = (/ 0.829_JPRB , 0.7819_JPRB /)
180	ZASWCF2 = (/ 2.482_JPRB , 5.0150_JPRB /)
181
182	!* Ice cloud properties - from Sun and Shine (1995)
183
184	ZSUSHE2 = (/ _ZERO_ , 8.6822_JPRB /)
185	ZSUSHF2 = (/ _ZERO_ , 9.6277_JPRB /)
186	ZSUSHH2 = (/ 0.8522_JPRB , 0.8819_JPRB /)
187	ZSUSHK2 = (/ 0.1620_JPRB , 0.1630_JPRB /)
188	ZSUSHA2 = (/ _ZERO_ , 23.204_JPRB /)
189	ZSUSHG2 = (/ 0.3270_JPRB , 0.4180_JPRB /)
190
191	!* Ice cloud properties - from Fu and Liou (1993)
192
193	ZFLAA02 = (/-6.656E-3_JPRB ,-6.656E-3_JPRB /)
194	ZFLAA12 = (/ 3.686_JPRB , 3.686_JPRB /)
195	ZFLBB02 = (/ .10998E-5_JPRB , .21136E-1_JPRB /)
196	ZFLBB12 = (/-.26101E-7_JPRB , .39150E-3_JPRB /)
197	ZFLBB22 = (/ .18096E-8_JPRB ,-.20740E-6_JPRB /)
198	ZFLBB32 = (/-.47387E-11_JPRB,-.28829E-8_JPRB /)
199	ZFLCC02 = (/ .22110E+1_JPRB , .22498E+1_JPRB /)
200	ZFLCC12 = (/-.10398E-2_JPRB , .23656E-3_JPRB /)
201	ZFLCC22 = (/ .65199E-4_JPRB , .51948E-4_JPRB /)
202	ZFLCC32 = (/-.34498E-6_JPRB ,-.29768E-6_JPRB /)
203	ZFLDD02 = (/ .12495_JPRB , .11716_JPRB /)
204	ZFLDD12 = (/-.43582E-3_JPRB ,-.45208E-3_JPRB /)
205	ZFLDD22 = (/ .14092E-4_JPRB , .12772E-4_JPRB /)
206	ZFLDD32 = (/-.69565E-7_JPRB ,-.62779E-7_JPRB /)
207
208	!* Ice cloud properties - from Fu (1996)
209
210	!ZFUAA02 = (/-.291721E-04_JPRB ,
211	!ZFUAA12 = (/ .251925E+01_JPRB ,
212	!ZFUBB02 = (/ .135403E-06_JPRB ,
213	!ZFUBB12 = (/ .992822E-07_JPRB ,
214	!ZFUBB22 = (/-.738432E-10_JPRB ,
215	!ZFUBB32 = (/ .331119E-12_JPRB ,
216	!ZFUCC02 = (/ .748127E+00_JPRB ,
217	!ZFUCC12 = (/ .956845E-03_JPRB ,
218	!ZFUCC22 = (/-.111517E-05_JPRB ,
219	!ZFUCC32 = (/-.815573E-08_JPRB ,
220	!ZFUDD02 = (/ .115730E+00_JPRB ,
221	!ZFUDD12 = (/ .256481E-03_JPRB ,
222	!ZFUDD22 = (/ .191313E-05_JPRB ,
223	!ZFUDD32 = (/-.124603E-07_JPRB ,
224	! ------------------------------------------------------------------
225
226	!* 1.2 FOUR SPECTRAL INTERVALS
227	! -----------------------
228
229	! SW : 4 spectral intervals (0.25 - 0.69 - 1.19 - 2.38 - 4.00)
230
231	!* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
232
233	! ZEBCUA4 = (/ 3.448E-03 , 3.448E-03 , 3.448E-03 , 3.448E-03 /)
234	! ZEBCUB4 = (/ 2.431 , 2.431 , 2.431 , 2.431 /)
235	! ZEBCUC4 = (/ 0.00001 , 0.00011 , 0.01861 , 0.46658 /)
236	! ZEBCUD4 = (/ 0. , 1.405E-05 , 8.328E-4 , 2.050E-05 /)
237	! ZEBCUE4 = (/ 0.7661 , 0.7730 , 0.7940 , 0.9595 /)
238	! ZEBCUF4 = (/ 5.851E-04 , 5.665E-04 , 7.267E-04 , 1.076E-04 /)
239
240	ZEBCUA4 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /)
241	ZEBCUB4 = (/ 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB /)
242	ZEBCUC4 = (/ 0.00001_JPRB , 0.00011_JPRB , 0.0197796_JPRB , 0.46658_JPRB /)
243	ZEBCUD4 = (/ _ZERO_ , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /)
244	ZEBCUE4 = (/ 0.7661_JPRB , 0.7730_JPRB , 0.795653_JPRB , 0.9595_JPRB /)
245	ZEBCUF4 = (/ 5.851E-04_JPRB , 5.665E-04_JPRB , 7.267E-04_JPRB , 1.076E-04_JPRB /)
246
247	!* Water cloud properties - from Fouquart (1987)
248
249	ZYFWCA4 = (/ _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ /)
250	ZYFWCB4 = (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /)
251	ZYFWCC4 = (/ 0.9999_JPRB , 0.9988_JPRB , 0.9988_JPRB , 0.9988_JPRB /)
252	ZYFWCD4 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB /)
253	ZYFWCE4 = (/ _HALF_ , 0.05_JPRB , 0.05_JPRB , 0.05_JPRB /)
254	ZYFWCF4 = (/ 0.865_JPRB , 0.910_JPRB , 0.910_JPRB , 0.910_JPRB /)
255
256	!* Water cloud properties - from Slingo (1989)
257
258	ZASWCA4 = (/ 2.817_JPRB , 2.682_JPRB , 2.264_JPRB , 1.281_JPRB /)
259	ZASWCB4 = (/ 1.305_JPRB , 1.346_JPRB , 1.455_JPRB , 1.641_JPRB /)
260	ZASWCC4 = (/-5.62E-08_JPRB ,-6.94E-06_JPRB , 4.64E-04_JPRB , 2.01E-01_JPRB /)
261	ZASWCD4 = (/ 1.63E-07_JPRB , 2.35E-05_JPRB , 1.24E-03_JPRB , 7.56E-03_JPRB /)
262	ZASWCE4 = (/ 0.829_JPRB , 0.794_JPRB , 0.754_JPRB , 0.826_JPRB /)
263	ZASWCF4 = (/ 2.482_JPRB , 4.226_JPRB , 6.560_JPRB , 4.353_JPRB /)
264
265	!* Ice cloud properties - from Sun and Shine (1995)
266
267	ZSUSHE4 = (/ _ZERO_ , 7.2471E-02_JPRB , 17.5933_JPRB , 48.7166_JPRB /)
268	ZSUSHF4 = (/ _ZERO_ , 4.01511_JPRB , 21.1249_JPRB , 1.20890_JPRB /)
269	ZSUSHH4 = (/ 0.8522_JPRB , 0.85841_JPRB , 0.90778_JPRB ,0.982046_JPRB /)
270	ZSUSHK4 = (/ 0.1620_JPRB , 0.160048_JPRB , 0.188521_JPRB ,0.0411446_JPRB /)
271	ZSUSHA4 = (/ _ZERO_ , 0.273455_JPRB , 41.7675_JPRB , 161.104_JPRB /)
272	ZSUSHG4 = (/ 0.3270_JPRB , 0.343668_JPRB , 0.526192_JPRB , 0.574040_JPRB /)
273
274	!* Ice cloud properties - from Fu and Liou (1993)
275
276	ZFLAA04 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /)
277	ZFLAA14 = (/ 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB /)
278	ZFLBB04 = (/ .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /)
279	ZFLBB14 = (/-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /)
280	ZFLBB24 = (/ .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /)
281	ZFLBB34 = (/-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /)
282	ZFLCC04 = (/ .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /)
283	ZFLCC14 = (/-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /)
284	ZFLCC24 = (/ .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /)
285	ZFLCC34 = (/-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /)
286	ZFLDD04 = (/ .12495_JPRB , .12363_JPRB , .11983_JPRB , .21834E-1_JPRB /)
287	ZFLDD14 = (/-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /)
288	ZFLDD24 = (/ .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /)
289	ZFLDD34 = (/-.69565E-7_JPRB ,-.68851E-7_JPRB,-.59367E-7_JPRB,-.46257E-8_JPRB /)
290
291	!* Ice cloud properties - from Fu (1996) as tabulated in Sun & Rikus (1999)
292
293	ZFUAA04 = (/-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /)
294	ZFUAA14 = (/ 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /)
295	ZFUBB04 = (/-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /)
296	ZFUBB14 = (/ 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /)
297	ZFUBB24 = (/-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /)
298	ZFUBB34 = (/ 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /)
299	ZFUCC04 = (/ 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /)
300	ZFUCC14 = (/ 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /)
301	ZFUCC24 = (/-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /)
302	ZFUCC34 = (/-1.05910E-08_JPRB,-3.03712E-09_JPRB, 1.97757E-08_JPRB, 8.23543E-08_JPRB /)
303
304	! ------------------------------------------------------------------
305
306	!* 1.3 SIX SPECTRAL INTERVALS
307	! ----------------------
308
309	! SW : 6 spectral intervals (0.185-0.25-0.44-0.69-1.19-2.38-4.00)
310
311	!* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
312
313	ZEBCUA6 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /)
314	ZEBCUB6 = (/ 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB , 2.431_JPRB /)
315	ZEBCUC6 = (/ 0.00001_JPRB , 0.00001_JPRB , 0.00001_JPRB , 0.00011_JPRB , 0.0197796_JPRB , 0.46658_JPRB /)
316	ZEBCUD6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /)
317	ZEBCUE6 = (/ 0.7661_JPRB , 0.7661_JPRB , 0.7661_JPRB , 0.7730_JPRB , 0.795653_JPRB , 0.9595_JPRB /)
318	ZEBCUF6 = (/ 5.851E-04_JPRB , 5.851E-04_JPRB , 5.851E-04_JPRB , 5.665E-04_JPRB , 7.267E-04_JPRB , 1.076E-04_JPRB /)
319
320	!* Water cloud properties - from Fouquart (1987)
321
322	ZYFWCA6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ , _ZERO_ /)
323	ZYFWCB6 = (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /)
324	ZYFWCC6 = (/ 0.9999_JPRB , 0.9999_JPRB , 0.9999_JPRB , 0.9988_JPRB , 0.9988_JPRB , 0.9988_JPRB /)
325	ZYFWCD6 = (/ 5.000E-04_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB /)
326	ZYFWCE6 = (/ _HALF_ , _HALF_ , _HALF_ , 0.05_JPRB , 0.05_JPRB , 0.05_JPRB /)
327	ZYFWCF6 = (/ 0.865_JPRB , 0.865_JPRB , 0.865_JPRB , 0.910_JPRB , 0.910_JPRB , 0.910_JPRB /)
328
329	!* Water cloud properties - from Slingo (1989)
330
331	ZASWCA6 = (/ 2.817_JPRB , 2.817_JPRB , 2.817_JPRB , 2.682_JPRB , 2.264_JPRB , 1.281_JPRB /)
332	ZASWCB6 = (/ 1.305_JPRB , 1.305_JPRB , 1.305_JPRB , 1.346_JPRB , 1.455_JPRB , 1.641_JPRB /)
333	ZASWCC6 = (/-5.62E-08_JPRB ,-5.62E-08_JPRB ,-5.62E-08_JPRB ,-6.94E-06_JPRB , 4.64E-04_JPRB , 2.01E-01_JPRB /)
334	ZASWCD6 = (/ 1.63E-07_JPRB , 1.63E-07_JPRB , 1.63E-07_JPRB , 2.35E-05_JPRB , 1.24E-03_JPRB , 7.56E-03_JPRB /)
335	ZASWCE6 = (/ 0.829_JPRB , 0.829_JPRB , 0.829_JPRB , 0.794_JPRB , 0.754_JPRB , 0.826_JPRB /)
336	ZASWCF6 = (/ 2.482_JPRB , 2.482_JPRB , 2.482_JPRB , 4.226_JPRB , 6.560_JPRB , 4.353_JPRB /)
337
338	!* Ice cloud properties - from Sun and Shine (1995)
339
340	ZSUSHE6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 7.2471E-02_JPRB , 17.5933_JPRB , 48.7166_JPRB /)
341	ZSUSHF6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 4.01511_JPRB , 21.1249_JPRB , 1.20890_JPRB /)
342	ZSUSHH6 = (/ 0.8522_JPRB , 0.8522_JPRB , 0.8522_JPRB , 0.85841_JPRB , 0.90778_JPRB ,0.982046_JPRB /)
343	ZSUSHK6 = (/ 0.1620_JPRB , 0.1620_JPRB , 0.1620_JPRB , 0.160048_JPRB , 0.188521_JPRB ,0.0411446_JPRB /)
344	ZSUSHA6 = (/ _ZERO_ , _ZERO_ , _ZERO_ , 0.273455_JPRB , 41.7675_JPRB , 161.104_JPRB /)
345	ZSUSHG6 = (/ 0.3270_JPRB , 0.3270_JPRB , 0.3270_JPRB , 0.343668_JPRB , 0.526192_JPRB , 0.574040_JPRB /)
346
347	!* Ice cloud properties - from Fu and Liou (1993)
348
349	ZFLAA06 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /)
350	ZFLAA16 = (/ 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB , 3.686_JPRB /)
351	ZFLBB06 = (/ .10998E-5_JPRB , .10998E-5_JPRB , .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /)
352	ZFLBB16 = (/-.26101E-7_JPRB ,-.26101E-7_JPRB ,-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /)
353	ZFLBB26 = (/ .18096E-8_JPRB , .18096E-8_JPRB , .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /)
354	ZFLBB36 = (/-.47387E-11_JPRB,-.47387E-11_JPRB,-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /)
355	ZFLCC06 = (/ .22110E+1_JPRB , .22110E+1_JPRB , .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /)
356	ZFLCC16 = (/-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /)
357	ZFLCC26 = (/ .65199E-4_JPRB , .65199E-4_JPRB , .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /)
358	ZFLCC36 = (/-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /)
359	ZFLDD06 = (/ .12495_JPRB , .12495_JPRB , .12495_JPRB , .12363_JPRB , .11983_JPRB , .21834E-1_JPRB /)
360	ZFLDD16 = (/-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /)
361	ZFLDD26 = (/ .14092E-4_JPRB , .14092E-4_JPRB , .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /)
362	ZFLDD36 = (/-.69565E-7_JPRB ,-.69565E-7_JPRB ,-.69565E-7_JPRB ,-.68851E-7_JPRB,-.59367E-7_JPRB,-.46257E-8_JPRB /)
363
364	!* Ice cloud properties - from Fu (1996) as tabulated in Sun & Rikus (1999)
365
366	ZFUAA06 = (/-1.30817E-04_JPRB,-1.30817E-04_JPRB,-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /)
367	ZFUAA16 = (/ 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /)
368	ZFUBB06 = (/-1.55357E-08_JPRB,-1.55357E-08_JPRB,-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /)
369	ZFUBB16 = (/ 1.95793E-07_JPRB, 1.95793E-07_JPRB, 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /)
370	ZFUBB26 = (/-2.31234E-10_JPRB,-2.31234E-10_JPRB,-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /)
371	ZFUBB36 = (/ 1.12247E-12_JPRB, 1.12247E-12_JPRB, 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /)
372	ZFUCC06 = (/ 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /)
373	ZFUCC16 = (/ 9.10564E-04_JPRB, 9.10564E-04_JPRB, 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /)
374	ZFUCC26 = (/-4.62479E-07_JPRB,-4.62479E-07_JPRB,-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /)
375	ZFUCC36 = (/-1.05910E-08_JPRB,-1.05910E-08_JPRB,-1.05910E-08_JPRB,-3.03712E-09_JPRB, 1.97757E-08_JPRB, 8.23543E-08_JPRB /)
376
377	! ------------------------------------------------------------------
378
379	! LW : absorption coefficient as a function of effective radius in RRTM
380
381	! water clouds from Savijarvi
382
383	RHSAVI( 1, :) = (/&
384	&0.1651082_JPRB , -0.003494839_JPRB , _ZERO_ /)
385	RHSAVI( 2, :) = (/&
386	&0.327820597_JPRB, -0.0219634383_JPRB , 0.000506783898_JPRB /)
387	RHSAVI( 3, :) = (/&
388	&0.504805453_JPRB, -0.0478602354_JPRB , 0.00141521102_JPRB /)
389	RHSAVI( 4, :) = (/&
390	&0.513169093_JPRB, -0.0508960145_JPRB , 0.00155498711_JPRB /)
391	RHSAVI( 5, :) = (/&
392	&0.448042082_JPRB, -0.0431857592_JPRB , 0.00130848978_JPRB /)
393	RHSAVI( 6, :) = (/&
394	&0.249547237_JPRB, -0.0185273835_JPRB , 0.00050361258_JPRB /)
395	RHSAVI( 7, :) = (/&
396	&0.135486796_JPRB, -0.00585852322_JPRB, 0.00011873119_JPRB /)
397	RHSAVI( 8, :) = (/&
398	&0.126668819_JPRB, -0.00497949082_JPRB, 0.0000927679172_JPRB/)
399	RHSAVI( 9, :) = (/&
400	&0.130938752_JPRB, -0.00531814674_JPRB, 0.0000984953029_JPRB/)
401	RHSAVI(10, :) = (/&
402	&0.147024519_JPRB, -0.00671655774_JPRB, 0.000130885091_JPRB /)
403	RHSAVI(11, :) = (/&
404	&0.241463914_JPRB, -0.0174486461_JPRB , 0.000457756556_JPRB /)
405	RHSAVI(12, :) = (/&
406	&0.07576579_JPRB , -0.001695588_JPRB , _ZERO_ /)
407	RHSAVI(13, :) = (/&
408	&0.1032178_JPRB , -0.00293412_JPRB , _ZERO_ /)
409	RHSAVI(14, :) = (/&
410	&0.07342832_JPRB , -0.001775135_JPRB , _ZERO_ /)
411	RHSAVI(15, :) = (/&
412	&0.04649514_JPRB , -0.0009165462_JPRB , _ZERO_ /)
413	RHSAVI(16, :) = (/&
414	&0.06893519_JPRB , -0.0001245402_JPRB , _ZERO_ /)
415
416	! water clouds from Lindner & Li (2000)
417
418	RLINLI( 1, :) = (/&
419	& 0.88116E-01_JPRB,-0.12857E-02_JPRB, 0.81658E+00_JPRB,-0.39428E+01_JPRB, 0.46652E+01_JPRB /)
420	RLINLI( 2, :) = (/&
421	& 0.41307E-03_JPRB,-0.59631E-04_JPRB, 0.24275E+01_JPRB,-0.90838E+01_JPRB, 0.96069E+01_JPRB /)
422	RLINLI( 3, :) = (/&
423	&-0.57709E-01_JPRB, 0.99071E-03_JPRB, 0.31118E+01_JPRB,-0.95540E+01_JPRB, 0.90189E+01_JPRB /)
424	RLINLI( 4, :) = (/&
425	&-0.53069E-01_JPRB, 0.99992E-03_JPRB, 0.28045E+01_JPRB,-0.72836E+01_JPRB, 0.62573E+01_JPRB /)
426	RLINLI( 5, :) = (/&
427	&-0.23627E-01_JPRB, 0.55291E-03_JPRB, 0.21785E+01_JPRB,-0.54664E+01_JPRB, 0.47379E+01_JPRB /)
428	RLINLI( 6, :) = (/&
429	& 0.29022E-01_JPRB,-0.39657E-03_JPRB, 0.14902E+01_JPRB,-0.50777E+01_JPRB, 0.52170E+01_JPRB /)
430	RLINLI( 7, :) = (/&
431	&-0.24901E-01_JPRB, 0.16195E-03_JPRB, 0.29375E+02_JPRB,-0.11437E+02_JPRB, 0.12273E+02_JPRB /)
432	RLINLI( 8, :) = (/&
433	&-0.14269E+00_JPRB, 0.22282E-02_JPRB, 0.46478E+01_JPRB,-0.16369E+02_JPRB, 0.16533E+02_JPRB /)
434	RLINLI( 9, :) = (/&
435	&-0.20398E+00_JPRB, 0.34708E-02_JPRB, 0.52858E+01_JPRB,-0.16603E+02_JPRB, 0.15392E+02_JPRB /)
436	RLINLI(10, :) = (/&
437	&-0.18318E+00_JPRB, 0.33080E-02_JPRB, 0.46120E+01_JPRB,-0.11550E+02_JPRB, 0.87086E+01_JPRB /)
438	RLINLI(11, :) = (/&
439	&-0.20420E+00_JPRB, 0.37167E-02_JPRB, 0.48566E+01_JPRB,-0.11972E+02_JPRB, 0.86344E+01_JPRB /)
440	RLINLI(12, :) = (/&
441	&-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
442	RLINLI(13, :) = (/&
443	&-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
444	RLINLI(14, :) = (/&
445	&-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
446	RLINLI(15, :) = (/&
447	&-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
448	RLINLI(16, :) = (/&
449	&-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
450
451	! ice clouds From Fu & Liou
452
453	RFULIO( 1, :) = (/&
454	& -7.752E-03_JPRB, 4.624_JPRB, -42.01_JPRB /)
455	RFULIO( 2, :) = (/&
456	& -1.741E-02_JPRB, 5.541_JPRB, -58.42_JPRB /)
457	RFULIO( 3, :) = (/&
458	& -1.704E-02_JPRB, 4.830_JPRB, 16.27_JPRB /)
459	RFULIO( 4, :) = (/&
460	& -1.151E-02_JPRB, 4.182_JPRB, 31.13_JPRB /)
461	RFULIO( 5, :) = (/&
462	& -1.026E-02_JPRB, 4.105_JPRB, 16.36_JPRB /)
463	RFULIO( 6, :) = (/&
464	& -8.294E-03_JPRB, 3.925_JPRB, 1.315_JPRB /)
465	RFULIO( 7, :) = (/&
466	& -1.153E-02_JPRB, 4.109_JPRB, 17.32_JPRB /)
467	RFULIO( 8, :) = (/&
468	& -9.609E-03_JPRB, 3.768_JPRB, 34.11_JPRB /)
469	RFULIO( 9, :) = (/&
470	& -9.061E-03_JPRB, 3.741_JPRB, 26.48_JPRB /)
471	RFULIO(10, :) = (/&
472	& -8.441E-03_JPRB, 3.715_JPRB, 19.48_JPRB /)
473	RFULIO(11, :) = (/&
474	& -8.088E-03_JPRB, 3.717_JPRB, 17.17_JPRB /)
475	RFULIO(12, :) = (/&
476	& -8.088E-03_JPRB, 3.717_JPRB, 17.17_JPRB /)
477	RFULIO(13, :) = (/&
478	& -7.770E-03_JPRB, 3.734_JPRB, 11.85_JPRB /)
479	RFULIO(14, :) = (/&
480	& -6.656E-03_JPRB, 3.686_JPRB, _ZERO_ /)
481	RFULIO(15, :) = (/&
482	& -6.656E-03_JPRB, 3.686_JPRB, _ZERO_ /)
483	RFULIO(16, :) = (/&
484	& -6.656E-03_JPRB, 3.686_JPRB, _ZERO_ /)
485
486	! ice clouds from Fu et al. (1998)
487
488	RFUETA( 1, :) = (/&
489	& 4.919685E-03_JPRB, 2.327741E+00_JPRB,-1.390858E+01_JPRB /)
490	RFUETA( 2, :) = (/&
491	& 3.325756E-03_JPRB, 2.601360E+00_JPRB,-1.909602E+01_JPRB /)
492	RFUETA( 3, :) = (/&
493	&-1.334860E-02_JPRB, 4.043808E+00_JPRB,-2.171029E+01_JPRB /)
494	RFUETA( 4, :) = (/&
495	&-9.524174E-03_JPRB, 3.587742E+00_JPRB,-1.068895E+01_JPRB /)
496	RFUETA( 5, :) = (/&
497	&-4.159424E-03_JPRB, 3.047325E+00_JPRB,-5.061568E+00_JPRB /)
498	RFUETA( 6, :) = (/&
499	&-1.691632E-03_JPRB, 2.765756E+00_JPRB,-8.331033E+00_JPRB /)
500	RFUETA( 7, :) = (/&
501	&-8.372696E-03_JPRB, 3.455018E+00_JPRB,-1.516692E+01_JPRB /)
502	RFUETA( 8, :) = (/&
503	&-8.178608E-03_JPRB, 3.401245E+00_JPRB,-8.812820E+00_JPRB /)
504	RFUETA( 9, :) = (/&
505	&-4.936610E-03_JPRB, 3.087764E+00_JPRB,-3.884262E+00_JPRB /)
506	RFUETA(10, :) = (/&
507	&-3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /)
508	RFUETA(11, :) = (/&
509	&-3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /)
510	RFUETA(12, :) = (/&
511	&-2.465236E-03_JPRB, 2.833187E+00_JPRB,-4.227573E-01_JPRB /)
512	RFUETA(13, :) = (/&
513	&-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
514	RFUETA(14, :) = (/&
515	&-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
516	RFUETA(15, :) = (/&
517	&-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
518	RFUETA(16, :) = (/&
519	&-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
520
521
522	! ----------------------------------------------------------------
523	! Ebert-Curry
524
525	! LW : spectrally defined for EC-OPE
526
527	REBCUI = (/ 1.136_JPRB, 1.338_JPRB, 1.166_JPRB, 1.166_JPRB, 1.118_JPRB, &
528	&0.600_JPRB /)
529	REBCUJ = (/ 0.0036_JPRB, 0.0003_JPRB, 0.0016_JPRB, 0.0016_JPRB, 0.0020_JPRB, &
530	&0.0068_JPRB /)
531
532	! LW : spectrally defined for RRTM
533	! mass-absorption coefficients for vertical path: no diffusivity factor
534
535	REBCUG = (/ 0.718_JPRB, 0.726_JPRB, 1.136_JPRB, 1.320_JPRB, 1.505_JPRB, &
536	& 1.290_JPRB, 0.911_JPRB, 0.949_JPRB, 1.021_JPRB, 1.193_JPRB, &
537	& 1.279_JPRB, 0.626_JPRB, 0.647_JPRB, 0.668_JPRB, 0.690_JPRB, &
538	& 0.690_JPRB /)
539
540
541	REBCUH = (/ 0.0069_JPRB, 0.0060_JPRB, 0.0024_JPRB, 0.0004_JPRB,-0.0016_JPRB, &
542	& 0.0003_JPRB, 0.0043_JPRB, 0.0038_JPRB, 0.0030_JPRB, 0.0013_JPRB, &
543	& 0.0005_JPRB, 0.0054_JPRB, 0.0052_JPRB, 0.0050_JPRB, 0.0048_JPRB, &
544	& 0.0048_JPRB /)
545
546
547	! Sun-Shine
548
549	RSUSHFA = (/ 1.047_JPRB, -0.913E-04_JPRB, 0.203E-03_JPRB, -0.106E-04_JPRB /)
550
551	! ------------------------------------------------------------------
552
553	!* 2.
554
555	!* Liquid/Solid water transition
556
557	RTIW= 263._JPRB
558	RRIW= 20._JPRB
559
560	! Ice particle Effective Radius as a function of LWC
561
562	REFFIA= 40._JPRB
563	REFFIB= 0._JPRB
564
565	! Sun-Shine
566
567	RSUSHC= 0.0306_JPRB
568	RSUSHD= 0.2548_JPRB
569
570	! ------------------------------------------------------------------
571
572	! SW : absorption coefficients
573
574	IF (KSW == 2) THEN
575	DO JNU=1,KSW
576	RASWCA(JNU)=ZASWCA2(JNU)*1.E-02_JPRB
577	RASWCB(JNU)=ZASWCB2(JNU)
578	RASWCC(JNU)=ZASWCC2(JNU)
579	RASWCD(JNU)=ZASWCD2(JNU)
580	RASWCE(JNU)=ZASWCE2(JNU)
581	RASWCF(JNU)=ZASWCF2(JNU)*1.E-03_JPRB
582
583	REBCUA(JNU)=ZEBCUA2(JNU)
584	REBCUB(JNU)=ZEBCUB2(JNU)
585	REBCUC(JNU)=ZEBCUC2(JNU)
586	REBCUD(JNU)=ZEBCUD2(JNU)
587	REBCUE(JNU)=ZEBCUE2(JNU)
588	REBCUF(JNU)=ZEBCUF2(JNU)
589
590	RYFWCA(JNU)=ZYFWCA2(JNU)
591	RYFWCB(JNU)=ZYFWCB2(JNU)
592	RYFWCC(JNU)=ZYFWCC2(JNU)
593	RYFWCD(JNU)=ZYFWCD2(JNU)
594	RYFWCE(JNU)=ZYFWCE2(JNU)
595	RYFWCF(JNU)=ZYFWCF2(JNU)
596
597	RSUSHE(JNU)=ZSUSHE2(JNU)*1.E-02_JPRB
598	RSUSHF(JNU)=ZSUSHF2(JNU)*1.E-02_JPRB
599	RSUSHH(JNU)=ZSUSHH2(JNU)
600	RSUSHK(JNU)=ZSUSHK2(JNU)*1.E-01_JPRB
601	RSUSHA(JNU)=ZSUSHA2(JNU)*1.E-03_JPRB
602	RSUSHG(JNU)=ZSUSHG2(JNU)*1.E-01_JPRB
603
604	RFLAA0(JNU)=ZFLAA02(JNU)
605	RFLAA1(JNU)=ZFLAA12(JNU)
606	RFLBB0(JNU)=ZFLBB02(JNU)
607	RFLBB1(JNU)=ZFLBB12(JNU)
608	RFLBB2(JNU)=ZFLBB22(JNU)
609	RFLBB3(JNU)=ZFLBB32(JNU)
610	RFLCC0(JNU)=ZFLCC02(JNU)
611	RFLCC1(JNU)=ZFLCC12(JNU)
612	RFLCC2(JNU)=ZFLCC22(JNU)
613	RFLCC3(JNU)=ZFLCC32(JNU)
614	RFLDD0(JNU)=ZFLDD02(JNU)
615	RFLDD1(JNU)=ZFLDD12(JNU)
616	RFLDD2(JNU)=ZFLDD22(JNU)
617	RFLDD3(JNU)=ZFLDD32(JNU)
618
619	! RFUAA0(JNU)=ZFUAA02(JNU)
620	! RFUAA1(JNU)=ZFUAA12(JNU)
621	! RFUBB0(JNU)=ZFUBB02(JNU)
622	! RFUBB1(JNU)=ZFUBB12(JNU)
623	! RFUBB2(JNU)=ZFUBB22(JNU)
624	! RFUBB3(JNU)=ZFUBB32(JNU)
625	! RFUCC0(JNU)=ZFUCC02(JNU)
626	! RFUCC1(JNU)=ZFUCC12(JNU)
627	! RFUCC2(JNU)=ZFUCC22(JNU)
628	! RFUCC3(JNU)=ZFUCC32(JNU)
629	! RFUDD0(JNU)=ZFUDD02(JNU)
630	! RFUDD1(JNU)=ZFUDD12(JNU)
631	! RFUDD2(JNU)=ZFUDD22(JNU)
632	! RFUDD3(JNU)=ZFUDD32(JNU)
633
634	ENDDO
635	ELSEIF (KSW == 4) THEN
636	DO JNU=1,KSW
637	RASWCA(JNU)=ZASWCA4(JNU)*1.E-02_JPRB
638	RASWCB(JNU)=ZASWCB4(JNU)
639	RASWCC(JNU)=ZASWCC4(JNU)
640	RASWCD(JNU)=ZASWCD4(JNU)
641	RASWCE(JNU)=ZASWCE4(JNU)
642	RASWCF(JNU)=ZASWCF4(JNU)*1.E-03_JPRB
643
644	REBCUA(JNU)=ZEBCUA4(JNU)
645	REBCUB(JNU)=ZEBCUB4(JNU)
646	REBCUC(JNU)=ZEBCUC4(JNU)
647	REBCUD(JNU)=ZEBCUD4(JNU)
648	REBCUE(JNU)=ZEBCUE4(JNU)
649	REBCUF(JNU)=ZEBCUF4(JNU)
650
651	RYFWCA(JNU)=ZYFWCA4(JNU)
652	RYFWCB(JNU)=ZYFWCB4(JNU)
653	RYFWCC(JNU)=ZYFWCC4(JNU)
654	RYFWCD(JNU)=ZYFWCD4(JNU)
655	RYFWCE(JNU)=ZYFWCE4(JNU)
656	RYFWCF(JNU)=ZYFWCF4(JNU)
657
658	RSUSHE(JNU)=ZSUSHE4(JNU)*1.E-02_JPRB
659	RSUSHF(JNU)=ZSUSHF4(JNU)*1.E-02_JPRB
660	RSUSHH(JNU)=ZSUSHH4(JNU)
661	RSUSHK(JNU)=ZSUSHK4(JNU)*1.E-01_JPRB
662	RSUSHA(JNU)=ZSUSHA4(JNU)*1.E-03_JPRB
663	RSUSHG(JNU)=ZSUSHG4(JNU)*1.E-01_JPRB
664
665	RFLAA0(JNU)=ZFLAA04(JNU)
666	RFLAA1(JNU)=ZFLAA14(JNU)
667	RFLBB0(JNU)=ZFLBB04(JNU)
668	RFLBB1(JNU)=ZFLBB14(JNU)
669	RFLBB2(JNU)=ZFLBB24(JNU)
670	RFLBB3(JNU)=ZFLBB34(JNU)
671	RFLCC0(JNU)=ZFLCC04(JNU)
672	RFLCC1(JNU)=ZFLCC14(JNU)
673	RFLCC2(JNU)=ZFLCC24(JNU)
674	RFLCC3(JNU)=ZFLCC34(JNU)
675	RFLDD0(JNU)=ZFLDD04(JNU)
676	RFLDD1(JNU)=ZFLDD14(JNU)
677	RFLDD2(JNU)=ZFLDD24(JNU)
678	RFLDD3(JNU)=ZFLDD34(JNU)
679
680	RFUAA0(JNU)=ZFUAA04(JNU)
681	RFUAA1(JNU)=ZFUAA14(JNU)
682	RFUBB0(JNU)=ZFUBB04(JNU)
683	RFUBB1(JNU)=ZFUBB14(JNU)
684	RFUBB2(JNU)=ZFUBB24(JNU)
685	RFUBB3(JNU)=ZFUBB34(JNU)
686	RFUCC0(JNU)=ZFUCC04(JNU)
687	RFUCC1(JNU)=ZFUCC14(JNU)
688	RFUCC2(JNU)=ZFUCC24(JNU)
689	RFUCC3(JNU)=ZFUCC34(JNU)
690	! RFUDD0(JNU)=ZFUDD04(JNU)
691	! RFUDD1(JNU)=ZFUDD14(JNU)
692	! RFUDD2(JNU)=ZFUDD24(JNU)
693	! RFUDD3(JNU)=ZFUDD34(JNU)
694
695	ENDDO
696	ELSEIF (KSW == 6) THEN
697	DO JNU=1,KSW
698	RASWCA(JNU)=ZASWCA6(JNU)*1.E-02_JPRB
699	RASWCB(JNU)=ZASWCB6(JNU)
700	RASWCC(JNU)=ZASWCC6(JNU)
701	RASWCD(JNU)=ZASWCD6(JNU)
702	RASWCE(JNU)=ZASWCE6(JNU)
703	RASWCF(JNU)=ZASWCF6(JNU)*1.E-03_JPRB
704
705	REBCUA(JNU)=ZEBCUA6(JNU)
706	REBCUB(JNU)=ZEBCUB6(JNU)
707	REBCUC(JNU)=ZEBCUC6(JNU)
708	REBCUD(JNU)=ZEBCUD6(JNU)
709	REBCUE(JNU)=ZEBCUE6(JNU)
710	REBCUF(JNU)=ZEBCUF6(JNU)
711
712	RYFWCA(JNU)=ZYFWCA6(JNU)
713	RYFWCB(JNU)=ZYFWCB6(JNU)
714	RYFWCC(JNU)=ZYFWCC6(JNU)
715	RYFWCD(JNU)=ZYFWCD6(JNU)
716	RYFWCE(JNU)=ZYFWCE6(JNU)
717	RYFWCF(JNU)=ZYFWCF6(JNU)
718
719	RSUSHE(JNU)=ZSUSHE6(JNU)*1.E-02_JPRB
720	RSUSHF(JNU)=ZSUSHF6(JNU)*1.E-02_JPRB
721	RSUSHH(JNU)=ZSUSHH6(JNU)
722	RSUSHK(JNU)=ZSUSHK6(JNU)*1.E-01_JPRB
723	RSUSHA(JNU)=ZSUSHA6(JNU)*1.E-03_JPRB
724	RSUSHG(JNU)=ZSUSHG6(JNU)*1.E-01_JPRB
725
726	RFLAA0(JNU)=ZFLAA06(JNU)
727	RFLAA1(JNU)=ZFLAA16(JNU)
728	RFLBB0(JNU)=ZFLBB06(JNU)
729	RFLBB1(JNU)=ZFLBB16(JNU)
730	RFLBB2(JNU)=ZFLBB26(JNU)
731	RFLBB3(JNU)=ZFLBB36(JNU)
732	RFLCC0(JNU)=ZFLCC06(JNU)
733	RFLCC1(JNU)=ZFLCC16(JNU)
734	RFLCC2(JNU)=ZFLCC26(JNU)
735	RFLCC3(JNU)=ZFLCC36(JNU)
736	RFLDD0(JNU)=ZFLDD06(JNU)
737	RFLDD1(JNU)=ZFLDD16(JNU)
738	RFLDD2(JNU)=ZFLDD26(JNU)
739	RFLDD3(JNU)=ZFLDD36(JNU)
740
741	RFUAA0(JNU)=ZFUAA06(JNU)
742	RFUAA1(JNU)=ZFUAA16(JNU)
743	RFUBB0(JNU)=ZFUBB06(JNU)
744	RFUBB1(JNU)=ZFUBB16(JNU)
745	RFUBB2(JNU)=ZFUBB26(JNU)
746	RFUBB3(JNU)=ZFUBB36(JNU)
747	RFUCC0(JNU)=ZFUCC06(JNU)
748	RFUCC1(JNU)=ZFUCC16(JNU)
749	RFUCC2(JNU)=ZFUCC26(JNU)
750	RFUCC3(JNU)=ZFUCC36(JNU)
751	! RFUDD0(JNU)=ZFUDD06(JNU)
752	! RFUDD1(JNU)=ZFUDD16(JNU)
753	! RFUDD2(JNU)=ZFUDD26(JNU)
754	! RFUDD3(JNU)=ZFUDD36(JNU)
755
756	ENDDO
757	ELSE
758	! CALL ABOR1('SUCLOPN: WRONG SW SPECTRAL RESOLUTION')
759	STOP 'SUCLOPN: WRONG SW SPECTRAL RESOLUTION'
760	ENDIF
761
762	! ------------------------------------------------------------------
763
764	!* 2. CLOUD OVERLAP PARAMETERS
765	! ------------------------
766
767	ZAOVLP = (/ -2.250E-05_JPRB,-7.316E-06_JPRB,-1.966E-05_JPRB /)
768	ZBOVLP = (/ 0.7865_JPRB , 0.8186_JPRB , 0.8900_JPRB /)
769
770	IF (KLEV == 19) THEN
771	RAOVLP=ZAOVLP(1)
772	RBOVLP=ZBOVLP(1)
773	ELSE IF (KLEV == 31) THEN
774	RAOVLP=ZAOVLP(2)
775	RBOVLP=ZBOVLP(2)
776	ELSE IF (KLEV == 60) THEN
777	RAOVLP=ZAOVLP(3)
778	RBOVLP=ZBOVLP(3)
779	ELSE
780	RAOVLP=ZAOVLP(3)
781	RBOVLP=ZBOVLP(3)
782	END IF
783
784	! ------------------------------------------------------------------
785
786	RETURN
787	END SUBROUTINE SUCLOPN

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

Download in other formats:

Original Format