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