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 ! Lindner & Li, 2000: J. Climate, 13, 1797-1805
43 
44 ! AUTHOR.
45 ! -------
46 ! JEAN-JACQUES MORCRETTE *ECMWF*
47 
48 ! MODIFICATIONS.
49 ! --------------
50 ! ORIGINAL : 92-02-29
51 ! 97-04-16 JJ Morcrette 2 and 4 spectral intervals
52 ! 01-01-16 JJ Morcrette ice cloud from Fu 96, Fu et al 98
53 ! 01-03-01 JJ Morcrette water clouds from Lindner, Li 2000
54 ! 00-10-25 JJMorcrette 6 spectral intervals
55 ! M.Hamrud 01-Oct-2003 CY28 Cleaning
56 
57 ! ------------------------------------------------------------------
58 
59 USE parkind1 ,ONLY : jpim ,jprb
60 USE yomhook ,ONLY : lhook, dr_hook
61 
62 USE yoesw , ONLY : ryfwca ,ryfwcb ,ryfwcc ,ryfwcd ,&
63  & ryfwce ,ryfwcf ,rebcua ,rebcub ,rebcuc ,&
64  & rebcud ,rebcue ,rebcuf ,rebcug ,rebcuh ,&
65  & rebcui ,rebcuj ,raswca ,raswcb ,raswcc ,&
66  & raswcd ,raswce ,raswcf ,rsushe ,rsushf ,&
67  & rsushh ,rsushk ,rsusha ,rsushg ,rsushfa ,&
68  & rsushc ,rsushd ,reffia ,reffib ,rhsavi ,&
69  & rfuaa0 ,rfuaa1 ,rfubb0 ,rfubb1 ,rfubb2 ,&
70  & rfubb3 ,rfucc0 ,rfucc1 ,rfucc2 ,rfucc3 ,&
71  & rfueta ,rfuetb ,rfuetc ,rfulio ,rflaa0 ,rflaa1 ,rflbb0 ,&
72  & rflbb1 ,rflbb2 ,rflbb3 ,rflcc0 ,rflcc1 ,&
73  & rflcc2 ,rflcc3 ,rfldd0 ,rfldd1 ,rfldd2 ,&
74  & rfldd3 ,rlilia ,rlilib ,rtiw ,rriw
75 USE yoerad , ONLY : raovlp ,rbovlp ,lsrtm
76 
77 IMPLICIT NONE
78 
79 INTEGER(KIND=JPIM),INTENT(IN) :: KTSW
80 INTEGER(KIND=JPIM),INTENT(IN) :: KSW
81 INTEGER(KIND=JPIM),INTENT(IN) :: KLEV
82 INTEGER(KIND=JPIM) :: JNU
83 
84 ! -----------------------------------------------------------------
85 REAL(KIND=JPRB) :: ZEBCUA2(2) ,ZEBCUB2(2) ,ZEBCUC2(2) ,ZEBCUD2(2)&
86  & , ZEBCUE2(2) ,ZEBCUF2(2) ,ZYFWCA2(2) ,ZYFWCB2(2)&
87  & , ZYFWCC2(2) ,ZYFWCD2(2) ,ZYFWCE2(2) ,ZYFWCF2(2)&
88  & , ZASWCA2(2) ,ZASWCB2(2) ,ZASWCC2(2) ,ZASWCD2(2)&
89  & , ZASWCE2(2) ,ZASWCF2(2) ,ZSUSHE2(2) ,ZSUSHF2(2)&
90  & , ZSUSHH2(2) ,ZSUSHK2(2) ,ZSUSHA2(2) ,ZSUSHG2(2)&
91  & , ZFLAA02(2) ,ZFLAA12(2) ,ZFLBB02(2) ,ZFLBB12(2)&
92  & , ZFLBB22(2) ,ZFLBB32(2) ,ZFLCC02(2) ,ZFLCC12(2)&
93  & , ZFLCC22(2) ,ZFLCC32(2) ,ZFLDD02(2) ,ZFLDD12(2)&
94  & , ZFLDD22(2) ,ZFLDD32(2)
95 
96 REAL(KIND=JPRB) :: ZEBCUA4(4) ,ZEBCUB4(4) ,ZEBCUC4(4) ,ZEBCUD4(4)&
97  & , ZEBCUE4(4) ,ZEBCUF4(4) ,ZYFWCA4(4) ,ZYFWCB4(4)&
98  & , ZYFWCC4(4) ,ZYFWCD4(4) ,ZYFWCE4(4) ,ZYFWCF4(4)&
99  & , ZASWCA4(4) ,ZASWCB4(4) ,ZASWCC4(4) ,ZASWCD4(4)&
100  & , ZASWCE4(4) ,ZASWCF4(4) ,ZSUSHE4(4) ,ZSUSHF4(4)&
101  & , ZSUSHH4(4) ,ZSUSHK4(4) ,ZSUSHA4(4) ,ZSUSHG4(4)&
102  & , ZFLAA04(4) ,ZFLAA14(4) ,ZFLBB04(4) ,ZFLBB14(4)&
103  & , ZFLBB24(4) ,ZFLBB34(4) ,ZFLCC04(4) ,ZFLCC14(4)&
104  & , ZFLCC24(4) ,ZFLCC34(4) ,ZFLDD04(4) ,ZFLDD14(4)&
105  & , ZFLDD24(4) ,ZFLDD34(4)&
106  & , ZFUAA04(4) ,ZFUAA14(4) ,ZFUBB04(4) ,ZFUBB14(4)&
107  & , ZFUBB24(4) ,ZFUBB34(4) ,ZFUCC04(4) ,ZFUCC14(4)&
108  & , ZFUCC24(4) ,ZFUCC34(4)
109 
110 REAL(KIND=JPRB) :: ZAOVLP(3), ZBOVLP(3)
111 
112 REAL(KIND=JPRB) :: ZEBCUA6(6) ,ZEBCUB6(6) ,ZEBCUC6(6) ,ZEBCUD6(6)&
113  & , ZEBCUE6(6) ,ZEBCUF6(6) ,ZYFWCA6(6) ,ZYFWCB6(6)&
114  & , ZYFWCC6(6) ,ZYFWCD6(6) ,ZYFWCE6(6) ,ZYFWCF6(6)&
115  & , ZASWCA6(6) ,ZASWCB6(6) ,ZASWCC6(6) ,ZASWCD6(6)&
116  & , ZASWCE6(6) ,ZASWCF6(6) ,ZSUSHE6(6) ,ZSUSHF6(6)&
117  & , ZSUSHH6(6) ,ZSUSHK6(6) ,ZSUSHA6(6) ,ZSUSHG6(6)&
118  & , ZFLAA06(6) ,ZFLAA16(6) ,ZFLBB06(6) ,ZFLBB16(6)&
119  & , ZFLBB26(6) ,ZFLBB36(6) ,ZFLCC06(6) ,ZFLCC16(6)&
120  & , ZFLCC26(6) ,ZFLCC36(6) ,ZFLDD06(6) ,ZFLDD16(6)&
121  & , ZFLDD26(6) ,ZFLDD36(6)&
122  & , ZFUAA06(6) ,ZFUAA16(6) ,ZFUBB06(6) ,ZFUBB16(6)&
123  & , ZFUBB26(6) ,ZFUBB36(6) ,ZFUCC06(6) ,ZFUCC16(6)&
124  & , ZFUCC26(6) ,ZFUCC36(6)
125 
126 REAL(KIND=JPRB) :: ZEBCUA14(14) ,ZEBCUB14(14) ,ZEBCUC14(14) ,ZEBCUD14(14)&
127  & , ZEBCUE14(14) ,ZEBCUF14(14) ,ZYFWCA14(14) ,ZYFWCB14(14)&
128  & , ZYFWCC14(14) ,ZYFWCD14(14) ,ZYFWCE14(14) ,ZYFWCF14(14)&
129  & , ZASWCA14(14) ,ZASWCB14(14) ,ZASWCC14(14) ,ZASWCD14(14)&
130  & , ZASWCE14(14) ,ZASWCF14(14) ,ZSUSHE14(14) ,ZSUSHF14(14)&
131  & , ZSUSHH14(14) ,ZSUSHK14(14) ,ZSUSHA14(14) ,ZSUSHG14(14)&
132  & , ZFLAA014(14) ,ZFLAA114(14) ,ZFLBB014(14) ,ZFLBB114(14)&
133  & , ZFLBB214(14) ,ZFLBB314(14) ,ZFLCC014(14) ,ZFLCC114(14)&
134  & , ZFLCC214(14) ,ZFLCC314(14) ,ZFLDD014(14) ,ZFLDD114(14)&
135  & , ZFLDD214(14) ,ZFLDD314(14)&
136  & , ZFUAA014(14) ,ZFUAA114(14) ,ZFUBB014(14) ,ZFUBB114(14)&
137  & , ZFUBB214(14) ,ZFUBB314(14) ,ZFUCC014(14) ,ZFUCC114(14)&
138  & , ZFUCC214(14) ,ZFUCC314(14)
139 REAL(KIND=JPRB) :: ZHOOK_HANDLE
140 
141 #include "abor1.intfb.h"
142 
143 ! -----------------------------------------------------------------
144 
145 !* 1. SHORTWAVE CLOUD OPTICAL PROPERTIES
146 ! ----------------------------------
147 
148 ! ------------------------------------------------------------------
149 
150 !* 1.1 TWO SPECTRAL INTERVALS
151 ! ----------------------
152 
153 ! SW : 0.25 - 0.68 - 4.00 microns
154 
155 !* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
156 
157 ! optical properties
158 IF (lhook) CALL dr_hook('SUCLOPN',0,zhook_handle)
159 zebcua2 = (/ 3.448e-04_jprb , 3.448e-04_jprb /)
160 zebcub2 = (/ 2.431_jprb , 2.431_jprb /)
161 zebcuc2 = (/ 0.00001_jprb , 0.024366_jprb /)
162 zebcud2 = (/ 0.0_jprb , 2.487e-04_jprb /)
163 zebcue2 = (/ 0.7661_jprb , 0.7866_jprb /)
164 zebcuf2 = (/ 5.851e-04_jprb , 5.937e-04_jprb /)
165 
166 ! optical properties
167 ! ZEBCUA2 = (/ 3.448E-04 , 3.448E-04 /)
168 ! ZEBCUB2 = (/ 2.431 , 2.431 /)
169 ! ZEBCUC2 = (/ 0.00001 , 0.035589 /)
170 ! ZEBCUD2 = (/ 0. , 2.757E-04 /)
171 ! ZEBCUE2 = (/ 0.7661 , 0.7921 /)
172 ! ZEBCUF2 = (/ 5.851E-04 , 5.893E-04 /)
173 
174 !* Water cloud properties - from Fouquart (1987)
175 
176 zyfwca2 = (/ 0.0_jprb , 0.0_jprb /)
177 zyfwcb2 = (/ 1.5_jprb , 1.5_jprb /)
178 zyfwcc2 = (/ 0.9999_jprb , 0.9988_jprb /)
179 zyfwcd2 = (/ 5.000e-04_jprb , 2.500e-03_jprb /)
180 zyfwce2 = (/ 0.5_jprb , 0.05_jprb /)
181 zyfwcf2 = (/ 0.865_jprb , 0.910_jprb /)
182 
183 !* Water cloud properties - from Slingo (1989)
184 
185 zaswca2 = (/ 2.817_jprb , 2.455_jprb /)
186 zaswcb2 = (/ 1.305_jprb , 1.400_jprb /)
187 zaswcc2 = (/-5.62e-08_jprb , 1.28e-02_jprb /)
188 zaswcd2 = (/ 1.63e-07_jprb , 8.97e-04_jprb /)
189 zaswce2 = (/ 0.829_jprb , 0.783_jprb /)
190 zaswcf2 = (/ 2.482_jprb , 4.988_jprb /)
191 
192 !* Ice cloud properties - from Sun and Shine (1995)
193 
194 zsushe2 = (/ 0.0_jprb , 8.6822_jprb /)
195 zsushf2 = (/ 0.0_jprb , 9.6277_jprb /)
196 zsushh2 = (/ 0.8522_jprb , 0.8819_jprb /)
197 zsushk2 = (/ 0.1620_jprb , 0.1630_jprb /)
198 zsusha2 = (/ 0.0_jprb , 23.204_jprb /)
199 zsushg2 = (/ 0.3270_jprb , 0.4180_jprb /)
200 
201 !* Ice cloud properties - from Fu and Liou (1993)
202 
203 zflaa02 = (/-6.656e-3_jprb ,-6.656e-3_jprb /)
204 zflaa12 = (/ 3.686_jprb , 3.686_jprb /)
205 zflbb02 = (/ .10998e-5_jprb , .21136e-1_jprb /)
206 zflbb12 = (/-.26101e-7_jprb , .39150e-3_jprb /)
207 zflbb22 = (/ .18096e-8_jprb ,-.20740e-6_jprb /)
208 zflbb32 = (/-.47387e-11_jprb,-.28829e-8_jprb /)
209 zflcc02 = (/ .22110e+1_jprb , .22498e+1_jprb /)
210 zflcc12 = (/-.10398e-2_jprb , .23656e-3_jprb /)
211 zflcc22 = (/ .65199e-4_jprb , .51948e-4_jprb /)
212 zflcc32 = (/-.34498e-6_jprb ,-.29768e-6_jprb /)
213 zfldd02 = (/ .12495_jprb , .11716_jprb /)
214 zfldd12 = (/-.43582e-3_jprb ,-.45208e-3_jprb /)
215 zfldd22 = (/ .14092e-4_jprb , .12772e-4_jprb /)
216 zfldd32 = (/-.69565e-7_jprb ,-.62779e-7_jprb /)
217 
218 !* Ice cloud properties - from Fu (1996)
219 
220 !ZFUAA02 = (/-.291721E-04_JPRB ,
221 !ZFUAA12 = (/ .251925E+01_JPRB ,
222 !ZFUBB02 = (/ .135403E-06_JPRB ,
223 !ZFUBB12 = (/ .992822E-07_JPRB ,
224 !ZFUBB22 = (/-.738432E-10_JPRB ,
225 !ZFUBB32 = (/ .331119E-12_JPRB ,
226 !ZFUCC02 = (/ .748127E+00_JPRB ,
227 !ZFUCC12 = (/ .956845E-03_JPRB ,
228 !ZFUCC22 = (/-.111517E-05_JPRB ,
229 !ZFUCC32 = (/-.815573E-08_JPRB ,
230 !ZFUDD02 = (/ .115730E+00_JPRB ,
231 !ZFUDD12 = (/ .256481E-03_JPRB ,
232 !ZFUDD22 = (/ .191313E-05_JPRB ,
233 !ZFUDD32 = (/-.124603E-07_JPRB ,
234 ! ------------------------------------------------------------------
235 
236 !* 1.2 FOUR SPECTRAL INTERVALS
237 ! -----------------------
238 
239 ! SW : 4 spectral intervals (0.25 - 0.69 - 1.19 - 2.38 - 4.00)
240 
241 !* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
242 
243 ! ZEBCUA4 = (/ 3.448E-03 , 3.448E-03 , 3.448E-03 , 3.448E-03 /)
244 ! ZEBCUB4 = (/ 2.431 , 2.431 , 2.431 , 2.431 /)
245 ! ZEBCUC4 = (/ 0.00001 , 0.00011 , 0.01861 , 0.46658 /)
246 ! ZEBCUD4 = (/ 0. , 1.405E-05 , 8.328E-4 , 2.050E-05 /)
247 ! ZEBCUE4 = (/ 0.7661 , 0.7730 , 0.7940 , 0.9595 /)
248 ! ZEBCUF4 = (/ 5.851E-04 , 5.665E-04 , 7.267E-04 , 1.076E-04 /)
249 
250 zebcua4 = (/ 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb /)
251 zebcub4 = (/ 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb /)
252 zebcuc4 = (/ 0.00001_jprb , 0.00011_jprb , 0.0197796_jprb , 0.46658_jprb /)
253 zebcud4 = (/ 0.0_jprb , 1.405e-05_jprb , 7.95513e-4_jprb, 2.050e-05_jprb /)
254 zebcue4 = (/ 0.7661_jprb , 0.7730_jprb , 0.795653_jprb , 0.9595_jprb /)
255 zebcuf4 = (/ 5.851e-04_jprb , 5.665e-04_jprb , 7.267e-04_jprb , 1.076e-04_jprb /)
256 
257 !* Water cloud properties - from Fouquart (1987)
258 
259 zyfwca4 = (/ 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb /)
260 zyfwcb4 = (/ 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb /)
261 zyfwcc4 = (/ 0.9999_jprb , 0.9988_jprb , 0.9988_jprb , 0.9988_jprb /)
262 zyfwcd4 = (/ 5.000e-04_jprb , 2.500e-03_jprb , 2.500e-03_jprb, 2.500e-03_jprb /)
263 zyfwce4 = (/ 0.5_jprb , 0.05_jprb , 0.05_jprb , 0.05_jprb /)
264 zyfwcf4 = (/ 0.865_jprb , 0.910_jprb , 0.910_jprb , 0.910_jprb /)
265 
266 !* Water cloud properties - from Slingo (1989)
267 
268 zaswca4 = (/ 2.817_jprb , 2.682_jprb , 2.264_jprb , 1.281_jprb /)
269 zaswcb4 = (/ 1.305_jprb , 1.346_jprb , 1.454_jprb , 1.641_jprb /)
270 zaswcc4 = (/-5.62e-08_jprb ,-6.94e-06_jprb , 4.64e-04_jprb , 2.01e-01_jprb /)
271 zaswcd4 = (/ 1.63e-07_jprb , 2.35e-05_jprb , 1.24e-03_jprb , 7.56e-03_jprb /)
272 zaswce4 = (/ 0.829_jprb , 0.794_jprb , 0.754_jprb , 0.826_jprb /)
273 zaswcf4 = (/ 2.482_jprb , 4.226_jprb , 6.560_jprb , 4.353_jprb /)
274 
275 !* Ice cloud properties - from Sun and Shine (1995)
276 
277 zsushe4 = (/ 0.0_jprb , 7.2471e-02_jprb , 17.5933_jprb , 48.7166_jprb /)
278 zsushf4 = (/ 0.0_jprb , 4.01511_jprb , 21.1249_jprb , 1.20890_jprb /)
279 zsushh4 = (/ 0.8522_jprb , 0.85841_jprb , 0.90778_jprb ,0.982046_jprb /)
280 zsushk4 = (/ 0.1620_jprb , 0.160048_jprb , 0.188521_jprb ,0.0411446_jprb /)
281 zsusha4 = (/ 0.0_jprb , 0.273455_jprb , 41.7675_jprb , 161.104_jprb /)
282 zsushg4 = (/ 0.3270_jprb , 0.343668_jprb , 0.526192_jprb , 0.574040_jprb /)
283 
284 !* Ice cloud properties - from Fu and Liou (1993)
285 
286 zflaa04 = (/-6.656e-03_jprb ,-6.656e-03_jprb,-6.656e-03_jprb,-6.656e-03_jprb /)
287 zflaa14 = (/ 3.686_jprb , 3.686_jprb , 3.686_jprb , 3.686_jprb /)
288 zflbb04 = (/ .10998e-5_jprb , .20208e-4_jprb, .51557e-3_jprb, .39517e+0_jprb /)
289 zflbb14 = (/-.26101e-7_jprb , .96483e-5_jprb, .10731e-2_jprb, .15787e-2_jprb /)
290 zflbb24 = (/ .18096e-8_jprb , .83009e-7_jprb, .17753e-5_jprb,-.14337e-4_jprb /)
291 zflbb34 = (/-.47387e-11_jprb,-.32217e-9_jprb,-.18379e-7_jprb, .46942e-7_jprb /)
292 zflcc04 = (/ .22110e+1_jprb , .22151e+1_jprb, .22534e+1_jprb, .26653e+1_jprb /)
293 zflcc14 = (/-.10398e-2_jprb ,-.77982e-3_jprb, .16163e-2_jprb, .56935e-2_jprb /)
294 zflcc24 = (/ .65199e-4_jprb , .63750e-4_jprb, .44037e-4_jprb,-.54316e-4_jprb /)
295 zflcc34 = (/-.34498e-6_jprb ,-.34466e-6_jprb,-.27627e-6_jprb, .17858e-6_jprb /)
296 zfldd04 = (/ .12495_jprb , .12363_jprb , .11983_jprb , .21834e-1_jprb /)
297 zfldd14 = (/-.43582e-3_jprb ,-.44419e-3_jprb,-.50108e-3_jprb,-.29204e-3_jprb /)
298 zfldd24 = (/ .14092e-4_jprb , .14038e-4_jprb, .11843e-4_jprb, .18060e-5_jprb /)
299 zfldd34 = (/-.69565e-7_jprb ,-.68851e-7_jprb,-.59367e-7_jprb,-.46257e-8_jprb /)
300 
301 !* Ice cloud properties - from Fu (1996) as tabulated in Sun & Rikus (1999)
302 
303 zfuaa04 = (/-1.30817e-04_jprb,-6.39479e-05_jprb,-6.74730e-06_jprb, 1.62674e-04_jprb /)
304 zfuaa14 = (/ 2.52883e+00_jprb, 2.52393e+00_jprb, 2.52056e+00_jprb, 2.49823e+00_jprb /)
305 zfubb04 = (/-1.55357e-08_jprb,-7.90657e-07_jprb, 9.64842e-04_jprb, 2.25112e-01_jprb /)
306 zfubb14 = (/ 1.95793e-07_jprb, 7.79991e-06_jprb, 9.09809e-04_jprb, 3.05017e-03_jprb /)
307 zfubb24 = (/-2.31234e-10_jprb, 2.90894e-10_jprb,-3.57557e-06_jprb,-2.54236e-05_jprb /)
308 zfubb34 = (/ 1.12247e-12_jprb,-2.02818e-12_jprb, 1.00197e-08_jprb, 8.49116e-08_jprb /)
309 zfucc04 = (/ 7.39781e-01_jprb, 7.52335e-01_jprb, 7.56307e-01_jprb, 8.30812e-01_jprb /)
310 zfucc14 = (/ 9.10564e-04_jprb, 1.06211e-03_jprb, 1.73364e-03_jprb, 2.62788e-03_jprb /)
311 zfucc24 = (/-4.62479e-07_jprb,-2.45770e-06_jprb,-8.92191e-06_jprb,-2.43196e-05_jprb /)
312 zfucc34 = (/-1.05910e-08_jprb,-3.03712e-09_jprb, 1.97757e-08_jprb, 8.23543e-08_jprb /)
313 
314 ! ------------------------------------------------------------------
315 
316 !* 1.3 SIX SPECTRAL INTERVALS
317 ! ----------------------
318 
319 ! SW : 6 spectral intervals (0.185-0.25-0.44-0.69-1.19-2.38-4.00)
320 
321 !* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
322 
323 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 /)
324 zebcub6 = (/ 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb /)
325 zebcuc6 = (/ 0.00001_jprb , 0.00001_jprb , 0.00001_jprb , 0.00011_jprb , 0.0197796_jprb , 0.46658_jprb /)
326 zebcud6 = (/ 0.0_jprb , 0.0_jprb , 0.0_jprb , 1.405e-05_jprb , 7.95513e-4_jprb, 2.050e-05_jprb /)
327 zebcue6 = (/ 0.7661_jprb , 0.7661_jprb , 0.7661_jprb , 0.7730_jprb , 0.795653_jprb , 0.9595_jprb /)
328 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 /)
329 
330 !* Water cloud properties - from Fouquart (1987)
331 
332 zyfwca6 = (/ 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb /)
333 zyfwcb6 = (/ 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb /)
334 zyfwcc6 = (/ 0.9999_jprb , 0.9999_jprb , 0.9999_jprb , 0.9988_jprb , 0.9988_jprb , 0.9988_jprb /)
335 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 /)
336 zyfwce6 = (/ 0.5_jprb , 0.5_jprb , 0.5_jprb , 0.05_jprb , 0.05_jprb , 0.05_jprb /)
337 zyfwcf6 = (/ 0.865_jprb , 0.865_jprb , 0.865_jprb , 0.910_jprb , 0.910_jprb , 0.910_jprb /)
338 
339 !* Water cloud properties - from Slingo (1989)
340 
341 zaswca6 = (/ 3.094_jprb , 2.881_jprb , 2.790_jprb , 2.682_jprb , 2.254_jprb , 1.281_jprb /)
342 zaswcb6 = (/ 1.252_jprb , 1.284_jprb , 1.313_jprb , 1.346_jprb , 1.456_jprb , 1.641_jprb /)
343 zaswcc6 = (/ 7.90e-07_jprb , 1.77e-07_jprb ,-1.55e-07_jprb ,-6.94e-06_jprb , 4.75e-04_jprb , 2.01e-01_jprb /)
344 zaswcd6 = (/ 3.69e-07_jprb , 1.23e-07_jprb , 1.80e-07_jprb , 2.35e-05_jprb , 1.26e-03_jprb , 7.56e-03_jprb /)
345 zaswce6 = (/ 0.844_jprb , 0.839_jprb , 0.824_jprb , 0.794_jprb , 0.754_jprb , 0.826_jprb /)
346 zaswcf6 = (/ 1.558_jprb , 1.914_jprb , 2.723_jprb , 4.226_jprb , 6.561_jprb , 4.353_jprb /)
347 
348 !* Ice cloud properties - from Sun and Shine (1995)
349 
350 zsushe6 = (/ 0.0_jprb , 0.0_jprb , 0.0_jprb , 7.2471e-02_jprb , 17.5933_jprb , 48.7166_jprb /)
351 zsushf6 = (/ 0.0_jprb , 0.0_jprb , 0.0_jprb , 4.01511_jprb , 21.1249_jprb , 1.20890_jprb /)
352 zsushh6 = (/ 0.8522_jprb , 0.8522_jprb , 0.8522_jprb , 0.85841_jprb , 0.90778_jprb ,0.982046_jprb /)
353 zsushk6 = (/ 0.1620_jprb , 0.1620_jprb , 0.1620_jprb , 0.160048_jprb , 0.188521_jprb ,0.0411446_jprb /)
354 zsusha6 = (/ 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.273455_jprb , 41.7675_jprb , 161.104_jprb /)
355 zsushg6 = (/ 0.3270_jprb , 0.3270_jprb , 0.3270_jprb , 0.343668_jprb , 0.526192_jprb , 0.574040_jprb /)
356 
357 !* Ice cloud properties - from Fu and Liou (1993)
358 
359 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 /)
360 zflaa16 = (/ 3.686_jprb , 3.686_jprb , 3.686_jprb , 3.686_jprb , 3.686_jprb , 3.686_jprb /)
361 zflbb06 = (/ .10998e-5_jprb , .10998e-5_jprb , .10998e-5_jprb , .20208e-4_jprb, .51557e-3_jprb, .39517e+0_jprb /)
362 zflbb16 = (/-.26101e-7_jprb ,-.26101e-7_jprb ,-.26101e-7_jprb , .96483e-5_jprb, .10731e-2_jprb, .15787e-2_jprb /)
363 zflbb26 = (/ .18096e-8_jprb , .18096e-8_jprb , .18096e-8_jprb , .83009e-7_jprb, .17753e-5_jprb,-.14337e-4_jprb /)
364 zflbb36 = (/-.47387e-11_jprb,-.47387e-11_jprb,-.47387e-11_jprb,-.32217e-9_jprb,-.18379e-7_jprb, .46942e-7_jprb /)
365 zflcc06 = (/ .22110e+1_jprb , .22110e+1_jprb , .22110e+1_jprb , .22151e+1_jprb, .22534e+1_jprb, .26653e+1_jprb /)
366 zflcc16 = (/-.10398e-2_jprb ,-.10398e-2_jprb ,-.10398e-2_jprb ,-.77982e-3_jprb, .16163e-2_jprb, .56935e-2_jprb /)
367 zflcc26 = (/ .65199e-4_jprb , .65199e-4_jprb , .65199e-4_jprb , .63750e-4_jprb, .44037e-4_jprb,-.54316e-4_jprb /)
368 zflcc36 = (/-.34498e-6_jprb ,-.34498e-6_jprb ,-.34498e-6_jprb ,-.34466e-6_jprb,-.27627e-6_jprb, .17858e-6_jprb /)
369 zfldd06 = (/ .12495_jprb , .12495_jprb , .12495_jprb , .12363_jprb , .11983_jprb , .21834e-1_jprb /)
370 zfldd16 = (/-.43582e-3_jprb ,-.43582e-3_jprb ,-.43582e-3_jprb ,-.44419e-3_jprb,-.50108e-3_jprb,-.29204e-3_jprb /)
371 zfldd26 = (/ .14092e-4_jprb , .14092e-4_jprb , .14092e-4_jprb , .14038e-4_jprb, .11843e-4_jprb, .18060e-5_jprb /)
372 zfldd36 = (/-.69565e-7_jprb ,-.69565e-7_jprb ,-.69565e-7_jprb ,-.68851e-7_jprb,-.59367e-7_jprb,-.46257e-8_jprb /)
373 
374 !* Ice cloud properties - from Fu (1996) as tabulated in Sun & Rikus (1999)
375 
376 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 /)
377 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 /)
378 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 /)
379 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 /)
380 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 /)
381 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 /)
382 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 /)
383 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 /)
384 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 /)
385 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 /)
386 
387 ! ------------------------------------------------------------------
388 
389 !* 1.4 FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW
390 ! ---------------------------------------
391 
392 ! SW : 14 spectral intervals
393 ! 3.846 - 3.077
394 ! 3.077 - 2.500
395 ! 2.500 - 2.150
396 ! 2.150 - 1.942
397 ! 1.942 - 1.626
398 ! 1.626 - 1.299
399 ! 1.299 - 1.242
400 ! 1.242 - 0.7782
401 ! 0.7782- 0.6250
402 ! 0.6250- 0.4415
403 ! 0.4415- 0.3448
404 ! 0.3448- 0.2632
405 ! 0.2632- 0.2000
406 ! 12.195 - 3.846
407 
408 !* Ice cloud properties - crystal: remapped from Ebert and Curry, 1992
409 
410 zebcua14 = (/ &
411  & 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb &
412  & , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb , 3.448e-03_jprb &
413  & , 3.448e-03_jprb , 3.448e-03_jprb /)
414 
415 zebcub14 = (/ &
416  & 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb &
417  & , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb , 2.431_jprb &
418  & , 2.431_jprb , 2.431_jprb /)
419 
420 zebcuc14 = (/ &
421  & 0.46658_jprb , 0.46658_jprb , 0.03779_jprb , 0.03779_jprb , 0.01240_jprb , 0.01240_jprb &
422  & , 0.00011_jprb , 0.00011_jprb , 0.00001_jprb , 0.00001_jprb , 0.00001_jprb , 0.00001_jprb &
423  & , 0.00001_jprb , 0.46658_jprb /)
424 
425 zebcud14 = (/ &
426  & 2.050e-05_jprb , 2.050e-05_jprb , 1.284e-03_jprb , 1.284e-03_jprb , 6.867e-04_jprb , 6.867e-04_jprb &
427  & , 1.405e-05_jprb , 1.405e-05_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
428  & , 0.0_jprb , 2.050e-05_jprb /)
429 
430 zebcue14 = (/ &
431  & 0.9595_jprb , 0.9595_jprb , 0.8172_jprb , 0.8172_jprb , 0.7865_jprb , 0.7865_jprb &
432  & , 0.7730_jprb , 0.7730_jprb , 0.7661_jprb , 0.7661_jprb , 0.7661_jprb , 0.7661_jprb &
433  & , 0.7661_jprb , 0.9595_jprb /)
434 
435 zebcuf14 = (/ &
436  & 1.076e-04_jprb , 1.076e-04_jprb , 7.463e-04_jprb , 7.463e-04_jprb , 7.204e-04_jprb , 7.204e-04_jprb &
437  & , 5.665e-04_jprb , 5.665e-04_jprb , 5.851e-04_jprb , 5.851e-04_jprb , 5.851e-04_jprb , 5.851e-04_jprb &
438  & , 5.851e-04_jprb , 1.076e-04_jprb /)
439 
440 !* Water cloud properties - from Fouquart (1987)
441 
442 zyfwca14 = (/ &
443  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
444  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
445  & , 0.0_jprb , 0.0_jprb /)
446 
447 zyfwcb14 = (/ &
448  & 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb &
449  & , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb , 1.5_jprb &
450  & , 1.5_jprb , 1.5_jprb /)
451 
452 zyfwcc14 = (/ &
453  & 0.9988_jprb , 0.9988_jprb , 0.9988_jprb , 0.9988_jprb , 0.9988_jprb , 0.9988_jprb &
454  & , 0.9988_jprb , 0.9988_jprb , 0.9999_jprb , 0.9999_jprb , 0.9999_jprb , 0.9999_jprb &
455  & , 0.9999_jprb , 0.9988_jprb /)
456 
457 zyfwcd14 = (/ &
458  & 2.500e-03_jprb , 2.500e-03_jprb , 2.500e-03_jprb , 2.500e-03_jprb , 2.500e-03_jprb , 2.500e-03_jprb &
459  & , 2.500e-03_jprb , 2.500e-03_jprb , 5.000e-04_jprb , 5.000e-04_jprb , 5.000e-04_jprb , 5.000e-04_jprb &
460  & , 5.000e-04_jprb , 2.500e-03_jprb /)
461 
462 zyfwce14 = (/ &
463  & 0.05_jprb , 0.05_jprb , 0.05_jprb , 0.05_jprb , 0.05_jprb , 0.05_jprb &
464  & , 0.05_jprb , 0.05_jprb , 0.5_jprb , 0.5_jprb , 0.5_jprb , 0.5_jprb &
465  & , 0.5_jprb , 0.05_jprb /)
466 
467 zyfwcf14 = (/ &
468  & 0.910_jprb , 0.910_jprb , 0.910_jprb , 0.910_jprb , 0.910_jprb , 0.910_jprb &
469  & , 0.910_jprb , 0.910_jprb , 0.865_jprb , 0.865_jprb , 0.865_jprb , 0.865_jprb &
470  & , 0.865_jprb , 0.910_jprb /)
471 
472 !* Water cloud properties - from Slingo (1989)
473 
474 zaswca14 = (/ &
475  & -1.023_jprb , 1.950_jprb , 1.850_jprb , 1.970_jprb , 1.970_jprb , 2.463_jprb &
476  & , 2.551_jprb , 2.622_jprb , 2.895_jprb , 2.672_jprb , 2.801_jprb , 2.944_jprb &
477  & , 3.094_jprb ,-1.023_jprb /)
478 
479 zaswcb14 = (/ &
480  & 1.933_jprb , 1.540_jprb , 1.556_jprb , 1.501_jprb , 1.501_jprb , 1.420_jprb &
481  & , 1.401_jprb , 1.362_jprb , 1.315_jprb , 1.320_jprb , 1.293_jprb , 1.270_jprb &
482  & , 1.252_jprb , 1.933_jprb /)
483 
484 zaswcc14 = (/ &
485  & 2.50e-02_jprb , 4.49e-01_jprb , 1.90e-04_jprb , 1.20e-03_jprb , 1.20e-03_jprb , 2.40e-04_jprb &
486  & , 6.20e-05_jprb , 3.30e-06_jprb ,-1.20e-07_jprb , 0.0_jprb , 1.00e-06_jprb ,-6.50e-07_jprb &
487  & , 7.90e-07_jprb , 2.50e-02_jprb /)
488 
489 zaswcd14 = (/ &
490  & 1.22e-02_jprb , 1.54e-03_jprb , 2.54e-03_jprb , 2.16e-03_jprb , 2.16e-03_jprb , 8.56e-04_jprb &
491  & , 2.60e-04_jprb , 2.80e-06_jprb , 4.40e-07_jprb , 0.0_jprb , 0.0_jprb , 4.33e-07_jprb &
492  & , 3.69e-07_jprb , 1.22e-02_jprb /)
493 
494 zaswce14 = (/ &
495  & 0.726_jprb , 0.831_jprb , 0.769_jprb , 0.740_jprb , 0.740_jprb , 0.754_jprb &
496  & , 0.773_jprb , 0.806_jprb , 0.818_jprb , 0.828_jprb , 0.836_jprb , 0.841_jprb &
497  & , 0.844_jprb , 0.726_jprb /)
498 
499 zaswcf14 = (/ &
500  & 6.652_jprb , 6.102_jprb , 5.171_jprb , 7.469_jprb , 7.469_jprb , 6.555_jprb &
501  & , 5.405_jprb , 3.355_jprb , 2.989_jprb , 2.467_jprb , 2.153_jprb , 1.680_jprb &
502  & , 1.558_jprb , 6.652_jprb /)
503 
504 !* Ice cloud properties - from Sun and Shine (1995)
505 
506 zsushe14 = (/ &
507  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
508  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
509  & , 0.0_jprb , 0.0_jprb /)
510 zsushf14 = (/ &
511  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
512  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
513  & , 0.0_jprb , 0.0_jprb /)
514 zsushh14 = (/ &
515  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
516  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
517  & , 0.0_jprb , 0.0_jprb /)
518 zsushk14 = (/ &
519  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
520  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
521  & , 0.0_jprb , 0.0_jprb /)
522 zsusha14 = (/ &
523  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
524  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
525  & , 0.0_jprb , 0.0_jprb /)
526 zsushg14 = (/ &
527  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
528  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
529  & , 0.0_jprb , 0.0_jprb /)
530 
531 !* Ice cloud properties - from Fu and Liou (1993)
532 
533 zflaa014 = (/ &
534  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
535  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
536  & , 0.0_jprb , 0.0_jprb /)
537 zflaa114 = (/ &
538  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
539  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
540  & , 0.0_jprb , 0.0_jprb /)
541 zflbb014 = (/ &
542  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
543  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
544  & , 0.0_jprb , 0.0_jprb /)
545 zflbb114 = (/ &
546  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
547  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
548  & , 0.0_jprb , 0.0_jprb /)
549 zflbb214 = (/ &
550  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
551  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
552  & , 0.0_jprb , 0.0_jprb /)
553 zflbb314 = (/ &
554  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
555  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
556  & , 0.0_jprb , 0.0_jprb /)
557 zflcc014 = (/ &
558  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
559  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
560  & , 0.0_jprb , 0.0_jprb /)
561 zflcc114 = (/ &
562  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
563  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
564  & , 0.0_jprb , 0.0_jprb /)
565 zflcc214 = (/ &
566  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
567  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
568  & , 0.0_jprb , 0.0_jprb /)
569 zflcc314 = (/ &
570  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
571  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
572  & , 0.0_jprb , 0.0_jprb /)
573 zfldd014 = (/ &
574  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
575  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
576  & , 0.0_jprb , 0.0_jprb /)
577 zfldd114 = (/ &
578  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
579  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
580  & , 0.0_jprb , 0.0_jprb /)
581 zfldd214 = (/ &
582  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
583  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
584  & , 0.0_jprb , 0.0_jprb /)
585 zfldd314 = (/ &
586  & 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
587  & , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb , 0.0_jprb &
588  & , 0.0_jprb , 0.0_jprb /)
589 
590 !* Ice cloud properties - from Fu (1996)
591 
592 zfuaa014 = (/ &
593  & 1.87598e-04_jprb, 2.97295e-04_jprb, 4.89477e-04_jprb,-8.37325e-06_jprb,-8.37325e-06_jprb,-8.05155e-04_jprb &
594  & , 6.51659e-05_jprb, 8.10443e-05_jprb, 1.61983e-04_jprb,-9.45458e-05_jprb,-2.58858e-04_jprb,-2.66955e-04_jprb &
595  & ,-2.36447e-04_jprb,-2.54823e-04_jprb /)
596 
597 zfuaa114 = (/ &
598  & 2.51396e+00_jprb, 2.48895e+00_jprb, 2.48776e+00_jprb, 2.52504e+00_jprb, 2.52504e+00_jprb, 2.57600e+00_jprb &
599  & , 2.51660e+00_jprb, 2.51619e+00_jprb, 2.50746e+00_jprb, 2.52061e+00_jprb, 2.53815e+00_jprb, 2.54179e+00_jprb &
600  & , 2.53817e+00_jprb, 2.52909e+00_jprb /)
601 
602 zfubb014 = (/ &
603  & 1.96793e-01_jprb, 4.64416e-01_jprb, 5.83469e-04_jprb, 2.53234e-03_jprb, 2.53234e-03_jprb,-2.85518e-05_jprb &
604  & ,-1.48012e-07_jprb,-1.57963e-07_jprb,-7.78001e-08_jprb, 5.08447e-07_jprb,-1.98529e-07_jprb,-1.00570e-07_jprb &
605  & ,-2.69916e-07_jprb, 2.60155e-01_jprb /)
606 
607 zfubb114 = (/ &
608  & 5.75235e-03_jprb, 2.04716e-05_jprb, 1.18127e-03_jprb, 1.75078e-03_jprb, 1.75078e-03_jprb, 1.71993e-03_jprb &
609  & , 9.02355e-05_jprb, 1.72475e-06_jprb, 2.53360e-07_jprb, 2.73206e-08_jprb, 9.39480e-08_jprb, 1.60441e-07_jprb &
610  & , 2.12909e-07_jprb, 5.45547e-03_jprb /)
611 
612 zfubb214 = (/ &
613  & -5.29220e-05_jprb,-4.60375e-07_jprb,-3.40011e-06_jprb,-8.00994e-06_jprb,-8.00994e-06_jprb,-7.43697e-06_jprb &
614  & ,-1.98190e-08_jprb, 9.02156e-11_jprb,-1.15489e-10_jprb, 4.96553e-11_jprb,-2.54540e-10_jprb,-2.05663e-10_jprb &
615  & ,-2.65397e-10_jprb,-5.58760e-05_jprb /)
616 
617 zfubb314 = (/ &
618  & 1.76618e-07_jprb, 2.03198e-09_jprb, 8.78549e-09_jprb, 2.31309e-08_jprb, 2.31309e-08_jprb, 2.09647e-08_jprb &
619  & , 4.01914e-11_jprb,-3.79423e-13_jprb, 4.65084e-13_jprb,-1.86001e-13_jprb, 1.10876e-12_jprb, 8.88595e-13_jprb &
620  & , 1.12983e-12_jprb, 1.97086e-07_jprb /)
621 
622 zfucc014 = (/ &
623  & 0.759183_jprb , 0.919599_jprb , 0.775916_jprb , 0.758748_jprb , 0.758748_jprb , 0.752528_jprb &
624  & , 0.751277_jprb , 0.752318_jprb , 0.749693_jprb , 0.749856_jprb , 0.743546_jprb , 0.737809_jprb &
625  & , 0.733260_jprb , 0.799084_jprb /)
626 
627 zfucc114 = (/ &
628  & 4.93765e-03_jprb, 5.03025e-04_jprb, 1.74517e-03_jprb, 2.02709e-03_jprb, 2.02709e-03_jprb, 1.95748e-03_jprb &
629  & , 1.29824e-03_jprb, 1.04224e-03_jprb, 1.05446e-03_jprb, 8.89161e-04_jprb, 9.08674e-04_jprb, 8.97515e-04_jprb &
630  & , 9.18317e-04_jprb, 4.81706e-03_jprb /)
631 
632 zfucc214 = (/ &
633  & -4.84059e-05_jprb,-5.74771e-06_jprb,-9.21314e-06_jprb,-1.17029e-05_jprb,-1.17029e-05_jprb,-1.02495e-05_jprb &
634  & ,-4.99075e-06_jprb,-2.26618e-06_jprb,-2.32576e-06_jprb,-3.49578e-07_jprb,-4.65326e-07_jprb,-2.17099e-07_jprb &
635  & ,-4.22974e-07_jprb,-5.13220e-05_jprb /)
636 
637 zfucc314 = (/ &
638  & 1.65801e-07_jprb, 2.01731e-08_jprb, 2.15003e-08_jprb, 2.95195e-08_jprb, 2.95195e-08_jprb, 2.35479e-08_jprb &
639  & , 6.33757e-09_jprb,-3.68283e-09_jprb,-3.58307e-09_jprb, 1.09913e-08_jprb,-1.05786e-08_jprb,-1.16090e-08_jprb &
640  & ,-1.07976e-08_jprb, 1.84420e-07_jprb /)
641 
642 ! ------------------------------------------------------------------
643 
644 ! LW : absorption coefficient as a function of effective radius in RRTM
645 
646 ! water clouds from Savijarvi
647 
648 rhsavi( 1, :) = (/ 0.1651082_jprb , -0.003494839_jprb , 0.0_jprb /)
649 rhsavi( 2, :) = (/ 0.327820597_jprb, -0.0219634383_jprb , 0.000506783898_jprb /)
650 rhsavi( 3, :) = (/ 0.504805453_jprb, -0.0478602354_jprb , 0.00141521102_jprb /)
651 rhsavi( 4, :) = (/ 0.513169093_jprb, -0.0508960145_jprb , 0.00155498711_jprb /)
652 rhsavi( 5, :) = (/ 0.448042082_jprb, -0.0431857592_jprb , 0.00130848978_jprb /)
653 rhsavi( 6, :) = (/ 0.249547237_jprb, -0.0185273835_jprb , 0.00050361258_jprb /)
654 rhsavi( 7, :) = (/ 0.135486796_jprb, -0.00585852322_jprb, 0.00011873119_jprb /)
655 rhsavi( 8, :) = (/ 0.126668819_jprb, -0.00497949082_jprb, 0.0000927679172_jprb/)
656 rhsavi( 9, :) = (/ 0.130938752_jprb, -0.00531814674_jprb, 0.0000984953029_jprb/)
657 rhsavi(10, :) = (/ 0.147024519_jprb, -0.00671655774_jprb, 0.000130885091_jprb /)
658 rhsavi(11, :) = (/ 0.241463914_jprb, -0.0174486461_jprb , 0.000457756556_jprb /)
659 rhsavi(12, :) = (/ 0.07576579_jprb , -0.001695588_jprb , 0.0_jprb /)
660 rhsavi(13, :) = (/ 0.1032178_jprb , -0.00293412_jprb , 0.0_jprb /)
661 rhsavi(14, :) = (/ 0.07342832_jprb , -0.001775135_jprb , 0.0_jprb /)
662 rhsavi(15, :) = (/ 0.04649514_jprb , -0.0009165462_jprb , 0.0_jprb /)
663 rhsavi(16, :) = (/ 0.06893519_jprb , -0.0001245402_jprb , 0.0_jprb /)
664 
665 ! water clouds from Lindner & Li (2000)
666 
667 rlilia( 1, :) = (/ 0.88116e-01_jprb,-0.12857e-02_jprb, 0.81658e+00_jprb,-0.39428e+01_jprb, 0.46652e+01_jprb /)
668 rlilia( 2, :) = (/ 0.41307e-03_jprb,-0.59631e-04_jprb, 0.24275e+01_jprb,-0.90838e+01_jprb, 0.96069e+01_jprb /)
669 rlilia( 3, :) = (/-0.57709e-01_jprb, 0.99071e-03_jprb, 0.31118e+01_jprb,-0.95540e+01_jprb, 0.90189e+01_jprb /)
670 rlilia( 4, :) = (/-0.53069e-01_jprb, 0.99992e-03_jprb, 0.28045e+01_jprb,-0.72836e+01_jprb, 0.62573e+01_jprb /)
671 rlilia( 5, :) = (/-0.23627e-01_jprb, 0.55291e-03_jprb, 0.21785e+01_jprb,-0.54664e+01_jprb, 0.47379e+01_jprb /)
672 rlilia( 6, :) = (/ 0.29022e-01_jprb,-0.39657e-03_jprb, 0.14902e+01_jprb,-0.50777e+01_jprb, 0.52170e+01_jprb /)
673 rlilia( 7, :) = (/-0.24901e-01_jprb, 0.16195e-03_jprb, 0.29375e+02_jprb,-0.11437e+02_jprb, 0.12273e+02_jprb /)
674 rlilia( 8, :) = (/-0.14269e+00_jprb, 0.22282e-02_jprb, 0.46478e+01_jprb,-0.16369e+02_jprb, 0.16533e+02_jprb /)
675 rlilia( 9, :) = (/-0.20398e+00_jprb, 0.34708e-02_jprb, 0.52858e+01_jprb,-0.16603e+02_jprb, 0.15392e+02_jprb /)
676 rlilia(10, :) = (/-0.18318e+00_jprb, 0.33080e-02_jprb, 0.46120e+01_jprb,-0.11550e+02_jprb, 0.87086e+01_jprb /)
677 rlilia(11, :) = (/-0.20420e+00_jprb, 0.37167e-02_jprb, 0.48566e+01_jprb,-0.11972e+02_jprb, 0.86344e+01_jprb /)
678 rlilia(12, :) = (/-0.14037e+00_jprb, 0.28058e-02_jprb, 0.34969e+01_jprb,-0.33770e+01_jprb,-0.23541e+01_jprb /)
679 rlilia(13, :) = (/-0.14037e+00_jprb, 0.28058e-02_jprb, 0.34969e+01_jprb,-0.33770e+01_jprb,-0.23541e+01_jprb /)
680 rlilia(14, :) = (/-0.14037e+00_jprb, 0.28058e-02_jprb, 0.34969e+01_jprb,-0.33770e+01_jprb,-0.23541e+01_jprb /)
681 rlilia(15, :) = (/-0.14037e+00_jprb, 0.28058e-02_jprb, 0.34969e+01_jprb,-0.33770e+01_jprb,-0.23541e+01_jprb /)
682 rlilia(16, :) = (/-0.14037e+00_jprb, 0.28058e-02_jprb, 0.34969e+01_jprb,-0.33770e+01_jprb,-0.23541e+01_jprb /)
683 
684 rlilib( 1, :) = (/ 0.76110e+00_jprb, 0.57048e+00_jprb,-0.16721e-01_jprb, 0.27886e-03_jprb /)
685 rlilib( 2, :) = (/ 0.51897e+00_jprb, 0.86352e+00_jprb,-0.29099e-02_jprb, 0.63126e-04_jprb /)
686 rlilib( 3, :) = (/ 0.50964e+00_jprb, 0.76057e+00_jprb,-0.15630e-02_jprb, 0.26245e-04_jprb /)
687 rlilib( 4, :) = (/ 0.55180e+00_jprb, 0.65812e+00_jprb,-0.41768e-02_jprb, 0.64243e-04_jprb /)
688 rlilib( 5, :) = (/ 0.60890e+00_jprb, 0.59525e+00_jprb,-0.74557e-02_jprb, 0.11757e-03_jprb /)
689 rlilib( 6, :) = (/ 0.49787e+00_jprb, 0.74581e+00_jprb,-0.37379e-02_jprb, 0.90555e-04_jprb /)
690 rlilib( 7, :) = (/ 0.46323e-01_jprb, 0.10131e+01_jprb, 0.18692e-01_jprb,-0.21898e-03_jprb /)
691 rlilib( 8, :) = (/-0.43354e-01_jprb, 0.78454e+00_jprb, 0.26449e-01_jprb,-0.36385e-03_jprb /)
692 rlilib( 9, :) = (/-0.27033e-01_jprb, 0.54717e+00_jprb, 0.27903e-01_jprb,-0.41074e-03_jprb /)
693 rlilib(10, :) = (/ 0.11595e+00_jprb, 0.30559e+00_jprb, 0.22727e-01_jprb,-0.36032e-03_jprb /)
694 rlilib(11, :) = (/-0.14764e-01_jprb, 0.28331e+00_jprb, 0.23515e-01_jprb,-0.33662e-03_jprb /)
695 rlilib(12, :) = (/ 0.68056e-02_jprb, 0.59837e-01_jprb, 0.21125e-01_jprb,-0.29740e-03_jprb /)
696 rlilib(13, :) = (/ 0.68056e-02_jprb, 0.59837e-01_jprb, 0.21125e-01_jprb,-0.29740e-03_jprb /)
697 rlilib(14, :) = (/ 0.68056e-02_jprb, 0.59837e-01_jprb, 0.21125e-01_jprb,-0.29740e-03_jprb /)
698 rlilib(15, :) = (/ 0.68056e-02_jprb, 0.59837e-01_jprb, 0.21125e-01_jprb,-0.29740e-03_jprb /)
699 rlilib(16, :) = (/ 0.68056e-02_jprb, 0.59837e-01_jprb, 0.21125e-01_jprb,-0.29740e-03_jprb /)
700 
701 ! ice clouds from Fu & Liou (1993)
702 
703 rfulio( 1, :) = (/ -7.752e-03_jprb, 4.624_jprb, -42.01_jprb /)
704 rfulio( 2, :) = (/ -1.741e-02_jprb, 5.541_jprb, -58.42_jprb /)
705 rfulio( 3, :) = (/ -1.704e-02_jprb, 4.830_jprb, 16.27_jprb /)
706 rfulio( 4, :) = (/ -1.151e-02_jprb, 4.182_jprb, 31.13_jprb /)
707 rfulio( 5, :) = (/ -1.026e-02_jprb, 4.105_jprb, 16.36_jprb /)
708 rfulio( 6, :) = (/ -8.294e-03_jprb, 3.925_jprb, 1.315_jprb /)
709 rfulio( 7, :) = (/ -1.153e-02_jprb, 4.109_jprb, 17.32_jprb /)
710 rfulio( 8, :) = (/ -9.609e-03_jprb, 3.768_jprb, 34.11_jprb /)
711 rfulio( 9, :) = (/ -9.061e-03_jprb, 3.741_jprb, 26.48_jprb /)
712 rfulio(10, :) = (/ -8.441e-03_jprb, 3.715_jprb, 19.48_jprb /)
713 rfulio(11, :) = (/ -8.088e-03_jprb, 3.717_jprb, 17.17_jprb /)
714 rfulio(12, :) = (/ -8.088e-03_jprb, 3.717_jprb, 17.17_jprb /)
715 rfulio(13, :) = (/ -7.770e-03_jprb, 3.734_jprb, 11.85_jprb /)
716 rfulio(14, :) = (/ -6.656e-03_jprb, 3.686_jprb, 0.0_jprb /)
717 rfulio(15, :) = (/ -6.656e-03_jprb, 3.686_jprb, 0.0_jprb /)
718 rfulio(16, :) = (/ -6.656e-03_jprb, 3.686_jprb, 0.0_jprb /)
719 
720 ! ice clouds from Fu et al. (1998)
721 
722 !-- extinction coefficient
723 rfueta( 1, :) = (/ 4.919685e-03_jprb, 2.327741e+00_jprb,-1.390858e+01_jprb /)
724 rfueta( 2, :) = (/ 3.325756e-03_jprb, 2.601360e+00_jprb,-1.909602e+01_jprb /)
725 rfueta( 3, :) = (/ -1.334860e-02_jprb, 4.043808e+00_jprb,-2.171029e+01_jprb /)
726 rfueta( 4, :) = (/ -9.524174e-03_jprb, 3.587742e+00_jprb,-1.068895e+01_jprb /)
727 rfueta( 5, :) = (/ -4.159424e-03_jprb, 3.047325e+00_jprb,-5.061568e+00_jprb /)
728 rfueta( 6, :) = (/ -1.691632e-03_jprb, 2.765756e+00_jprb,-8.331033e+00_jprb /)
729 rfueta( 7, :) = (/ -8.372696e-03_jprb, 3.455018e+00_jprb,-1.516692e+01_jprb /)
730 rfueta( 8, :) = (/ -8.178608e-03_jprb, 3.401245e+00_jprb,-8.812820e+00_jprb /)
731 rfueta( 9, :) = (/ -4.936610e-03_jprb, 3.087764e+00_jprb,-3.884262e+00_jprb /)
732 rfueta(10, :) = (/ -3.034573e-03_jprb, 2.900043e+00_jprb,-1.849911e+00_jprb /)
733 rfueta(11, :) = (/ -3.034573e-03_jprb, 2.900043e+00_jprb,-1.849911e+00_jprb /)
734 rfueta(12, :) = (/ -2.465236e-03_jprb, 2.833187e+00_jprb,-4.227573e-01_jprb /)
735 rfueta(13, :) = (/ -2.308881e-03_jprb, 2.814002e+00_jprb, 1.072211e+00_jprb /)
736 rfueta(14, :) = (/ -2.308881e-03_jprb, 2.814002e+00_jprb, 1.072211e+00_jprb /)
737 rfueta(15, :) = (/ -2.308881e-03_jprb, 2.814002e+00_jprb, 1.072211e+00_jprb /)
738 rfueta(16, :) = (/ -2.308881e-03_jprb, 2.814002e+00_jprb, 1.072211e+00_jprb /)
739 
740 !-- absorption coefficient
741 rfuetb( 1, :) = (/ 8.869787e-01_jprb, 2.118409e-02_jprb,-2.781429e-04_jprb, 1.094562e-06_jprb /)
742 rfuetb( 2, :) = (/ 2.005578e-01_jprb, 2.132614e-02_jprb,-1.751052e-04_jprb, 5.355885e-07_jprb /)
743 rfuetb( 3, :) = (/ 3.003701e-01_jprb, 2.051529e-02_jprb,-1.931684e-04_jprb, 6.583031e-07_jprb /)
744 rfuetb( 4, :) = (/ 9.551440e-01_jprb, 1.309792e-02_jprb,-1.793694e-04_jprb, 7.313392e-07_jprb /)
745 rfuetb( 5, :) = (/ 1.466481e+00_jprb,-2.129226e-03_jprb,-1.361630e-05_jprb, 1.193649e-07_jprb /)
746 rfuetb( 6, :) = (/ 1.195515e+00_jprb, 3.350616e-03_jprb,-5.266996e-05_jprb, 2.233377e-07_jprb /)
747 rfuetb( 7, :) = (/ 5.409536e-01_jprb, 1.949649e-02_jprb,-2.050908e-04_jprb, 7.364680e-07_jprb /)
748 rfuetb( 8, :) = (/ 5.874323e-01_jprb, 1.876628e-02_jprb,-2.045834e-04_jprb, 7.510080e-07_jprb /)
749 rfuetb( 9, :) = (/ 7.152274e-01_jprb, 1.621734e-02_jprb,-1.868544e-04_jprb, 7.078738e-07_jprb /)
750 rfuetb(10, :) = (/ 8.862434e-01_jprb, 1.226538e-02_jprb,-1.523076e-04_jprb, 6.000892e-07_jprb /)
751 rfuetb(11, :) = (/ 8.862434e-01_jprb, 1.226538e-02_jprb,-1.523076e-04_jprb, 6.000892e-07_jprb /)
752 rfuetb(12, :) = (/ 7.428957e-01_jprb, 1.279601e-02_jprb,-1.391803e-04_jprb, 5.180104e-07_jprb /)
753 rfuetb(13, :) = (/ 4.346482e-01_jprb, 1.721457e-02_jprb,-1.623227e-04_jprb, 5.561523e-07_jprb /)
754 rfuetb(14, :) = (/ 4.346482e-01_jprb, 1.721457e-02_jprb,-1.623227e-04_jprb, 5.561523e-07_jprb /)
755 rfuetb(15, :) = (/ 4.346482e-01_jprb, 1.721457e-02_jprb,-1.623227e-04_jprb, 5.561523e-07_jprb /)
756 rfuetb(16, :) = (/ 4.346482e-01_jprb, 1.721457e-02_jprb,-1.623227e-04_jprb, 5.561523e-07_jprb /)
757 
758 !-- asymmetry factor
759 rfuetc( 1, :) = (/ 4.949276e-01_jprb, 1.186174e-02_jprb,-1.267629e-04_jprb, 4.603574e-07_jprb /)
760 rfuetc( 2, :) = (/ 6.891414e-01_jprb, 6.192281e-03_jprb,-6.459514e-05_jprb, 2.436963e-07_jprb /)
761 rfuetc( 3, :) = (/ 7.260484e-01_jprb, 2.664334e-03_jprb,-1.251136e-05_jprb, 2.243377e-08_jprb /)
762 rfuetc( 4, :) = (/ 7.363466e-01_jprb, 4.798266e-03_jprb,-4.413293e-05_jprb, 1.525774e-07_jprb /)
763 rfuetc( 5, :) = (/ 7.984021e-01_jprb, 3.977117e-03_jprb,-4.471984e-05_jprb, 1.694919e-07_jprb /)
764 rfuetc( 6, :) = (/ 8.663385e-01_jprb, 2.797934e-03_jprb,-3.187011e-05_jprb, 1.217209e-07_jprb /)
765 rfuetc( 7, :) = (/ 8.906280e-01_jprb, 1.903269e-03_jprb,-1.733552e-05_jprb, 5.855071e-08_jprb /)
766 rfuetc( 8, :) = (/ 8.609604e-01_jprb, 2.200445e-03_jprb,-1.748105e-05_jprb, 5.176616e-08_jprb /)
767 rfuetc( 9, :) = (/ 8.522816e-01_jprb, 2.523627e-03_jprb,-2.149196e-05_jprb, 6.685067e-08_jprb /)
768 rfuetc(10, :) = (/ 8.741665e-01_jprb, 2.455409e-03_jprb,-2.456935e-05_jprb, 8.641223e-08_jprb /)
769 rfuetc(11, :) = (/ 8.741665e-01_jprb, 2.455409e-03_jprb,-2.456935e-05_jprb, 8.641223e-08_jprb /)
770 rfuetc(12, :) = (/ 8.472918e-01_jprb, 2.559953e-03_jprb,-2.182660e-05_jprb, 6.879977e-08_jprb /)
771 rfuetc(13, :) = (/ 7.962716e-01_jprb, 3.003488e-03_jprb,-2.082376e-05_jprb, 5.366545e-08_jprb /)
772 
773 ! ----------------------------------------------------------------
774 ! Ebert-Curry
775 
776 ! LW : spectrally defined for EC-OPE
777 
778 rebcui = (/ 1.136_jprb, 1.338_jprb, 1.166_jprb, 1.166_jprb, 1.118_jprb, &
779  & 0.600_jprb /)
780 rebcuj = (/ 0.0036_jprb, 0.0003_jprb, 0.0016_jprb, 0.0016_jprb, 0.0020_jprb, &
781  & 0.0068_jprb /)
782 
783 ! LW : spectrally defined for RRTM
784 ! mass-absorption coefficients for vertical path: no diffusivity factor
785 
786 rebcug = (/ 0.718_jprb, 0.726_jprb, 1.136_jprb, 1.320_jprb, 1.505_jprb, &
787  & 1.290_jprb, 0.911_jprb, 0.949_jprb, 1.021_jprb, 1.193_jprb, &
788  & 1.279_jprb, 0.626_jprb, 0.647_jprb, 0.668_jprb, 0.690_jprb, &
789  & 0.690_jprb /)
790 
791 rebcuh = (/ 0.0069_jprb, 0.0060_jprb, 0.0024_jprb, 0.0004_jprb,-0.0016_jprb, &
792  & 0.0003_jprb, 0.0043_jprb, 0.0038_jprb, 0.0030_jprb, 0.0013_jprb, &
793  & 0.0005_jprb, 0.0054_jprb, 0.0052_jprb, 0.0050_jprb, 0.0048_jprb, &
794  & 0.0048_jprb /)
795 
796 ! Sun-Shine
797 
798 rsushfa = (/ 1.047_jprb, -0.913e-04_jprb, 0.203e-03_jprb, -0.106e-04_jprb /)
799 
800 ! ------------------------------------------------------------------
801 
802 !* 2.
803 
804 !* Liquid/Solid water transition
805 
806 rtiw= 263._jprb
807 rriw= 20._jprb
808 
809 ! Ice particle Effective Radius as a function of LWC
810 
811 reffia= 40._jprb
812 reffib= 0._jprb
813 
814 ! Sun-Shine
815 
816 rsushc= 0.0306_jprb
817 rsushd= 0.2548_jprb
818 
819 ! ------------------------------------------------------------------
820 
821 ! SW : absorption coefficients
822 
823 IF (ksw == 2) THEN
824  DO jnu=1,ksw
825  raswca(jnu)=zaswca2(jnu)*1.e-02_jprb
826  raswcb(jnu)=zaswcb2(jnu)
827  raswcc(jnu)=zaswcc2(jnu)
828  raswcd(jnu)=zaswcd2(jnu)
829  raswce(jnu)=zaswce2(jnu)
830  raswcf(jnu)=zaswcf2(jnu)*1.e-03_jprb
831 
832  rebcua(jnu)=zebcua2(jnu)
833  rebcub(jnu)=zebcub2(jnu)
834  rebcuc(jnu)=zebcuc2(jnu)
835  rebcud(jnu)=zebcud2(jnu)
836  rebcue(jnu)=zebcue2(jnu)
837  rebcuf(jnu)=zebcuf2(jnu)
838 
839  ryfwca(jnu)=zyfwca2(jnu)
840  ryfwcb(jnu)=zyfwcb2(jnu)
841  ryfwcc(jnu)=zyfwcc2(jnu)
842  ryfwcd(jnu)=zyfwcd2(jnu)
843  ryfwce(jnu)=zyfwce2(jnu)
844  ryfwcf(jnu)=zyfwcf2(jnu)
845 
846  rsushe(jnu)=zsushe2(jnu)*1.e-02_jprb
847  rsushf(jnu)=zsushf2(jnu)*1.e-02_jprb
848  rsushh(jnu)=zsushh2(jnu)
849  rsushk(jnu)=zsushk2(jnu)*1.e-01_jprb
850  rsusha(jnu)=zsusha2(jnu)*1.e-03_jprb
851  rsushg(jnu)=zsushg2(jnu)*1.e-01_jprb
852 
853  rflaa0(jnu)=zflaa02(jnu)
854  rflaa1(jnu)=zflaa12(jnu)
855  rflbb0(jnu)=zflbb02(jnu)
856  rflbb1(jnu)=zflbb12(jnu)
857  rflbb2(jnu)=zflbb22(jnu)
858  rflbb3(jnu)=zflbb32(jnu)
859  rflcc0(jnu)=zflcc02(jnu)
860  rflcc1(jnu)=zflcc12(jnu)
861  rflcc2(jnu)=zflcc22(jnu)
862  rflcc3(jnu)=zflcc32(jnu)
863  rfldd0(jnu)=zfldd02(jnu)
864  rfldd1(jnu)=zfldd12(jnu)
865  rfldd2(jnu)=zfldd22(jnu)
866  rfldd3(jnu)=zfldd32(jnu)
867 
868 ! Fu (1996)'s optical properties are not available for 2 spectral intervals
869 ! RFUAA0(JNU)=ZFUAA02(JNU)
870 ! RFUAA1(JNU)=ZFUAA12(JNU)
871 ! RFUBB0(JNU)=ZFUBB02(JNU)
872 ! RFUBB1(JNU)=ZFUBB12(JNU)
873 ! RFUBB2(JNU)=ZFUBB22(JNU)
874 ! RFUBB3(JNU)=ZFUBB32(JNU)
875 ! RFUCC0(JNU)=ZFUCC02(JNU)
876 ! RFUCC1(JNU)=ZFUCC12(JNU)
877 ! RFUCC2(JNU)=ZFUCC22(JNU)
878 ! RFUCC3(JNU)=ZFUCC32(JNU)
879 
880  ENDDO
881 ELSEIF (ksw == 4) THEN
882  DO jnu=1,ksw
883  raswca(jnu)=zaswca4(jnu)*1.e-02_jprb
884  raswcb(jnu)=zaswcb4(jnu)
885  raswcc(jnu)=zaswcc4(jnu)
886  raswcd(jnu)=zaswcd4(jnu)
887  raswce(jnu)=zaswce4(jnu)
888  raswcf(jnu)=zaswcf4(jnu)*1.e-03_jprb
889 
890  rebcua(jnu)=zebcua4(jnu)
891  rebcub(jnu)=zebcub4(jnu)
892  rebcuc(jnu)=zebcuc4(jnu)
893  rebcud(jnu)=zebcud4(jnu)
894  rebcue(jnu)=zebcue4(jnu)
895  rebcuf(jnu)=zebcuf4(jnu)
896 
897  ryfwca(jnu)=zyfwca4(jnu)
898  ryfwcb(jnu)=zyfwcb4(jnu)
899  ryfwcc(jnu)=zyfwcc4(jnu)
900  ryfwcd(jnu)=zyfwcd4(jnu)
901  ryfwce(jnu)=zyfwce4(jnu)
902  ryfwcf(jnu)=zyfwcf4(jnu)
903 
904  rsushe(jnu)=zsushe4(jnu)*1.e-02_jprb
905  rsushf(jnu)=zsushf4(jnu)*1.e-02_jprb
906  rsushh(jnu)=zsushh4(jnu)
907  rsushk(jnu)=zsushk4(jnu)*1.e-01_jprb
908  rsusha(jnu)=zsusha4(jnu)*1.e-03_jprb
909  rsushg(jnu)=zsushg4(jnu)*1.e-01_jprb
910 
911  rflaa0(jnu)=zflaa04(jnu)
912  rflaa1(jnu)=zflaa14(jnu)
913  rflbb0(jnu)=zflbb04(jnu)
914  rflbb1(jnu)=zflbb14(jnu)
915  rflbb2(jnu)=zflbb24(jnu)
916  rflbb3(jnu)=zflbb34(jnu)
917  rflcc0(jnu)=zflcc04(jnu)
918  rflcc1(jnu)=zflcc14(jnu)
919  rflcc2(jnu)=zflcc24(jnu)
920  rflcc3(jnu)=zflcc34(jnu)
921  rfldd0(jnu)=zfldd04(jnu)
922  rfldd1(jnu)=zfldd14(jnu)
923  rfldd2(jnu)=zfldd24(jnu)
924  rfldd3(jnu)=zfldd34(jnu)
925 
926  rfuaa0(jnu)=zfuaa04(jnu)
927  rfuaa1(jnu)=zfuaa14(jnu)
928  rfubb0(jnu)=zfubb04(jnu)
929  rfubb1(jnu)=zfubb14(jnu)
930  rfubb2(jnu)=zfubb24(jnu)
931  rfubb3(jnu)=zfubb34(jnu)
932  rfucc0(jnu)=zfucc04(jnu)
933  rfucc1(jnu)=zfucc14(jnu)
934  rfucc2(jnu)=zfucc24(jnu)
935  rfucc3(jnu)=zfucc34(jnu)
936 
937  ENDDO
938 ELSEIF (ksw == 6) THEN
939  DO jnu=1,ksw
940  raswca(jnu)=zaswca6(jnu)*1.e-02_jprb
941  raswcb(jnu)=zaswcb6(jnu)
942  raswcc(jnu)=zaswcc6(jnu)
943  raswcd(jnu)=zaswcd6(jnu)
944  raswce(jnu)=zaswce6(jnu)
945  raswcf(jnu)=zaswcf6(jnu)*1.e-03_jprb
946 
947  rebcua(jnu)=zebcua6(jnu)
948  rebcub(jnu)=zebcub6(jnu)
949  rebcuc(jnu)=zebcuc6(jnu)
950  rebcud(jnu)=zebcud6(jnu)
951  rebcue(jnu)=zebcue6(jnu)
952  rebcuf(jnu)=zebcuf6(jnu)
953 
954  ryfwca(jnu)=zyfwca6(jnu)
955  ryfwcb(jnu)=zyfwcb6(jnu)
956  ryfwcc(jnu)=zyfwcc6(jnu)
957  ryfwcd(jnu)=zyfwcd6(jnu)
958  ryfwce(jnu)=zyfwce6(jnu)
959  ryfwcf(jnu)=zyfwcf6(jnu)
960 
961  rsushe(jnu)=zsushe6(jnu)*1.e-02_jprb
962  rsushf(jnu)=zsushf6(jnu)*1.e-02_jprb
963  rsushh(jnu)=zsushh6(jnu)
964  rsushk(jnu)=zsushk6(jnu)*1.e-01_jprb
965  rsusha(jnu)=zsusha6(jnu)*1.e-03_jprb
966  rsushg(jnu)=zsushg6(jnu)*1.e-01_jprb
967 
968  rflaa0(jnu)=zflaa06(jnu)
969  rflaa1(jnu)=zflaa16(jnu)
970  rflbb0(jnu)=zflbb06(jnu)
971  rflbb1(jnu)=zflbb16(jnu)
972  rflbb2(jnu)=zflbb26(jnu)
973  rflbb3(jnu)=zflbb36(jnu)
974  rflcc0(jnu)=zflcc06(jnu)
975  rflcc1(jnu)=zflcc16(jnu)
976  rflcc2(jnu)=zflcc26(jnu)
977  rflcc3(jnu)=zflcc36(jnu)
978  rfldd0(jnu)=zfldd06(jnu)
979  rfldd1(jnu)=zfldd16(jnu)
980  rfldd2(jnu)=zfldd26(jnu)
981  rfldd3(jnu)=zfldd36(jnu)
982 
983  rfuaa0(jnu)=zfuaa06(jnu)
984  rfuaa1(jnu)=zfuaa16(jnu)
985  rfubb0(jnu)=zfubb06(jnu)
986  rfubb1(jnu)=zfubb16(jnu)
987  rfubb2(jnu)=zfubb26(jnu)
988  rfubb3(jnu)=zfubb36(jnu)
989  rfucc0(jnu)=zfucc06(jnu)
990  rfucc1(jnu)=zfucc16(jnu)
991  rfucc2(jnu)=zfucc26(jnu)
992  rfucc3(jnu)=zfucc36(jnu)
993 
994  ENDDO
995 ELSEIF (ktsw == 14 .AND. lsrtm) THEN
996  print *,'SUCLOPN: 14-SPECTRAL INTERVALS --> RRTM_SW'
997 ! DO JNU=1,KSW
998 ! RASWCA(JNU)=ZASWCA14(JNU)*1.E-02_JPRB
999 ! RASWCB(JNU)=ZASWCB14(JNU)
1000 ! RASWCC(JNU)=ZASWCC14(JNU)
1001 ! RASWCD(JNU)=ZASWCD14(JNU)
1002 ! RASWCE(JNU)=ZASWCE14(JNU)
1003 ! RASWCF(JNU)=ZASWCF14(JNU)*1.E-03_JPRB
1004 !
1005 ! REBCUA(JNU)=ZEBCUA14(JNU)
1006 ! REBCUB(JNU)=ZEBCUB14(JNU)
1007 ! REBCUC(JNU)=ZEBCUC14(JNU)
1008 ! REBCUD(JNU)=ZEBCUD14(JNU)
1009 ! REBCUE(JNU)=ZEBCUE14(JNU)
1010 ! REBCUF(JNU)=ZEBCUF14(JNU)
1011 !
1012 ! RYFWCA(JNU)=ZYFWCA14(JNU)
1013 ! RYFWCB(JNU)=ZYFWCB14(JNU)
1014 ! RYFWCC(JNU)=ZYFWCC14(JNU)
1015 ! RYFWCD(JNU)=ZYFWCD14(JNU)
1016 ! RYFWCE(JNU)=ZYFWCE14(JNU)
1017 ! RYFWCF(JNU)=ZYFWCF14(JNU)
1018 !
1019 ! RSUSHE(JNU)=ZSUSHE14(JNU)*1.E-02_JPRB
1020 ! RSUSHF(JNU)=ZSUSHF14(JNU)*1.E-02_JPRB
1021 ! RSUSHH(JNU)=ZSUSHH14(JNU)
1022 ! RSUSHK(JNU)=ZSUSHK14(JNU)*1.E-01_JPRB
1023 ! RSUSHA(JNU)=ZSUSHA14(JNU)*1.E-03_JPRB
1024 ! RSUSHG(JNU)=ZSUSHG14(JNU)*1.E-01_JPRB
1025 !
1026 ! RFLAA0(JNU)=ZFLAA014(JNU)
1027 ! RFLAA1(JNU)=ZFLAA114(JNU)
1028 ! RFLBB0(JNU)=ZFLBB014(JNU)
1029 ! RFLBB1(JNU)=ZFLBB114(JNU)
1030 ! RFLBB2(JNU)=ZFLBB214(JNU)
1031 ! RFLBB3(JNU)=ZFLBB314(JNU)
1032 ! RFLCC0(JNU)=ZFLCC014(JNU)
1033 ! RFLCC1(JNU)=ZFLCC114(JNU)
1034 ! RFLCC2(JNU)=ZFLCC214(JNU)
1035 ! RFLCC3(JNU)=ZFLCC314(JNU)
1036 ! RFLDD0(JNU)=ZFLDD014(JNU)
1037 ! RFLDD1(JNU)=ZFLDD114(JNU)
1038 ! RFLDD2(JNU)=ZFLDD214(JNU)
1039 ! RFLDD3(JNU)=ZFLDD314(JNU)
1040 !
1041 ! RFUAA0(JNU)=ZFUAA014(JNU)
1042 ! RFUAA1(JNU)=ZFUAA114(JNU)
1043 ! RFUBB0(JNU)=ZFUBB014(JNU)
1044 ! RFUBB1(JNU)=ZFUBB114(JNU)
1045 ! RFUBB2(JNU)=ZFUBB214(JNU)
1046 ! RFUBB3(JNU)=ZFUBB314(JNU)
1047 ! RFUCC0(JNU)=ZFUCC014(JNU)
1048 ! RFUCC1(JNU)=ZFUCC114(JNU)
1049 ! RFUCC2(JNU)=ZFUCC214(JNU)
1050 ! RFUCC3(JNU)=ZFUCC314(JNU)
1051 !! RFUDD0(JNU)=ZFUDD014(JNU)
1052 !! RFUDD1(JNU)=ZFUDD114(JNU)
1053 !! RFUDD2(JNU)=ZFUDD214(JNU)
1054 !! RFUDD3(JNU)=ZFUDD314(JNU)
1055 !
1056 ! PRINT *,'SUCLOPN: 14-SPECTRAL INTERVALS --> RRTM_SW'
1057 ! ENDDO
1058 ELSE
1059  CALL abor1('SUCLOPN: WRONG SW SPECTRAL RESOLUTION')
1060 ENDIF
1061 ! ------------------------------------------------------------------
1062 
1063 !* 2. CLOUD OVERLAP PARAMETERS
1064 ! ------------------------
1065 
1066 zaovlp = (/ -2.250e-05_jprb,-7.316e-06_jprb,-1.966e-05_jprb /)
1067 zbovlp = (/ 0.7865_jprb , 0.8186_jprb , 0.8900_jprb /)
1068 
1069 IF (klev == 19) THEN
1070  raovlp=zaovlp(1)
1071  rbovlp=zbovlp(1)
1072 ELSEIF (klev == 31) THEN
1073  raovlp=zaovlp(2)
1074  rbovlp=zbovlp(2)
1075 ELSEIF (klev == 60) THEN
1076  raovlp=zaovlp(3)
1077  rbovlp=zbovlp(3)
1078 ELSE
1079  raovlp=zaovlp(3)
1080  rbovlp=zbovlp(3)
1081 ENDIF
1082 
1083 ! ------------------------------------------------------------------
1084 
1085 IF (lhook) CALL dr_hook('SUCLOPN',1,zhook_handle)
1086 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::rfuetc
real(kind=jprb), dimension(16, 4) rfuetc
Definition: yoesw.F90:57
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
abor1
subroutine abor1(CDTEXT)
Definition: abor1.F90:2
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::rfuetb
real(kind=jprb), dimension(16, 4) rfuetb
Definition: yoesw.F90:57
yoerad::lsrtm
logical lsrtm
Definition: yoerad.F90:53
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
parkind1::jprb
integer, parameter jprb
Definition: parkind1.F90:31
yoesw::rlilia
real(kind=jprb), dimension(16, 5) rlilia
Definition: yoesw.F90:60
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
yomhook::lhook
logical lhook
Definition: yomhook.F90:12
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
yoesw::rlilib
real(kind=jprb), dimension(16, 4) rlilib
Definition: yoesw.F90:60
yoerad::raovlp
real(kind=jprb) raovlp
Definition: yoerad.F90:68
yoesw::rflaa1
real(kind=jprb), dimension(6) rflaa1
Definition: yoesw.F90:63
yomhook::dr_hook
subroutine dr_hook(CDNAME, KSWITCH, PKEY)
Definition: yomhook.F90:17
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
parkind1::jpim
integer, parameter jpim
Definition: parkind1.F90:13
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