LMDZ
susrtcop.F90
Go to the documentation of this file.
1 SUBROUTINE 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 
54 USE parkind1 ,ONLY : jpim ,jprb
55 USE yomhook ,ONLY : lhook, dr_hook
56 
57 USE 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 
68 IMPLICIT NONE
69 
70 INTEGER(KIND=JPIM) :: JNU
71 
72 ! -----------------------------------------------------------------
73 
74 REAL(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)
87 REAL(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 
116 IF (lhook) CALL dr_hook('SUSRTCOP',0,zhook_handle)
117 zebcua14 = (/ &
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 
122 zebcub14 = (/ &
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 
127 zebcuc14 = (/ &
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 
132 zebcud14 = (/ &
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 
137 zebcue14 = (/ &
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 
142 zebcuf14 = (/ &
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 
149 zyfwca14 = (/ &
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 
154 zyfwcb14 = (/ &
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 
159 zyfwcc14 = (/ &
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 
164 zyfwcd14 = (/ &
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 
169 zyfwce14 = (/ &
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 
174 zyfwcf14 = (/ &
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 
181 zaswca14 = (/ &
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 
186 zaswcb14 = (/ &
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 
191 zaswcc14 = (/ &
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 
196 zaswcd14 = (/ &
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 
201 zaswce14 = (/ &
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 
206 zaswcf14 = (/ &
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 
213 zsushe14 = (/ &
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 /)
217 zsushf14 = (/ &
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 /)
221 zsushh14 = (/ &
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 /)
225 zsushk14 = (/ &
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 /)
229 zsusha14 = (/ &
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 /)
233 zsushg14 = (/ &
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 
240 zflaa014 = (/ &
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 /)
244 zflaa114 = (/ &
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 /)
248 zflbb014 = (/ &
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 /)
252 zflbb114 = (/ &
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 /)
256 zflbb214 = (/ &
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 /)
260 zflbb314 = (/ &
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 /)
264 zflcc014 = (/ &
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 /)
268 zflcc114 = (/ &
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 /)
272 zflcc214 = (/ &
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 /)
276 zflcc314 = (/ &
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 /)
280 zfldd014 = (/ &
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 /)
284 zfldd114 = (/ &
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 /)
288 zfldd214 = (/ &
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 /)
292 zfldd314 = (/ &
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 
299 zfuaa014 = (/ &
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 
304 zfuaa114 = (/ &
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 
309 zfubb014 = (/ &
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 
314 zfubb114 = (/ &
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 
319 zfubb214 = (/ &
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 
324 zfubb314 = (/ &
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 
329 zfucc014 = (/ &
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 
334 zfucc114 = (/ &
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 
339 zfucc214 = (/ &
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 
344 zfucc314 = (/ &
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 
354 DO 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 
409 ENDDO
410 
411 ! ------------------------------------------------------------------
412 
413 IF (lhook) CALL dr_hook('SUSRTCOP',1,zhook_handle)
414 END SUBROUTINE susrtcop
415 
yoesrtcop::rsyfwe
real(kind=jprb), dimension(14) rsyfwe
Definition: yoesrtcop.F90:21
yoesrtcop::rsaswa
real(kind=jprb), dimension(14) rsaswa
Definition: yoesrtcop.F90:31
yoesrtcop::rsfld3
real(kind=jprb), dimension(14) rsfld3
Definition: yoesrtcop.F90:62
yoesrtcop::rsssif
real(kind=jprb), dimension(14) rsssif
Definition: yoesrtcop.F90:65
yoesrtcop::rsfub0
real(kind=jprb), dimension(14) rsfub0
Definition: yoesrtcop.F90:40
yoesrtcop::rsfub2
real(kind=jprb), dimension(14) rsfub2
Definition: yoesrtcop.F90:42
yoesrtcop::rsfld1
real(kind=jprb), dimension(14) rsfld1
Definition: yoesrtcop.F90:60
yoesrtcop::rsfub1
real(kind=jprb), dimension(14) rsfub1
Definition: yoesrtcop.F90:41
yoesrtcop::rsfld0
real(kind=jprb), dimension(14) rsfld0
Definition: yoesrtcop.F90:59
yoesrtcop::rsflb3
real(kind=jprb), dimension(14) rsflb3
Definition: yoesrtcop.F90:54
yoesrtcop::rsssig
real(kind=jprb), dimension(14) rsssig
Definition: yoesrtcop.F90:69
yoesrtcop::rsflb0
real(kind=jprb), dimension(14) rsflb0
Definition: yoesrtcop.F90:51
yoesrtcop::rsecid
real(kind=jprb), dimension(14) rsecid
Definition: yoesrtcop.F90:27
yoesrtcop::rsflc0
real(kind=jprb), dimension(14) rsflc0
Definition: yoesrtcop.F90:55
yoesrtcop::rsaswb
real(kind=jprb), dimension(14) rsaswb
Definition: yoesrtcop.F90:32
yoesrtcop::rsaswe
real(kind=jprb), dimension(14) rsaswe
Definition: yoesrtcop.F90:35
yoesrtcop::rsfuc0
real(kind=jprb), dimension(14) rsfuc0
Definition: yoesrtcop.F90:44
yoesrtcop::rsfub3
real(kind=jprb), dimension(14) rsfub3
Definition: yoesrtcop.F90:43
susrtcop
subroutine susrtcop
Definition: susrtcop.F90:2
parkind1::jprb
integer, parameter jprb
Definition: parkind1.F90:31
yoesrtcop::rsfua0
real(kind=jprb), dimension(14) rsfua0
Definition: yoesrtcop.F90:38
yoesrtcop::rsssik
real(kind=jprb), dimension(14) rsssik
Definition: yoesrtcop.F90:67
yoesrtcop::rsecic
real(kind=jprb), dimension(14) rsecic
Definition: yoesrtcop.F90:26
yoesrtcop::rsfuc1
real(kind=jprb), dimension(14) rsfuc1
Definition: yoesrtcop.F90:45
yoesrtcop::rsyfwc
real(kind=jprb), dimension(14) rsyfwc
Definition: yoesrtcop.F90:19
yoesrtcop::rsssie
real(kind=jprb), dimension(14) rsssie
Definition: yoesrtcop.F90:64
yoesrtcop::rsyfwf
real(kind=jprb), dimension(14) rsyfwf
Definition: yoesrtcop.F90:22
yoesrtcop::rsflc3
real(kind=jprb), dimension(14) rsflc3
Definition: yoesrtcop.F90:58
yoesrtcop::rsflb2
real(kind=jprb), dimension(14) rsflb2
Definition: yoesrtcop.F90:53
yoesrtcop::rsecib
real(kind=jprb), dimension(14) rsecib
Definition: yoesrtcop.F90:25
yomhook::lhook
logical lhook
Definition: yomhook.F90:12
yoesrtcop::rsflc1
real(kind=jprb), dimension(14) rsflc1
Definition: yoesrtcop.F90:56
yoesrtcop::rsecia
real(kind=jprb), dimension(14) rsecia
Definition: yoesrtcop.F90:24
yoesrtcop::rsflc2
real(kind=jprb), dimension(14) rsflc2
Definition: yoesrtcop.F90:57
yoesrtcop::rsaswd
real(kind=jprb), dimension(14) rsaswd
Definition: yoesrtcop.F90:34
yomhook::dr_hook
subroutine dr_hook(CDNAME, KSWITCH, PKEY)
Definition: yomhook.F90:17
yoesrtcop::rsflb1
real(kind=jprb), dimension(14) rsflb1
Definition: yoesrtcop.F90:52
yoesrtcop::rsfuc3
real(kind=jprb), dimension(14) rsfuc3
Definition: yoesrtcop.F90:47
yoesrtcop::rsaswc
real(kind=jprb), dimension(14) rsaswc
Definition: yoesrtcop.F90:33
yoesrtcop::rsfla0
real(kind=jprb), dimension(14) rsfla0
Definition: yoesrtcop.F90:49
yoesrtcop::rsfld2
real(kind=jprb), dimension(14) rsfld2
Definition: yoesrtcop.F90:61
parkind1::jpim
integer, parameter jpim
Definition: parkind1.F90:13
yoesrtcop::rsssih
real(kind=jprb), dimension(14) rsssih
Definition: yoesrtcop.F90:66
yoesrtcop::rsyfwb
real(kind=jprb), dimension(14) rsyfwb
Definition: yoesrtcop.F90:18
yoesrtcop::rsyfwd
real(kind=jprb), dimension(14) rsyfwd
Definition: yoesrtcop.F90:20
yoesrtcop::rsaswf
real(kind=jprb), dimension(14) rsaswf
Definition: yoesrtcop.F90:36
yoesrtcop::rsfuc2
real(kind=jprb), dimension(14) rsfuc2
Definition: yoesrtcop.F90:46
yoesrtcop::rsyfwa
real(kind=jprb), dimension(14) rsyfwa
Definition: yoesrtcop.F90:17
yoesrtcop::rsecie
real(kind=jprb), dimension(14) rsecie
Definition: yoesrtcop.F90:28
yoesrtcop::rsssia
real(kind=jprb), dimension(14) rsssia
Definition: yoesrtcop.F90:68
yoesrtcop::rsecif
real(kind=jprb), dimension(14) rsecif
Definition: yoesrtcop.F90:29
yoesrtcop::rsfla1
real(kind=jprb), dimension(14) rsfla1
Definition: yoesrtcop.F90:50
yoesrtcop::rsfua1
real(kind=jprb), dimension(14) rsfua1
Definition: yoesrtcop.F90:39