LMDZ
suclopn.F90
Go to the documentation of this file.
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
yoesw::rebcua
real(kind=jprb), dimension(6) rebcua
Definition: yoesw.F90:39
yoesw::ryfwcb
real(kind=jprb), dimension(6) ryfwcb
Definition: yoesw.F90:33
yoesw::reffib
real(kind=jprb) reffib
Definition: yoesw.F90:99
yoesw::ryfwcc
real(kind=jprb), dimension(6) ryfwcc
Definition: yoesw.F90:34
yoesw
Definition: yoesw.F90:1
yoesw::reffia
real(kind=jprb) reffia
Definition: yoesw.F90:98
yoesw::rfubb1
real(kind=jprb), dimension(6) rfubb1
Definition: yoesw.F90:76
yoesw::rfueta
real(kind=jprb), dimension(16, 3) rfueta
Definition: yoesw.F90:57
yoesw::ryfwcd
real(kind=jprb), dimension(6) ryfwcd
Definition: yoesw.F90:35
yoesw::rflbb2
real(kind=jprb), dimension(6) rflbb2
Definition: yoesw.F90:66
yoesw::rflbb1
real(kind=jprb), dimension(6) rflbb1
Definition: yoesw.F90:65
yoesw::raswce
real(kind=jprb), dimension(6) raswce
Definition: yoesw.F90:54
yoesw::rflcc3
real(kind=jprb), dimension(6) rflcc3
Definition: yoesw.F90:71
yoesw::rsushfa
real(kind=jprb), dimension(4) rsushfa
Definition: yoesw.F90:94
yoesw::rflbb3
real(kind=jprb), dimension(6) rflbb3
Definition: yoesw.F90:67
yoesw::rebcuc
real(kind=jprb), dimension(6) rebcuc
Definition: yoesw.F90:41
yoesw::raswcb
real(kind=jprb), dimension(6) raswcb
Definition: yoesw.F90:51
yoesw::rebcuf
real(kind=jprb), dimension(6) rebcuf
Definition: yoesw.F90:44
dimphy::klev
integer, save klev
Definition: dimphy.F90:7
yoesw::rfucc2
real(kind=jprb), dimension(6) rfucc2
Definition: yoesw.F90:81
yoesw::rtiw
real(kind=jprb) rtiw
Definition: yoesw.F90:100
yoesw::ryfwcf
real(kind=jprb), dimension(6) ryfwcf
Definition: yoesw.F90:37
yoesw::rflcc1
real(kind=jprb), dimension(6) rflcc1
Definition: yoesw.F90:69
yoesw::rsusha
real(kind=jprb), dimension(6) rsusha
Definition: yoesw.F90:92
yoesw::rsushd
real(kind=jprb) rsushd
Definition: yoesw.F90:96
yoesw::rebcud
real(kind=jprb), dimension(6) rebcud
Definition: yoesw.F90:42
yoesw::rfucc1
real(kind=jprb), dimension(6) rfucc1
Definition: yoesw.F90:80
yoesw::rfldd3
real(kind=jprb), dimension(6) rfldd3
Definition: yoesw.F90:86
yoesw::rebcuh
real(kind=jprb), dimension(16) rebcuh
Definition: yoesw.F90:46
yoesw::rebcug
real(kind=jprb), dimension(16) rebcug
Definition: yoesw.F90:45
yoesw::ryfwca
real(kind=jprb), dimension(6) ryfwca
Definition: yoesw.F90:32
yoesw::rhsavi
real(kind=jprb), dimension(16, 3) rhsavi
Definition: yoesw.F90:59
yoesw::rriw
real(kind=jprb) rriw
Definition: yoesw.F90:101
yoesw::rebcub
real(kind=jprb), dimension(6) rebcub
Definition: yoesw.F90:40
yoerad
Definition: yoerad.F90:1
yoesw::rflcc0
real(kind=jprb), dimension(6) rflcc0
Definition: yoesw.F90:68
yoesw::ryfwce
real(kind=jprb), dimension(6) ryfwce
Definition: yoesw.F90:36
yoesw::rfuaa1
real(kind=jprb), dimension(6) rfuaa1
Definition: yoesw.F90:74
yoesw::rfucc0
real(kind=jprb), dimension(6) rfucc0
Definition: yoesw.F90:79
yoesw::rfldd2
real(kind=jprb), dimension(6) rfldd2
Definition: yoesw.F90:85
yoesw::rebcuj
real(kind=jprb), dimension(6) rebcuj
Definition: yoesw.F90:48
yoesw::rfubb2
real(kind=jprb), dimension(6) rfubb2
Definition: yoesw.F90:77
yoesw::rsushh
real(kind=jprb), dimension(6) rsushh
Definition: yoesw.F90:90
yoesw::rsushe
real(kind=jprb), dimension(6) rsushe
Definition: yoesw.F90:88
yoesw::rfuaa0
real(kind=jprb), dimension(6) rfuaa0
Definition: yoesw.F90:73
yoesw::rfldd1
real(kind=jprb), dimension(6) rfldd1
Definition: yoesw.F90:84
yoerad::raovlp
real(kind=jprb) raovlp
Definition: yoerad.F90:68
yoesw::rflaa1
real(kind=jprb), dimension(6) rflaa1
Definition: yoesw.F90:63
yoesw::rfucc3
real(kind=jprb), dimension(6) rfucc3
Definition: yoesw.F90:82
yoesw::raswcc
real(kind=jprb), dimension(6) raswcc
Definition: yoesw.F90:52
yoesw::rfulio
real(kind=jprb), dimension(16, 3) rfulio
Definition: yoesw.F90:58
yoesw::rflbb0
real(kind=jprb), dimension(6) rflbb0
Definition: yoesw.F90:64
yoesw::rfubb3
real(kind=jprb), dimension(6) rfubb3
Definition: yoesw.F90:78
suclopn
subroutine suclopn(KTSW, KSW, KLEV)
Definition: suclopn.F90:2
yoesw::rsushf
real(kind=jprb), dimension(6) rsushf
Definition: yoesw.F90:89
yoesw::rsushc
real(kind=jprb) rsushc
Definition: yoesw.F90:95
yoesw::rsushg
real(kind=jprb), dimension(6) rsushg
Definition: yoesw.F90:93
yoerad::rbovlp
real(kind=jprb) rbovlp
Definition: yoerad.F90:68
yoesw::raswcf
real(kind=jprb), dimension(6) raswcf
Definition: yoesw.F90:55
yoesw::rfldd0
real(kind=jprb), dimension(6) rfldd0
Definition: yoesw.F90:83
yoesw::rebcue
real(kind=jprb), dimension(6) rebcue
Definition: yoesw.F90:43
yoesw::rsushk
real(kind=jprb), dimension(6) rsushk
Definition: yoesw.F90:91
yoesw::rflcc2
real(kind=jprb), dimension(6) rflcc2
Definition: yoesw.F90:70
yoesw::rflaa0
real(kind=jprb), dimension(6) rflaa0
Definition: yoesw.F90:62
yoesw::raswcd
real(kind=jprb), dimension(6) raswcd
Definition: yoesw.F90:53
yoesw::raswca
real(kind=jprb), dimension(6) raswca
Definition: yoesw.F90:50
yoesw::rebcui
real(kind=jprb), dimension(6) rebcui
Definition: yoesw.F90:47
yoesw::rfubb0
real(kind=jprb), dimension(6) rfubb0
Definition: yoesw.F90:75