source: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suclopn.F90 @ 5209

Last change on this file since 5209 was 5159, checked in by abarral, 3 months ago

Put dimensions.h and paramet.h into modules
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 52.3 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!     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
59USE PARKIND1  ,ONLY : JPIM     ,JPRB
60USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
61
62USE 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 
75USE YOERAD   , ONLY : RAOVLP   ,RBOVLP   ,LSRTM
76
77IMPLICIT NONE
78
79INTEGER(KIND=JPIM),INTENT(IN)    :: KTSW
80INTEGER(KIND=JPIM),INTENT(IN)    :: KSW
81INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV
82INTEGER(KIND=JPIM) ::  JNU
83
84!     -----------------------------------------------------------------
85REAL(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
96REAL(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 
110REAL(KIND=JPRB) :: ZAOVLP(3), ZBOVLP(3) 
111
112REAL(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
126REAL(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)   
139REAL(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
158IF (LHOOK) CALL DR_HOOK('SUCLOPN',0,ZHOOK_HANDLE)
159ZEBCUA2 = (/ 3.448E-04_JPRB , 3.448E-04_JPRB /)
160ZEBCUB2 = (/ 2.431_JPRB     , 2.431_JPRB     /)
161ZEBCUC2 = (/ 0.00001_JPRB   , 0.024366_JPRB  /)
162ZEBCUD2 = (/ 0.0_JPRB         , 2.487E-04_JPRB /)
163ZEBCUE2 = (/ 0.7661_JPRB    , 0.7866_JPRB    /)
164ZEBCUF2 = (/ 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
176ZYFWCA2 = (/ 0.0_JPRB       , 0.0_JPRB       /)
177ZYFWCB2 = (/ 1.5_JPRB       , 1.5_JPRB       /)
178ZYFWCC2 = (/ 0.9999_JPRB    , 0.9988_JPRB    /)
179ZYFWCD2 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB /)
180ZYFWCE2 = (/ 0.5_JPRB       , 0.05_JPRB      /)
181ZYFWCF2 = (/ 0.865_JPRB     , 0.910_JPRB     /)
182
183!* Water cloud properties - from Slingo (1989)
184
185ZASWCA2 = (/ 2.817_JPRB     , 2.455_JPRB    /)
186ZASWCB2 = (/ 1.305_JPRB     , 1.400_JPRB    /)
187ZASWCC2 = (/-5.62E-08_JPRB  , 1.28E-02_JPRB /)
188ZASWCD2 = (/ 1.63E-07_JPRB  , 8.97E-04_JPRB /)
189ZASWCE2 = (/ 0.829_JPRB     , 0.783_JPRB    /)
190ZASWCF2 = (/ 2.482_JPRB     , 4.988_JPRB    /)
191
192!* Ice cloud properties - from Sun and Shine (1995)
193
194ZSUSHE2 = (/ 0.0_JPRB       , 8.6822_JPRB    /)
195ZSUSHF2 = (/ 0.0_JPRB       , 9.6277_JPRB    /)
196ZSUSHH2 = (/ 0.8522_JPRB    , 0.8819_JPRB    /)
197ZSUSHK2 = (/ 0.1620_JPRB    , 0.1630_JPRB    /)
198ZSUSHA2 = (/ 0.0_JPRB       , 23.204_JPRB    /)
199ZSUSHG2 = (/ 0.3270_JPRB    , 0.4180_JPRB    /)
200
201!* Ice cloud properties - from Fu and Liou (1993)
202
203ZFLAA02 = (/-6.656E-3_JPRB  ,-6.656E-3_JPRB  /)
204ZFLAA12 = (/ 3.686_JPRB     , 3.686_JPRB     /)
205ZFLBB02 = (/ .10998E-5_JPRB , .21136E-1_JPRB /)
206ZFLBB12 = (/-.26101E-7_JPRB , .39150E-3_JPRB /)
207ZFLBB22 = (/ .18096E-8_JPRB ,-.20740E-6_JPRB /)
208ZFLBB32 = (/-.47387E-11_JPRB,-.28829E-8_JPRB /)
209ZFLCC02 = (/ .22110E+1_JPRB , .22498E+1_JPRB /)
210ZFLCC12 = (/-.10398E-2_JPRB , .23656E-3_JPRB /)
211ZFLCC22 = (/ .65199E-4_JPRB , .51948E-4_JPRB /)
212ZFLCC32 = (/-.34498E-6_JPRB ,-.29768E-6_JPRB /)
213ZFLDD02 = (/ .12495_JPRB    , .11716_JPRB    /)
214ZFLDD12 = (/-.43582E-3_JPRB ,-.45208E-3_JPRB /)
215ZFLDD22 = (/ .14092E-4_JPRB , .12772E-4_JPRB /)
216ZFLDD32 = (/-.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
250ZEBCUA4 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /)
251ZEBCUB4 = (/ 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     /)
252ZEBCUC4 = (/ 0.00001_JPRB   , 0.00011_JPRB   , 0.0197796_JPRB , 0.46658_JPRB   /)
253ZEBCUD4 = (/ 0.0_JPRB       , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /)
254ZEBCUE4 = (/ 0.7661_JPRB    , 0.7730_JPRB    , 0.795653_JPRB  , 0.9595_JPRB    /)
255ZEBCUF4 = (/ 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
259ZYFWCA4 = (/ 0.0_JPRB       , 0.0_JPRB       , 0.0_JPRB      , 0.0_JPRB        /)
260ZYFWCB4 = (/ 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB      , 1.5_JPRB        /)
261ZYFWCC4 = (/ 0.9999_JPRB    , 0.9988_JPRB    , 0.9988_JPRB   , 0.9988_JPRB     /)
262ZYFWCD4 = (/ 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB  /)
263ZYFWCE4 = (/ 0.5_JPRB       , 0.05_JPRB      , 0.05_JPRB     , 0.05_JPRB       /)
264ZYFWCF4 = (/ 0.865_JPRB     , 0.910_JPRB     , 0.910_JPRB    , 0.910_JPRB      /)
265
266!* Water cloud properties - from Slingo (1989)
267
268ZASWCA4 = (/ 2.817_JPRB     , 2.682_JPRB     , 2.264_JPRB    , 1.281_JPRB      /)
269ZASWCB4 = (/ 1.305_JPRB     , 1.346_JPRB     , 1.454_JPRB    , 1.641_JPRB      /)
270ZASWCC4 = (/-5.62E-08_JPRB  ,-6.94E-06_JPRB  , 4.64E-04_JPRB , 2.01E-01_JPRB   /)
271ZASWCD4 = (/ 1.63E-07_JPRB  , 2.35E-05_JPRB  , 1.24E-03_JPRB , 7.56E-03_JPRB   /)
272ZASWCE4 = (/ 0.829_JPRB     , 0.794_JPRB     , 0.754_JPRB    , 0.826_JPRB      /)
273ZASWCF4 = (/ 2.482_JPRB     , 4.226_JPRB     , 6.560_JPRB    , 4.353_JPRB      /)
274
275!* Ice cloud properties - from Sun and Shine (1995)
276
277ZSUSHE4 = (/ 0.0_JPRB     , 7.2471E-02_JPRB , 17.5933_JPRB  , 48.7166_JPRB  /)
278ZSUSHF4 = (/ 0.0_JPRB     , 4.01511_JPRB    , 21.1249_JPRB  , 1.20890_JPRB  /)
279ZSUSHH4 = (/ 0.8522_JPRB  , 0.85841_JPRB    , 0.90778_JPRB  ,0.982046_JPRB  /)
280ZSUSHK4 = (/ 0.1620_JPRB  , 0.160048_JPRB   , 0.188521_JPRB ,0.0411446_JPRB /)
281ZSUSHA4 = (/ 0.0_JPRB     , 0.273455_JPRB   , 41.7675_JPRB  , 161.104_JPRB  /)
282ZSUSHG4 = (/ 0.3270_JPRB  , 0.343668_JPRB   , 0.526192_JPRB , 0.574040_JPRB /)
283
284!* Ice cloud properties - from Fu and Liou (1993)
285
286ZFLAA04 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /)
287ZFLAA14 = (/ 3.686_JPRB     , 3.686_JPRB    , 3.686_JPRB    , 3.686_JPRB     /)
288ZFLBB04 = (/ .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /)
289ZFLBB14 = (/-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /)
290ZFLBB24 = (/ .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /)
291ZFLBB34 = (/-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /)
292ZFLCC04 = (/ .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /)
293ZFLCC14 = (/-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /)
294ZFLCC24 = (/ .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /)
295ZFLCC34 = (/-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /)
296ZFLDD04 = (/ .12495_JPRB    , .12363_JPRB   , .11983_JPRB   , .21834E-1_JPRB /)
297ZFLDD14 = (/-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /)
298ZFLDD24 = (/ .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /)
299ZFLDD34 = (/-.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
303ZFUAA04 = (/-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /)
304ZFUAA14 = (/ 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /)
305ZFUBB04 = (/-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /)
306ZFUBB14 = (/ 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /)
307ZFUBB24 = (/-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /)
308ZFUBB34 = (/ 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /)
309ZFUCC04 = (/ 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /)
310ZFUCC14 = (/ 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /)
311ZFUCC24 = (/-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /)
312ZFUCC34 = (/-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
323ZEBCUA6 = (/ 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB /)
324ZEBCUB6 = (/ 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     /)
325ZEBCUC6 = (/ 0.00001_JPRB   , 0.00001_JPRB   , 0.00001_JPRB   , 0.00011_JPRB   , 0.0197796_JPRB , 0.46658_JPRB   /)
326ZEBCUD6 = (/ 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 1.405E-05_JPRB , 7.95513E-4_JPRB, 2.050E-05_JPRB /)
327ZEBCUE6 = (/ 0.7661_JPRB    , 0.7661_JPRB    , 0.7661_JPRB    , 0.7730_JPRB    , 0.795653_JPRB  , 0.9595_JPRB    /)
328ZEBCUF6 = (/ 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
332ZYFWCA6 = (/ 0.0_JPRB       , 0.0_JPRB       , 0.0_JPRB       , 0.0_JPRB       , 0.0_JPRB      , 0.0_JPRB        /)
333ZYFWCB6 = (/ 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB      , 1.5_JPRB        /)
334ZYFWCC6 = (/ 0.9999_JPRB    , 0.9999_JPRB    , 0.9999_JPRB    , 0.9988_JPRB    , 0.9988_JPRB   , 0.9988_JPRB     /)
335ZYFWCD6 = (/ 5.000E-04_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB, 2.500E-03_JPRB  /)
336ZYFWCE6 = (/ 0.5_JPRB       , 0.5_JPRB       , 0.5_JPRB       , 0.05_JPRB      , 0.05_JPRB     , 0.05_JPRB       /)
337ZYFWCF6 = (/ 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
341ZASWCA6 = (/ 3.094_JPRB     , 2.881_JPRB     , 2.790_JPRB     , 2.682_JPRB     , 2.254_JPRB    , 1.281_JPRB      /)
342ZASWCB6 = (/ 1.252_JPRB     , 1.284_JPRB     , 1.313_JPRB     , 1.346_JPRB     , 1.456_JPRB    , 1.641_JPRB      /)
343ZASWCC6 = (/ 7.90E-07_JPRB  , 1.77E-07_JPRB  ,-1.55E-07_JPRB  ,-6.94E-06_JPRB  , 4.75E-04_JPRB , 2.01E-01_JPRB   /)
344ZASWCD6 = (/ 3.69E-07_JPRB  , 1.23E-07_JPRB  , 1.80E-07_JPRB  , 2.35E-05_JPRB  , 1.26E-03_JPRB , 7.56E-03_JPRB   /)
345ZASWCE6 = (/ 0.844_JPRB     , 0.839_JPRB     , 0.824_JPRB     , 0.794_JPRB     , 0.754_JPRB    , 0.826_JPRB      /)
346ZASWCF6 = (/ 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
350ZSUSHE6 = (/ 0.0_JPRB     , 0.0_JPRB     , 0.0_JPRB     , 7.2471E-02_JPRB , 17.5933_JPRB  , 48.7166_JPRB  /)
351ZSUSHF6 = (/ 0.0_JPRB     , 0.0_JPRB     , 0.0_JPRB     , 4.01511_JPRB    , 21.1249_JPRB  , 1.20890_JPRB  /)
352ZSUSHH6 = (/ 0.8522_JPRB  , 0.8522_JPRB  , 0.8522_JPRB  , 0.85841_JPRB    , 0.90778_JPRB  ,0.982046_JPRB  /)
353ZSUSHK6 = (/ 0.1620_JPRB  , 0.1620_JPRB  , 0.1620_JPRB  , 0.160048_JPRB   , 0.188521_JPRB ,0.0411446_JPRB /)
354ZSUSHA6 = (/ 0.0_JPRB     , 0.0_JPRB     , 0.0_JPRB     , 0.273455_JPRB   , 41.7675_JPRB  , 161.104_JPRB  /)
355ZSUSHG6 = (/ 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
359ZFLAA06 = (/-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB ,-6.656E-03_JPRB,-6.656E-03_JPRB,-6.656E-03_JPRB /)
360ZFLAA16 = (/ 3.686_JPRB     , 3.686_JPRB     , 3.686_JPRB     , 3.686_JPRB    , 3.686_JPRB    , 3.686_JPRB     /)
361ZFLBB06 = (/ .10998E-5_JPRB , .10998E-5_JPRB , .10998E-5_JPRB , .20208E-4_JPRB, .51557E-3_JPRB, .39517E+0_JPRB /)
362ZFLBB16 = (/-.26101E-7_JPRB ,-.26101E-7_JPRB ,-.26101E-7_JPRB , .96483E-5_JPRB, .10731E-2_JPRB, .15787E-2_JPRB /)
363ZFLBB26 = (/ .18096E-8_JPRB , .18096E-8_JPRB , .18096E-8_JPRB , .83009E-7_JPRB, .17753E-5_JPRB,-.14337E-4_JPRB /)
364ZFLBB36 = (/-.47387E-11_JPRB,-.47387E-11_JPRB,-.47387E-11_JPRB,-.32217E-9_JPRB,-.18379E-7_JPRB, .46942E-7_JPRB /)
365ZFLCC06 = (/ .22110E+1_JPRB , .22110E+1_JPRB , .22110E+1_JPRB , .22151E+1_JPRB, .22534E+1_JPRB, .26653E+1_JPRB /)
366ZFLCC16 = (/-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.10398E-2_JPRB ,-.77982E-3_JPRB, .16163E-2_JPRB, .56935E-2_JPRB /)
367ZFLCC26 = (/ .65199E-4_JPRB , .65199E-4_JPRB , .65199E-4_JPRB , .63750E-4_JPRB, .44037E-4_JPRB,-.54316E-4_JPRB /)
368ZFLCC36 = (/-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34498E-6_JPRB ,-.34466E-6_JPRB,-.27627E-6_JPRB, .17858E-6_JPRB /)
369ZFLDD06 = (/ .12495_JPRB    , .12495_JPRB    , .12495_JPRB    , .12363_JPRB   , .11983_JPRB   , .21834E-1_JPRB /)
370ZFLDD16 = (/-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.43582E-3_JPRB ,-.44419E-3_JPRB,-.50108E-3_JPRB,-.29204E-3_JPRB /)
371ZFLDD26 = (/ .14092E-4_JPRB , .14092E-4_JPRB , .14092E-4_JPRB , .14038E-4_JPRB, .11843E-4_JPRB, .18060E-5_JPRB /)
372ZFLDD36 = (/-.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
376ZFUAA06 = (/-1.30817E-04_JPRB,-1.30817E-04_JPRB,-1.30817E-04_JPRB,-6.39479E-05_JPRB,-6.74730E-06_JPRB, 1.62674E-04_JPRB /)
377ZFUAA16 = (/ 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52883E+00_JPRB, 2.52393E+00_JPRB, 2.52056E+00_JPRB, 2.49823E+00_JPRB /)
378ZFUBB06 = (/-1.55357E-08_JPRB,-1.55357E-08_JPRB,-1.55357E-08_JPRB,-7.90657E-07_JPRB, 9.64842E-04_JPRB, 2.25112E-01_JPRB /)
379ZFUBB16 = (/ 1.95793E-07_JPRB, 1.95793E-07_JPRB, 1.95793E-07_JPRB, 7.79991E-06_JPRB, 9.09809E-04_JPRB, 3.05017E-03_JPRB /)
380ZFUBB26 = (/-2.31234E-10_JPRB,-2.31234E-10_JPRB,-2.31234E-10_JPRB, 2.90894E-10_JPRB,-3.57557E-06_JPRB,-2.54236E-05_JPRB /)
381ZFUBB36 = (/ 1.12247E-12_JPRB, 1.12247E-12_JPRB, 1.12247E-12_JPRB,-2.02818E-12_JPRB, 1.00197E-08_JPRB, 8.49116E-08_JPRB /)
382ZFUCC06 = (/ 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.39781E-01_JPRB, 7.52335E-01_JPRB, 7.56307E-01_JPRB, 8.30812E-01_JPRB /)
383ZFUCC16 = (/ 9.10564E-04_JPRB, 9.10564E-04_JPRB, 9.10564E-04_JPRB, 1.06211E-03_JPRB, 1.73364E-03_JPRB, 2.62788E-03_JPRB /)
384ZFUCC26 = (/-4.62479E-07_JPRB,-4.62479E-07_JPRB,-4.62479E-07_JPRB,-2.45770E-06_JPRB,-8.92191E-06_JPRB,-2.43196E-05_JPRB /)
385ZFUCC36 = (/-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
410ZEBCUA14 = (/ &
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
415ZEBCUB14 = (/ &
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
420ZEBCUC14 = (/ &
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
425ZEBCUD14 = (/ &
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
430ZEBCUE14 = (/ &
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
435ZEBCUF14 = (/ &
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
442ZYFWCA14 = (/ &
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 
447ZYFWCB14 = (/ &
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 
452ZYFWCC14 = (/ &
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
457ZYFWCD14 = (/ &
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
462ZYFWCE14 = (/ &
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
467ZYFWCF14 = (/ &
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
474ZASWCA14 = (/ &
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
479ZASWCB14 = (/ &
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   
484ZASWCC14 = (/ &
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   
489ZASWCD14 = (/ &
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
494ZASWCE14 = (/ &
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
499ZASWCF14 = (/ &
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
506ZSUSHE14 = (/ &
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         /) 
510ZSUSHF14 = (/ &
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         /) 
514ZSUSHH14 = (/ &
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         /) 
518ZSUSHK14 = (/ &
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         /) 
522ZSUSHA14 = (/ &
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         /) 
526ZSUSHG14 = (/ &
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
533ZFLAA014 = (/ &
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         /) 
537ZFLAA114 = (/ &
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         /) 
541ZFLBB014 = (/ &
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         /) 
545ZFLBB114 = (/ &
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         /) 
549ZFLBB214 = (/ &
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         /) 
553ZFLBB314 = (/ &
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         /) 
557ZFLCC014 = (/ &
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         /) 
561ZFLCC114 = (/ &
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         /) 
565ZFLCC214 = (/ &
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         /) 
569ZFLCC314 = (/ &
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         /) 
573ZFLDD014 = (/ &
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         /) 
577ZFLDD114 = (/ &
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         /) 
581ZFLDD214 = (/ &
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         /) 
585ZFLDD314 = (/ &
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
592ZFUAA014 = (/ &
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
597ZFUAA114 = (/ &
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
602ZFUBB014 = (/ &
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
607ZFUBB114 = (/ &
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
612ZFUBB214 = (/ &
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
617ZFUBB314 = (/ &
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
622ZFUCC014 = (/ &
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   
627ZFUCC114 = (/ &
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
632ZFUCC214 = (/ &
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
637ZFUCC314 = (/ &
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
648RHSAVI( 1, :) = (/  0.1651082_JPRB  , -0.003494839_JPRB  , 0.0_JPRB            /) 
649RHSAVI( 2, :) = (/  0.327820597_JPRB, -0.0219634383_JPRB , 0.000506783898_JPRB /)   
650RHSAVI( 3, :) = (/  0.504805453_JPRB, -0.0478602354_JPRB , 0.00141521102_JPRB  /) 
651RHSAVI( 4, :) = (/  0.513169093_JPRB, -0.0508960145_JPRB , 0.00155498711_JPRB  /) 
652RHSAVI( 5, :) = (/  0.448042082_JPRB, -0.0431857592_JPRB , 0.00130848978_JPRB  /) 
653RHSAVI( 6, :) = (/  0.249547237_JPRB, -0.0185273835_JPRB , 0.00050361258_JPRB  /) 
654RHSAVI( 7, :) = (/  0.135486796_JPRB, -0.00585852322_JPRB, 0.00011873119_JPRB  /) 
655RHSAVI( 8, :) = (/  0.126668819_JPRB, -0.00497949082_JPRB, 0.0000927679172_JPRB/) 
656RHSAVI( 9, :) = (/  0.130938752_JPRB, -0.00531814674_JPRB, 0.0000984953029_JPRB/) 
657RHSAVI(10, :) = (/  0.147024519_JPRB, -0.00671655774_JPRB, 0.000130885091_JPRB /) 
658RHSAVI(11, :) = (/  0.241463914_JPRB, -0.0174486461_JPRB , 0.000457756556_JPRB /) 
659RHSAVI(12, :) = (/  0.07576579_JPRB , -0.001695588_JPRB  , 0.0_JPRB            /) 
660RHSAVI(13, :) = (/  0.1032178_JPRB  , -0.00293412_JPRB   , 0.0_JPRB            /) 
661RHSAVI(14, :) = (/  0.07342832_JPRB , -0.001775135_JPRB  , 0.0_JPRB            /) 
662RHSAVI(15, :) = (/  0.04649514_JPRB , -0.0009165462_JPRB , 0.0_JPRB            /) 
663RHSAVI(16, :) = (/  0.06893519_JPRB , -0.0001245402_JPRB , 0.0_JPRB            /) 
664   
665! water clouds from Lindner & Li (2000)
666
667RLILIA( 1, :) = (/ 0.88116E-01_JPRB,-0.12857E-02_JPRB, 0.81658E+00_JPRB,-0.39428E+01_JPRB, 0.46652E+01_JPRB /)       
668RLILIA( 2, :) = (/ 0.41307E-03_JPRB,-0.59631E-04_JPRB, 0.24275E+01_JPRB,-0.90838E+01_JPRB, 0.96069E+01_JPRB /)     
669RLILIA( 3, :) = (/-0.57709E-01_JPRB, 0.99071E-03_JPRB, 0.31118E+01_JPRB,-0.95540E+01_JPRB, 0.90189E+01_JPRB /)     
670RLILIA( 4, :) = (/-0.53069E-01_JPRB, 0.99992E-03_JPRB, 0.28045E+01_JPRB,-0.72836E+01_JPRB, 0.62573E+01_JPRB /) 
671RLILIA( 5, :) = (/-0.23627E-01_JPRB, 0.55291E-03_JPRB, 0.21785E+01_JPRB,-0.54664E+01_JPRB, 0.47379E+01_JPRB /)       
672RLILIA( 6, :) = (/ 0.29022E-01_JPRB,-0.39657E-03_JPRB, 0.14902E+01_JPRB,-0.50777E+01_JPRB, 0.52170E+01_JPRB /) 
673RLILIA( 7, :) = (/-0.24901E-01_JPRB, 0.16195E-03_JPRB, 0.29375E+02_JPRB,-0.11437E+02_JPRB, 0.12273E+02_JPRB /) 
674RLILIA( 8, :) = (/-0.14269E+00_JPRB, 0.22282E-02_JPRB, 0.46478E+01_JPRB,-0.16369E+02_JPRB, 0.16533E+02_JPRB /) 
675RLILIA( 9, :) = (/-0.20398E+00_JPRB, 0.34708E-02_JPRB, 0.52858E+01_JPRB,-0.16603E+02_JPRB, 0.15392E+02_JPRB /)       
676RLILIA(10, :) = (/-0.18318E+00_JPRB, 0.33080E-02_JPRB, 0.46120E+01_JPRB,-0.11550E+02_JPRB, 0.87086E+01_JPRB /) 
677RLILIA(11, :) = (/-0.20420E+00_JPRB, 0.37167E-02_JPRB, 0.48566E+01_JPRB,-0.11972E+02_JPRB, 0.86344E+01_JPRB /)       
678RLILIA(12, :) = (/-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) 
679RLILIA(13, :) = (/-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) 
680RLILIA(14, :) = (/-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) 
681RLILIA(15, :) = (/-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) 
682RLILIA(16, :) = (/-0.14037E+00_JPRB, 0.28058E-02_JPRB, 0.34969E+01_JPRB,-0.33770E+01_JPRB,-0.23541E+01_JPRB /) 
683
684RLILIB( 1, :) = (/ 0.76110E+00_JPRB, 0.57048E+00_JPRB,-0.16721E-01_JPRB, 0.27886E-03_JPRB /)
685RLILIB( 2, :) = (/ 0.51897E+00_JPRB, 0.86352E+00_JPRB,-0.29099E-02_JPRB, 0.63126E-04_JPRB /)
686RLILIB( 3, :) = (/ 0.50964E+00_JPRB, 0.76057E+00_JPRB,-0.15630E-02_JPRB, 0.26245E-04_JPRB /)
687RLILIB( 4, :) = (/ 0.55180E+00_JPRB, 0.65812E+00_JPRB,-0.41768E-02_JPRB, 0.64243E-04_JPRB /)
688RLILIB( 5, :) = (/ 0.60890E+00_JPRB, 0.59525E+00_JPRB,-0.74557E-02_JPRB, 0.11757E-03_JPRB /)
689RLILIB( 6, :) = (/ 0.49787E+00_JPRB, 0.74581E+00_JPRB,-0.37379E-02_JPRB, 0.90555E-04_JPRB /)
690RLILIB( 7, :) = (/ 0.46323E-01_JPRB, 0.10131E+01_JPRB, 0.18692E-01_JPRB,-0.21898E-03_JPRB /)
691RLILIB( 8, :) = (/-0.43354E-01_JPRB, 0.78454E+00_JPRB, 0.26449E-01_JPRB,-0.36385E-03_JPRB /)
692RLILIB( 9, :) = (/-0.27033E-01_JPRB, 0.54717E+00_JPRB, 0.27903E-01_JPRB,-0.41074E-03_JPRB /)
693RLILIB(10, :) = (/ 0.11595E+00_JPRB, 0.30559E+00_JPRB, 0.22727E-01_JPRB,-0.36032E-03_JPRB /)
694RLILIB(11, :) = (/-0.14764E-01_JPRB, 0.28331E+00_JPRB, 0.23515E-01_JPRB,-0.33662E-03_JPRB /)
695RLILIB(12, :) = (/ 0.68056E-02_JPRB, 0.59837E-01_JPRB, 0.21125E-01_JPRB,-0.29740E-03_JPRB /)
696RLILIB(13, :) = (/ 0.68056E-02_JPRB, 0.59837E-01_JPRB, 0.21125E-01_JPRB,-0.29740E-03_JPRB /)
697RLILIB(14, :) = (/ 0.68056E-02_JPRB, 0.59837E-01_JPRB, 0.21125E-01_JPRB,-0.29740E-03_JPRB /)
698RLILIB(15, :) = (/ 0.68056E-02_JPRB, 0.59837E-01_JPRB, 0.21125E-01_JPRB,-0.29740E-03_JPRB /)
699RLILIB(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
703RFULIO( 1, :) = (/ -7.752E-03_JPRB, 4.624_JPRB, -42.01_JPRB /) 
704RFULIO( 2, :) = (/ -1.741E-02_JPRB, 5.541_JPRB, -58.42_JPRB /) 
705RFULIO( 3, :) = (/ -1.704E-02_JPRB, 4.830_JPRB,  16.27_JPRB /) 
706RFULIO( 4, :) = (/ -1.151E-02_JPRB, 4.182_JPRB,  31.13_JPRB /) 
707RFULIO( 5, :) = (/ -1.026E-02_JPRB, 4.105_JPRB,  16.36_JPRB /) 
708RFULIO( 6, :) = (/ -8.294E-03_JPRB, 3.925_JPRB,  1.315_JPRB /) 
709RFULIO( 7, :) = (/ -1.153E-02_JPRB, 4.109_JPRB,  17.32_JPRB /) 
710RFULIO( 8, :) = (/ -9.609E-03_JPRB, 3.768_JPRB,  34.11_JPRB /) 
711RFULIO( 9, :) = (/ -9.061E-03_JPRB, 3.741_JPRB,  26.48_JPRB /) 
712RFULIO(10, :) = (/ -8.441E-03_JPRB, 3.715_JPRB,  19.48_JPRB /) 
713RFULIO(11, :) = (/ -8.088E-03_JPRB, 3.717_JPRB,  17.17_JPRB /) 
714RFULIO(12, :) = (/ -8.088E-03_JPRB, 3.717_JPRB,  17.17_JPRB /) 
715RFULIO(13, :) = (/ -7.770E-03_JPRB, 3.734_JPRB,  11.85_JPRB /) 
716RFULIO(14, :) = (/ -6.656E-03_JPRB, 3.686_JPRB,  0.0_JPRB   /) 
717RFULIO(15, :) = (/ -6.656E-03_JPRB, 3.686_JPRB,  0.0_JPRB   /) 
718RFULIO(16, :) = (/ -6.656E-03_JPRB, 3.686_JPRB,  0.0_JPRB   /) 
719 
720! ice clouds from Fu et al. (1998)
721
722!-- extinction coefficient
723RFUETA( 1, :) = (/  4.919685E-03_JPRB, 2.327741E+00_JPRB,-1.390858E+01_JPRB /) 
724RFUETA( 2, :) = (/  3.325756E-03_JPRB, 2.601360E+00_JPRB,-1.909602E+01_JPRB /) 
725RFUETA( 3, :) = (/ -1.334860E-02_JPRB, 4.043808E+00_JPRB,-2.171029E+01_JPRB /) 
726RFUETA( 4, :) = (/ -9.524174E-03_JPRB, 3.587742E+00_JPRB,-1.068895E+01_JPRB /) 
727RFUETA( 5, :) = (/ -4.159424E-03_JPRB, 3.047325E+00_JPRB,-5.061568E+00_JPRB /) 
728RFUETA( 6, :) = (/ -1.691632E-03_JPRB, 2.765756E+00_JPRB,-8.331033E+00_JPRB /) 
729RFUETA( 7, :) = (/ -8.372696E-03_JPRB, 3.455018E+00_JPRB,-1.516692E+01_JPRB /) 
730RFUETA( 8, :) = (/ -8.178608E-03_JPRB, 3.401245E+00_JPRB,-8.812820E+00_JPRB /) 
731RFUETA( 9, :) = (/ -4.936610E-03_JPRB, 3.087764E+00_JPRB,-3.884262E+00_JPRB /) 
732RFUETA(10, :) = (/ -3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /) 
733RFUETA(11, :) = (/ -3.034573E-03_JPRB, 2.900043E+00_JPRB,-1.849911E+00_JPRB /) 
734RFUETA(12, :) = (/ -2.465236E-03_JPRB, 2.833187E+00_JPRB,-4.227573E-01_JPRB /) 
735RFUETA(13, :) = (/ -2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) 
736RFUETA(14, :) = (/ -2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) 
737RFUETA(15, :) = (/ -2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) 
738RFUETA(16, :) = (/ -2.308881E-03_JPRB, 2.814002E+00_JPRB, 1.072211E+00_JPRB /) 
739
740!-- absorption coefficient
741RFUETB( 1, :) = (/  8.869787E-01_JPRB, 2.118409E-02_JPRB,-2.781429E-04_JPRB, 1.094562E-06_JPRB /)
742RFUETB( 2, :) = (/  2.005578E-01_JPRB, 2.132614E-02_JPRB,-1.751052E-04_JPRB, 5.355885E-07_JPRB /)
743RFUETB( 3, :) = (/  3.003701E-01_JPRB, 2.051529E-02_JPRB,-1.931684E-04_JPRB, 6.583031E-07_JPRB /)
744RFUETB( 4, :) = (/  9.551440E-01_JPRB, 1.309792E-02_JPRB,-1.793694E-04_JPRB, 7.313392E-07_JPRB /)
745RFUETB( 5, :) = (/  1.466481E+00_JPRB,-2.129226E-03_JPRB,-1.361630E-05_JPRB, 1.193649E-07_JPRB /)
746RFUETB( 6, :) = (/  1.195515E+00_JPRB, 3.350616E-03_JPRB,-5.266996E-05_JPRB, 2.233377E-07_JPRB /)
747RFUETB( 7, :) = (/  5.409536E-01_JPRB, 1.949649E-02_JPRB,-2.050908E-04_JPRB, 7.364680E-07_JPRB /)
748RFUETB( 8, :) = (/  5.874323E-01_JPRB, 1.876628E-02_JPRB,-2.045834E-04_JPRB, 7.510080E-07_JPRB /)
749RFUETB( 9, :) = (/  7.152274E-01_JPRB, 1.621734E-02_JPRB,-1.868544E-04_JPRB, 7.078738E-07_JPRB /)
750RFUETB(10, :) = (/  8.862434E-01_JPRB, 1.226538E-02_JPRB,-1.523076E-04_JPRB, 6.000892E-07_JPRB /)
751RFUETB(11, :) = (/  8.862434E-01_JPRB, 1.226538E-02_JPRB,-1.523076E-04_JPRB, 6.000892E-07_JPRB /)
752RFUETB(12, :) = (/  7.428957E-01_JPRB, 1.279601E-02_JPRB,-1.391803E-04_JPRB, 5.180104E-07_JPRB /)
753RFUETB(13, :) = (/  4.346482E-01_JPRB, 1.721457E-02_JPRB,-1.623227E-04_JPRB, 5.561523E-07_JPRB /)
754RFUETB(14, :) = (/  4.346482E-01_JPRB, 1.721457E-02_JPRB,-1.623227E-04_JPRB, 5.561523E-07_JPRB /)
755RFUETB(15, :) = (/  4.346482E-01_JPRB, 1.721457E-02_JPRB,-1.623227E-04_JPRB, 5.561523E-07_JPRB /)
756RFUETB(16, :) = (/  4.346482E-01_JPRB, 1.721457E-02_JPRB,-1.623227E-04_JPRB, 5.561523E-07_JPRB /)
757
758!-- asymmetry factor
759RFUETC( 1, :) = (/  4.949276E-01_JPRB, 1.186174E-02_JPRB,-1.267629E-04_JPRB, 4.603574E-07_JPRB /)
760RFUETC( 2, :) = (/  6.891414E-01_JPRB, 6.192281E-03_JPRB,-6.459514E-05_JPRB, 2.436963E-07_JPRB /)
761RFUETC( 3, :) = (/  7.260484E-01_JPRB, 2.664334E-03_JPRB,-1.251136E-05_JPRB, 2.243377E-08_JPRB /)
762RFUETC( 4, :) = (/  7.363466E-01_JPRB, 4.798266E-03_JPRB,-4.413293E-05_JPRB, 1.525774E-07_JPRB /)
763RFUETC( 5, :) = (/  7.984021E-01_JPRB, 3.977117E-03_JPRB,-4.471984E-05_JPRB, 1.694919E-07_JPRB /)
764RFUETC( 6, :) = (/  8.663385E-01_JPRB, 2.797934E-03_JPRB,-3.187011E-05_JPRB, 1.217209E-07_JPRB /)
765RFUETC( 7, :) = (/  8.906280E-01_JPRB, 1.903269E-03_JPRB,-1.733552E-05_JPRB, 5.855071E-08_JPRB /)
766RFUETC( 8, :) = (/  8.609604E-01_JPRB, 2.200445E-03_JPRB,-1.748105E-05_JPRB, 5.176616E-08_JPRB /)
767RFUETC( 9, :) = (/  8.522816E-01_JPRB, 2.523627E-03_JPRB,-2.149196E-05_JPRB, 6.685067E-08_JPRB /)
768RFUETC(10, :) = (/  8.741665E-01_JPRB, 2.455409E-03_JPRB,-2.456935E-05_JPRB, 8.641223E-08_JPRB /)
769RFUETC(11, :) = (/  8.741665E-01_JPRB, 2.455409E-03_JPRB,-2.456935E-05_JPRB, 8.641223E-08_JPRB /)
770RFUETC(12, :) = (/  8.472918E-01_JPRB, 2.559953E-03_JPRB,-2.182660E-05_JPRB, 6.879977E-08_JPRB /)
771RFUETC(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
778REBCUI = (/  1.136_JPRB,  1.338_JPRB,  1.166_JPRB,  1.166_JPRB,  1.118_JPRB,  &
779 & 0.600_JPRB /) 
780REBCUJ = (/ 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
786REBCUG = (/ 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
791REBCUH = (/ 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
798RSUSHFA = (/ 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
806RTIW= 263._JPRB
807RRIW= 20._JPRB
808
809! Ice particle Effective Radius as a function of LWC
810
811REFFIA= 40._JPRB
812REFFIB= 0._JPRB
813
814! Sun-Shine
815
816RSUSHC= 0.0306_JPRB
817RSUSHD= 0.2548_JPRB
818
819!     ------------------------------------------------------------------
820
821! SW : absorption coefficients
822
823IF (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
881ELSEIF (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
938ELSEIF (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
995ELSEIF (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
1058ELSE
1059  CALL ABOR1('SUCLOPN: WRONG SW SPECTRAL RESOLUTION')
1060ENDIF
1061!     ------------------------------------------------------------------
1062
1063!*          2.    CLOUD OVERLAP PARAMETERS
1064!                 ------------------------
1065
1066ZAOVLP = (/ -2.250E-05_JPRB,-7.316E-06_JPRB,-1.966E-05_JPRB /)
1067ZBOVLP = (/  0.7865_JPRB   , 0.8186_JPRB   , 0.8900_JPRB    /)
1068
1069IF (KLEV == 19) THEN
1070  RAOVLP=ZAOVLP(1)
1071  RBOVLP=ZBOVLP(1)
1072ELSEIF (KLEV == 31) THEN
1073  RAOVLP=ZAOVLP(2)
1074  RBOVLP=ZBOVLP(2)
1075ELSEIF (KLEV == 60) THEN
1076  RAOVLP=ZAOVLP(3)
1077  RBOVLP=ZBOVLP(3)
1078ELSE 
1079  RAOVLP=ZAOVLP(3)
1080  RBOVLP=ZBOVLP(3)
1081ENDIF 
1082
1083!     ------------------------------------------------------------------
1084
1085IF (LHOOK) CALL DR_HOOK('SUCLOPN',1,ZHOOK_HANDLE)
1086END SUBROUTINE SUCLOPN
Note: See TracBrowser for help on using the repository browser.