GCC Code Coverage Report

Directory: ./
File: rad/suswn.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 244 306 79.7%
Branches: 259 294 88.1%
Line Branch Exec Source
1 2 SUBROUTINE SUSWN (KTSW, KSW)
2
3 !**** *SUSW* - INITIALIZE COMMON YOESW
4
5 ! PURPOSE.
6 ! --------
7 ! INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS
8 ! NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES
9
10 !** INTERFACE.
11 ! ----------
12 ! *CALL* *SUSW
13
14 ! EXPLICIT ARGUMENTS :
15 ! --------------------
16 ! NONE
17
18 ! IMPLICIT ARGUMENTS :
19 ! --------------------
20 ! COMMON YOESW
21
22 ! METHOD.
23 ! -------
24 ! SEE DOCUMENTATION
25
26 ! EXTERNALS.
27 ! ----------
28
29 ! REFERENCE.
30 ! ----------
31 ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
32
33 ! AUTHOR.
34 ! -------
35 ! JEAN-JACQUES MORCRETTE *ECMWF*
36
37 ! MODIFICATIONS.
38 ! --------------
39 ! ORIGINAL : 88-12-15
40 ! 97-04-16 JJ Morcrette 2 and 4 interval spectral resolution
41 ! 00-10-24 JJ Morcrette sea-ice albedo revisited
42 ! 00-12-14 JJ Morcrette
43 ! and Ph.Dubuisson B.Bonnel 6 spectral interval resolution
44 ! 01-04-17 Ph.Dubuisson, B.Bonnel, JJ.Morcrette 6 sp.int.resolu'n
45 ! 01-06-28 B.Bonnel, JJ.Morcrette, Ph.Dubuisson Rayleigh (2/4/6)
46 ! 01-11-05 Ph.Dubuisson, JJMorcrette (new 2 intervals for TL/AD)
47 ! 03-04-24 JJ Morcrette SRTM with 14 spectral intervals
48 ! M.Hamrud 01-Oct-2003 CY28 Cleaning
49
50 ! ------------------------------------------------------------------
51
52 USE PARKIND1 ,ONLY : JPIM ,JPRB
53 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK
54
55 USE YOMCST , ONLY : RG
56 USE YOERAD , ONLY : NTSW ,LONEWSW ,LRAYL
57 USE YOESW , ONLY : APAD ,BPAD ,RRAY ,RSUN ,&
58 & RPDH1 ,RPDU1 ,RPNH ,RPNU ,RSWCE ,&
59 & RSWCP ,RTDH2O ,RTDUMG ,RTH2O ,RTUMG ,&
60 & D ,RROMA ,RROMB ,RRASY ,RHSRA ,&
61 & RHSRB ,RHSRC ,RHSRD ,RHSRE ,RHSRF ,&
62 & RHSRTA ,RHSRTB ,REXPO3 ,&
63 & RADJUST ,NEXPO3
64
65 IMPLICIT NONE
66
67 INTEGER(KIND=JPIM),INTENT(IN) :: KTSW
68 INTEGER(KIND=JPIM),INTENT(IN) :: KSW
69 INTEGER(KIND=JPIM) :: IMPSRTM(14)
70
71 ! ----------------------------------------------------------------
72 REAL(KIND=JPRB) :: ZAPAD2(2,3,7) , ZBPAD2(2,3,7) , ZD2(2,3)&
73 & , ZRAY2(2,6) , ZSUN2(2) , ZSWCE2(2) , ZSWCP2(2)&
74 & , ZROMA2(2) , ZROMB2(2) , ZRASY2(2)&
75 & , ZRA2(2) , ZRB2(2) , ZRC2(2)&
76 & , ZRD2(2) , ZRE2(2) , ZRF2(2)
77
78 REAL(KIND=JPRB) :: ZAPAD4(4,3,7) , ZBPAD4(4,3,7) , ZD4(4,3)&
79 & , ZRAY4(4,6) , ZSUN4(4) , ZSWCE4(4) , ZSWCP4(4)&
80 & , ZROMA4(4) , ZROMB4(4) , ZRASY4(4)&
81 & , ZRA4(4) , ZRB4(4) , ZRC4(4)&
82 & , ZRD4(4) , ZRE4(4) , ZRF4(4)
83
84 REAL(KIND=JPRB) :: ZAPAD6(6,3,7) , ZBPAD6(6,3,7) , ZD6(6,3)&
85 & , ZRAY6(6,6) , ZSUN6(6) , ZSWCE6(6) , ZSWCP6(6)&
86 & , ZROMA6(6) , ZROMB6(6) , ZRASY6(6)&
87 & , ZRA6(6) , ZRB6(6) , ZRC6(6)&
88 & , ZRD6(6) , ZRE6(6) , ZRF6(6)&
89 & , ZRBY6(6,6)
90
91
92 INTEGER(KIND=JPIM) :: JC3, JC6, JI, JJ, JW, I_K, ISW6
93
94 REAL(KIND=JPRB) :: ZH2O, ZPDHNIS, ZPDH2O, ZPDUNIS,&
95 & ZPDUMG, ZPRHNIS, ZPRH2O, &
96 & ZPRUNIS, ZPRUMG, ZRTO1, ZRTO2, ZTDHNIS, &
97 & ZTDUNIS, ZTHNIS, &
98 & ZTUNIS, ZUMG, ZADJUS2, ZADJUS4, ZADJUS6
99
100 REAL(KIND=JPRB) :: ZTDH2IS, ZTDU2IS, ZTH2IS, ZTU2IS, ZPDH2IS, ZPDU2IS, ZPRH2IS, ZPRU2IS
101 REAL(KIND=JPRB) :: ZTDH4IS, ZTDU4IS, ZTH4IS, ZTU4IS, ZPDH4IS, ZPDU4IS, ZPRH4IS, ZPRU4IS
102 REAL(KIND=JPRB) :: ZTDH6IS, ZTDU6IS, ZTH6IS, ZTU6IS, ZPDH6IS, ZPDU6IS, ZPRH6IS, ZPRU6IS
103 REAL(KIND=JPRB) :: ZHOOK_HANDLE
104
105 INTERFACE
106 SUBROUTINE ABOR1(CDTEXT)
107 CHARACTER(LEN=*) :: CDTEXT
108 END SUBROUTINE ABOR1
109 END INTERFACE
110
111 ! ----------------------------------------------------------------
112
113 !* 1. CLEAR-SKY ABSORPTION COEFFICIENTS FOR N SPECTRAL INTERVALS
114 ! --------------------------------------------------------
115
116 !* spectrally independent
117
118 !-- New coefficients
119
120
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2 IF (LHOOK) CALL DR_HOOK('SUSWN',0,ZHOOK_HANDLE)
121 ZTDHNIS = 0.450_JPRB
122 ZTDUNIS = 0.375_JPRB
123 ZTHNIS = 273._JPRB
124 ZTUNIS = 273._JPRB
125 ZPDHNIS = 0.90_JPRB
126 ZPDUNIS = 0.75_JPRB
127 ZPRHNIS = 101300._JPRB
128 ZPRUNIS = 101300._JPRB
129 ISW6=6
130
131 ! ----------------------------------------------------------------
132
133 !* 1.1 COEFFICIENTS FOR TWO SPECTRAL INTERVALS
134 ! ---------------------------------------
135
136 !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
137 ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
138 ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
139 ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
140
141 ZADJUS2 = 1.0_JPRB
142 ZTDH2IS = 0.450_JPRB
143 ZTDU2IS = 0.375_JPRB
144 ZTH2IS = 273._JPRB
145 ZTU2IS = 273._JPRB
146 ZPDH2IS = 0.90_JPRB
147 ZPDU2IS = 0.75_JPRB
148 ZPRH2IS = 101300._JPRB
149 ZPRU2IS = 101300._JPRB
150
151 !* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
152
153 2 ZSUN2(1) = 0.459760_JPRB
154
155 !-- New coefficients
156
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD2(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
157
158 ZAPAD2(1, 1, :) = (/&
159 & 0.1095121E+03_JPRB,&
160 & 0.8906647E+03_JPRB,&
161 & 0.1030601E+03_JPRB,&
162 & 0.000000000E-00_JPRB,&
163 & 0.000000000E-00_JPRB,&
164 & 0.000000000E-00_JPRB,&
165
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
166 ZAPAD2(1, 2, :) = (/&
167 & 0.5124541E+03_JPRB,&
168 & 0.3478800E+04_JPRB,&
169 & 0.1244272E+04_JPRB,&
170 & 0.000000000E-00_JPRB,&
171 & 0.000000000E-00_JPRB,&
172 & 0.000000000E-00_JPRB,&
173
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
174 ZAPAD2(1, 3, :) = (/&
175 & 0.4766857E-04_JPRB,&
176 & 0.4019806E-02_JPRB,&
177 & 0.1936247E+00_JPRB,&
178 & 0.1886522E+02_JPRB,&
179 & 0.6498215E+02_JPRB,&
180 & 0.1799387E+02_JPRB,&
181
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
182
183 ZBPAD2(1, 1, :) = (/&
184 & 0.1095121E+03_JPRB,&
185 & 0.8919103E+03_JPRB,&
186 & 0.1142218E+03_JPRB,&
187 & 0.100000000E+01_JPRB,&
188 & 0.000000000E-00_JPRB,&
189 & 0.000000000E-00_JPRB,&
190
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
191 ZBPAD2(1, 2, :) = (/&
192 & 0.5124541E+03_JPRB,&
193 & 0.3499516E+04_JPRB,&
194 & 0.1278186E+04_JPRB,&
195 & 0.100000000E+01_JPRB,&
196 & 0.000000000E-00_JPRB,&
197 & 0.000000000E-00_JPRB,&
198
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
199 ZBPAD2(1, 3, :) = (/&
200 & 0.4766857E-04_JPRB,&
201 & 0.4153343E-02_JPRB,&
202 & 0.1925810E+00_JPRB,&
203 & 0.1938938E+02_JPRB,&
204 & 0.6918191E+02_JPRB,&
205 & 0.2278795E+02_JPRB,&
206
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.100000000E+01_JPRB/)
207
208 ZRAY2(1,:)= (/&
209 & .389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,&
210
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & .659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/)
211
212 !-- Original coefficients of ECMWF model (before 27 June 2000)
213
214 !ZTDH2IS = 0.450_JPRB
215 !ZTDU2IS = 0.375_JPRB
216 !ZTH2IS = 296._JPRB
217 !ZTU2IS = 296._JPRB
218 !ZPDH2IS = 0.90_JPRB
219 !ZPDU2IS = 0.75_JPRB
220 !ZPRH2IS = 100000._JPRB
221 !ZPRU2IS = 100000._JPRB
222
223 !ZD2(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
224
225 !ZAPAD2(1, 1, :) = (/&
226 ! &0.184678379E+06_JPRB,&
227 ! &0.553080884E+05_JPRB,&
228 ! &0.248143712E+04_JPRB,&
229 ! &0.000000000E-00_JPRB,&
230 ! &0.000000000E-00_JPRB,&
231 ! &0.000000000E-00_JPRB,&
232 ! &0.000000000E-00_JPRB/)
233 !ZAPAD2(1, 2, :) = (/&
234 ! &0.715303869E+01_JPRB,&
235 ! &0.219386847E+03_JPRB,&
236 ! &0.830001089E+03_JPRB,&
237 ! &0.000000000E-00_JPRB,&
238 ! &0.000000000E-00_JPRB,&
239 ! &0.000000000E-00_JPRB,&
240 ! &0.000000000E-00_JPRB/)
241 !ZAPAD2(1, 3, :) = (/&
242 ! &0.925887084E-04_JPRB,&
243 ! &0.129353723E-01_JPRB,&
244 ! &0.800821928E+00_JPRB,&
245 ! &0.242715973E+02_JPRB,&
246 ! &0.878331486E+02_JPRB,&
247 ! &0.191559725E+02_JPRB,&
248 ! &0.000000000E+00_JPRB/)
249
250 !ZBPAD2(1, 1, :) = (/&
251 ! &0.184678379E+06_JPRB,&
252 ! &0.555188347E+05_JPRB,&
253 ! &0.253257443E+04_JPRB,&
254 ! &0.100000000E+01_JPRB,&
255 ! &0.000000000E-00_JPRB,&
256 ! &0.000000000E-00_JPRB,&
257 ! &0.000000000E-00_JPRB/)
258 !ZBPAD2(1, 2, :) = (/&
259 ! &0.715303869E+01_JPRB,&
260 ! &0.219441875E+03_JPRB,&
261 ! &0.831119997E+03_JPRB,&
262 ! &0.100000000E+01_JPRB,&
263 ! &0.000000000E-00_JPRB,&
264 ! &0.000000000E-00_JPRB,&
265 ! &0.000000000E-00_JPRB/)
266 !ZBPAD2(1, 3, :) = (/&
267 ! &0.925887084E-04_JPRB,&
268 ! &0.131812683E-01_JPRB,&
269 ! &0.812706117E+00_JPRB,&
270 ! &0.249863591E+02_JPRB,&
271 ! &0.931071925E+02_JPRB,&
272 ! &0.252233437E+02_JPRB,&
273 ! &0.100000000E+01_JPRB/)
274
275 ZRAY2(1,:)= (/&
276 & .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
277
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
278
279 2 ZSWCE2(1) = 0.0_JPRB
280 2 ZSWCP2(1) = 0.0_JPRB
281
282 ! ----------------------------------------------------------------
283
284 !* Near-Infrared (0.69 - 4.0 Microns) is not sub-divided
285
286 ! ----------------------------------------------------------------
287
288 !* 0.68 - 4.00 Micron
289
290 2 ZSUN2(2) = 0.540240_JPRB
291
292 !-- New coefficients
293
294
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD2(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB/)
295
296 ZAPAD2(2, 1, :) = (/&
297 & 0.1345066E-03_JPRB,&
298 & 0.9276888E+00_JPRB,&
299 & 0.3449105E+03_JPRB,&
300 & 0.7143025E+04_JPRB,&
301 & 0.6531248E+04_JPRB,&
302 & 0.2210815E+03_JPRB,&
303
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
304 ZAPAD2(2, 2, :) = (/&
305 & 0.7776703E-05_JPRB,&
306 & 0.6851592E-02_JPRB,&
307 & 0.5533947E+00_JPRB,&
308 & 0.2538585E+01_JPRB,&
309 & 0.000000000E+00_JPRB,&
310 & 0.000000000E+00_JPRB,&
311
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
312 ZAPAD2(2, 3, :) = (/&
313 & 0.932012720E+02_JPRB,&
314 & 0.437866983E+03_JPRB,&
315 & 0.000000000E-00_JPRB,&
316 & 0.000000000E-00_JPRB,&
317 & 0.000000000E-00_JPRB,&
318 & 0.000000000E-00_JPRB,&
319
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
320
321 ZBPAD2(2, 1, :) = (/&
322 & 0.1345066E-03_JPRB,&
323 & 0.9297728E+00_JPRB,&
324 & 0.3523821E+03_JPRB,&
325 & 0.7931590E+04_JPRB,&
326 & 0.8979086E+04_JPRB,&
327 & 0.5050864E+03_JPRB,&
328
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.1000000E+01_JPRB/)
329 ZBPAD2(2, 2, :) = (/&
330 & 0.7776703E-05_JPRB,&
331 & 0.6976301E-02_JPRB,&
332 & 0.6074521E+00_JPRB,&
333 & 0.3517521E+01_JPRB,&
334 & 0.100000000E+01_JPRB,&
335 & 0.000000000E+00_JPRB,&
336
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
337 ZBPAD2(2, 3, :) = (/&
338 & 0.932012720E+02_JPRB,&
339 & 0.439078599E+03_JPRB,&
340 & 0.100000000E+01_JPRB,&
341 & 0.000000000E+00_JPRB,&
342 & 0.000000000E+00_JPRB,&
343 & 0.000000000E+00_JPRB,&
344
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
345
346 ZRAY2(2,:)= (/&
347 & .104635E-01_JPRB, .102535E+00_JPRB,-.536771E+00_JPRB,&
348
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .123640E+01_JPRB,-.128930E+01_JPRB, .502426E+00_JPRB/)
349
350 !-- Original coefficients of ECMWF model (before 27 June 2000)
351
352 !ZD2(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB/)
353
354 !ZAPAD2(2, 1, :) = (/&
355 ! &0.577160808E-06_JPRB,&
356 ! &0.681760991E-02_JPRB,&
357 ! &0.502416370E+01_JPRB,&
358 ! &0.209771270E+03_JPRB,&
359 ! &0.644493139E+03_JPRB,&
360 ! &0.120654163E+03_JPRB,&
361 ! &0.000000000E+00_JPRB/)
362 !ZAPAD2(2, 2, :) = (/&
363 ! &0.422874202E-09_JPRB,&
364 ! &0.257499477E-05_JPRB,&
365 ! &0.349570458E-02_JPRB,&
366 ! &0.917758876E+00_JPRB,&
367 ! &0.327669989E+02_JPRB,&
368 ! &0.874155115E+02_JPRB,&
369 ! &0.000000000E+00_JPRB/)
370 !ZAPAD2(2, 3, :) = (/&
371 ! &0.932012720E+02_JPRB,&
372 ! &0.437866983E+03_JPRB,&
373 ! &0.000000000E-00_JPRB,&
374 ! &0.000000000E-00_JPRB,&
375 ! &0.000000000E-00_JPRB,&
376 ! &0.000000000E-00_JPRB,&
377 ! &0.000000000E+00_JPRB/)
378
379 !ZBPAD2(2, 1, :) = (/&
380 ! &0.577160808E-06_JPRB,&
381 ! &0.682595589E-02_JPRB,&
382 ! &0.508797283E+01_JPRB,&
383 ! &0.223901000E+03_JPRB,&
384 ! &0.770640908E+03_JPRB,&
385 ! &0.186812480E+03_JPRB,&
386 ! &0.100000000E+01_JPRB/)
387 !ZBPAD2(2, 2, :) = (/&
388 ! &0.422874202E-09_JPRB,&
389 ! &0.257567242E-05_JPRB,&
390 ! &0.349912408E-02_JPRB,&
391 ! &0.920770481E+00_JPRB,&
392 ! &0.330990960E+02_JPRB,&
393 ! &0.903249595E+02_JPRB,&
394 ! &0.100000000E+01_JPRB/)
395 !ZBPAD2(2, 3, :) = (/&
396 ! &0.932012720E+02_JPRB,&
397 ! &0.439078599E+03_JPRB,&
398 ! &0.100000000E+01_JPRB,&
399 ! &0.000000000E+00_JPRB,&
400 ! &0.000000000E+00_JPRB,&
401 ! &0.000000000E+00_JPRB,&
402 ! &0.000000000E+00_JPRB/)
403
404 !ZRAY2(2,:)= (/&
405 ! &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,&
406 ! &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/)
407
408 2 ZSWCE2(2) = 0.0_JPRB
409 2 ZSWCP2(2) = 0.0_JPRB
410
411 !=====================================================================
412
413 !* 1.2 COEFFICIENTS FOR FOUR SPECTRAL INTERVALS
414 ! ----------------------------------------
415
416 !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
417 ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
418 ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
419 ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
420
421 ZADJUS4 = 1.0_JPRB
422
423 !* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
424
425 2 ZSUN4(1) = 0.459760_JPRB
426
427 !-- New coefficients
428
429 !ZD4(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
430
431 !ZAPAD4(1, 1, :) = (/&
432 ! &0.1095121E+03_JPRB,&
433 ! &0.8906647E+03_JPRB,&
434 ! &0.1030601E+03_JPRB,&
435 ! &0.0000000E-00_JPRB,&
436 ! &0.0000000E-00_JPRB,&
437 ! &0.0000000E-00_JPRB,&
438 ! &0.0000000E-00_JPRB/)
439 !ZAPAD4(1, 2, :) = (/&
440 ! &0.5124541E+03_JPRB,&
441 ! &0.3478800E+04_JPRB,&
442 ! &0.1244272E+04_JPRB,&
443 ! &0.0000000E-00_JPRB,&
444 ! &0.0000000E-00_JPRB,&
445 ! &0.0000000E-00_JPRB,&
446 ! &0.0000000E-00_JPRB/)
447 !ZAPAD4(1, 3, :) = (/&
448 ! &0.4766857E-04_JPRB,&
449 ! &0.4019806E-02_JPRB,&
450 ! &0.1936247E+00_JPRB,&
451 ! &0.1886522E+02_JPRB,&
452 ! &0.6498215E+02_JPRB,&
453 ! &0.1799387E+02_JPRB,&
454 ! &0.0000000E+00_JPRB/)
455
456 !ZBPAD4(1, 1, :) = (/&
457 ! &0.1095121E+03_JPRB,&
458 ! &0.8919103E+03_JPRB,&
459 ! &0.1142218E+03_JPRB,&
460 ! &0.1000000E+01_JPRB,&
461 ! &0.0000000E-00_JPRB,&
462 ! &0.0000000E-00_JPRB,&
463 ! &0.0000000E-00_JPRB/)
464 !ZBPAD4(1, 2, :) = (/&
465 ! &0.5124541E+03_JPRB,&
466 ! &0.3499516E+04_JPRB,&
467 ! &0.1278186E+04_JPRB,&
468 ! &0.1000000E+01_JPRB,&
469 ! &0.0000000E-00_JPRB,&
470 ! &0.0000000E-00_JPRB,&
471 ! &0.0000000E-00_JPRB/)
472 !ZBPAD4(1, 3, :) = (/&
473 ! &0.4766857E-04_JPRB,&
474 ! &0.4153343E-02_JPRB,&
475 ! &0.1925810E+00_JPRB,&
476 ! &0.1938938E+02_JPRB,&
477 ! &0.6918191E+02_JPRB,&
478 ! &0.2278795E+02_JPRB,&
479 ! &0.1000000E+01_JPRB/)
480
481 !ZRAY4(1,:)= (/&
482 ! &.389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,&
483 ! &.659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/)
484
485 !-- Original coefficients in ECMWF scheme
486
487 ZTDH4IS = 0.450_JPRB
488 ZTDU4IS = 0.375_JPRB
489 ZTH4IS = 296._JPRB
490 ZTU4IS = 296._JPRB
491 ZPDH4IS = 0.90_JPRB
492 ZPDU4IS = 0.75_JPRB
493 ZPRH4IS = 100000._JPRB
494 ZPRU4IS = 100000._JPRB
495
496
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD4(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
497
498 ZAPAD4(1, 1, :) = (/&
499 & 0.184678379E+06_JPRB,&
500 & 0.553080884E+05_JPRB,&
501 & 0.248143712E+04_JPRB,&
502 & 0.000000000E-00_JPRB,&
503 & 0.000000000E-00_JPRB,&
504 & 0.000000000E-00_JPRB,&
505
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
506 ZAPAD4(1, 2, :) = (/&
507 & 0.715303869E+01_JPRB,&
508 & 0.219386847E+03_JPRB,&
509 & 0.830001089E+03_JPRB,&
510 & 0.000000000E-00_JPRB,&
511 & 0.000000000E-00_JPRB,&
512 & 0.000000000E-00_JPRB,&
513
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
514 ZAPAD4(1, 3, :) = (/&
515 & 0.925887084E-04_JPRB,&
516 & 0.129353723E-01_JPRB,&
517 & 0.800821928E+00_JPRB,&
518 & 0.242715973E+02_JPRB,&
519 & 0.878331486E+02_JPRB,&
520 & 0.191559725E+02_JPRB,&
521
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
522
523 ZBPAD4(1, 1, :) = (/&
524 & 0.184678379E+06_JPRB,&
525 & 0.555188347E+05_JPRB,&
526 & 0.253257443E+04_JPRB,&
527 & 0.100000000E+01_JPRB,&
528 & 0.000000000E-00_JPRB,&
529 & 0.000000000E-00_JPRB,&
530
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
531 ZBPAD4(1, 2, :) = (/&
532 & 0.715303869E+01_JPRB,&
533 & 0.219441875E+03_JPRB,&
534 & 0.831119997E+03_JPRB,&
535 & 0.100000000E+01_JPRB,&
536 & 0.000000000E-00_JPRB,&
537 & 0.000000000E-00_JPRB,&
538
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
539 ZBPAD4(1, 3, :) = (/&
540 & 0.925887084E-04_JPRB,&
541 & 0.131812683E-01_JPRB,&
542 & 0.812706117E+00_JPRB,&
543 & 0.249863591E+02_JPRB,&
544 & 0.931071925E+02_JPRB,&
545 & 0.252233437E+02_JPRB,&
546
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.100000000E+01_JPRB/)
547
548 ZRAY4(1,:)= (/&
549 & .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
550
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
551
552 2 ZSWCE4(1) = 0.0_JPRB
553 2 ZSWCP4(1) = 0.0_JPRB
554
555 ! ----------------------------------------------------------------
556
557 !* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
558
559 ! ----------------------------------------------------------------
560
561 !* 0.69 - 1.19 Micron
562
563 2 ZSUN4(2) = 0.326158_JPRB
564
565 !-- New coefficients
566
567 !ZD4(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /)
568
569 !ZAPAD4(2, 1, :) = (/&
570 ! &0.1335726E+02_JPRB,&
571 ! &0.2939136E+04_JPRB,&
572 ! &0.4010585E+05_JPRB,&
573 ! &0.7195030E+05_JPRB,&
574 ! &0.1648338E+05_JPRB,&
575 ! &0.3373738E+03_JPRB,&
576 ! &0.0000000E+00_JPRB/)
577 !ZAPAD4(2, 2, :) = (/&
578 ! &0.2001271E-01_JPRB,&
579 ! &0.2480831E+01_JPRB,&
580 ! &0.3444162E+02_JPRB,&
581 ! &0.4788946E+02_JPRB,&
582 ! &0.0000000E+00_JPRB,&
583 ! &0.0000000E+00_JPRB,&
584 ! &0.0000000E+00_JPRB/)
585 !ZAPAD4(2, 3, :) = (/&
586 ! &0.000000000E+00_JPRB,&
587 ! &0.000000000E+00_JPRB,&
588 ! &0.000000000E+00_JPRB,&
589 ! &0.000000000E+00_JPRB,&
590 ! &0.000000000E+00_JPRB,&
591 ! &0.000000000E+00_JPRB,&
592 ! &0.000000000E+00_JPRB/)
593
594 !ZBPAD4(2, 1, :) = (/&
595 ! &0.1335726E+02_JPRB,&
596 ! &0.2942327E+04_JPRB,&
597 ! &0.4077237E+05_JPRB,&
598 ! &0.7749017E+05_JPRB,&
599 ! &0.2123132E+05_JPRB,&
600 ! &0.6659687E+03_JPRB,&
601 ! &0.1000000E+01_JPRB/)
602 !ZBPAD4(2, 2, :) = (/&
603 ! &0.2001271E-01_JPRB,&
604 ! &0.2549067E+01_JPRB,&
605 ! &0.3752433E+02_JPRB,&
606 ! &0.6276637E+02_JPRB,&
607 ! &0.1000000E+01_JPRB,&
608 ! &0.0000000E+00_JPRB,&
609 ! &0.0000000E+00_JPRB/)
610 !ZBPAD4(2, 3, :) = (/&
611 ! &1.000000000E+00_JPRB,&
612 ! &0.000000000E+00_JPRB,&
613 ! &0.000000000E+00_JPRB,&
614 ! &0.000000000E+00_JPRB,&
615 ! &0.000000000E+00_JPRB,&
616 ! &0.000000000E+00_JPRB,&
617 ! &0.000000000E+00_JPRB/)
618
619 !DATA (ZRAY4(2,K),K=1,6) /&
620 ! & .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,&
621 ! &-.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/
622
623 !-- Original coefficients in the ECMWF model
624
625
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD4(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /)
626
627 ZAPAD4(2, 1, :) = (/&
628 & 0.690730834E-02_JPRB,&
629 & 0.151704275E+01_JPRB,&
630 & 0.751477543E+02_JPRB,&
631 & 0.759770236E+03_JPRB,&
632 & 0.109800326E+04_JPRB,&
633 & 0.148407574E+03_JPRB,&
634
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
635 ZAPAD4(2, 2, :) = (/&
636 & 0.863790752E-03_JPRB,&
637 & 0.448762291E+00_JPRB,&
638 & 0.332530367E+02_JPRB,&
639 & 0.190914146E+03_JPRB,&
640 & 0.000000000E+00_JPRB,&
641 & 0.000000000E+00_JPRB,&
642
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
643 ZAPAD4(2, 3, :) = (/&
644 & 0.000000000E+00_JPRB,&
645 & 0.000000000E+00_JPRB,&
646 & 0.000000000E+00_JPRB,&
647 & 0.000000000E+00_JPRB,&
648 & 0.000000000E+00_JPRB,&
649 & 0.000000000E+00_JPRB,&
650
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
651
652 ZBPAD4(2, 1, :) = (/&
653 & 0.690730834E-02_JPRB,&
654 & 0.151954406E+01_JPRB,&
655 & 0.756512527E+02_JPRB,&
656 & 0.779384997E+03_JPRB,&
657 & 0.121113108E+04_JPRB,&
658 & 0.207678436E+03_JPRB,&
659
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.100000000E+01_JPRB/)
660 ZBPAD4(2, 2, :) = (/&
661 & 0.863790752E-03_JPRB,&
662 & 0.448948107E+00_JPRB,&
663 & 0.333186750E+02_JPRB,&
664 & 0.192727216E+03_JPRB,&
665 & 0.100000000E+01_JPRB,&
666 & 0.000000000E+00_JPRB,&
667
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
668 ZBPAD4(2, 3, :) = (/&
669 & 1.000000000E+00_JPRB,&
670 & 0.000000000E+00_JPRB,&
671 & 0.000000000E+00_JPRB,&
672 & 0.000000000E+00_JPRB,&
673 & 0.000000000E+00_JPRB,&
674 & 0.000000000E+00_JPRB,&
675
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.000000000E+00_JPRB/)
676
677 DATA (ZRAY4(2,I_K),I_K=1,6) /&
678 & .164261E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
679 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/
680
681 2 ZSWCE4(2) = 0.0_JPRB
682 2 ZSWCP4(2) = 0.0_JPRB
683
684 ! ----------------------------------------------------------------
685
686 !* 1.19 - 2.38 Microns
687
688 2 ZSUN4(3) = 0.180608_JPRB
689
690 !-- New coefficients
691
692 !ZD4(3,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
693
694 !ZAPAD4(3, 1, :) = (/&
695 ! &0.3325841E-03_JPRB,&
696 ! &0.6194496E+00_JPRB,&
697 ! &0.1497138E+03_JPRB,&
698 ! &0.2314864E+04_JPRB,&
699 ! &0.2380109E+04_JPRB,&
700 ! &0.9553823E+02_JPRB,&
701 ! &0.0000000E+00_JPRB/)
702 !ZAPAD4(3, 2, :) = (/&
703 ! &0.4552471E-03_JPRB,&
704 ! &0.4084154E+00_JPRB,&
705 ! &0.6114905E+01_JPRB,&
706 ! &0.7102540E+01_JPRB,&
707 ! &0.0000000E+00_JPRB,&
708 ! &0.0000000E+00_JPRB,&
709 ! &0.0000000E+00_JPRB/)
710 !ZAPAD4(3, 3, :) = (/&
711 ! &0.000000000E+00_JPRB,&
712 ! &0.000000000E+00_JPRB,&
713 ! &0.000000000E+00_JPRB,&
714 ! &0.000000000E+00_JPRB,&
715 ! &0.000000000E+00_JPRB,&
716 ! &0.000000000E+00_JPRB,&
717 ! &0.000000000E+00_JPRB/)
718
719 !ZBPAD4(3, 1, :) = (/&
720 ! &0.3325841E-03_JPRB,&
721 ! &0.6231947E+00_JPRB,&
722 ! &0.1553098E+03_JPRB,&
723 ! &0.2822458E+04_JPRB,&
724 ! &0.3885194E+04_JPRB,&
725 ! &0.2700235E+03_JPRB,&
726 ! &0.1000000E+01_JPRB/)
727 !ZBPAD4(3, 2, :) = (/&
728 ! &0.4552471E-03_JPRB,&
729 ! &0.4088242E+00_JPRB,&
730 ! &0.6411905E+01_JPRB,&
731 ! &0.9444439E+01_JPRB,&
732 ! &0.1000000E+01_JPRB,&
733 ! &0.0000000E+00_JPRB,&
734 ! &0.0000000E+00_JPRB/)
735 !ZBPAD4(3, 3, :) = (/&
736 ! &1.000000000E+00_JPRB,&
737 ! &0.000000000E+00_JPRB,&
738 ! &0.000000000E+00_JPRB,&
739 ! &0.000000000E+00_JPRB,&
740 ! &0.000000000E+00_JPRB,&
741 ! &0.000000000E+00_JPRB,&
742 ! &0.000000000E+00_JPRB/)
743
744 !ZRAY4(3,:)= (/&
745 ! &.288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,&
746 ! &.540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/)
747
748
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD4(3,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /)
749
750 ZAPAD4(3, 1, :) = (/&
751 & 0.837531303E-05_JPRB,&
752 & 0.173886341E-01_JPRB,&
753 & 0.518852799E+01_JPRB,&
754 & 0.159078416E+03_JPRB,&
755 & 0.493273523E+03_JPRB,&
756 & 0.102567293E+03_JPRB,&
757
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
758 ZAPAD4(3, 2, :) = (/&
759 & 0.657978575E-02_JPRB,&
760 & 0.752617872E+00_JPRB,&
761 & 0.158209734E+02_JPRB,&
762 & 0.410274915E+02_JPRB,&
763 & 0.000000000E+00_JPRB,&
764 & 0.000000000E+00_JPRB,&
765
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
766 ZAPAD4(3, 3, :) = (/&
767 & 0.000000000E+00_JPRB,&
768 & 0.000000000E+00_JPRB,&
769 & 0.000000000E+00_JPRB,&
770 & 0.000000000E+00_JPRB,&
771 & 0.000000000E+00_JPRB,&
772 & 0.000000000E+00_JPRB,&
773
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
774
775 ZBPAD4(3, 1, :) = (/&
776 & 0.837531303E-05_JPRB,&
777 & 0.174882536E-01_JPRB,&
778 & 0.534536580E+01_JPRB,&
779 & 0.180351767E+03_JPRB,&
780 & 0.673126838E+03_JPRB,&
781 & 0.182718543E+03_JPRB,&
782
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.100000000E+01_JPRB/)
783 ZBPAD4(3, 2, :) = (/&
784 & 0.657978575E-02_JPRB,&
785 & 0.753752065E+00_JPRB,&
786 & 0.159286262E+02_JPRB,&
787 & 0.424278450E+02_JPRB,&
788 & 0.100000000E+01_JPRB,&
789 & 0.000000000E+00_JPRB,&
790
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
791 ZBPAD4(3, 3, :) = (/&
792 & 1.000000000E+00_JPRB,&
793 & 0.000000000E+00_JPRB,&
794 & 0.000000000E+00_JPRB,&
795 & 0.000000000E+00_JPRB,&
796 & 0.000000000E+00_JPRB,&
797 & 0.000000000E+00_JPRB,&
798
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
799
800 ZRAY4(3,:)= (/&
801 & .180438E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
802
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
803
804 2 ZSWCE4(3) = 0.0_JPRB
805 2 ZSWCP4(3) = 0.0_JPRB
806
807 ! ----------------------------------------------------------------
808
809 !* 2.38 - 4.00 Microns
810
811 2 ZSUN4(4) = 0.033474_JPRB
812
813 !-- New coefficients
814
815 !ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
816
817 !ZAPAD4(4, 1, :) = (/&
818 ! &0.2122889E-06_JPRB,&
819 ! &0.9030576E-03_JPRB,&
820 ! &0.2431282E+00_JPRB,&
821 ! &0.4901345E+01_JPRB,&
822 ! &0.3996347E+01_JPRB,&
823 ! &0.3910227E+01_JPRB,&
824 ! &0.0000000E+00_JPRB/)
825 !ZAPAD4(4, 2, :) = (/&
826 ! &0.1215163E-03_JPRB,&
827 ! &0.1222574E+00_JPRB,&
828 ! &0.9382420E+01_JPRB,&
829 ! &0.6875727E+02_JPRB,&
830 ! &0.2746421E+02_JPRB,&
831 ! &0.0000000E+00_JPRB,&
832 ! &0.0000000E+00_JPRB/)
833 !ZAPAD4(4, 3, :) = (/&
834 ! &0.263068898E+02_JPRB,&
835 ! &0.146425875E+03_JPRB,&
836 ! &0.860137809E+02_JPRB,&
837 ! &0.000000000E+00_JPRB,&
838 ! &0.000000000E+00_JPRB,&
839 ! &0.000000000E+00_JPRB,&
840 ! &0.000000000E+00_JPRB/)
841
842 !ZBPAD4(4, 1, :) = (/&
843 ! &0.2122889E-06_JPRB,&
844 ! &0.9379083E-03_JPRB,&
845 ! &0.2957335E+00_JPRB,&
846 ! &0.8747190E+01_JPRB,&
847 ! &0.1015794E+02_JPRB,&
848 ! &0.1361277E+02_JPRB,&
849 ! &0.1000000E+01_JPRB/)
850 !ZBPAD4(4, 2, :) = (/&
851 ! &0.1215163E-03_JPRB,&
852 ! &0.1255648E+00_JPRB,&
853 ! &0.1060119E+02_JPRB,&
854 ! &0.8414439E+02_JPRB,&
855 ! &0.4299438E+02_JPRB,&
856 ! &0.1000000E+01_JPRB,&
857 ! &0.0000000E+00_JPRB/)
858 !ZBPAD4(4, 3, :) = (/&
859 ! &0.263068898E+02_JPRB,&
860 ! &0.152569217E+03_JPRB,&
861 ! &0.976791971E+02_JPRB,&
862 ! &0.100000000E+01_JPRB,&
863 ! &0.000000000E+00_JPRB,&
864 ! &0.000000000E+00_JPRB,&
865 ! &0.000000000E+00_JPRB/)
866
867 !ZRAY4(4,:)= (/&
868 ! &.493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,&
869 ! &.741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/)
870
871
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
872
873 ZAPAD4(4, 1, :) = (/&
874 & 0.122118185E-06_JPRB,&
875 & 0.154042531E-02_JPRB,&
876 & 0.141152193E+01_JPRB,&
877 & 0.685368761E+02_JPRB,&
878 & 0.216522281E+03_JPRB,&
879 & 0.421228746E+02_JPRB,&
880
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
881 ZAPAD4(4, 2, :) = (/&
882 & 0.364233560E-10_JPRB,&
883 & 0.217340835E-06_JPRB,&
884 & 0.292623386E-03_JPRB,&
885 & 0.797100631E-01_JPRB,&
886 & 0.319103672E+01_JPRB,&
887 & 0.110530283E+02_JPRB,&
888
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
889 ZAPAD4(4, 3, :) = (/&
890 & 0.263068898E+02_JPRB,&
891 & 0.146425875E+03_JPRB,&
892 & 0.860137809E+02_JPRB,&
893 & 0.000000000E+00_JPRB,&
894 & 0.000000000E+00_JPRB,&
895 & 0.000000000E+00_JPRB,&
896
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
897
898 ZBPAD4(4, 1, :) = (/&
899 & 0.122118185E-06_JPRB,&
900 & 0.156287582E-02_JPRB,&
901 & 0.156994562E+01_JPRB,&
902 & 0.102304103E+03_JPRB,&
903 & 0.475479878E+03_JPRB,&
904 & 0.188714799E+03_JPRB,&
905
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.100000000E+01_JPRB/)
906 ZBPAD4(4, 2, :) = (/&
907 & 0.364233560E-10_JPRB,&
908 & 0.218265806E-06_JPRB,&
909 & 0.297085917E-03_JPRB,&
910 & 0.834253403E-01_JPRB,&
911 & 0.358290002E+01_JPRB,&
912 & 0.139206164E+02_JPRB,&
913
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.100000000E+01_JPRB/)
914 ZBPAD4(4, 3, :) = (/&
915 & 0.263068898E+02_JPRB,&
916 & 0.152569217E+03_JPRB,&
917 & 0.976791971E+02_JPRB,&
918 & 0.100000000E+01_JPRB,&
919 & 0.000000000E+00_JPRB,&
920 & 0.000000000E+00_JPRB,&
921
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
922
923 ZRAY4(4,:)= (/&
924 & .136020E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
925
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
926
927 2 ZSWCE4(4) = 0.0_JPRB
928 2 ZSWCP4(4) = 0.0_JPRB
929
930 !=====================================================================
931
932 !* 1.3 COEFFICIENTS FOR SIX SPECTRAL INTERVALS
933 ! ---------------------------------------
934
935 !* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
936 ! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
937 ! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
938 ! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
939
940 ZADJUS6 = 0.9905709_JPRB
941 ZADJUS6 = 1.0_JPRB
942
943 ZTDH6IS = 0.450_JPRB
944 ZTDU6IS = 0.375_JPRB
945 ZTH6IS = 273._JPRB
946 ZTU6IS = 273._JPRB
947 ZPDH6IS = 0.90_JPRB
948 ZPDU6IS = 0.75_JPRB
949 ZPRH6IS = 101300._JPRB
950 ZPRU6IS = 101300._JPRB
951
952 !* 1st spectral interval: U.V. (0.18 - 0.25 Micron)
953
954 2 ZSUN6(1) = 0.001917_JPRB
955
956
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD6(1,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
957
958 ZAPAD6(1, 1, :) = (/&
959 & 0.000000000E-00_JPRB,&
960 & 0.000000000E-00_JPRB,&
961 & 0.000000000E-00_JPRB,&
962 & 0.000000000E-00_JPRB,&
963 & 0.000000000E-00_JPRB,&
964 & 0.000000000E-00_JPRB,&
965
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
966 ZAPAD6(1, 2, :) = (/&
967 & 0.000000000E-00_JPRB,&
968 & 0.000000000E-00_JPRB,&
969 & 0.000000000E-00_JPRB,&
970 & 0.000000000E-00_JPRB,&
971 & 0.000000000E-00_JPRB,&
972 & 0.000000000E-00_JPRB,&
973
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
974 ZAPAD6(1, 3, :) = (/&
975 & 0.000000000E-00_JPRB,&
976 & 0.000000000E-00_JPRB,&
977 & 0.000000000E-00_JPRB,&
978 & 0.000000000E-00_JPRB,&
979 & 0.000000000E-00_JPRB,&
980 & 0.000000000E-00_JPRB,&
981
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
982
983 ZBPAD6(1, 1, :) = (/&
984 & 0.100000000E+01_JPRB,&
985 & 0.000000000E-00_JPRB,&
986 & 0.000000000E-00_JPRB,&
987 & 0.000000000E-00_JPRB,&
988 & 0.000000000E-00_JPRB,&
989 & 0.000000000E-00_JPRB,&
990
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
991 ZBPAD6(1, 2, :) = (/&
992 & 0.100000000E+01_JPRB,&
993 & 0.000000000E-00_JPRB,&
994 & 0.000000000E-00_JPRB,&
995 & 0.000000000E-00_JPRB,&
996 & 0.000000000E-00_JPRB,&
997 & 0.000000000E-00_JPRB,&
998
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
999 ZBPAD6(1, 3, :) = (/&
1000 & 0.100000000E+01_JPRB,&
1001 & 0.000000000E-00_JPRB,&
1002 & 0.000000000E-00_JPRB,&
1003 & 0.000000000E-00_JPRB,&
1004 & 0.000000000E-00_JPRB,&
1005 & 0.000000000E-00_JPRB,&
1006
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1007
1008 ZRAY6(1,:)= (/&
1009 & .334018E+00_JPRB,-.298776E+01_JPRB, .613911E+02_JPRB,&
1010
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & -.164331E+03_JPRB, .197946E+03_JPRB,-.652206E+02_JPRB/)
1011 ZRBY6(1,:)= (/&
1012 & .43959002E+01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
1013
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
1014
1015 2 ZSWCE6(1) = 0.0_JPRB
1016 2 ZSWCP6(1) = 0.0_JPRB
1017
1018 2 NEXPO3(1) = 7
1019 REXPO3(1, 1, :) =(/&
1020 & 0.051395E+00_JPRB, 0.048250E+00_JPRB, 0.112339E+00_JPRB,&
1021 & 0.101426E+00_JPRB, 0.007700E+00_JPRB, 0.441320E+00_JPRB,&
1022
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.237571E+00_JPRB /)
1023 REXPO3(1, 2, :) =(/&
1024 & 0.100022E+02_JPRB, 0.851159E+02_JPRB, 0.346737E+03_JPRB,&
1025 & 0.158501E+02_JPRB, 0.724223E+01_JPRB, 0.177828E+03_JPRB,&
1026
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.467708E+02_JPRB /)
1027
1028 !* 2nd spectral interval: U.V. (0.25 - 0.44 Micron)
1029
1030 2 ZSUN6(2) = 0.135708_JPRB
1031
1032
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD6(2,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
1033
1034 ZAPAD6(2, 1, :) = (/&
1035 & 0.000000000E-00_JPRB,&
1036 & 0.000000000E-00_JPRB,&
1037 & 0.000000000E-00_JPRB,&
1038 & 0.000000000E-00_JPRB,&
1039 & 0.000000000E-00_JPRB,&
1040 & 0.000000000E-00_JPRB,&
1041
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1042 ZAPAD6(2, 2, :) = (/&
1043 & 0.000000000E-00_JPRB,&
1044 & 0.000000000E-00_JPRB,&
1045 & 0.000000000E-00_JPRB,&
1046 & 0.000000000E-00_JPRB,&
1047 & 0.000000000E-00_JPRB,&
1048 & 0.000000000E-00_JPRB,&
1049
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1050 ZAPAD6(2, 3, :) = (/&
1051 & 0.000000000E-00_JPRB,&
1052 & 0.000000000E-00_JPRB,&
1053 & 0.000000000E-00_JPRB,&
1054 & 0.000000000E-00_JPRB,&
1055 & 0.000000000E-00_JPRB,&
1056 & 0.000000000E-00_JPRB,&
1057
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1058
1059 ZBPAD6(2, 1, :) = (/&
1060 & 0.100000000E+01_JPRB,&
1061 & 0.000000000E-00_JPRB,&
1062 & 0.000000000E-00_JPRB,&
1063 & 0.000000000E-00_JPRB,&
1064 & 0.000000000E-00_JPRB,&
1065 & 0.000000000E-00_JPRB,&
1066
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1067 ZBPAD6(2, 2, :) = (/&
1068 & 0.100000000E+01_JPRB,&
1069 & 0.000000000E-00_JPRB,&
1070 & 0.000000000E-00_JPRB,&
1071 & 0.000000000E-00_JPRB,&
1072 & 0.000000000E-00_JPRB,&
1073 & 0.000000000E-00_JPRB,&
1074
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1075 ZBPAD6(2, 3, :) = (/&
1076 & 0.100000000E+01_JPRB,&
1077 & 0.000000000E-00_JPRB,&
1078 & 0.000000000E-00_JPRB,&
1079 & 0.000000000E-00_JPRB,&
1080 & 0.000000000E-00_JPRB,&
1081 & 0.000000000E-00_JPRB,&
1082
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1083
1084 ZRAY6(2,:)= (/&
1085 & -.895408E-03_JPRB, .283130E+01_JPRB,-.867917E+01_JPRB,&
1086
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & .148672E+02_JPRB,-.127282E+02_JPRB, .422212E+01_JPRB/)
1087 ZRBY6(2,:)= (/&
1088 & .55503070E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
1089
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
1090
1091 2 ZSWCE6(2) = 0.0_JPRB
1092 2 ZSWCP6(2) = 0.0_JPRB
1093
1094 2 NEXPO3(2) = 7
1095 REXPO3(2, 1, :) =(/&
1096 & 0.043801E+00_JPRB, 0.078893E+00_JPRB, 0.036839E+00_JPRB,&
1097 & 0.022503E+00_JPRB, 0.042333E+00_JPRB, 0.037870E+00_JPRB,&
1098
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.737762E+00_JPRB /)
1099 REXPO3(2, 2, :) =(/&
1100 & 0.234249E+01_JPRB, 0.125170E+00_JPRB, 0.549527E+02_JPRB,&
1101 & 0.257041E+03_JPRB, 0.476838E+00_JPRB, 0.911993E+01_JPRB,&
1102
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.000000E+00_JPRB /)
1103
1104 !* 3rd spectral interval: Visible (0.44 - 0.69 Micron)
1105
1106 2 ZSUN6(3) = 0.322135_JPRB
1107
1108
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD6(3,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
1109
1110 ZAPAD6(3, 1, :) = (/&
1111 & 0.1762097E+03_JPRB,&
1112 & 0.1641762E+03_JPRB,&
1113 & 0.8687919E+02_JPRB,&
1114 & 0.0000000E-00_JPRB,&
1115 & 0.0000000E-00_JPRB,&
1116 & 0.0000000E-00_JPRB,&
1117
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E-00_JPRB/)
1118 ZAPAD6(3, 2, :) = (/&
1119 & 0.5581224E+00_JPRB,&
1120 & 0.1748430E+03_JPRB,&
1121 & 0.1134123E+04_JPRB,&
1122 & 0.3490429E+03_JPRB,&
1123 & 0.0000000E-00_JPRB,&
1124 & 0.0000000E-00_JPRB,&
1125
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E-00_JPRB/)
1126 ZAPAD6(3, 3, :) = (/&
1127 & 0.000000000E-00_JPRB,&
1128 & 0.000000000E-00_JPRB,&
1129 & 0.000000000E-00_JPRB,&
1130 & 0.000000000E-00_JPRB,&
1131 & 0.000000000E-00_JPRB,&
1132 & 0.000000000E-00_JPRB,&
1133
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1134
1135 ZBPAD6(3, 1, :) = (/&
1136 & 0.1762097E+03_JPRB,&
1137 & 0.1663950E+03_JPRB,&
1138 & 0.8939724E+02_JPRB,&
1139 & 0.1000000E+01_JPRB,&
1140 & 0.0000000E-00_JPRB,&
1141 & 0.0000000E-00_JPRB,&
1142
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E-00_JPRB/)
1143 ZBPAD6(3, 2, :) = (/&
1144 & 0.5581224E+00_JPRB,&
1145 & 0.1749251E+03_JPRB,&
1146 & 0.1159910E+04_JPRB,&
1147 & 0.3893268E+03_JPRB,&
1148 & 0.1000000E+01_JPRB,&
1149 & 0.0000000E-00_JPRB,&
1150
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E-00_JPRB/)
1151 ZBPAD6(3, 3, :) = (/&
1152 & 0.100000000E+01_JPRB,&
1153 & 0.000000000E-00_JPRB,&
1154 & 0.000000000E-00_JPRB,&
1155 & 0.000000000E-00_JPRB,&
1156 & 0.000000000E-00_JPRB,&
1157 & 0.000000000E-00_JPRB,&
1158
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E-00_JPRB/)
1159
1160 ZRAY6(3,:)= (/&
1161 & .474852E-01_JPRB, .605110E+00_JPRB,-.248070E+01_JPRB,&
1162
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & .485580E+01_JPRB,-.444043E+01_JPRB, .152369E+01_JPRB/)
1163 ZRBY6(3,:)= (/&
1164 & .10528199E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
1165
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
1166
1167 2 ZSWCE6(3) = 0.0_JPRB
1168 2 ZSWCP6(3) = 0.0_JPRB
1169
1170 2 NEXPO3(3) = 6
1171 REXPO3(3, 1, :) =(/&
1172 & 0.063442E+00_JPRB, 0.058550E+00_JPRB, 0.237534E+00_JPRB,&
1173 & 0.412292E+00_JPRB, 0.126141E+00_JPRB, 0.102041E+00_JPRB,&
1174
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000E+00_JPRB /)
1175 REXPO3(3, 2, :) =(/&
1176 & 0.125170E+00_JPRB, 0.119209E-01_JPRB, 0.119209E+00_JPRB,&
1177 & 0.417233E-01_JPRB, 0.894070E-01_JPRB, 0.000000E+00_JPRB,&
1178
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.000000E+00_JPRB /)
1179
1180 ! ----------------------------------------------------------------
1181
1182 !* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
1183
1184 ! ----------------------------------------------------------------
1185
1186 !* 0.69 - 1.19 Micron
1187
1188 !* UMG is O2 only
1189
1190 2 ZSUN6(4) = 0.326158_JPRB
1191
1192
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD6(4,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /)
1193
1194 ZAPAD6(4, 1, :) = (/&
1195 & 0.1335726E+02_JPRB,&
1196 & 0.2939136E+04_JPRB,&
1197 & 0.4010585E+05_JPRB,&
1198 & 0.7195030E+05_JPRB,&
1199 & 0.1648338E+05_JPRB,&
1200 & 0.3373738E+03_JPRB,&
1201
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1202 ZAPAD6(4, 2, :) = (/&
1203 & 0.2001271E-01_JPRB,&
1204 & 0.2480831E+01_JPRB,&
1205 & 0.3444162E+02_JPRB,&
1206 & 0.4788946E+02_JPRB,&
1207 & 0.0000000E+00_JPRB,&
1208 & 0.0000000E+00_JPRB,&
1209
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1210 ZAPAD6(4, 3, :) = (/&
1211 & 0.000000000E+00_JPRB,&
1212 & 0.000000000E+00_JPRB,&
1213 & 0.000000000E+00_JPRB,&
1214 & 0.000000000E+00_JPRB,&
1215 & 0.000000000E+00_JPRB,&
1216 & 0.000000000E+00_JPRB,&
1217
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
1218
1219 ZBPAD6(4, 1, :) = (/&
1220 & 0.1335726E+02_JPRB,&
1221 & 0.2942327E+04_JPRB,&
1222 & 0.4077237E+05_JPRB,&
1223 & 0.7749017E+05_JPRB,&
1224 & 0.2123132E+05_JPRB,&
1225 & 0.6659687E+03_JPRB,&
1226
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.1000000E+01_JPRB/)
1227 ZBPAD6(4, 2, :) = (/&
1228 & 0.2001271E-01_JPRB,&
1229 & 0.2549067E+01_JPRB,&
1230 & 0.3752433E+02_JPRB,&
1231 & 0.6276637E+02_JPRB,&
1232 & 0.1000000E+01_JPRB,&
1233 & 0.0000000E+00_JPRB,&
1234
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1235 ZBPAD6(4, 3, :) = (/&
1236 & 1.000000000E+00_JPRB,&
1237 & 0.000000000E+00_JPRB,&
1238 & 0.000000000E+00_JPRB,&
1239 & 0.000000000E+00_JPRB,&
1240 & 0.000000000E+00_JPRB,&
1241 & 0.000000000E+00_JPRB,&
1242
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
1243
1244 ZRAY6(4,:)= (/&
1245 & .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,&
1246
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & -.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/)
1247 ZRBY6(4,:)= (/&
1248 & .16436996E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
1249
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
1250
1251 2 ZSWCE6(4) = 0.0_JPRB
1252 2 ZSWCP6(4) = 0.0_JPRB
1253
1254 2 NEXPO3(4) = 4
1255 REXPO3(4, 1, :) =(/&
1256 & 0.000074E+00_JPRB, 0.320194E+00_JPRB, 0.082915E+00_JPRB,&
1257 & 0.596816E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1258
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000E+00_JPRB /)
1259 REXPO3(4, 2, :) =(/&
1260 & 0.232458E+00_JPRB, 0.119209E-01_JPRB, 0.178814E-01_JPRB,&
1261 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1262
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.000000E+00_JPRB /)
1263
1264 ! ----------------------------------------------------------------
1265
1266 !* 1.19 - 2.38 Microns
1267
1268 !* UMG is CO2 only
1269
1270 2 ZSUN6(5) = 0.180608_JPRB
1271
1272
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD6(5,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
1273
1274 ZAPAD6(5, 1, :) = (/&
1275 & 0.3325841E-03_JPRB,&
1276 & 0.6194496E+00_JPRB,&
1277 & 0.1497138E+03_JPRB,&
1278 & 0.2314864E+04_JPRB,&
1279 & 0.2380109E+04_JPRB,&
1280 & 0.9553823E+02_JPRB,&
1281
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1282 ZAPAD6(5, 2, :) = (/&
1283 & 0.4552471E-03_JPRB,&
1284 & 0.4084154E+00_JPRB,&
1285 & 0.6114905E+01_JPRB,&
1286 & 0.7102540E+01_JPRB,&
1287 & 0.0000000E+00_JPRB,&
1288 & 0.0000000E+00_JPRB,&
1289
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1290 ZAPAD6(5, 3, :) = (/&
1291 & 0.000000000E+00_JPRB,&
1292 & 0.000000000E+00_JPRB,&
1293 & 0.000000000E+00_JPRB,&
1294 & 0.000000000E+00_JPRB,&
1295 & 0.000000000E+00_JPRB,&
1296 & 0.000000000E+00_JPRB,&
1297
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
1298
1299 ZBPAD6(5, 1, :) = (/&
1300 & 0.3325841E-03_JPRB,&
1301 & 0.6231947E+00_JPRB,&
1302 & 0.1553098E+03_JPRB,&
1303 & 0.2822458E+04_JPRB,&
1304 & 0.3885194E+04_JPRB,&
1305 & 0.2700235E+03_JPRB,&
1306
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.1000000E+01_JPRB/)
1307 ZBPAD6(5, 2, :) = (/&
1308 & 0.4552471E-03_JPRB,&
1309 & 0.4088242E+00_JPRB,&
1310 & 0.6411905E+01_JPRB,&
1311 & 0.9444439E+01_JPRB,&
1312 & 0.1000000E+01_JPRB,&
1313 & 0.0000000E+00_JPRB,&
1314
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1315 ZBPAD6(5, 3, :) = (/&
1316 & 1.000000000E+00_JPRB,&
1317 & 0.000000000E+00_JPRB,&
1318 & 0.000000000E+00_JPRB,&
1319 & 0.000000000E+00_JPRB,&
1320 & 0.000000000E+00_JPRB,&
1321 & 0.000000000E+00_JPRB,&
1322
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
1323
1324 ZRAY6(5,:)= (/&
1325 & .288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,&
1326
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & .540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/)
1327 ZRBY6(5,:)= (/&
1328 & .18073079E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
1329
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
1330
1331 2 ZSWCE6(5) = 0.0_JPRB
1332 2 ZSWCP6(5) = 0.0_JPRB
1333
1334 2 NEXPO3(5) = 0
1335 REXPO3(5, 1, :) =(/&
1336 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1337 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1338
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000E+00_JPRB /)
1339 REXPO3(5, 2, :) =(/&
1340 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1341 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1342
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.000000E+00_JPRB /)
1343
1344 ! ----------------------------------------------------------------
1345
1346 !* 2.38 - 4.00 Microns
1347
1348 2 ZSUN6(6) = 0.033474_JPRB
1349
1350
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
 8 ZD6(6,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
1351
1352 ZAPAD6(6, 1, :) = (/&
1353 & 0.2122889E-06_JPRB,&
1354 & 0.9030576E-03_JPRB,&
1355 & 0.2431282E+00_JPRB,&
1356 & 0.4901345E+01_JPRB,&
1357 & 0.3996347E+01_JPRB,&
1358 & 0.3910227E+01_JPRB,&
1359
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1360 ZAPAD6(6, 2, :) = (/&
1361 & 0.1215163E-03_JPRB,&
1362 & 0.1222574E+00_JPRB,&
1363 & 0.9382420E+01_JPRB,&
1364 & 0.6875727E+02_JPRB,&
1365 & 0.2746421E+02_JPRB,&
1366 & 0.0000000E+00_JPRB,&
1367
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1368 ZAPAD6(6, 3, :) = (/&
1369 & 0.263068898E+02_JPRB,&
1370 & 0.146425875E+03_JPRB,&
1371 & 0.860137809E+02_JPRB,&
1372 & 0.000000000E+00_JPRB,&
1373 & 0.000000000E+00_JPRB,&
1374 & 0.000000000E+00_JPRB,&
1375
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
1376
1377 ZBPAD6(6, 1, :) = (/&
1378 & 0.2122889E-06_JPRB,&
1379 & 0.9379083E-03_JPRB,&
1380 & 0.2957335E+00_JPRB,&
1381 & 0.8747190E+01_JPRB,&
1382 & 0.1015794E+02_JPRB,&
1383 & 0.1361277E+02_JPRB,&
1384
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.1000000E+01_JPRB/)
1385 ZBPAD6(6, 2, :) = (/&
1386 & 0.1215163E-03_JPRB,&
1387 & 0.1255648E+00_JPRB,&
1388 & 0.1060119E+02_JPRB,&
1389 & 0.8414439E+02_JPRB,&
1390 & 0.4299438E+02_JPRB,&
1391 & 0.1000000E+01_JPRB,&
1392
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.0000000E+00_JPRB/)
1393 ZBPAD6(6, 3, :) = (/&
1394 & 0.263068898E+02_JPRB,&
1395 & 0.152569217E+03_JPRB,&
1396 & 0.976791971E+02_JPRB,&
1397 & 0.100000000E+01_JPRB,&
1398 & 0.000000000E+00_JPRB,&
1399 & 0.000000000E+00_JPRB,&
1400
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000000E+00_JPRB/)
1401
1402 ZRAY6(6,:)= (/&
1403 & .493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,&
1404
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 & .741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/)
1405 ZRBY6(6,:)= (/&
1406 & .13618247E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
1407
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 & .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
1408
1409 2 ZSWCE6(6) = 0.0_JPRB
1410 2 ZSWCP6(6) = 0.0_JPRB
1411
1412 2 NEXPO3(6) = 0
1413 REXPO3(6, 1, :) =(/&
1414 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1415 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1416
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 14 times.
 16 & 0.000000E+00_JPRB /)
1417 REXPO3(6, 2, :) =(/&
1418 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1419 & 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1420
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
 16 & 0.000000E+00_JPRB /)
1421
1422 !=====================================================================
1423 !* 2.4 OPTICAL PARAMETERS FOR RAIN DROPS
1424 ! Savijarvi et al. (1996)
1425
1426 ZRTO1 = 0.003_JPRB
1427 ZRTO2 = -0.22_JPRB
1428 ! CAUTION JUST TEMPORARY PARAMETERS
1429 2 ZROMA2(1:2)= (/ 0.00008_JPRB , 0.264_JPRB /)
1430 2 ZROMB2(1:2)= (/ 0.23_JPRB , 0.09_JPRB /)
1431 2 ZRASY2(1:2)= (/ 0.88_JPRB , 0.94_JPRB /)
1432
1433 2 ZRA2(1:2)= (/ 1.5_JPRB , 1.5_JPRB /)
1434 2 ZRB2(1:2)= (/ 0.50_JPRB , 1.13_JPRB /)
1435 2 ZRC2(1:2)= (/ 5.58E-7_JPRB , 8.55E-4_JPRB/)
1436 2 ZRD2(1:2)= (/ 1.25E-7_JPRB , 1.28E-3_JPRB/)
1437 2 ZRE2(1:2)= (/ 0.841_JPRB , 0.786_JPRB /)
1438 2 ZRF2(1:2)= (/ 2.08E-3_JPRB , 5.32E-3_JPRB/)
1439
1440 2 ZROMA4(1:4)= (/ 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB , 0.465_JPRB /)
1441 2 ZROMB4(1:4)= (/ 0.23_JPRB , 0.22_JPRB , 0.09_JPRB , 0.001_JPRB /)
1442 2 ZRASY4(1:4)= (/ 0.88_JPRB , 0.89_JPRB , 0.94_JPRB , 0.97_JPRB /)
1443
1444 2 ZRA4(1:4)= (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /)
1445 2 ZRB4(1:4)= (/ 0.50_JPRB , 0.78_JPRB , 1.13_JPRB , 2.00_JPRB /)
1446 2 ZRC4(1:4)= (/ 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /)
1447 2 ZRD4(1:4)= (/ 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /)
1448 2 ZRE4(1:4)= (/ 0.841_JPRB , 0.821_JPRB , 0.786_JPRB , 0.820_JPRB /)
1449 2 ZRF4(1:4)= (/ 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /)
1450
1451 2 ZROMA6(1:6)= (/ 0.00008_JPRB , 0.00008_JPRB , 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB , 0.465_JPRB /)
1452 2 ZROMB6(1:6)= (/ 0.23_JPRB , 0.23_JPRB , 0.23_JPRB , 0.22_JPRB , 0.09_JPRB , 0.001_JPRB /)
1453 2 ZRASY6(1:6)= (/ 0.88_JPRB , 0.88_JPRB , 0.88_JPRB , 0.89_JPRB , 0.94_JPRB , 0.97_JPRB /)
1454
1455 2 ZRA6(1:6)= (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /)
1456 2 ZRB6(1:6)= (/ 0.50_JPRB , 0.50_JPRB , 0.50_JPRB , 0.78_JPRB , 1.13_JPRB , 2.00_JPRB /)
1457 2 ZRC6(1:6)= (/ 5.58E-7_JPRB , 5.58E-7_JPRB , 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /)
1458 2 ZRD6(1:6)= (/ 1.25E-7_JPRB , 1.25E-7_JPRB , 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /)
1459 2 ZRE6(1:6)= (/ 0.841_JPRB , 0.841_JPRB , 0.841_JPRB , 0.821_JPRB , 0.786_JPRB , 0.820_JPRB /)
1460 2 ZRF6(1:6)= (/ 2.08E-3_JPRB , 2.08E-3_JPRB , 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /)
1461
1462 !=====================================================================
1463 ! 3. FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW
1464
1465 ! SW : 14 spectral intervals : mapping to SW6 spectral intervals
1466 ! 0.185-0.25-0.44-0.69-1.19-2.38-4.00
1467 ! 3.846 - 3.077 6
1468 ! 3.077 - 2.500 6
1469 ! 2.500 - 2.150 5
1470 ! 2.150 - 1.942 5
1471 ! 1.942 - 1.626 5
1472 ! 1.626 - 1.299 5
1473 ! 1.299 - 1.242 5
1474 ! 1.242 - 0.7782 4
1475 ! 0.7782- 0.6250 4
1476 ! 0.6250- 0.4415 3
1477 ! 0.4415- 0.3448 2
1478 ! 0.3448- 0.2632 2
1479 ! 0.2632- 0.2000 1
1480 ! 12.195 - 3.846 6
1481
1482 2 IMPSRTM(1:14)= (/ 6,6,5,5,5,5,5,4,4,3,2,2,1,6 /)
1483 !=====================================================================
1484
1485 !* 4. SET VALUES.
1486 ! -----------
1487
1488
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2 IF (LONEWSW) THEN
1489 ZPDH2O = ZPDHNIS
1490 ZPDUMG = ZPDUNIS
1491 ZPRH2O = ZPRHNIS
1492 ZPRUMG = ZPRUNIS
1493 2 RTDH2O = ZTDHNIS
1494 2 RTDUMG = ZTDUNIS
1495 2 RTH2O = ZTHNIS
1496 2 RTUMG = ZTUNIS
1497 ENDIF
1498
1499 2 RPDH1=ZPDH2O+1.0_JPRB
1500 2 RPDU1=ZPDUMG+1.0_JPRB
1501 2 ZH2O=1.0_JPRB/( 10._JPRB* RG * RPDH1 )
1502 2 ZUMG=1.0_JPRB/( 10._JPRB* RG * RPDU1 )
1503 2 RPNU = ZUMG/(ZPRUMG**ZPDUMG)
1504 2 RPNH = ZH2O/(ZPRH2O**ZPDH2O)
1505
1506
4/6
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
 2 IF(.NOT.ALLOCATED(RSUN)) ALLOCATE (RSUN(NTSW))
1507
1508
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 2 times.
 14 DO JW=1,ISW6
1509 12 RSUN(JW)=0.0_JPRB
1510 12 RSWCE(JW)=0.0_JPRB
1511 12 RSWCP(JW)=0.0_JPRB
1512
1513
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 36 times.
 48 DO JC3=1,3
1514 48 D(JW,JC3)=0.0_JPRB
1515 ENDDO
1516
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 72 times.
 84 DO JC6=1,6
1517 84 RRAY(JW,JC6)=0.0_JPRB
1518 ENDDO
1519
2/2
✓ Branch 0 taken 36 times.
✓ Branch 1 taken 12 times.
 50 DO JI=1,3
1520
2/2
✓ Branch 0 taken 252 times.
✓ Branch 1 taken 36 times.
 300 DO JJ=1,7
1521 252 APAD(JW,JI,JJ)=0.0_JPRB
1522 288 BPAD(JW,JI,JJ)=0.0_JPRB
1523 ENDDO
1524 ENDDO
1525 ENDDO
1526
1527
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 IF (KSW == 2) THEN
1528 ZPDH2O = ZPDH2IS
1529 ZPDUMG = ZPDU2IS
1530 ZPRH2O = ZPRH2IS
1531 ZPRUMG = ZPRU2IS
1532 ✗ RTDH2O = ZTDH2IS
1533 ✗ RTDUMG = ZTDU2IS
1534 ✗ RTH2O = ZTH2IS
1535 ✗ RTUMG = ZTU2IS
1536
1537 ✗ RADJUST=ZADJUS2
1538 ✗ RHSRTA=ZRTO1
1539 ✗ RHSRTB=ZRTO2
1540 ✗ DO JW=1,KSW
1541 ✗ RSUN (JW)=ZSUN2(JW)
1542
1543 ✗ RSWCE(JW)=ZSWCE2(JW)
1544 ✗ RSWCP(JW)=ZSWCP2(JW)
1545
1546 ✗ RROMA(JW)=ZROMA2(JW)
1547 ✗ RROMB(JW)=ZROMB2(JW)
1548 ✗ RRASY(JW)=ZRASY2(JW)
1549 ✗ RHSRA(JW)=ZRA2(JW)
1550 ✗ RHSRB(JW)=ZRB2(JW)
1551 ✗ RHSRC(JW)=ZRC2(JW)
1552 ✗ RHSRD(JW)=ZRD2(JW)
1553 ✗ RHSRE(JW)=ZRE2(JW)
1554 ✗ RHSRF(JW)=ZRF2(JW)
1555
1556 ✗ DO JC3=1,3
1557 ✗ D(JW,JC3)=ZD2(JW,JC3)
1558 ENDDO
1559 ✗ DO JC6=1,6
1560 ✗ RRAY(JW,JC6)=ZRAY2(JW,JC6)
1561 ENDDO
1562 ✗ DO JI=1,3
1563 ✗ DO JJ=1,7
1564 ✗ APAD(JW,JI,JJ)=ZAPAD2(JW,JI,JJ)
1565 ✗ BPAD(JW,JI,JJ)=ZBPAD2(JW,JI,JJ)
1566 ENDDO
1567 ENDDO
1568 ENDDO
1569
1570
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 ELSEIF (KSW == 4) THEN
1571 ZPDH2O = ZPDH4IS
1572 ZPDUMG = ZPDU4IS
1573 ZPRH2O = ZPRH4IS
1574 ZPRUMG = ZPRU4IS
1575 ✗ RTDH2O = ZTDH4IS
1576 ✗ RTDUMG = ZTDU4IS
1577 ✗ RTH2O = ZTH4IS
1578 ✗ RTUMG = ZTU4IS
1579
1580 ✗ RADJUST=ZADJUS4
1581 ✗ RHSRTA=ZRTO1
1582 ✗ RHSRTB=ZRTO2
1583 ✗ DO JW=1,KSW
1584 ✗ RSUN (JW)=ZSUN4(JW)
1585
1586 ✗ RSWCE(JW)=ZSWCE4(JW)
1587 ✗ RSWCP(JW)=ZSWCP4(JW)
1588
1589 ✗ RROMA(JW)=ZROMA4(JW)
1590 ✗ RROMB(JW)=ZROMB4(JW)
1591 ✗ RRASY(JW)=ZRASY4(JW)
1592 ✗ RHSRA(JW)=ZRA4(JW)
1593 ✗ RHSRB(JW)=ZRB4(JW)
1594 ✗ RHSRC(JW)=ZRC4(JW)
1595 ✗ RHSRD(JW)=ZRD4(JW)
1596 ✗ RHSRE(JW)=ZRE4(JW)
1597 ✗ RHSRF(JW)=ZRF4(JW)
1598
1599 ✗ DO JC3=1,3
1600 ✗ D(JW,JC3)=ZD4(JW,JC3)
1601 ENDDO
1602 ✗ DO JC6=1,6
1603 ✗ RRAY(JW,JC6)=ZRAY4(JW,JC6)
1604 ENDDO
1605 ✗ DO JI=1,3
1606 ✗ DO JJ=1,7
1607 ✗ APAD(JW,JI,JJ)=ZAPAD4(JW,JI,JJ)
1608 ✗ BPAD(JW,JI,JJ)=ZBPAD4(JW,JI,JJ)
1609 ENDDO
1610 ENDDO
1611 ENDDO
1612
1613
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2 ELSEIF (KSW == 6) THEN
1614 ZPDH2O = ZPDH6IS
1615 ZPDUMG = ZPDU6IS
1616 ZPRH2O = ZPRH6IS
1617 ZPRUMG = ZPRU6IS
1618 2 RTDH2O = ZTDH6IS
1619 2 RTDUMG = ZTDU6IS
1620 2 RTH2O = ZTH6IS
1621 2 RTUMG = ZTU6IS
1622
1623 2 RADJUST=ZADJUS6
1624 2 RHSRTA=ZRTO1
1625 2 RHSRTB=ZRTO2
1626
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 12 times.
 14 DO JW=1,KSW
1627 12 RSUN (JW)=ZSUN6(JW)
1628
1629 12 RSWCE(JW)=ZSWCE6(JW)
1630 12 RSWCP(JW)=ZSWCP6(JW)
1631
1632 12 RROMA(JW)=ZROMA6(JW)
1633 12 RROMB(JW)=ZROMB6(JW)
1634 12 RRASY(JW)=ZRASY6(JW)
1635 12 RHSRA(JW)=ZRA6(JW)
1636 12 RHSRB(JW)=ZRB6(JW)
1637 12 RHSRC(JW)=ZRC6(JW)
1638 12 RHSRD(JW)=ZRD6(JW)
1639 12 RHSRE(JW)=ZRE6(JW)
1640 12 RHSRF(JW)=ZRF6(JW)
1641
1642
2/2
✓ Branch 0 taken 36 times.
✓ Branch 1 taken 12 times.
 48 DO JC3=1,3
1643 48 D(JW,JC3)=ZD6(JW,JC3)
1644 ENDDO
1645
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 12 times.
 12 IF (LRAYL) THEN
1646
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 72 times.
 84 DO JC6=1,6
1647 84 RRAY(JW,JC6)=ZRAY6(JW,JC6)
1648 ENDDO
1649 ELSE
1650 ✗ DO JC6=1,6
1651 ✗ RRAY(JW,JC6)=ZRBY6(JW,JC6)
1652 ENDDO
1653 ENDIF
1654
2/2
✓ Branch 0 taken 36 times.
✓ Branch 1 taken 12 times.
 50 DO JI=1,3
1655
2/2
✓ Branch 0 taken 252 times.
✓ Branch 1 taken 36 times.
 300 DO JJ=1,7
1656 252 APAD(JW,JI,JJ)=ZAPAD6(JW,JI,JJ)
1657 288 BPAD(JW,JI,JJ)=ZBPAD6(JW,JI,JJ)
1658 ENDDO
1659 ENDDO
1660 ENDDO
1661
1662 ✗ ELSEIF (KTSW == 14 .AND. KSW == 14) THEN
1663 ✗ PRINT *,'SUSWN: 14-SPECTRAL INTERVALS --> RRTM_SW'
1664
1665 ELSE
1666 ✗ PRINT *,'SUSWN: NTSW, NSW: ',KTSW,KSW
1667 ✗ CALL ABOR1('SUSWN: WRONG NUMBER OF SW INTERVALS')
1668 ENDIF
1669
1670 ! ----------------------------------------------------------------
1671
1672 !* 5. PRINT VALUES.
1673 ! ------------
1674
1675 ! PRINT 9700,RPDH1,RPDU1,RPNU,RPNH
1676 !9700 FORMAT(1X,'RPDH1,RPDU1,RPNU,RPNH ',4F10.5)
1677 ! DO JW=1,KSW
1678 ! PRINT *,' SW INTERVAL ',JW
1679 ! PRINT 9701,JW,RSUN(JW),RSWCE(JW),RSWCP(JW)
1680 !9701 FORMAT(1X,I3,F10.7,2F10.5)
1681 ! PRINT 9702,JW,(D(JW,JC3),JC3=1,3)
1682 !9702 FORMAT(1X,I3,3F10.5)
1683 ! PRINT 9703,JW,(RRAY(JW,JC6),JC6=1,6)
1684 !9703 FORMAT(1X,I3,6E15.8)
1685 ! PRINT 9704,JW,(JABS,(APAD(JW,JABS,JJ),JJ=1,7),JABS=1,3)
1686 !9704 FORMAT(1X,I3,(I3,7E15.8))
1687 ! PRINT 9706,JW,RROMA(JW),RROMB(JW),RRASY(JW),RHSRTA,RHSRTB
1688 !9706 FORMAT(1X,I3,' Rain Optics from Rain Rate: ',F8.5,F6.3,3F6.3)
1689 ! PRINT 9707,JW,RHSRA(JW),RHSRB(JW),RHSRC(JW),RHSRD(JW)
1690 ! & ,RHSRE(JW),RHSRF(JW)
1691 !9707 FORMAT(1X,I3,' Rain Optics as function of Reff: ',6E11.4)
1692 ! ENDDO
1693
1694 ! ----------------------------------------------------------------
1695
1696
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2 IF (LHOOK) CALL DR_HOOK('SUSWN',1,ZHOOK_HANDLE)
1697 2 END SUBROUTINE SUSWN
1698