GCC Code Coverage Report


Directory: ./
File: rad/suclopn.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 363 460 78.9%
Branches: 228 244 93.4%

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