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 |