source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phymar/suclopn.F90 @ 4622

Last change on this file since 4622 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR


Integration of MAR physics

File size: 31.7 KB
Line 
1SUBROUTINE 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
58USE 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
71USE YOERAD   , ONLY : RAOVLP   ,RBOVLP
72
73IMPLICIT NONE
74
75
76!     DUMMY INTEGER SCALARS
77INTEGER_M :: KSW
78INTEGER_M :: KTSW
79INTEGER_M :: KLEV
80INTEGER_M :: JC, JNU
81
82
83
84
85!     -----------------------------------------------------------------
86REAL_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
101REAL_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
116REAL_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 
132REAL_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
149ZEBCUA2 = (/ 3.448E-04_JPRB , 3.448E-04_JPRB /)
150ZEBCUB2 = (/ 2.431_JPRB     , 2.431_JPRB     /)
151ZEBCUC2 = (/ 0.00001_JPRB   , 0.024366_JPRB  /)
152ZEBCUD2 = (/ _ZERO_         , 2.487E-04_JPRB /)
153ZEBCUE2 = (/ 0.7661_JPRB    , 0.7866_JPRB    /)
154ZEBCUF2 = (/ 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
166ZYFWCA2 = (/ _ZERO_         , _ZERO_         /)
167ZYFWCB2 = (/ 1.5_JPRB       , 1.5_JPRB       /)
168ZYFWCC2 = (/ 0.9999_JPRB    , 0.9988_JPRB    /)
169ZYFWCD2 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB /)
170ZYFWCE2 = (/ _HALF_         , 0.05_JPRB      /)
171ZYFWCF2 = (/ 0.865_JPRB     , 0.910_JPRB     /)
172
173!* Water cloud properties - from Slingo (1989)
174
175ZASWCA2 = (/ 2.817_JPRB     , 2.4584_JPRB    /)
176ZASWCB2 = (/ 1.305_JPRB     , 1.3994_JPRB    /)
177ZASWCC2 = (/-5.62E-08_JPRB  , 1.111E-02_JPRB /)
178ZASWCD2 = (/ 1.63E-07_JPRB  , 8.613E-04_JPRB /)
179ZASWCE2 = (/ 0.829_JPRB     , 0.7819_JPRB    /)
180ZASWCF2 = (/ 2.482_JPRB     , 5.0150_JPRB    /)
181
182!* Ice cloud properties - from Sun and Shine (1995)
183
184ZSUSHE2 = (/ _ZERO_         , 8.6822_JPRB    /)
185ZSUSHF2 = (/ _ZERO_         , 9.6277_JPRB    /)
186ZSUSHH2 = (/ 0.8522_JPRB    , 0.8819_JPRB    /)
187ZSUSHK2 = (/ 0.1620_JPRB    , 0.1630_JPRB    /)
188ZSUSHA2 = (/ _ZERO_         , 23.204_JPRB    /)
189ZSUSHG2 = (/ 0.3270_JPRB    , 0.4180_JPRB    /)
190
191!* Ice cloud properties - from Fu and Liou (1993)
192
193ZFLAA02 = (/-6.656E-3_JPRB  ,-6.656E-3_JPRB  /)
194ZFLAA12 = (/ 3.686_JPRB     , 3.686_JPRB     /)
195ZFLBB02 = (/ .10998E-5_JPRB , .21136E-1_JPRB /)
196ZFLBB12 = (/-.26101E-7_JPRB , .39150E-3_JPRB /)
197ZFLBB22 = (/ .18096E-8_JPRB ,-.20740E-6_JPRB /)
198ZFLBB32 = (/-.47387E-11_JPRB,-.28829E-8_JPRB /)
199ZFLCC02 = (/ .22110E+1_JPRB , .22498E+1_JPRB /)
200ZFLCC12 = (/-.10398E-2_JPRB , .23656E-3_JPRB /)
201ZFLCC22 = (/ .65199E-4_JPRB , .51948E-4_JPRB /)
202ZFLCC32 = (/-.34498E-6_JPRB ,-.29768E-6_JPRB /)
203ZFLDD02 = (/ .12495_JPRB    , .11716_JPRB    /)
204ZFLDD12 = (/-.43582E-3_JPRB ,-.45208E-3_JPRB /)
205ZFLDD22 = (/ .14092E-4_JPRB , .12772E-4_JPRB /)
206ZFLDD32 = (/-.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
240ZEBCUA4 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /)
241ZEBCUB4 = (/ 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     /)
242ZEBCUC4 = (/ 0.00001_JPRB   , 0.00011_JPRB   , 0.0197796_JPRB , 0.46658_JPRB   /)
243ZEBCUD4 = (/ _ZERO_         , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /)
244ZEBCUE4 = (/ 0.7661_JPRB    , 0.7730_JPRB    , 0.795653_JPRB  , 0.9595_JPRB    /)
245ZEBCUF4 = (/ 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
249ZYFWCA4 = (/ _ZERO_         , _ZERO_         , _ZERO_        , _ZERO_          /)
250ZYFWCB4 = (/ 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB      , 1.5_JPRB        /)
251ZYFWCC4 = (/ 0.9999_JPRB    , 0.9988_JPRB    , 0.9988_JPRB   , 0.9988_JPRB     /)
252ZYFWCD4 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB  /)
253ZYFWCE4 = (/ _HALF_         , 0.05_JPRB      , 0.05_JPRB     , 0.05_JPRB       /)
254ZYFWCF4 = (/ 0.865_JPRB     , 0.910_JPRB     , 0.910_JPRB    , 0.910_JPRB      /)
255
256!* Water cloud properties - from Slingo (1989)
257
258ZASWCA4 = (/ 2.817_JPRB     , 2.682_JPRB     , 2.264_JPRB    , 1.281_JPRB      /)
259ZASWCB4 = (/ 1.305_JPRB     , 1.346_JPRB     , 1.455_JPRB    , 1.641_JPRB      /)
260ZASWCC4 = (/-5.62E-08_JPRB  ,-6.94E-06_JPRB  , 4.64E-04_JPRB , 2.01E-01_JPRB   /)
261ZASWCD4 = (/ 1.63E-07_JPRB  , 2.35E-05_JPRB  , 1.24E-03_JPRB , 7.56E-03_JPRB   /)
262ZASWCE4 = (/ 0.829_JPRB     , 0.794_JPRB     , 0.754_JPRB    , 0.826_JPRB      /)
263ZASWCF4 = (/ 2.482_JPRB     , 4.226_JPRB     , 6.560_JPRB    , 4.353_JPRB      /)
264
265!* Ice cloud properties - from Sun and Shine (1995)
266
267ZSUSHE4 = (/ _ZERO_       , 7.2471E-02_JPRB , 17.5933_JPRB  , 48.7166_JPRB  /)
268ZSUSHF4 = (/ _ZERO_       , 4.01511_JPRB    , 21.1249_JPRB  , 1.20890_JPRB  /)
269ZSUSHH4 = (/ 0.8522_JPRB  , 0.85841_JPRB    , 0.90778_JPRB  ,0.982046_JPRB  /)
270ZSUSHK4 = (/ 0.1620_JPRB  , 0.160048_JPRB   , 0.188521_JPRB ,0.0411446_JPRB /)
271ZSUSHA4 = (/ _ZERO_       , 0.273455_JPRB   , 41.7675_JPRB  , 161.104_JPRB  /)
272ZSUSHG4 = (/ 0.3270_JPRB  , 0.343668_JPRB   , 0.526192_JPRB , 0.574040_JPRB /)
273
274!* Ice cloud properties - from Fu and Liou (1993)
275
276ZFLAA04 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /)
277ZFLAA14 = (/ 3.686_JPRB     , 3.686_JPRB    , 3.686_JPRB    , 3.686_JPRB     /)
278ZFLBB04 = (/ .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /)
279ZFLBB14 = (/-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /)
280ZFLBB24 = (/ .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /)
281ZFLBB34 = (/-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /)
282ZFLCC04 = (/ .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /)
283ZFLCC14 = (/-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /)
284ZFLCC24 = (/ .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /)
285ZFLCC34 = (/-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /)
286ZFLDD04 = (/ .12495_JPRB    , .12363_JPRB   , .11983_JPRB   , .21834E-1_JPRB /)
287ZFLDD14 = (/-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /)
288ZFLDD24 = (/ .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /)
289ZFLDD34 = (/-.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
293ZFUAA04 = (/-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /)
294ZFUAA14 = (/ 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /)
295ZFUBB04 = (/-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /)
296ZFUBB14 = (/ 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /)
297ZFUBB24 = (/-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /)
298ZFUBB34 = (/ 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /)
299ZFUCC04 = (/ 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /)
300ZFUCC14 = (/ 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /)
301ZFUCC24 = (/-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /)
302ZFUCC34 = (/-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
313ZEBCUA6 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /)
314ZEBCUB6 = (/ 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     /)
315ZEBCUC6 = (/ 0.00001_JPRB   , 0.00001_JPRB   , 0.00001_JPRB   , 0.00011_JPRB   , 0.0197796_JPRB , 0.46658_JPRB   /)
316ZEBCUD6 = (/ _ZERO_         , _ZERO_         , _ZERO_         , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /)
317ZEBCUE6 = (/ 0.7661_JPRB    , 0.7661_JPRB    , 0.7661_JPRB    , 0.7730_JPRB    , 0.795653_JPRB  , 0.9595_JPRB    /)
318ZEBCUF6 = (/ 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
322ZYFWCA6 = (/ _ZERO_         , _ZERO_         , _ZERO_         , _ZERO_         , _ZERO_        , _ZERO_          /)
323ZYFWCB6 = (/ 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB      , 1.5_JPRB        /)
324ZYFWCC6 = (/ 0.9999_JPRB    , 0.9999_JPRB    , 0.9999_JPRB    , 0.9988_JPRB    , 0.9988_JPRB   , 0.9988_JPRB     /)
325ZYFWCD6 = (/ 5.000E-04_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB  /)
326ZYFWCE6 = (/ _HALF_         , _HALF_         , _HALF_         , 0.05_JPRB      , 0.05_JPRB     , 0.05_JPRB       /)
327ZYFWCF6 = (/ 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
331ZASWCA6 = (/ 2.817_JPRB     , 2.817_JPRB     , 2.817_JPRB     , 2.682_JPRB     , 2.264_JPRB    , 1.281_JPRB      /)
332ZASWCB6 = (/ 1.305_JPRB     , 1.305_JPRB     , 1.305_JPRB     , 1.346_JPRB     , 1.455_JPRB    , 1.641_JPRB      /)
333ZASWCC6 = (/-5.62E-08_JPRB  ,-5.62E-08_JPRB  ,-5.62E-08_JPRB  ,-6.94E-06_JPRB  , 4.64E-04_JPRB , 2.01E-01_JPRB   /)
334ZASWCD6 = (/ 1.63E-07_JPRB  , 1.63E-07_JPRB  , 1.63E-07_JPRB  , 2.35E-05_JPRB  , 1.24E-03_JPRB , 7.56E-03_JPRB   /)
335ZASWCE6 = (/ 0.829_JPRB     , 0.829_JPRB     , 0.829_JPRB     , 0.794_JPRB     , 0.754_JPRB    , 0.826_JPRB      /)
336ZASWCF6 = (/ 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
340ZSUSHE6 = (/ _ZERO_       , _ZERO_       , _ZERO_       , 7.2471E-02_JPRB , 17.5933_JPRB  , 48.7166_JPRB  /)
341ZSUSHF6 = (/ _ZERO_       , _ZERO_       , _ZERO_       , 4.01511_JPRB    , 21.1249_JPRB  , 1.20890_JPRB  /)
342ZSUSHH6 = (/ 0.8522_JPRB  , 0.8522_JPRB  , 0.8522_JPRB  , 0.85841_JPRB    , 0.90778_JPRB  ,0.982046_JPRB  /)
343ZSUSHK6 = (/ 0.1620_JPRB  , 0.1620_JPRB  , 0.1620_JPRB  , 0.160048_JPRB   , 0.188521_JPRB ,0.0411446_JPRB /)
344ZSUSHA6 = (/ _ZERO_       , _ZERO_       , _ZERO_       , 0.273455_JPRB   , 41.7675_JPRB  , 161.104_JPRB  /)
345ZSUSHG6 = (/ 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
349ZFLAA06 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /)
350ZFLAA16 = (/ 3.686_JPRB     , 3.686_JPRB     , 3.686_JPRB     , 3.686_JPRB    , 3.686_JPRB    , 3.686_JPRB     /)
351ZFLBB06 = (/ .10998E-5_JPRB , .10998E-5_JPRB , .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /)
352ZFLBB16 = (/-.26101E-7_JPRB ,-.26101E-7_JPRB ,-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /)
353ZFLBB26 = (/ .18096E-8_JPRB , .18096E-8_JPRB , .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /)
354ZFLBB36 = (/-.47387E-11_JPRB,-.47387E-11_JPRB,-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /)
355ZFLCC06 = (/ .22110E+1_JPRB , .22110E+1_JPRB , .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /)
356ZFLCC16 = (/-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /)
357ZFLCC26 = (/ .65199E-4_JPRB , .65199E-4_JPRB , .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /)
358ZFLCC36 = (/-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /)
359ZFLDD06 = (/ .12495_JPRB    , .12495_JPRB    , .12495_JPRB    , .12363_JPRB   , .11983_JPRB   , .21834E-1_JPRB /)
360ZFLDD16 = (/-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /)
361ZFLDD26 = (/ .14092E-4_JPRB , .14092E-4_JPRB , .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /)
362ZFLDD36 = (/-.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
366ZFUAA06 = (/-1.30817E-04_JPRB,-1.30817E-04_JPRB,-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /)
367ZFUAA16 = (/ 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /)
368ZFUBB06 = (/-1.55357E-08_JPRB,-1.55357E-08_JPRB,-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /)
369ZFUBB16 = (/ 1.95793E-07_JPRB, 1.95793E-07_JPRB, 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /)
370ZFUBB26 = (/-2.31234E-10_JPRB,-2.31234E-10_JPRB,-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /)
371ZFUBB36 = (/ 1.12247E-12_JPRB, 1.12247E-12_JPRB, 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /)
372ZFUCC06 = (/ 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /)
373ZFUCC16 = (/ 9.10564E-04_JPRB, 9.10564E-04_JPRB, 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /)
374ZFUCC26 = (/-4.62479E-07_JPRB,-4.62479E-07_JPRB,-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /)
375ZFUCC36 = (/-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
383RHSAVI( 1, :) = (/&
384    &0.1651082_JPRB  , -0.003494839_JPRB  , _ZERO_              /)
385RHSAVI( 2, :) = (/&
386    &0.327820597_JPRB, -0.0219634383_JPRB , 0.000506783898_JPRB /)
387RHSAVI( 3, :) = (/&
388    &0.504805453_JPRB, -0.0478602354_JPRB , 0.00141521102_JPRB  /)
389RHSAVI( 4, :) = (/&
390    &0.513169093_JPRB, -0.0508960145_JPRB , 0.00155498711_JPRB  /)
391RHSAVI( 5, :) = (/&
392    &0.448042082_JPRB, -0.0431857592_JPRB , 0.00130848978_JPRB  /)
393RHSAVI( 6, :) = (/&
394    &0.249547237_JPRB, -0.0185273835_JPRB , 0.00050361258_JPRB  /)
395RHSAVI( 7, :) = (/&
396    &0.135486796_JPRB, -0.00585852322_JPRB, 0.00011873119_JPRB  /)
397RHSAVI( 8, :) = (/&
398    &0.126668819_JPRB, -0.00497949082_JPRB, 0.0000927679172_JPRB/)
399RHSAVI( 9, :) = (/&
400    &0.130938752_JPRB, -0.00531814674_JPRB, 0.0000984953029_JPRB/)
401RHSAVI(10, :) = (/&
402    &0.147024519_JPRB, -0.00671655774_JPRB, 0.000130885091_JPRB /)
403RHSAVI(11, :) = (/&
404    &0.241463914_JPRB, -0.0174486461_JPRB , 0.000457756556_JPRB /)
405RHSAVI(12, :) = (/&
406    &0.07576579_JPRB , -0.001695588_JPRB  , _ZERO_              /)
407RHSAVI(13, :) = (/&
408    &0.1032178_JPRB  , -0.00293412_JPRB   , _ZERO_              /)
409RHSAVI(14, :) = (/&
410    &0.07342832_JPRB , -0.001775135_JPRB  , _ZERO_              /)
411RHSAVI(15, :) = (/&
412    &0.04649514_JPRB , -0.0009165462_JPRB , _ZERO_              /)
413RHSAVI(16, :) = (/&
414    &0.06893519_JPRB , -0.0001245402_JPRB , _ZERO_              /)
415   
416! water clouds from Lindner & Li (2000)
417
418RLINLI( 1, :) = (/&   
419  & 0.88116E-01_JPRB,-0.12857E-02_JPRB, 0.81658E+00_JPRB,-0.39428E+01_JPRB, 0.46652E+01_JPRB /)     
420RLINLI( 2, :) = (/&     
421  & 0.41307E-03_JPRB,-0.59631E-04_JPRB, 0.24275E+01_JPRB,-0.90838E+01_JPRB, 0.96069E+01_JPRB /)   
422RLINLI( 3, :) = (/&     
423  &-0.57709E-01_JPRB, 0.99071E-03_JPRB, 0.31118E+01_JPRB,-0.95540E+01_JPRB, 0.90189E+01_JPRB /)   
424RLINLI( 4, :) = (/&         
425  &-0.53069E-01_JPRB, 0.99992E-03_JPRB, 0.28045E+01_JPRB,-0.72836E+01_JPRB, 0.62573E+01_JPRB /)
426RLINLI( 5, :) = (/&   
427  &-0.23627E-01_JPRB, 0.55291E-03_JPRB, 0.21785E+01_JPRB,-0.54664E+01_JPRB, 0.47379E+01_JPRB /)     
428RLINLI( 6, :) = (/&         
429  & 0.29022E-01_JPRB,-0.39657E-03_JPRB, 0.14902E+01_JPRB,-0.50777E+01_JPRB, 0.52170E+01_JPRB /)
430RLINLI( 7, :) = (/&         
431  &-0.24901E-01_JPRB, 0.16195E-03_JPRB, 0.29375E+02_JPRB,-0.11437E+02_JPRB, 0.12273E+02_JPRB /)
432RLINLI( 8, :) = (/&         
433  &-0.14269E+00_JPRB, 0.22282E-02_JPRB, 0.46478E+01_JPRB,-0.16369E+02_JPRB, 0.16533E+02_JPRB /)
434RLINLI( 9, :) = (/&   
435  &-0.20398E+00_JPRB, 0.34708E-02_JPRB, 0.52858E+01_JPRB,-0.16603E+02_JPRB, 0.15392E+02_JPRB /)     
436RLINLI(10, :) = (/&         
437  &-0.18318E+00_JPRB, 0.33080E-02_JPRB, 0.46120E+01_JPRB,-0.11550E+02_JPRB, 0.87086E+01_JPRB /)
438RLINLI(11, :) = (/&   
439  &-0.20420E+00_JPRB, 0.37167E-02_JPRB, 0.48566E+01_JPRB,-0.11972E+02_JPRB, 0.86344E+01_JPRB /)     
440RLINLI(12, :) = (/&         
441  &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
442RLINLI(13, :) = (/&
443  &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
444RLINLI(14, :) = (/&         
445  &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
446RLINLI(15, :) = (/&         
447  &-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /)
448RLINLI(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
453RFULIO( 1, :) = (/&
454 & -7.752E-03_JPRB, 4.624_JPRB, -42.01_JPRB /)
455RFULIO( 2, :) = (/&
456 & -1.741E-02_JPRB, 5.541_JPRB, -58.42_JPRB /)
457RFULIO( 3, :) = (/&
458 & -1.704E-02_JPRB, 4.830_JPRB,  16.27_JPRB /)
459RFULIO( 4, :) = (/&
460 & -1.151E-02_JPRB, 4.182_JPRB,  31.13_JPRB /)
461RFULIO( 5, :) = (/&
462 & -1.026E-02_JPRB, 4.105_JPRB,  16.36_JPRB /)
463RFULIO( 6, :) = (/&
464 & -8.294E-03_JPRB, 3.925_JPRB,  1.315_JPRB /)
465RFULIO( 7, :) = (/&
466 & -1.153E-02_JPRB, 4.109_JPRB,  17.32_JPRB /)
467RFULIO( 8, :) = (/&
468 & -9.609E-03_JPRB, 3.768_JPRB,  34.11_JPRB /)
469RFULIO( 9, :) = (/&
470 & -9.061E-03_JPRB, 3.741_JPRB,  26.48_JPRB /)
471RFULIO(10, :) = (/&
472 & -8.441E-03_JPRB, 3.715_JPRB,  19.48_JPRB /)
473RFULIO(11, :) = (/&
474 & -8.088E-03_JPRB, 3.717_JPRB,  17.17_JPRB /)
475RFULIO(12, :) = (/&
476 & -8.088E-03_JPRB, 3.717_JPRB,  17.17_JPRB /)
477RFULIO(13, :) = (/&
478 & -7.770E-03_JPRB, 3.734_JPRB,  11.85_JPRB /)
479RFULIO(14, :) = (/&
480 & -6.656E-03_JPRB, 3.686_JPRB,  _ZERO_     /)
481RFULIO(15, :) = (/&
482 & -6.656E-03_JPRB, 3.686_JPRB,  _ZERO_     /)
483RFULIO(16, :) = (/&
484 & -6.656E-03_JPRB, 3.686_JPRB,  _ZERO_     /)
485 
486! ice clouds from Fu et al. (1998)
487
488RFUETA( 1, :) = (/&
489 & 4.919685E-03_JPRB, 2.327741E+00_JPRB,-1.390858E+01_JPRB /)
490RFUETA( 2, :) = (/&
491 & 3.325756E-03_JPRB, 2.601360E+00_JPRB,-1.909602E+01_JPRB /)
492RFUETA( 3, :) = (/&
493 &-1.334860E-02_JPRB, 4.043808E+00_JPRB,-2.171029E+01_JPRB /)
494RFUETA( 4, :) = (/&
495 &-9.524174E-03_JPRB, 3.587742E+00_JPRB,-1.068895E+01_JPRB /)
496RFUETA( 5, :) = (/&
497 &-4.159424E-03_JPRB, 3.047325E+00_JPRB,-5.061568E+00_JPRB /)
498RFUETA( 6, :) = (/&
499 &-1.691632E-03_JPRB, 2.765756E+00_JPRB,-8.331033E+00_JPRB /)
500RFUETA( 7, :) = (/&
501 &-8.372696E-03_JPRB, 3.455018E+00_JPRB,-1.516692E+01_JPRB /)
502RFUETA( 8, :) = (/&
503 &-8.178608E-03_JPRB, 3.401245E+00_JPRB,-8.812820E+00_JPRB /)
504RFUETA( 9, :) = (/&
505 &-4.936610E-03_JPRB, 3.087764E+00_JPRB,-3.884262E+00_JPRB /)
506RFUETA(10, :) = (/&
507 &-3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /)
508RFUETA(11, :) = (/&
509 &-3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /)
510RFUETA(12, :) = (/&
511 &-2.465236E-03_JPRB, 2.833187E+00_JPRB,-4.227573E-01_JPRB /)
512RFUETA(13, :) = (/&
513 &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
514RFUETA(14, :) = (/&
515 &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
516RFUETA(15, :) = (/&
517 &-2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /)
518RFUETA(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
527REBCUI = (/  1.136_JPRB,  1.338_JPRB,  1.166_JPRB,  1.166_JPRB,  1.118_JPRB,  &
528 &0.600_JPRB /)
529REBCUJ = (/ 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
535REBCUG = (/ 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
541REBCUH = (/ 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
549RSUSHFA = (/ 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
557RTIW= 263._JPRB
558RRIW= 20._JPRB
559
560! Ice particle Effective Radius as a function of LWC
561
562REFFIA= 40._JPRB
563REFFIB= 0._JPRB
564
565! Sun-Shine
566
567RSUSHC= 0.0306_JPRB
568RSUSHD= 0.2548_JPRB
569
570!     ------------------------------------------------------------------
571
572! SW : absorption coefficients
573
574IF (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
635ELSEIF (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
696ELSEIF (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
757ELSE
758!  CALL ABOR1('SUCLOPN: WRONG SW SPECTRAL RESOLUTION')
759  STOP 'SUCLOPN: WRONG SW SPECTRAL RESOLUTION'
760ENDIF
761
762!     ------------------------------------------------------------------
763
764!*          2.    CLOUD OVERLAP PARAMETERS
765!                 ------------------------
766
767ZAOVLP = (/ -2.250E-05_JPRB,-7.316E-06_JPRB,-1.966E-05_JPRB /)
768ZBOVLP = (/  0.7865_JPRB   , 0.8186_JPRB   , 0.8900_JPRB    /)
769
770IF (KLEV == 19) THEN
771  RAOVLP=ZAOVLP(1)
772  RBOVLP=ZBOVLP(1)
773ELSE IF (KLEV == 31) THEN
774  RAOVLP=ZAOVLP(2)
775  RBOVLP=ZBOVLP(2)
776ELSE IF (KLEV == 60) THEN
777  RAOVLP=ZAOVLP(3)
778  RBOVLP=ZBOVLP(3)
779ELSE 
780  RAOVLP=ZAOVLP(3)
781  RBOVLP=ZBOVLP(3)
782END IF 
783
784!     ------------------------------------------------------------------
785
786RETURN
787END SUBROUTINE SUCLOPN
Note: See TracBrowser for help on using the repository browser.