source: LMDZ5/branches/IPSLCM5A2.1/libf/phylmd/rrtm/susrtcop.F90 @ 5443

Last change on this file since 5443 was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM

Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence
  • 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: 19.0 KB
Line 
1SUBROUTINE SUSRTCOP
2
3!**** *SUSRTCOP*  - INITIALIZE COMMON YOESRTCOP
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOESRTCOP, WITH CLOUD OPTICAL PARAMETERS
8
9!**   INTERFACE.
10!     ----------
11!        *CALL*  SUSRTCOP
12!        FROM *SUECRAD*
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        COMMON YOESRTCOP
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 : 03-03-06
50!        M.Hamrud      01-Oct-2003 CY28 Cleaning
51
52!     ------------------------------------------------------------------
53
54USE PARKIND1  ,ONLY : JPIM     ,JPRB
55USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
56
57USE YOESRTCOP , ONLY : &
58 & RSYFWA   ,RSYFWB   ,RSYFWC   ,RSYFWD   ,RSYFWE   ,RSYFWF  &
59 & , RSECIA   ,RSECIB   ,RSECIC   ,RSECID   ,RSECIE   ,RSECIF &
60 & , RSASWA   ,RSASWB   ,RSASWC   ,RSASWD   ,RSASWE   ,RSASWF &
61 & , RSSSIE   ,RSSSIF   ,RSSSIH   ,RSSSIK   ,RSSSIA   ,RSSSIG &
62 & , RSFUA0   ,RSFUA1   ,RSFUB0   ,RSFUB1   ,RSFUB2   ,RSFUB3 &
63 & , RSFUC0   ,RSFUC1   ,RSFUC2   ,RSFUC3   &
64 & , RSFLA0   ,RSFLA1   ,RSFLB0   ,RSFLB1   ,RSFLB2   ,RSFLB3 &
65 & , RSFLC0   ,RSFLC1   ,RSFLC2   ,RSFLC3   ,RSFLD0   ,RSFLD1 &
66 & , RSFLD2   ,RSFLD3 
67
68IMPLICIT NONE
69
70INTEGER(KIND=JPIM) ::  JNU
71
72!     -----------------------------------------------------------------
73
74REAL(KIND=JPRB) :: ZEBCUA14(14)  ,ZEBCUB14(14)  ,ZEBCUC14(14)  ,ZEBCUD14(14)&
75 & ,  ZEBCUE14(14)  ,ZEBCUF14(14)  ,ZYFWCA14(14)  ,ZYFWCB14(14)&
76 & ,  ZYFWCC14(14)  ,ZYFWCD14(14)  ,ZYFWCE14(14)  ,ZYFWCF14(14)&
77 & ,  ZASWCA14(14)  ,ZASWCB14(14)  ,ZASWCC14(14)  ,ZASWCD14(14)&
78 & ,  ZASWCE14(14)  ,ZASWCF14(14)  ,ZSUSHE14(14)  ,ZSUSHF14(14)&
79 & ,  ZSUSHH14(14)  ,ZSUSHK14(14)  ,ZSUSHA14(14)  ,ZSUSHG14(14)&
80 & ,  ZFLAA014(14)  ,ZFLAA114(14)  ,ZFLBB014(14)  ,ZFLBB114(14)&
81 & ,  ZFLBB214(14)  ,ZFLBB314(14)  ,ZFLCC014(14)  ,ZFLCC114(14)&
82 & ,  ZFLCC214(14)  ,ZFLCC314(14)  ,ZFLDD014(14)  ,ZFLDD114(14)&
83 & ,  ZFLDD214(14)  ,ZFLDD314(14)&
84 & ,  ZFUAA014(14)  ,ZFUAA114(14)  ,ZFUBB014(14)  ,ZFUBB114(14)&
85 & ,  ZFUBB214(14)  ,ZFUBB314(14)  ,ZFUCC014(14)  ,ZFUCC114(14)&
86 & ,  ZFUCC214(14)  ,ZFUCC314(14) 
87REAL(KIND=JPRB) :: ZHOOK_HANDLE
88!     -----------------------------------------------------------------
89
90!*          1.    SHORTWAVE CLOUD OPTICAL PROPERTIES
91!                 ----------------------------------
92
93!     ------------------------------------------------------------------
94
95!*          1.1    FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW
96!                  ---------------------------------------
97
98! SW : 14 spectral intervals
99!  3.846 -  3.077
100!  3.077 -  2.500
101!  2.500 -  2.150
102!  2.150 -  1.942
103!  1.942 -  1.626
104!  1.626 -  1.299
105!  1.299 -  1.242
106!  1.242 -  0.7782
107!  0.7782-  0.6250
108!  0.6250-  0.4415
109!  0.4415-  0.3448
110!  0.3448-  0.2632
111!  0.2632-  0.2000
112! 12.195 -  3.846
113
114!* Ice cloud properties - crystal: remapped from Ebert and Curry, 1992
115
116IF (LHOOK) CALL DR_HOOK('SUSRTCOP',0,ZHOOK_HANDLE)
117ZEBCUA14 = (/ &
118 & 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB &
119 & , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB , 3.448E-03_JPRB &
120 & , 3.448E-03_JPRB , 3.448E-03_JPRB /) 
121
122ZEBCUB14 = (/ &
123 & 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     &
124 & , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     , 2.431_JPRB     &
125 & , 2.431_JPRB     , 2.431_JPRB     /) 
126
127ZEBCUC14 = (/ &
128 & 0.46658_JPRB   , 0.46658_JPRB   , 0.03779_JPRB   , 0.03779_JPRB   , 0.01240_JPRB   , 0.01240_JPRB   &
129 & , 0.00011_JPRB   , 0.00011_JPRB   , 0.00001_JPRB   , 0.00001_JPRB   , 0.00001_JPRB   , 0.00001_JPRB   &
130 & , 0.00001_JPRB   , 0.46658_JPRB   /) 
131
132ZEBCUD14 = (/ &
133 & 2.050E-05_JPRB , 2.050E-05_JPRB , 1.284E-03_JPRB , 1.284E-03_JPRB , 6.867E-04_JPRB , 6.867E-04_JPRB &
134 & , 1.405E-05_JPRB , 1.405E-05_JPRB , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB          &
135 & , 0.0_JPRB         , 2.050E-05_JPRB /) 
136
137ZEBCUE14 = (/ &
138 & 0.9595_JPRB    , 0.9595_JPRB    , 0.8172_JPRB    , 0.8172_JPRB    , 0.7865_JPRB    , 0.7865_JPRB    &
139 & , 0.7730_JPRB    , 0.7730_JPRB    , 0.7661_JPRB    , 0.7661_JPRB    , 0.7661_JPRB    , 0.7661_JPRB    &
140 & , 0.7661_JPRB    , 0.9595_JPRB    /) 
141
142ZEBCUF14 = (/ &
143 & 1.076E-04_JPRB , 1.076E-04_JPRB , 7.463E-04_JPRB , 7.463E-04_JPRB , 7.204E-04_JPRB , 7.204E-04_JPRB &
144 & , 5.665E-04_JPRB , 5.665E-04_JPRB , 5.851E-04_JPRB , 5.851E-04_JPRB , 5.851E-04_JPRB , 5.851E-04_JPRB &
145 & , 5.851E-04_JPRB , 1.076E-04_JPRB /) 
146
147!* Water cloud properties - from Fouquart (1987)
148
149ZYFWCA14 = (/ &
150 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
151 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
152 & , 0.0_JPRB         , 0.0_JPRB         /) 
153 
154ZYFWCB14 = (/ &
155 & 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB      , 1.5_JPRB        &
156 & , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB       , 1.5_JPRB      , 1.5_JPRB        &
157 & , 1.5_JPRB       , 1.5_JPRB       /) 
158 
159ZYFWCC14 = (/ &
160 & 0.9988_JPRB    , 0.9988_JPRB    , 0.9988_JPRB    , 0.9988_JPRB    , 0.9988_JPRB    , 0.9988_JPRB    &
161 & , 0.9988_JPRB    , 0.9988_JPRB    , 0.9999_JPRB    , 0.9999_JPRB    , 0.9999_JPRB    , 0.9999_JPRB    &
162 & , 0.9999_JPRB    , 0.9988_JPRB    /) 
163
164ZYFWCD14 = (/ &
165 & 2.500E-03_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB , 2.500E-03_JPRB &
166 & , 2.500E-03_JPRB , 2.500E-03_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB , 5.000E-04_JPRB &
167 & , 5.000E-04_JPRB , 2.500E-03_JPRB /) 
168
169ZYFWCE14 = (/ &
170 & 0.05_JPRB      , 0.05_JPRB      , 0.05_JPRB      , 0.05_JPRB      , 0.05_JPRB      , 0.05_JPRB      &
171 & , 0.05_JPRB      , 0.05_JPRB      , 0.5_JPRB         , 0.5_JPRB         , 0.5_JPRB         , 0.5_JPRB               &
172 & , 0.5_JPRB         , 0.05_JPRB      /) 
173
174ZYFWCF14 = (/ &
175 & 0.910_JPRB     , 0.910_JPRB     , 0.910_JPRB     , 0.910_JPRB     , 0.910_JPRB     , 0.910_JPRB     &
176 & , 0.910_JPRB     , 0.910_JPRB     , 0.865_JPRB     , 0.865_JPRB     , 0.865_JPRB     , 0.865_JPRB      &
177 & , 0.865_JPRB     , 0.910_JPRB     /) 
178
179!* Water cloud properties - from Slingo (1989)
180
181ZASWCA14 = (/ &
182 & -1.023_JPRB     , 1.950_JPRB     , 1.850_JPRB     , 1.970_JPRB     , 1.970_JPRB     , 2.463_JPRB     &
183 & , 2.551_JPRB     , 2.622_JPRB     , 2.895_JPRB     , 2.672_JPRB     , 2.801_JPRB    , 2.944_JPRB     &
184 & , 3.094_JPRB     ,-1.023_JPRB     /) 
185
186ZASWCB14 = (/ &
187 & 1.933_JPRB     , 1.540_JPRB     , 1.556_JPRB     , 1.501_JPRB     , 1.501_JPRB     , 1.420_JPRB     &
188 & , 1.401_JPRB     , 1.362_JPRB     , 1.315_JPRB     , 1.320_JPRB     , 1.293_JPRB     , 1.270_JPRB     &
189 & , 1.252_JPRB     , 1.933_JPRB     /) 
190   
191ZASWCC14 = (/ &
192 & 2.50E-02_JPRB  , 4.49E-01_JPRB  , 1.90E-04_JPRB  , 1.20E-03_JPRB  , 1.20E-03_JPRB  , 2.40E-04_JPRB  &
193 & , 6.20E-05_JPRB  , 3.30E-06_JPRB  ,-1.20E-07_JPRB  , 0.0_JPRB         , 1.00E-06_JPRB  ,-6.50E-07_JPRB  &
194 & , 7.90E-07_JPRB  , 2.50E-02_JPRB  /) 
195   
196ZASWCD14 = (/ &
197 & 1.22E-02_JPRB  , 1.54E-03_JPRB  , 2.54E-03_JPRB  , 2.16E-03_JPRB  , 2.16E-03_JPRB  , 8.56E-04_JPRB  &
198 & , 2.60E-04_JPRB  , 2.80E-06_JPRB  , 4.40E-07_JPRB  , 0.0_JPRB         , 0.0_JPRB         , 4.33E-07_JPRB  &
199 & , 3.69E-07_JPRB  , 1.22E-02_JPRB  /) 
200
201ZASWCE14 = (/ &
202 & 0.726_JPRB     , 0.831_JPRB     , 0.769_JPRB     , 0.740_JPRB     , 0.740_JPRB     , 0.754_JPRB     &
203 & , 0.773_JPRB     , 0.806_JPRB     , 0.818_JPRB     , 0.828_JPRB     , 0.836_JPRB     , 0.841_JPRB     &
204 & , 0.844_JPRB     , 0.726_JPRB     /) 
205
206ZASWCF14 = (/ &
207 & 6.652_JPRB     , 6.102_JPRB     , 5.171_JPRB     , 7.469_JPRB     , 7.469_JPRB     , 6.555_JPRB     &
208 & , 5.405_JPRB     , 3.355_JPRB     , 2.989_JPRB     , 2.467_JPRB     , 2.153_JPRB     , 1.680_JPRB     &
209 & , 1.558_JPRB     , 6.652_JPRB     /) 
210
211!* Ice cloud properties - from Sun and Shine (1995)
212
213ZSUSHE14 = (/ &
214 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
215 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
216 & , 0.0_JPRB         , 0.0_JPRB         /) 
217ZSUSHF14 = (/ &
218 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
219 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
220 & , 0.0_JPRB         , 0.0_JPRB         /) 
221ZSUSHH14 = (/ &
222 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
223 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
224 & , 0.0_JPRB         , 0.0_JPRB         /) 
225ZSUSHK14 = (/ &
226 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
227 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
228 & , 0.0_JPRB         , 0.0_JPRB         /) 
229ZSUSHA14 = (/ &
230 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
231 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
232 & , 0.0_JPRB         , 0.0_JPRB         /) 
233ZSUSHG14 = (/ &
234 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
235 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
236 & , 0.0_JPRB         , 0.0_JPRB         /) 
237
238!* Ice cloud properties - from Fu and Liou (1993)
239
240ZFLAA014 = (/ &
241 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
242 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
243 & , 0.0_JPRB         , 0.0_JPRB         /) 
244ZFLAA114 = (/ &
245 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
246 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
247 & , 0.0_JPRB         , 0.0_JPRB         /) 
248ZFLBB014 = (/ &
249 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
250 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
251 & , 0.0_JPRB         , 0.0_JPRB         /) 
252ZFLBB114 = (/ &
253 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
254 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
255 & , 0.0_JPRB         , 0.0_JPRB         /) 
256ZFLBB214 = (/ &
257 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
258 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
259 & , 0.0_JPRB         , 0.0_JPRB         /) 
260ZFLBB314 = (/ &
261 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
262 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
263 & , 0.0_JPRB         , 0.0_JPRB         /) 
264ZFLCC014 = (/ &
265 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
266 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
267 & , 0.0_JPRB         , 0.0_JPRB         /) 
268ZFLCC114 = (/ &
269 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
270 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
271 & , 0.0_JPRB         , 0.0_JPRB         /) 
272ZFLCC214 = (/ &
273 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
274 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
275 & , 0.0_JPRB         , 0.0_JPRB         /) 
276ZFLCC314 = (/ &
277 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
278 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
279 & , 0.0_JPRB         , 0.0_JPRB         /) 
280ZFLDD014 = (/ &
281 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
282 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
283 & , 0.0_JPRB         , 0.0_JPRB         /) 
284ZFLDD114 = (/ &
285 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
286 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
287 & , 0.0_JPRB         , 0.0_JPRB         /) 
288ZFLDD214 = (/ &
289 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
290 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
291 & , 0.0_JPRB         , 0.0_JPRB         /) 
292ZFLDD314 = (/ &
293 & 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
294 & , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB         , 0.0_JPRB        , 0.0_JPRB          &
295 & , 0.0_JPRB         , 0.0_JPRB         /) 
296
297!* Ice cloud properties - from Fu (1996)
298
299ZFUAA014 = (/ &
300 & 1.87598E-04_JPRB, 2.97295E-04_JPRB, 4.89477E-04_JPRB,-8.37325E-06_JPRB,-8.37325E-06_JPRB,-8.05155E-04_JPRB  &
301 & , 6.51659E-05_JPRB, 8.10443E-05_JPRB, 1.61983E-04_JPRB,-9.45458E-05_JPRB,-2.58858E-04_JPRB,-2.66955E-04_JPRB &
302 & ,-2.36447E-04_JPRB,-2.54823E-04_JPRB /) 
303
304ZFUAA114 = (/ &
305 & 2.51396E+00_JPRB, 2.48895E+00_JPRB, 2.48776E+00_JPRB, 2.52504E+00_JPRB, 2.52504E+00_JPRB, 2.57600E+00_JPRB &
306 & , 2.51660E+00_JPRB, 2.51619E+00_JPRB, 2.50746E+00_JPRB, 2.52061E+00_JPRB, 2.53815E+00_JPRB, 2.54179E+00_JPRB &
307 & , 2.53817E+00_JPRB, 2.52909E+00_JPRB /) 
308
309ZFUBB014 = (/ &
310 & 1.96793E-01_JPRB, 4.64416E-01_JPRB, 5.83469E-04_JPRB, 2.53234E-03_JPRB, 2.53234E-03_JPRB,-2.85518E-05_JPRB &
311 & ,-1.48012E-07_JPRB,-1.57963E-07_JPRB,-7.78001E-08_JPRB, 5.08447E-07_JPRB,-1.98529E-07_JPRB,-1.00570E-07_JPRB &
312 & ,-2.69916E-07_JPRB, 2.60155E-01_JPRB /)   
313
314ZFUBB114 = (/ &
315 & 5.75235E-03_JPRB, 2.04716E-05_JPRB, 1.18127E-03_JPRB, 1.75078E-03_JPRB, 1.75078E-03_JPRB, 1.71993E-03_JPRB &
316 & , 9.02355E-05_JPRB, 1.72475E-06_JPRB, 2.53360E-07_JPRB, 2.73206E-08_JPRB, 9.39480E-08_JPRB, 1.60441E-07_JPRB &
317 & , 2.12909E-07_JPRB, 5.45547E-03_JPRB /) 
318
319ZFUBB214 = (/ &
320 & -5.29220E-05_JPRB,-4.60375E-07_JPRB,-3.40011E-06_JPRB,-8.00994E-06_JPRB,-8.00994E-06_JPRB,-7.43697E-06_JPRB &
321 & ,-1.98190E-08_JPRB, 9.02156E-11_JPRB,-1.15489E-10_JPRB, 4.96553E-11_JPRB,-2.54540E-10_JPRB,-2.05663E-10_JPRB &
322 & ,-2.65397E-10_JPRB,-5.58760E-05_JPRB /) 
323
324ZFUBB314 = (/ &
325 & 1.76618E-07_JPRB, 2.03198E-09_JPRB, 8.78549E-09_JPRB, 2.31309E-08_JPRB, 2.31309E-08_JPRB, 2.09647E-08_JPRB &
326 & , 4.01914E-11_JPRB,-3.79423E-13_JPRB, 4.65084E-13_JPRB,-1.86001E-13_JPRB, 1.10876E-12_JPRB, 8.88595E-13_JPRB &
327 & , 1.12983E-12_JPRB, 1.97086E-07_JPRB /) 
328
329ZFUCC014 = (/ &
330 & 0.759183_JPRB   , 0.919599_JPRB   , 0.775916_JPRB   , 0.758748_JPRB   , 0.758748_JPRB   , 0.752528_JPRB    &
331 & , 0.751277_JPRB   , 0.752318_JPRB   , 0.749693_JPRB   , 0.749856_JPRB   , 0.743546_JPRB   , 0.737809_JPRB    &
332 & , 0.733260_JPRB   , 0.799084_JPRB   /) 
333   
334ZFUCC114 = (/ &
335 & 4.93765E-03_JPRB, 5.03025E-04_JPRB, 1.74517E-03_JPRB, 2.02709E-03_JPRB, 2.02709E-03_JPRB, 1.95748E-03_JPRB &
336 & , 1.29824E-03_JPRB, 1.04224E-03_JPRB, 1.05446E-03_JPRB, 8.89161E-04_JPRB, 9.08674E-04_JPRB, 8.97515E-04_JPRB &
337 & , 9.18317E-04_JPRB, 4.81706E-03_JPRB /) 
338
339ZFUCC214 = (/ &
340 & -4.84059E-05_JPRB,-5.74771E-06_JPRB,-9.21314E-06_JPRB,-1.17029E-05_JPRB,-1.17029E-05_JPRB,-1.02495E-05_JPRB &
341 & ,-4.99075E-06_JPRB,-2.26618E-06_JPRB,-2.32576E-06_JPRB,-3.49578E-07_JPRB,-4.65326E-07_JPRB,-2.17099E-07_JPRB &
342 & ,-4.22974E-07_JPRB,-5.13220E-05_JPRB /) 
343
344ZFUCC314 = (/ &
345 & 1.65801E-07_JPRB, 2.01731E-08_JPRB, 2.15003E-08_JPRB, 2.95195E-08_JPRB, 2.95195E-08_JPRB, 2.35479E-08_JPRB &
346 & , 6.33757E-09_JPRB,-3.68283E-09_JPRB,-3.58307E-09_JPRB, 1.09913E-08_JPRB,-1.05786E-08_JPRB,-1.16090E-08_JPRB &
347 & ,-1.07976E-08_JPRB, 1.84420E-07_JPRB /) 
348
349!     ------------------------------------------------------------------
350!     ------------------------------------------------------------------
351
352! SW : absorption coefficients
353
354DO JNU=1,14
355  RSASWA(JNU)=ZASWCA14(JNU)*1.E-02_JPRB
356  RSASWB(JNU)=ZASWCB14(JNU)
357  RSASWC(JNU)=ZASWCC14(JNU)
358  RSASWD(JNU)=ZASWCD14(JNU)
359  RSASWE(JNU)=ZASWCE14(JNU)
360  RSASWF(JNU)=ZASWCF14(JNU)*1.E-03_JPRB
361
362  RSECIA(JNU)=ZEBCUA14(JNU)
363  RSECIB(JNU)=ZEBCUB14(JNU)
364  RSECIC(JNU)=ZEBCUC14(JNU)
365  RSECID(JNU)=ZEBCUD14(JNU)
366  RSECIE(JNU)=ZEBCUE14(JNU)
367  RSECIF(JNU)=ZEBCUF14(JNU)
368
369  RSYFWA(JNU)=ZYFWCA14(JNU)
370  RSYFWB(JNU)=ZYFWCB14(JNU)
371  RSYFWC(JNU)=ZYFWCC14(JNU)
372  RSYFWD(JNU)=ZYFWCD14(JNU)
373  RSYFWE(JNU)=ZYFWCE14(JNU)
374  RSYFWF(JNU)=ZYFWCF14(JNU)
375
376  RSSSIE(JNU)=ZSUSHE14(JNU)*1.E-02_JPRB
377  RSSSIF(JNU)=ZSUSHF14(JNU)*1.E-02_JPRB
378  RSSSIH(JNU)=ZSUSHH14(JNU)
379  RSSSIK(JNU)=ZSUSHK14(JNU)*1.E-01_JPRB
380  RSSSIA(JNU)=ZSUSHA14(JNU)*1.E-03_JPRB
381  RSSSIG(JNU)=ZSUSHG14(JNU)*1.E-01_JPRB
382   
383  RSFLA0(JNU)=ZFLAA014(JNU)
384  RSFLA1(JNU)=ZFLAA114(JNU)
385  RSFLB0(JNU)=ZFLBB014(JNU)
386  RSFLB1(JNU)=ZFLBB114(JNU)
387  RSFLB2(JNU)=ZFLBB214(JNU)
388  RSFLB3(JNU)=ZFLBB314(JNU)
389  RSFLC0(JNU)=ZFLCC014(JNU)
390  RSFLC1(JNU)=ZFLCC114(JNU)
391  RSFLC2(JNU)=ZFLCC214(JNU)
392  RSFLC3(JNU)=ZFLCC314(JNU)
393  RSFLD0(JNU)=ZFLDD014(JNU)
394  RSFLD1(JNU)=ZFLDD114(JNU)
395  RSFLD2(JNU)=ZFLDD214(JNU)
396  RSFLD3(JNU)=ZFLDD314(JNU)
397   
398  RSFUA0(JNU)=ZFUAA014(JNU)
399  RSFUA1(JNU)=ZFUAA114(JNU)
400  RSFUB0(JNU)=ZFUBB014(JNU)
401  RSFUB1(JNU)=ZFUBB114(JNU)
402  RSFUB2(JNU)=ZFUBB214(JNU)
403  RSFUB3(JNU)=ZFUBB314(JNU)
404  RSFUC0(JNU)=ZFUCC014(JNU)
405  RSFUC1(JNU)=ZFUCC114(JNU)
406  RSFUC2(JNU)=ZFUCC214(JNU)
407  RSFUC3(JNU)=ZFUCC314(JNU)
408   
409ENDDO
410
411!     ------------------------------------------------------------------
412
413IF (LHOOK) CALL DR_HOOK('SUSRTCOP',1,ZHOOK_HANDLE)
414END SUBROUTINE SUSRTCOP
415
Note: See TracBrowser for help on using the repository browser.