source: LMDZ5/branches/testing/libf/phymar/suswn.F90 @ 3990

Last change on this file since 3990 was 2160, checked in by Laurent Fairhead, 10 years ago

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand
File size: 37.8 KB
Line 
1SUBROUTINE 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
45!     ------------------------------------------------------------------
46
47#include "tsmbkind.h"
48
49USE YOMCST   , ONLY : RG
50USE YOERAD   , ONLY : NTSW     ,LONEWSW
51USE YOESW    , ONLY : APAD     ,BPAD     ,RRAY     ,RSUN     ,&
52           &RPDH1    ,RPDU1    ,RPNH     ,RPNU     ,RSWCE    ,&
53           &RSWCP    ,RTDH2O   ,RTDUMG   ,RTH2O    ,RTUMG    ,&
54           &D        ,RROMA    ,RROMB    ,RRASY    ,RHSRA    ,&
55           &RHSRB    ,RHSRC    ,RHSRD    ,RHSRE    ,RHSRF    ,&
56           &RHSRTA   ,RHSRTB   ,REXPO3   ,&
57           &RALBICE_AR         ,RALBICE_AN         ,RWEIGHT  ,&
58           &RADJUST  ,NTYPS    ,NEXPO3
59
60
61IMPLICIT NONE
62
63
64!     DUMMY INTEGER SCALARS
65INTEGER_M :: KSW
66INTEGER_M :: KTSW
67
68
69
70!     ----------------------------------------------------------------
71REAL_B :: ZAPAD2(2,3,7)  , ZBPAD2(2,3,7)  , ZD2(2,3)&
72  &,  ZRAY2(2,6)     , ZSUN2(2)       , ZSWCE2(2)  , ZSWCP2(2)&
73  &,  ZWEIGHT2(2,8)  , ZROMA2(2)      , ZROMB2(2)  , ZRASY2(2)&
74  &,  ZRA2(2)        , ZRB2(2)        , ZRC2(2)&
75  &,  ZRD2(2)        , ZRE2(2)        , ZRF2(2)    , ZALBICE2(12,2)
76
77REAL_B :: ZAPAD4(4,3,7)  , ZBPAD4(4,3,7)  , ZD4(4,3)&
78  &,  ZRAY4(4,6)     , ZSUN4(4)       , ZSWCE4(4)  , ZSWCP4(4)&
79  &,  ZWEIGHT4(4,8)  , ZROMA4(4)      , ZROMB4(4)  , ZRASY4(4)&
80  &,  ZRA4(4)        , ZRB4(4)        , ZRC4(4)&
81  &,  ZRD4(4)        , ZRE4(4)        , ZRF4(4)    , ZALBICE4(12,4)
82 
83REAL_B :: ZAPAD6(6,3,7)  , ZBPAD6(6,3,7)  , ZD6(6,3)&
84  &,  ZRAY6(6,6)     , ZSUN6(6)       , ZSWCE6(6)  , ZSWCP6(6)&
85  &,  ZWEIGHT6(6,8)  , ZROMA6(6)      , ZROMB6(6)  , ZRASY6(6)&
86  &,  ZRA6(6)        , ZRB6(6)        , ZRC6(6)&
87  &,  ZRD6(6)        , ZRE6(6)        , ZRF6(6)    , ZALBICE6(12,6)
88 
89REAL_B :: ZTWEIGHT(8) 
90
91!     LOCAL INTEGER SCALARS
92INTEGER_M :: JC3, JC6, JI, JJ, JW, K, JTYPS, IM, JM
93
94!     LOCAL REAL SCALARS
95REAL_B :: ZH2O, ZPDHNIS, ZPDH2O, ZPDUNIS,&
96          &ZPDUMG, ZPRHNIS, ZPRH2O, &
97          &ZPRUNIS, ZPRUMG, ZRTO1, ZRTO2, ZTDHNIS, &
98          &ZTDUNIS, ZTHNIS, &
99          &ZTUNIS, ZUMG, ZADJUS2, ZADJUS4, ZADJUS6
100
101
102!     ----------------------------------------------------------------
103
104!*        1.  CLEAR-SKY ABSORPTION COEFFICIENTS FOR N SPECTRAL INTERVALS
105!             --------------------------------------------------------
106
107!* spectrally independent
108
109ZTDHNIS = 0.450_JPRB
110ZTDUNIS = 0.375_JPRB
111ZTHNIS  = 296._JPRB
112ZTUNIS  = 296._JPRB
113ZPDHNIS = 0.90_JPRB
114ZPDUNIS = 0.75_JPRB
115ZPRHNIS = 100000._JPRB
116ZPRUNIS = 100000._JPRB
117
118!     ----------------------------------------------------------------
119
120!*        1.1  COEFFICIENTS FOR TWO SPECTRAL INTERVALS
121!              ---------------------------------------
122
123
124!* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
125!       H2O:  Pref=1000hPa, Tref=296K, Pdep=0.9
126!       UMG:  Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
127!       O3 :  unchanged in interval 1, from HITRAN 92 in interval 4
128
129ZADJUS2 = _ONE_
130
131!* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
132
133ZSUN2(1) = 0.459760_JPRB
134
135ZD2(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
136
137ZAPAD2(1, 1, :) = (/&
138 &0.184678379E+06_JPRB,&
139 &0.553080884E+05_JPRB,&
140 &0.248143712E+04_JPRB,&
141 &0.000000000E-00_JPRB,&
142 &0.000000000E-00_JPRB,&
143 &0.000000000E-00_JPRB,&
144 &0.000000000E-00_JPRB/)
145ZAPAD2(1, 2, :) = (/&
146 &0.715303869E+01_JPRB,&
147 &0.219386847E+03_JPRB,&
148 &0.830001089E+03_JPRB,&
149 &0.000000000E-00_JPRB,&
150 &0.000000000E-00_JPRB,&
151 &0.000000000E-00_JPRB,&
152 &0.000000000E-00_JPRB/)
153ZAPAD2(1, 3, :) = (/&
154 &0.925887084E-04_JPRB,&
155 &0.129353723E-01_JPRB,&
156 &0.800821928E+00_JPRB,&
157 &0.242715973E+02_JPRB,&
158 &0.878331486E+02_JPRB,&
159 &0.191559725E+02_JPRB,&
160 &0.000000000E+00_JPRB/)
161
162ZBPAD2(1, 1, :) = (/&
163 &0.184678379E+06_JPRB,&
164 &0.555188347E+05_JPRB,&
165 &0.253257443E+04_JPRB,&
166 &0.100000000E+01_JPRB,&
167 &0.000000000E-00_JPRB,&
168 &0.000000000E-00_JPRB,&
169 &0.000000000E-00_JPRB/)
170ZBPAD2(1, 2, :) = (/&
171 &0.715303869E+01_JPRB,&
172 &0.219441875E+03_JPRB,&
173 &0.831119997E+03_JPRB,&
174 &0.100000000E+01_JPRB,&
175 &0.000000000E-00_JPRB,&
176 &0.000000000E-00_JPRB,&
177 &0.000000000E-00_JPRB/)
178ZBPAD2(1, 3, :) = (/&
179 &0.925887084E-04_JPRB,&
180 &0.131812683E-01_JPRB,&
181 &0.812706117E+00_JPRB,&
182 &0.249863591E+02_JPRB,&
183 &0.931071925E+02_JPRB,&
184 &0.252233437E+02_JPRB,&
185 &0.100000000E+01_JPRB/)
186
187ZRAY2(1,:)= (/&
188 &.428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
189 &.522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
190
191ZSWCE2(1) = _ZERO_
192ZSWCP2(1) = _ZERO_
193
194!     ----------------------------------------------------------------
195
196!* Near-Infrared (0.69 - 4.0 Microns) is not sub-divided
197
198!     ----------------------------------------------------------------
199
200!* 0.68 - 4.00 Micron
201
202ZSUN2(2) = 0.540240_JPRB
203
204ZD2(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB/)
205
206ZAPAD2(2, 1, :) = (/&
207 &0.577160808E-06_JPRB,&
208 &0.681760991E-02_JPRB,&
209 &0.502416370E+01_JPRB,&
210 &0.209771270E+03_JPRB,&
211 &0.644493139E+03_JPRB,&
212 &0.120654163E+03_JPRB,&
213 &0.000000000E+00_JPRB/)
214ZAPAD2(2, 2, :) = (/&
215 &0.422874202E-09_JPRB,&
216 &0.257499477E-05_JPRB,&
217 &0.349570458E-02_JPRB,&
218 &0.917758876E+00_JPRB,&
219 &0.327669989E+02_JPRB,&
220 &0.874155115E+02_JPRB,&
221 &0.000000000E+00_JPRB/)
222ZAPAD2(2, 3, :) = (/&
223 &0.932012720E+02_JPRB,&
224 &0.437866983E+03_JPRB,&
225 &0.000000000E-00_JPRB,&
226 &0.000000000E-00_JPRB,&
227 &0.000000000E-00_JPRB,&
228 &0.000000000E-00_JPRB,&
229 &0.000000000E+00_JPRB/)
230
231ZBPAD2(2, 1, :) = (/&
232 &0.577160808E-06_JPRB,&
233 &0.682595589E-02_JPRB,&
234 &0.508797283E+01_JPRB,&
235 &0.223901000E+03_JPRB,&
236 &0.770640908E+03_JPRB,&
237 &0.186812480E+03_JPRB,&
238 &0.100000000E+01_JPRB/)
239ZBPAD2(2, 2, :) = (/&
240 &0.422874202E-09_JPRB,&
241 &0.257567242E-05_JPRB,&
242 &0.349912408E-02_JPRB,&
243 &0.920770481E+00_JPRB,&
244 &0.330990960E+02_JPRB,&
245 &0.903249595E+02_JPRB,&
246 &0.100000000E+01_JPRB/)
247ZBPAD2(2, 3, :) = (/&
248 &0.932012720E+02_JPRB,&
249 &0.439078599E+03_JPRB,&
250 &0.100000000E+01_JPRB,&
251 &0.000000000E+00_JPRB,&
252 &0.000000000E+00_JPRB,&
253 &0.000000000E+00_JPRB,&
254 &0.000000000E+00_JPRB/)
255
256ZRAY2(2,:)= (/&
257 &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,&
258 &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/)
259
260ZSWCE2(2) = _ZERO_
261ZSWCP2(2) = _ZERO_
262
263!=====================================================================
264
265!*        1.2  COEFFICIENTS FOR FOUR SPECTRAL INTERVALS
266!              ----------------------------------------
267
268
269!* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
270!       H2O:  Pref=1000hPa, Tref=296K, Pdep=0.9
271!       UMG:  Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
272!       O3 :  unchanged in interval 1, from HITRAN 92 in interval 4
273
274ZADJUS4 = _ONE_
275
276!* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
277
278ZSUN4(1) = 0.459760_JPRB
279
280ZD4(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
281
282ZAPAD4(1, 1, :) = (/&
283 &0.184678379E+06_JPRB,&
284 &0.553080884E+05_JPRB,&
285 &0.248143712E+04_JPRB,&
286 &0.000000000E-00_JPRB,&
287 &0.000000000E-00_JPRB,&
288 &0.000000000E-00_JPRB,&
289 &0.000000000E-00_JPRB/)
290ZAPAD4(1, 2, :) = (/&
291 &0.715303869E+01_JPRB,&
292 &0.219386847E+03_JPRB,&
293 &0.830001089E+03_JPRB,&
294 &0.000000000E-00_JPRB,&
295 &0.000000000E-00_JPRB,&
296 &0.000000000E-00_JPRB,&
297 &0.000000000E-00_JPRB/)
298ZAPAD4(1, 3, :) = (/&
299 &0.925887084E-04_JPRB,&
300 &0.129353723E-01_JPRB,&
301 &0.800821928E+00_JPRB,&
302 &0.242715973E+02_JPRB,&
303 &0.878331486E+02_JPRB,&
304 &0.191559725E+02_JPRB,&
305 &0.000000000E+00_JPRB/)
306
307ZBPAD4(1, 1, :) = (/&
308 &0.184678379E+06_JPRB,&
309 &0.555188347E+05_JPRB,&
310 &0.253257443E+04_JPRB,&
311 &0.100000000E+01_JPRB,&
312 &0.000000000E-00_JPRB,&
313 &0.000000000E-00_JPRB,&
314 &0.000000000E-00_JPRB/)
315ZBPAD4(1, 2, :) = (/&
316 &0.715303869E+01_JPRB,&
317 &0.219441875E+03_JPRB,&
318 &0.831119997E+03_JPRB,&
319 &0.100000000E+01_JPRB,&
320 &0.000000000E-00_JPRB,&
321 &0.000000000E-00_JPRB,&
322 &0.000000000E-00_JPRB/)
323ZBPAD4(1, 3, :) = (/&
324 &0.925887084E-04_JPRB,&
325 &0.131812683E-01_JPRB,&
326 &0.812706117E+00_JPRB,&
327 &0.249863591E+02_JPRB,&
328 &0.931071925E+02_JPRB,&
329 &0.252233437E+02_JPRB,&
330 &0.100000000E+01_JPRB/)
331
332ZRAY4(1,:)= (/&
333 &.428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
334 &.522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
335
336ZSWCE4(1) = _ZERO_
337ZSWCP4(1) = _ZERO_
338
339!     ----------------------------------------------------------------
340
341!* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
342
343!     ----------------------------------------------------------------
344
345!* 0.69 - 1.19 Micron
346
347ZSUN4(2) = 0.326158_JPRB
348
349ZD4(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /)
350
351ZAPAD4(2, 1, :) = (/&
352 &0.690730834E-02_JPRB,&
353 &0.151704275E+01_JPRB,&
354 &0.751477543E+02_JPRB,&
355 &0.759770236E+03_JPRB,&
356 &0.109800326E+04_JPRB,&
357 &0.148407574E+03_JPRB,&
358 &0.000000000E+00_JPRB/)
359ZAPAD4(2, 2, :) = (/&
360 &0.863790752E-03_JPRB,&
361 &0.448762291E+00_JPRB,&
362 &0.332530367E+02_JPRB,&
363 &0.190914146E+03_JPRB,&
364 &0.000000000E+00_JPRB,&
365 &0.000000000E+00_JPRB,&
366 &0.000000000E+00_JPRB/)
367ZAPAD4(2, 3, :) = (/&
368 &0.000000000E+00_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 &0.000000000E+00_JPRB/)
375
376ZBPAD4(2, 1, :) = (/&
377 &0.690730834E-02_JPRB,&
378 &0.151954406E+01_JPRB,&
379 &0.756512527E+02_JPRB,&
380 &0.779384997E+03_JPRB,&
381 &0.121113108E+04_JPRB,&
382 &0.207678436E+03_JPRB,&
383 &0.100000000E+01_JPRB/)
384ZBPAD4(2, 2, :) = (/&
385 &0.863790752E-03_JPRB,&
386 &0.448948107E+00_JPRB,&
387 &0.333186750E+02_JPRB,&
388 &0.192727216E+03_JPRB,&
389 &0.100000000E+01_JPRB,&
390 &0.000000000E+00_JPRB,&
391 &0.000000000E+00_JPRB/)
392ZBPAD4(2, 3, :) = (/&
393 &1.000000000E+00_JPRB,&
394 &0.000000000E+00_JPRB,&
395 &0.000000000E+00_JPRB,&
396 &0.000000000E+00_JPRB,&
397 &0.000000000E+00_JPRB,&
398 &0.000000000E+00_JPRB,&
399 &0.000000000E+00_JPRB/)
400 
401DATA (ZRAY4(2,K),K=1,6) /&
402 &.164261E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
403 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/
404
405ZSWCE4(2) = _ZERO_
406ZSWCP4(2) = _ZERO_
407
408!     ----------------------------------------------------------------
409
410!* 1.19 - 2.38 Microns
411
412ZSUN4(3) = 0.180608_JPRB
413
414ZD4(3,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /)
415
416ZAPAD4(3, 1, :) = (/&
417 &0.837531303E-05_JPRB,&
418 &0.173886341E-01_JPRB,&
419 &0.518852799E+01_JPRB,&
420 &0.159078416E+03_JPRB,&
421 &0.493273523E+03_JPRB,&
422 &0.102567293E+03_JPRB,&
423 &0.000000000E+00_JPRB/)
424ZAPAD4(3, 2, :) = (/&
425 &0.657978575E-02_JPRB,&
426 &0.752617872E+00_JPRB,&
427 &0.158209734E+02_JPRB,&
428 &0.410274915E+02_JPRB,&
429 &0.000000000E+00_JPRB,&
430 &0.000000000E+00_JPRB,&
431 &0.000000000E+00_JPRB/)
432ZAPAD4(3, 3, :) = (/&
433 &0.000000000E+00_JPRB,&
434 &0.000000000E+00_JPRB,&
435 &0.000000000E+00_JPRB,&
436 &0.000000000E+00_JPRB,&
437 &0.000000000E+00_JPRB,&
438 &0.000000000E+00_JPRB,&
439 &0.000000000E+00_JPRB/)
440
441ZBPAD4(3, 1, :) = (/&
442 &0.837531303E-05_JPRB,&
443 &0.174882536E-01_JPRB,&
444 &0.534536580E+01_JPRB,&
445 &0.180351767E+03_JPRB,&
446 &0.673126838E+03_JPRB,&
447 &0.182718543E+03_JPRB,&
448 &0.100000000E+01_JPRB/)
449ZBPAD4(3, 2, :) = (/&
450 &0.657978575E-02_JPRB,&
451 &0.753752065E+00_JPRB,&
452 &0.159286262E+02_JPRB,&
453 &0.424278450E+02_JPRB,&
454 &0.100000000E+01_JPRB,&
455 &0.000000000E+00_JPRB,&
456 &0.000000000E+00_JPRB/)
457ZBPAD4(3, 3, :) = (/&
458 &1.000000000E+00_JPRB,&
459 &0.000000000E+00_JPRB,&
460 &0.000000000E+00_JPRB,&
461 &0.000000000E+00_JPRB,&
462 &0.000000000E+00_JPRB,&
463 &0.000000000E+00_JPRB,&
464 &0.000000000E+00_JPRB/)
465
466ZRAY4(3,:)= (/&
467 &.180438E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
468 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
469
470ZSWCE4(3) = _ZERO_
471ZSWCP4(3) = _ZERO_
472
473!     ----------------------------------------------------------------
474
475!* 2.38 - 4.00 Microns
476
477ZSUN4(4) = 0.033474_JPRB
478
479ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
480
481ZAPAD4(4, 1, :) = (/&
482 &0.122118185E-06_JPRB,&
483 &0.154042531E-02_JPRB,&
484 &0.141152193E+01_JPRB,&
485 &0.685368761E+02_JPRB,&
486 &0.216522281E+03_JPRB,&
487 &0.421228746E+02_JPRB,&
488 &0.000000000E+00_JPRB/)
489ZAPAD4(4, 2, :) = (/&
490 &0.364233560E-10_JPRB,&
491 &0.217340835E-06_JPRB,&
492 &0.292623386E-03_JPRB,&
493 &0.797100631E-01_JPRB,&
494 &0.319103672E+01_JPRB,&
495 &0.110530283E+02_JPRB,&
496 &0.000000000E+00_JPRB/)
497ZAPAD4(4, 3, :) = (/&
498 &0.263068898E+02_JPRB,&
499 &0.146425875E+03_JPRB,&
500 &0.860137809E+02_JPRB,&
501 &0.000000000E+00_JPRB,&
502 &0.000000000E+00_JPRB,&
503 &0.000000000E+00_JPRB,&
504 &0.000000000E+00_JPRB/)
505
506ZBPAD4(4, 1, :) = (/&
507 &0.122118185E-06_JPRB,&
508 &0.156287582E-02_JPRB,&
509 &0.156994562E+01_JPRB,&
510 &0.102304103E+03_JPRB,&
511 &0.475479878E+03_JPRB,&
512 &0.188714799E+03_JPRB,&
513 &0.100000000E+01_JPRB/)
514ZBPAD4(4, 2, :) = (/&
515 &0.364233560E-10_JPRB,&
516 &0.218265806E-06_JPRB,&
517 &0.297085917E-03_JPRB,&
518 &0.834253403E-01_JPRB,&
519 &0.358290002E+01_JPRB,&
520 &0.139206164E+02_JPRB,&
521 &0.100000000E+01_JPRB/)
522ZBPAD4(4, 3, :) = (/&
523 &0.263068898E+02_JPRB,&
524 &0.152569217E+03_JPRB,&
525 &0.976791971E+02_JPRB,&
526 &0.100000000E+01_JPRB,&
527 &0.000000000E+00_JPRB,&
528 &0.000000000E+00_JPRB,&
529 &0.000000000E+00_JPRB/)
530
531ZRAY4(4,:)= (/&
532 &.136020E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
533 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
534
535ZSWCE4(4) = _ZERO_
536ZSWCP4(4) = _ZERO_
537
538!=====================================================================
539
540!*        1.3  COEFFICIENTS FOR SIX SPECTRAL INTERVALS
541!              ---------------------------------------
542
543
544!* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
545!       H2O:  Pref=1000hPa, Tref=296K, Pdep=0.9
546!       UMG:  Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
547!       O3 :  unchanged in interval 1, from HITRAN 92 in interval 4
548
549ZADJUS6 = 0.9905709_JPRB
550
551!* 1st spectral interval: U.V.  (0.18 - 0.25 Micron)
552
553ZSUN6(1) = 0.0019124_JPRB
554
555ZD6(1,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
556
557ZAPAD6(1, 1, :) = (/&
558 &0.000000000E-00_JPRB,&
559 &0.000000000E-00_JPRB,&
560 &0.000000000E-00_JPRB,&
561 &0.000000000E-00_JPRB,&
562 &0.000000000E-00_JPRB,&
563 &0.000000000E-00_JPRB,&
564 &0.000000000E-00_JPRB/)
565ZAPAD6(1, 2, :) = (/&
566 &0.000000000E-00_JPRB,&
567 &0.000000000E-00_JPRB,&
568 &0.000000000E-00_JPRB,&
569 &0.000000000E-00_JPRB,&
570 &0.000000000E-00_JPRB,&
571 &0.000000000E-00_JPRB,&
572 &0.000000000E-00_JPRB/)
573ZAPAD6(1, 3, :) = (/&
574 &0.000000000E-00_JPRB,&
575 &0.000000000E-00_JPRB,&
576 &0.000000000E-00_JPRB,&
577 &0.000000000E-00_JPRB,&
578 &0.000000000E-00_JPRB,&
579 &0.000000000E-00_JPRB,&
580 &0.000000000E-00_JPRB/)
581 
582ZBPAD6(1, 1, :) = (/&
583 &0.100000000E+01_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 &0.000000000E-00_JPRB/)
590ZBPAD6(1, 2, :) = (/&
591 &0.100000000E+01_JPRB,&
592 &0.000000000E-00_JPRB,&
593 &0.000000000E-00_JPRB,&
594 &0.000000000E-00_JPRB,&
595 &0.000000000E-00_JPRB,&
596 &0.000000000E-00_JPRB,&
597 &0.000000000E-00_JPRB/)
598ZBPAD6(1, 3, :) = (/&
599 &0.100000000E+01_JPRB,&
600 &0.000000000E-00_JPRB,&
601 &0.000000000E-00_JPRB,&
602 &0.000000000E-00_JPRB,&
603 &0.000000000E-00_JPRB,&
604 &0.000000000E-00_JPRB,&
605 &0.000000000E-00_JPRB/)
606 
607ZRAY6(1,:)= (/&
608 &.43959002E+01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
609 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
610
611ZSWCE6(1) = _ZERO_
612ZSWCP6(1) = _ZERO_
613
614NEXPO3(1) = 6
615REXPO3(1, 1, :) =(/&
616  &0.000334E+00_JPRB, 0.011131E+00_JPRB, 0.000064E+00_JPRB,&
617  &0.218900E+00_JPRB, 0.134118E+00_JPRB, 0.635453E+00_JPRB,&
618  &0.000000E+00_JPRB /)
619REXPO3(1, 2, :) =(/&
620  &0.467737E+03_JPRB, 0.389056E+02_JPRB, 0.223872E+04_JPRB,&
621  &0.102346E+02_JPRB, 0.794561E+01_JPRB, 0.138054E+02_JPRB,&
622  &0.000000E+00_JPRB /)
623 
624
625
626!* 2nd spectral interval: U.V.  (0.25 - 0.44 Micron)
627
628ZSUN6(2) = 0.1353560_JPRB
629
630ZD6(2,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
631
632ZAPAD6(2, 1, :) = (/&
633 &0.000000000E-00_JPRB,&
634 &0.000000000E-00_JPRB,&
635 &0.000000000E-00_JPRB,&
636 &0.000000000E-00_JPRB,&
637 &0.000000000E-00_JPRB,&
638 &0.000000000E-00_JPRB,&
639 &0.000000000E-00_JPRB/)
640ZAPAD6(2, 2, :) = (/&
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 &0.000000000E-00_JPRB,&
647 &0.000000000E-00_JPRB/)
648ZAPAD6(2, 3, :) = (/&
649 &0.000000000E-00_JPRB,&
650 &0.000000000E-00_JPRB,&
651 &0.000000000E-00_JPRB,&
652 &0.000000000E-00_JPRB,&
653 &0.000000000E-00_JPRB,&
654 &0.000000000E-00_JPRB,&
655 &0.000000000E-00_JPRB/)
656 
657ZBPAD6(2, 1, :) = (/&
658 &0.100000000E+01_JPRB,&
659 &0.000000000E-00_JPRB,&
660 &0.000000000E-00_JPRB,&
661 &0.000000000E-00_JPRB,&
662 &0.000000000E-00_JPRB,&
663 &0.000000000E-00_JPRB,&
664 &0.000000000E-00_JPRB/)
665ZBPAD6(2, 2, :) = (/&
666 &0.100000000E+01_JPRB,&
667 &0.000000000E-00_JPRB,&
668 &0.000000000E-00_JPRB,&
669 &0.000000000E-00_JPRB,&
670 &0.000000000E-00_JPRB,&
671 &0.000000000E-00_JPRB,&
672 &0.000000000E-00_JPRB/)
673ZBPAD6(2, 3, :) = (/&
674 &0.100000000E+01_JPRB,&
675 &0.000000000E-00_JPRB,&
676 &0.000000000E-00_JPRB,&
677 &0.000000000E-00_JPRB,&
678 &0.000000000E-00_JPRB,&
679 &0.000000000E-00_JPRB,&
680 &0.000000000E-00_JPRB/)
681 
682ZRAY6(2,:)= (/&
683 &.55503070E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
684 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
685
686ZSWCE6(2) = _ZERO_
687ZSWCP6(2) = _ZERO_
688
689NEXPO3(2) = 7
690REXPO3(2, 1, :) =(/&
691  &0.043801E+00_JPRB, 0.078893E+00_JPRB, 0.036839E+00_JPRB,&
692  &0.022503E+00_JPRB, 0.042333E+00_JPRB, 0.037870E+00_JPRB,&
693  &0.737762E+00_JPRB /)
694REXPO3(2, 2, :) =(/&
695  &0.234249E+01_JPRB, 0.125170E+00_JPRB, 0.549527E+02_JPRB,&
696  &0.257041E+03_JPRB, 0.476838E+00_JPRB, 0.911993E+01_JPRB,&
697  &0.000000E+00_JPRB /)
698
699
700!* 3rd spectral interval: Visible  (0.44 - 0.69 Micron)
701
702ZSUN6(3) = 0.3213000_JPRB
703
704ZD6(3,:)= (/ 0.800000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
705
706ZAPAD6(3, 1, :) = (/&
707 &0.3682644E+02_JPRB,&
708 &0.6491415E+03_JPRB,&
709 &0.8720916E+02_JPRB,&
710 &0.0000000E-00_JPRB,&
711 &0.0000000E-00_JPRB,&
712 &0.0000000E-00_JPRB,&
713 &0.0000000E-00_JPRB/)
714ZAPAD6(3, 2, :) = (/&
715 &0.0000000E-00_JPRB,&
716 &0.0000000E-00_JPRB,&
717 &0.0000000E-00_JPRB,&
718 &0.0000000E-00_JPRB,&
719 &0.0000000E-00_JPRB,&
720 &0.0000000E-00_JPRB,&
721 &0.0000000E-00_JPRB/)
722ZAPAD6(3, 3, :) = (/&
723 &0.000000000E-00_JPRB,&
724 &0.000000000E-00_JPRB,&
725 &0.000000000E-00_JPRB,&
726 &0.000000000E-00_JPRB,&
727 &0.000000000E-00_JPRB,&
728 &0.000000000E-00_JPRB,&
729 &0.000000000E-00_JPRB/)
730 
731ZBPAD6(3, 1, :) = (/&
732 &0.3682644E+02_JPRB,&
733 &0.6495968E+03_JPRB,&
734 &0.9637657E+02_JPRB,&
735 &0.1000000E+01_JPRB,&
736 &0.0000000E-00_JPRB,&
737 &0.0000000E-00_JPRB,&
738 &0.0000000E-00_JPRB/)
739ZBPAD6(3, 2, :) = (/&
740 &0.1000000E+01_JPRB,&
741 &0.0000000E-00_JPRB,&
742 &0.0000000E-00_JPRB,&
743 &0.0000000E-00_JPRB,&
744 &0.0000000E-00_JPRB,&
745 &0.0000000E-00_JPRB,&
746 &0.0000000E-00_JPRB/)
747ZBPAD6(3, 3, :) = (/&
748 &0.100000000E+01_JPRB,&
749 &0.000000000E-00_JPRB,&
750 &0.000000000E-00_JPRB,&
751 &0.000000000E-00_JPRB,&
752 &0.000000000E-00_JPRB,&
753 &0.000000000E-00_JPRB,&
754 &0.000000000E-00_JPRB/)
755 
756ZRAY6(3,:)= (/&
757 &.10528199E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
758 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
759
760ZSWCE6(3) = _ZERO_
761ZSWCP6(3) = _ZERO_
762
763NEXPO3(3) = 6
764REXPO3(3, 1, :) =(/&
765  &0.063442E+00_JPRB, 0.058550E+00_JPRB, 0.237534E+00_JPRB,&
766  &0.412292E+00_JPRB, 0.126141E+00_JPRB, 0.102041E+00_JPRB,&
767  &0.000000E+00_JPRB /)
768REXPO3(3, 2, :) =(/&
769  &0.125170E+00_JPRB, 0.119209E-01_JPRB, 0.119209E+00_JPRB,&
770  &0.417233E-01_JPRB, 0.894070E-01_JPRB, 0.000000E+00_JPRB,&
771  &0.000000E+00_JPRB /)
772
773!     ----------------------------------------------------------------
774
775!* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
776
777!     ----------------------------------------------------------------
778
779!* 0.69 - 1.19 Micron
780
781!* UMG is O2 only
782
783ZSUN6(4) = 0.3211833_JPRB
784
785ZD6(4,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
786
787ZAPAD6(4, 1, :) = (/&
788 &0.9969941E+02_JPRB,&
789 &0.7496307E+04_JPRB,&
790 &0.7003009E+05_JPRB,&
791 &0.9381446E+05_JPRB,&
792 &0.1827735E+05_JPRB,&
793 &0.3459400E+03_JPRB,&
794 &0.0000000E+00_JPRB/)
795ZAPAD6(4, 2, :) = (/&
796 &0.1567941E+10_JPRB,&
797 &0.3094160E+09_JPRB,&
798 &0.3650192E+07_JPRB,&
799 &0.5026130E+04_JPRB,&
800 &0.0000000E+00_JPRB,&
801 &0.0000000E+00_JPRB,&
802 &0.0000000E+00_JPRB/)
803ZAPAD6(4, 3, :) = (/&
804 &0.000000000E+00_JPRB,&
805 &0.000000000E+00_JPRB,&
806 &0.000000000E+00_JPRB,&
807 &0.000000000E+00_JPRB,&
808 &0.000000000E+00_JPRB,&
809 &0.000000000E+00_JPRB,&
810 &0.000000000E+00_JPRB/)
811
812ZBPAD6(4, 1, :) = (/&
813 &0.9969941E+02_JPRB,&
814 &0.7517916E+04_JPRB,&
815 &0.7149485E+05_JPRB,&
816 &0.1020564E+06_JPRB,&
817 &0.2337706E+05_JPRB,&
818 &0.6874490E+03_JPRB,&
819 &0.1000000E+01_JPRB/)
820ZBPAD6(4, 2, :) = (/&
821 &0.1567941E+10_JPRB,&
822 &0.3332380E+09_JPRB,&
823 &0.4940093E+07_JPRB,&
824 &0.1026382E+05_JPRB,&
825 &0.1000000E+01_JPRB,&
826 &0.0000000E+00_JPRB,&
827 &0.0000000E+00_JPRB/)
828ZBPAD6(4, 3, :) = (/&
829 &1.000000000E+00_JPRB,&
830 &0.000000000E+00_JPRB,&
831 &0.000000000E+00_JPRB,&
832 &0.000000000E+00_JPRB,&
833 &0.000000000E+00_JPRB,&
834 &0.000000000E+00_JPRB,&
835 &0.000000000E+00_JPRB/)
836 
837DATA (ZRAY6(4,K),K=1,6) /&
838 &.16436996E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
839 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/
840
841ZSWCE6(4) = _ZERO_
842ZSWCP6(4) = _ZERO_
843
844NEXPO3(4) = 4
845REXPO3(4, 1, :) =(/&
846  &0.000074E+00_JPRB, 0.320194E+00_JPRB, 0.082915E+00_JPRB,&
847  &0.596816E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
848  &0.000000E+00_JPRB /)
849REXPO3(4, 2, :) =(/&
850  &0.232458E+00_JPRB, 0.119209E-01_JPRB, 0.178814E-01_JPRB,&
851  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
852  &0.000000E+00_JPRB /)
853
854
855!     ----------------------------------------------------------------
856
857!* 1.19 - 2.38 Microns
858
859!* UMG is CO2 only
860
861ZSUN6(5) = 0.1778533_JPRB
862
863ZD6(5,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
864
865ZAPAD6(5, 1, :) = (/&
866 &0.3325841E-03_JPRB,&
867 &0.6194496E+00_JPRB,&
868 &0.1497138E+03_JPRB,&
869 &0.2314864E+04_JPRB,&
870 &0.2380109E+04_JPRB,&
871 &0.9553823E+02_JPRB,&
872 &0.0000000E+00_JPRB/)
873ZAPAD6(5, 2, :) = (/&
874 &0.2419549E-02_JPRB,&
875 &0.1504652E+01_JPRB,&
876 &0.1339335E+02_JPRB,&
877 &0.8932533E+01_JPRB,&
878 &0.0000000E+00_JPRB,&
879 &0.0000000E+00_JPRB,&
880 &0.0000000E+00_JPRB/)
881ZAPAD6(5, 3, :) = (/&
882 &0.000000000E+00_JPRB,&
883 &0.000000000E+00_JPRB,&
884 &0.000000000E+00_JPRB,&
885 &0.000000000E+00_JPRB,&
886 &0.000000000E+00_JPRB,&
887 &0.000000000E+00_JPRB,&
888 &0.000000000E+00_JPRB/)
889
890ZBPAD6(5, 1, :) = (/&
891 &0.3325841E-03_JPRB,&
892 &0.6231947E+00_JPRB,&
893 &0.1553098E+03_JPRB,&
894 &0.2822458E+04_JPRB,&
895 &0.3885194E+04_JPRB,&
896 &0.2700235E+03_JPRB,&
897 &0.1000000E+01_JPRB/)
898ZBPAD6(5, 2, :) = (/&
899 &0.2419549E-02_JPRB,&
900 &0.1507850E+01_JPRB,&
901 &0.1489144E+03_JPRB,&
902 &0.1425322E+02_JPRB,&
903 &0.1000000E+01_JPRB,&
904 &0.0000000E+00_JPRB,&
905 &0.0000000E+00_JPRB/)
906ZBPAD6(5, 3, :) = (/&
907 &1.000000000E+00_JPRB,&
908 &0.000000000E+00_JPRB,&
909 &0.000000000E+00_JPRB,&
910 &0.000000000E+00_JPRB,&
911 &0.000000000E+00_JPRB,&
912 &0.000000000E+00_JPRB,&
913 &0.000000000E+00_JPRB/)
914
915ZRAY6(5,:)= (/&
916 &.18073079E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
917 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
918
919ZSWCE6(5) = _ZERO_
920ZSWCP6(5) = _ZERO_
921
922NEXPO3(5) = 0
923REXPO3(5, 1, :) =(/&
924  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
925  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
926  &0.000000E+00_JPRB /)
927REXPO3(5, 2, :) =(/&
928  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
929  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
930  &0.000000E+00_JPRB /)
931
932!     ----------------------------------------------------------------
933
934!* 2.38 - 4.00 Microns
935
936ZSUN6(6) = 0.0329659_JPRB
937
938ZD6(6,:)= (/ 0.000000000_JPRB, 0.600000000_JPRB, 0.000000000_JPRB /)
939
940ZAPAD6(6, 1, :) = (/&
941 &0.2122889E-06_JPRB,&
942 &0.9030576E-03_JPRB,&
943 &0.2431282E+00_JPRB,&
944 &0.4901345E+01_JPRB,&
945 &0.3996347E+01_JPRB,&
946 &0.3910227E+01_JPRB,&
947 &0.0000000E+00_JPRB/)
948ZAPAD6(6, 2, :) = (/&
949 &0.2051629E-02_JPRB,&
950 &0.8053325E+00_JPRB,&
951 &0.1826847E+02_JPRB,&
952 &0.1991453E+02_JPRB,&
953 &0.0000000E+00_JPRB,&
954 &0.0000000E+00_JPRB,&
955 &0.0000000E+00_JPRB/)
956ZAPAD6(6, 3, :) = (/&
957 &0.263068898E+02_JPRB,&
958 &0.146425875E+03_JPRB,&
959 &0.860137809E+02_JPRB,&
960 &0.000000000E+00_JPRB,&
961 &0.000000000E+00_JPRB,&
962 &0.000000000E+00_JPRB,&
963 &0.000000000E+00_JPRB/)
964
965ZBPAD6(6, 1, :) = (/&
966 &0.2122889E-06_JPRB,&
967 &0.9379083E-03_JPRB,&
968 &0.2957335E+00_JPRB,&
969 &0.8747190E+01_JPRB,&
970 &0.1015794E+02_JPRB,&
971 &0.1361277E+02_JPRB,&
972 &0.1000000E+01_JPRB/)
973ZBPAD6(6, 2, :) = (/&
974 &0.2051629E-02_JPRB,&
975 &0.9455932E+00_JPRB,&
976 &0.3030076E+02_JPRB,&
977 &0.4307387E+02_JPRB,&
978 &0.1000000E+01_JPRB,&
979 &0.0000000E+00_JPRB,&
980 &0.0000000E+00_JPRB/)
981ZBPAD6(6, 3, :) = (/&
982 &0.263068898E+02_JPRB,&
983 &0.152569217E+03_JPRB,&
984 &0.976791971E+02_JPRB,&
985 &0.100000000E+01_JPRB,&
986 &0.000000000E+00_JPRB,&
987 &0.000000000E+00_JPRB,&
988 &0.000000000E+00_JPRB/)
989
990ZRAY6(6,:)= (/&
991 &.13618247E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
992 &.000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
993
994ZSWCE6(6) = _ZERO_
995ZSWCP6(6) = _ZERO_
996
997NEXPO3(6) = 0
998REXPO3(6, 1, :) =(/&
999  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1000  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1001  &0.000000E+00_JPRB /)
1002REXPO3(6, 2, :) =(/&
1003  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1004  &0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
1005  &0.000000E+00_JPRB /)
1006
1007!=====================================================================
1008!     ----------------------------------------------------------------
1009
1010!*    2.3   WEIGHTS FOR SPECTRAL SURFACE ALBEDO
1011!           Water
1012!           Sea-Ice (Ebert, Curry, 1993)
1013!           Wet skin
1014!           Low vegetation (snow free) (BR, 1982)
1015!           Snow on low vegetation (Warren, 1982)
1016!           High vegetation (snow free) (BR, 1982)
1017!           Snow under high vegetation  (Warren, 1982)
1018!           Bare soil (Briegleb & Ramanathan, 1982)
1019
1020NTYPS = 8
1021ZWEIGHT2(:,1)= (/ 1._JPRB, 1._JPRB /)
1022ZWEIGHT2(:,2)= (/ 1._JPRB, 1._JPRB /)
1023ZWEIGHT2(:,3)= (/ 1._JPRB, 1._JPRB /)
1024ZWEIGHT2(:,4)= (/ 1._JPRB, 3._JPRB /)
1025ZWEIGHT2(:,5)= (/ 0.920_JPRB, 0.537_JPRB /)
1026ZWEIGHT2(:,6)= (/ 1._JPRB, 3._JPRB /)
1027ZWEIGHT2(:,7)= (/ 0.860_JPRB, 0.432_JPRB /)
1028ZWEIGHT2(:,8)= (/ 1._JPRB, 2._JPRB /)
1029
1030ZWEIGHT4(:,1)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB /)
1031ZWEIGHT4(:,2)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB /)
1032ZWEIGHT4(:,3)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB /)
1033ZWEIGHT4(:,4)= (/ 1._JPRB, 4._JPRB, 2._JPRB, 1._JPRB /)
1034ZWEIGHT4(:,5)= (/ 0.920_JPRB, 0.798_JPRB , 0.159_JPRB, 0.010_JPRB /)
1035ZWEIGHT4(:,6)= (/ 1._JPRB, 4._JPRB, 2._JPRB, 1._JPRB /)
1036ZWEIGHT4(:,7)= (/ 0.860_JPRB, 0.664_JPRB , 0.092_JPRB, 0.010_JPRB /)
1037ZWEIGHT4(:,8)= (/ 1._JPRB, 2._JPRB, 2._JPRB, 1._JPRB /)
1038
1039ZWEIGHT6(:,1)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB /)
1040ZWEIGHT6(:,2)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB /)
1041ZWEIGHT6(:,3)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB, 1._JPRB /)
1042ZWEIGHT6(:,4)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 4._JPRB, 2._JPRB, 1._JPRB /)
1043ZWEIGHT6(:,5)= (/ 0.920_JPRB, 0.920_JPRB, 0.920_JPRB, 0.798_JPRB , 0.159_JPRB, 0.010_JPRB /)
1044ZWEIGHT6(:,6)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 4._JPRB, 2._JPRB, 1._JPRB /)
1045ZWEIGHT6(:,7)= (/ 0.860_JPRB, 0.860_JPRB, 0.860_JPRB, 0.664_JPRB , 0.092_JPRB, 0.010_JPRB /)
1046ZWEIGHT6(:,8)= (/ 1._JPRB, 1._JPRB, 1._JPRB, 2._JPRB, 2._JPRB, 1._JPRB /)
1047
1048!*  For sea ice, monthly values are based on Ebert and Curry, 1993, Table 2.
1049!   We take dry snow albedo as the representative value for non-summer
1050!   months, and bare sea-ice as the representative value for summer
1051!   months. The values for Antarctic are shifted six-months.
1052
1053!- 2-spectral intervals
1054
1055!*  Sea ice surf. albedo for 0.25-0.69 micron (snow covered; Ebert and Curry, 1993)
1056ZALBICE2(1:12,1) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
1057                   &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
1058                   &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
1059!*  Sea ice surf. albedo for 0.69-4.00 microns (snow covered; Ebert and Curry, 1993)
1060!ZALBICE2(1:12,2) = (/0.664_JPRB,0.664_JPRB,0.664_JPRB,0.664_JPRB,&
1061!                    &0.664_JPRB,0.476_JPRB,0.288_JPRB,0.288_JPRB,&
1062!                    &0.664_JPRB,0.664_JPRB,0.664_JPRB,0.664_JPRB/)
1063ZALBICE2(1:12,2) = (/0.587_JPRB,0.587_JPRB,0.587_JPRB,0.587_JPRB,&
1064                    &0.587_JPRB,0.438_JPRB,0.288_JPRB,0.288_JPRB,&
1065                    &0.587_JPRB,0.587_JPRB,0.587_JPRB,0.587_JPRB/)
1066                     
1067!- 4-spectral intervals
1068
1069!*  Sea ice surf. albedo for 0.25-0.69 micron (snow covered; Ebert and Curry, 1993)
1070ZALBICE4(1:12,1) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
1071                    &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
1072                    &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
1073!*  Sea ice surf. albedo for 0.69-1.19 micron (snow covered; Ebert and Curry, 1993)
1074ZALBICE4(1:12,2) = (/0.832_JPRB,0.832_JPRB,0.832_JPRB,0.832_JPRB,&
1075                    &0.832_JPRB,0.638_JPRB,0.443_JPRB,0.443_JPRB,&
1076                    &0.832_JPRB,0.832_JPRB,0.832_JPRB,0.832_JPRB/)
1077!*  Sea ice surf. albedo for 1.19-2.38 micron (snow covered; Ebert and Curry, 1993)
1078ZALBICE4(1:12,3) = (/0.250_JPRB,0.250_JPRB,0.250_JPRB,0.250_JPRB,&
1079                    &0.250_JPRB,0.153_JPRB,0.055_JPRB,0.055_JPRB,&
1080                    &0.250_JPRB,0.250_JPRB,0.250_JPRB,0.250_JPRB/)
1081!*  Sea ice surf. albedo for 2.38-4.00 microns (snow covered; Ebert and Curry, 1993)
1082ZALBICE4(1:12,4) = (/0.025_JPRB,0.025_JPRB,0.025_JPRB,0.025_JPRB,&
1083                    &0.025_JPRB,0.030_JPRB,0.036_JPRB,0.036_JPRB,&
1084                    &0.025_JPRB,0.025_JPRB,0.025_JPRB,0.025_JPRB/)
1085                     
1086!- 6-spectral intervals
1087
1088!*  Sea ice surf. albedo for 0.185-0.25 micron (snow covered; Ebert and Curry, 1993)
1089ZALBICE6(1:12,1) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
1090                    &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
1091                    &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
1092!*  Sea ice surf. albedo for 0.25-0.44 micron (snow covered; Ebert and Curry, 1993)
1093ZALBICE6(1:12,2) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
1094                    &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
1095                    &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
1096!*  Sea ice surf. albedo for 0.44-0.69 micron (snow covered; Ebert and Curry, 1993)
1097ZALBICE6(1:12,3) = (/0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB,&
1098                    &0.975_JPRB,0.876_JPRB,0.778_JPRB,0.778_JPRB,&
1099                    &0.975_JPRB,0.975_JPRB,0.975_JPRB,0.975_JPRB/)
1100!*  Sea ice surf. albedo for 0.69-1.19 micron (snow covered; Ebert and Curry, 1993)
1101ZALBICE6(1:12,4) = (/0.832_JPRB,0.832_JPRB,0.832_JPRB,0.832_JPRB,&
1102                    &0.832_JPRB,0.638_JPRB,0.443_JPRB,0.443_JPRB,&
1103                    &0.832_JPRB,0.832_JPRB,0.832_JPRB,0.832_JPRB/)
1104!*  Sea ice surf. albedo for 1.19-2.38 micron (snow covered; Ebert and Curry, 1993)
1105ZALBICE6(1:12,5) = (/0.250_JPRB,0.250_JPRB,0.250_JPRB,0.250_JPRB,&
1106                    &0.250_JPRB,0.153_JPRB,0.055_JPRB,0.055_JPRB,&
1107                    &0.250_JPRB,0.250_JPRB,0.250_JPRB,0.250_JPRB/)
1108!*  Sea ice surf. albedo for 2.38-4.00 microns (snow covered; Ebert and Curry, 1993)
1109ZALBICE6(1:12,6) = (/0.025_JPRB,0.025_JPRB,0.025_JPRB,0.025_JPRB,&
1110                    &0.025_JPRB,0.030_JPRB,0.036_JPRB,0.036_JPRB,&
1111                    &0.025_JPRB,0.025_JPRB,0.025_JPRB,0.025_JPRB/)
1112                                           
1113!     ----------------------------------------------------------------
1114
1115!*    2.4   OPTICAL PARAMETERS FOR RAIN DROPS
1116!           Savijarvi et al. (1996)
1117
1118ZRTO1 =  0.003_JPRB
1119ZRTO2 = -0.22_JPRB
1120! CAUTION JUST TEMPORARY PARAMETERS     
1121ZROMA2(1:2)= (/ 0.00008_JPRB , 0.264_JPRB  /)
1122ZROMB2(1:2)= (/ 0.23_JPRB    , 0.09_JPRB   /)
1123ZRASY2(1:2)= (/ 0.88_JPRB    , 0.94_JPRB   /)
1124
1125ZRA2(1:2)= (/ 1.5_JPRB     , 1.5_JPRB    /)
1126ZRB2(1:2)= (/ 0.50_JPRB    , 1.13_JPRB   /)
1127ZRC2(1:2)= (/ 5.58E-7_JPRB , 8.55E-4_JPRB/)
1128ZRD2(1:2)= (/ 1.25E-7_JPRB , 1.28E-3_JPRB/)
1129ZRE2(1:2)= (/ 0.841_JPRB   , 0.786_JPRB  /)
1130ZRF2(1:2)= (/ 2.08E-3_JPRB , 5.32E-3_JPRB/)
1131
1132ZROMA4(1:4)= (/ 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB  , 0.465_JPRB   /)
1133ZROMB4(1:4)= (/ 0.23_JPRB    , 0.22_JPRB   , 0.09_JPRB   , 0.001_JPRB   /)
1134ZRASY4(1:4)= (/ 0.88_JPRB    , 0.89_JPRB   , 0.94_JPRB   , 0.97_JPRB    /)
1135
1136ZRA4(1:4)= (/ 1.5_JPRB     , 1.5_JPRB    , 1.5_JPRB    , 1.5_JPRB     /)
1137ZRB4(1:4)= (/ 0.50_JPRB    , 0.78_JPRB   , 1.13_JPRB   , 2.00_JPRB    /)
1138ZRC4(1:4)= (/ 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /)
1139ZRD4(1:4)= (/ 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /)
1140ZRE4(1:4)= (/ 0.841_JPRB   , 0.821_JPRB  , 0.786_JPRB  , 0.820_JPRB   /)
1141ZRF4(1:4)= (/ 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /)
1142
1143ZROMA6(1:6)= (/ 0.00008_JPRB , 0.00008_JPRB , 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB  , 0.465_JPRB   /)
1144ZROMB6(1:6)= (/ 0.23_JPRB    , 0.23_JPRB    , 0.23_JPRB    , 0.22_JPRB   , 0.09_JPRB   , 0.001_JPRB   /)
1145ZRASY6(1:6)= (/ 0.88_JPRB    , 0.88_JPRB    , 0.88_JPRB    , 0.89_JPRB   , 0.94_JPRB   , 0.97_JPRB    /)
1146
1147ZRA6(1:6)= (/ 1.5_JPRB     , 1.5_JPRB     , 1.5_JPRB     , 1.5_JPRB    , 1.5_JPRB    , 1.5_JPRB     /)
1148ZRB6(1:6)= (/ 0.50_JPRB    , 0.50_JPRB    , 0.50_JPRB    , 0.78_JPRB   , 1.13_JPRB   , 2.00_JPRB    /)
1149ZRC6(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 /)
1150ZRD6(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 /)
1151ZRE6(1:6)= (/ 0.841_JPRB   , 0.841_JPRB   , 0.841_JPRB   , 0.821_JPRB  , 0.786_JPRB  , 0.820_JPRB   /)
1152ZRF6(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 /)
1153
1154
1155!=====================================================================
1156
1157!*       2.    SET VALUES.
1158!              -----------
1159
1160
1161IF (LONEWSW) THEN
1162  ZPDH2O = ZPDHNIS
1163  ZPDUMG = ZPDUNIS
1164  ZPRH2O = ZPRHNIS
1165  ZPRUMG = ZPRUNIS
1166  RTDH2O = ZTDHNIS
1167  RTDUMG = ZTDUNIS
1168  RTH2O  = ZTHNIS
1169  RTUMG  = ZTUNIS
1170ENDIF
1171
1172RPDH1=ZPDH2O+_ONE_
1173RPDU1=ZPDUMG+_ONE_
1174ZH2O=_ONE_/( 10._JPRB* RG * RPDH1 )
1175ZUMG=_ONE_/( 10._JPRB* RG * RPDU1 )
1176RPNU = ZUMG/(ZPRUMG**ZPDUMG)
1177RPNH = ZH2O/(ZPRH2O**ZPDH2O)
1178
1179DO JTYPS=1,NTYPS
1180  ZTWEIGHT(JTYPS)=_ZERO_
1181END DO
1182 
1183DO JW=1,NTSW
1184  RSUN(JW)=_ZERO_
1185  RSWCE(JW)=_ZERO_
1186  RSWCP(JW)=_ZERO_
1187 
1188  DO JTYPS=1,NTYPS
1189    RWEIGHT(JW,JTYPS)=_ZERO_
1190  END DO
1191 
1192  DO JC3=1,3
1193    D(JW,JC3)=_ZERO_
1194  ENDDO
1195  DO JC6=1,6
1196    RRAY(JW,JC6)=_ZERO_
1197  ENDDO
1198  DO JI=1,3
1199    DO JJ=1,7
1200      APAD(JW,JI,JJ)=_ZERO_
1201      BPAD(JW,JI,JJ)=_ZERO_
1202    ENDDO
1203  ENDDO
1204ENDDO
1205
1206IF (KSW == 2) THEN
1207  RADJUST=ZADJUS2
1208  RHSRTA=ZRTO1
1209  RHSRTB=ZRTO2
1210  DO JW=1,KSW
1211    RSUN (JW)=ZSUN2(JW)
1212
1213    RSWCE(JW)=ZSWCE2(JW)
1214    RSWCP(JW)=ZSWCP2(JW)
1215
1216    DO JTYPS=1,NTYPS
1217      ZTWEIGHT(JTYPS)=ZTWEIGHT(JTYPS)+ZWEIGHT2(JW,JTYPS)*ZSUN2(JW)
1218    END DO
1219
1220    RROMA(JW)=ZROMA2(JW)
1221    RROMB(JW)=ZROMB2(JW)
1222    RRASY(JW)=ZRASY2(JW)
1223    RHSRA(JW)=ZRA2(JW)
1224    RHSRB(JW)=ZRB2(JW)
1225    RHSRC(JW)=ZRC2(JW)
1226    RHSRD(JW)=ZRD2(JW)
1227    RHSRE(JW)=ZRE2(JW)
1228    RHSRF(JW)=ZRF2(JW)
1229
1230    DO JC3=1,3
1231      D(JW,JC3)=ZD2(JW,JC3)
1232    ENDDO
1233    DO JC6=1,6
1234      RRAY(JW,JC6)=ZRAY2(JW,JC6)
1235    ENDDO
1236    DO JI=1,3
1237      DO JJ=1,7
1238        APAD(JW,JI,JJ)=ZAPAD2(JW,JI,JJ)
1239        BPAD(JW,JI,JJ)=ZBPAD2(JW,JI,JJ)
1240      ENDDO
1241    ENDDO
1242  ENDDO
1243 
1244  DO JW=1,KSW
1245    DO JTYPS=1,NTYPS
1246      RWEIGHT(JW,JTYPS)=ZWEIGHT2(JW,JTYPS)/ZTWEIGHT(JTYPS)
1247    END DO
1248 
1249    DO JM=1,12
1250      IM=MOD(JM+5,12)+1
1251      RALBICE_AR(IM,JW)=ZALBICE2(IM,JW)
1252      RALBICE_AN(JM,JW)=RALBICE_AR(IM,JW)
1253    END DO 
1254  END DO   
1255
1256ELSEIF (KSW == 4) THEN
1257  RADJUST=ZADJUS4
1258  RHSRTA=ZRTO1
1259  RHSRTB=ZRTO2
1260  DO JW=1,KSW
1261    RSUN (JW)=ZSUN4(JW)
1262
1263    RSWCE(JW)=ZSWCE4(JW)
1264    RSWCP(JW)=ZSWCP4(JW)
1265
1266    DO JTYPS=1,NTYPS
1267      ZTWEIGHT(JTYPS)=ZTWEIGHT(JTYPS)+ZWEIGHT4(JW,JTYPS)*ZSUN4(JW)
1268    END DO
1269
1270    RROMA(JW)=ZROMA4(JW)
1271    RROMB(JW)=ZROMB4(JW)
1272    RRASY(JW)=ZRASY4(JW)
1273    RHSRA(JW)=ZRA4(JW)
1274    RHSRB(JW)=ZRB4(JW)
1275    RHSRC(JW)=ZRC4(JW)
1276    RHSRD(JW)=ZRD4(JW)
1277    RHSRE(JW)=ZRE4(JW)
1278    RHSRF(JW)=ZRF4(JW)
1279
1280    DO JC3=1,3
1281      D(JW,JC3)=ZD4(JW,JC3)
1282    ENDDO
1283    DO JC6=1,6
1284      RRAY(JW,JC6)=ZRAY4(JW,JC6)
1285    ENDDO
1286    DO JI=1,3
1287      DO JJ=1,7
1288        APAD(JW,JI,JJ)=ZAPAD4(JW,JI,JJ)
1289        BPAD(JW,JI,JJ)=ZBPAD4(JW,JI,JJ)
1290      ENDDO
1291    ENDDO
1292  ENDDO
1293 
1294  DO JW=1,KSW
1295    DO JTYPS=1,NTYPS
1296      RWEIGHT(JW,JTYPS)=ZWEIGHT4(JW,JTYPS)/ZTWEIGHT(JTYPS)
1297    END DO
1298   
1299    DO JM=1,12
1300      IM=MOD(JM+5,12)+1
1301      RALBICE_AR(IM,JW)=ZALBICE4(IM,JW)
1302      RALBICE_AN(JM,JW)=RALBICE_AR(IM,JW)
1303    END DO 
1304  END DO   
1305 
1306ELSEIF (KSW == 6) THEN
1307  RADJUST=ZADJUS6
1308  RHSRTA=ZRTO1
1309  RHSRTB=ZRTO2
1310  DO JW=1,KSW
1311    RSUN (JW)=ZSUN6(JW)
1312
1313    RSWCE(JW)=ZSWCE6(JW)
1314    RSWCP(JW)=ZSWCP6(JW)
1315
1316    DO JTYPS=1,NTYPS
1317      ZTWEIGHT(JTYPS)=ZTWEIGHT(JTYPS)+ZWEIGHT6(JW,JTYPS)*ZSUN6(JW)
1318    END DO
1319
1320    RROMA(JW)=ZROMA6(JW)
1321    RROMB(JW)=ZROMB6(JW)
1322    RRASY(JW)=ZRASY6(JW)
1323    RHSRA(JW)=ZRA6(JW)
1324    RHSRB(JW)=ZRB6(JW)
1325    RHSRC(JW)=ZRC6(JW)
1326    RHSRD(JW)=ZRD6(JW)
1327    RHSRE(JW)=ZRE6(JW)
1328    RHSRF(JW)=ZRF6(JW)
1329
1330    DO JC3=1,3
1331      D(JW,JC3)=ZD6(JW,JC3)
1332    ENDDO
1333    DO JC6=1,6
1334      RRAY(JW,JC6)=ZRAY6(JW,JC6)
1335    ENDDO
1336    DO JI=1,3
1337      DO JJ=1,7
1338        APAD(JW,JI,JJ)=ZAPAD6(JW,JI,JJ)
1339        BPAD(JW,JI,JJ)=ZBPAD6(JW,JI,JJ)
1340      ENDDO
1341    ENDDO
1342  ENDDO
1343 
1344  DO JW=1,KSW
1345    DO JTYPS=1,NTYPS
1346      RWEIGHT(JW,JTYPS)=ZWEIGHT6(JW,JTYPS)/ZTWEIGHT(JTYPS)
1347    END DO
1348   
1349    DO JM=1,12
1350      IM=MOD(JM+5,12)+1
1351      RALBICE_AR(IM,JW)=ZALBICE6(IM,JW)
1352      RALBICE_AN(JM,JW)=RALBICE_AR(IM,JW)
1353    END DO 
1354  END DO   
1355 
1356ELSE
1357!  CALL ABOR1('SUSWN: WRONG NUMBER OF SW INTERVALS')
1358  STOP 'SUSWN: WRONG NUMBER OF SW INTERVALS'
1359ENDIF
1360
1361!     ----------------------------------------------------------------
1362
1363!*       7.    PRINT VALUES.
1364!              ------------
1365
1366
1367!     PRINT 9700,RPDH1,RPDU1,RPNU,RPNH
1368!9700 FORMAT(1X,'RPDH1,RPDU1,RPNU,RPNH ',4F10.5)     
1369!     DO JW=1,KSW     
1370!       PRINT *,' SW INTERVAL ',JW
1371!       PRINT 9701,JW,RSUN(JW),RSWCE(JW),RSWCP(JW)
1372!9701   FORMAT(1X,I3,F10.7,2F10.5)
1373!       PRINT 9702,JW,(D(JW,JC3),JC3=1,3)
1374!9702   FORMAT(1X,I3,3F10.5)
1375!       PRINT 9703,JW,(RRAY(JW,JC6),JC6=1,6)
1376!9703   FORMAT(1X,I3,6E15.8)
1377!       PRINT 9704,JW,(JABS,(APAD(JW,JABS,JJ),JJ=1,7),JABS=1,3)
1378!9704   FORMAT(1X,I3,(I3,7E15.8))
1379!       PRINT 9706,JW,RROMA(JW),RROMB(JW),RRASY(JW),RHSRTA,RHSRTB
1380!9706   FORMAT(1X,I3,' Rain Optics from Rain Rate: ',F8.5,F6.3,3F6.3)
1381!       PRINT 9707,JW,RHSRA(JW),RHSRB(JW),RHSRC(JW),RHSRD(JW)
1382!    &   ,RHSRE(JW),RHSRF(JW) 
1383!9707   FORMAT(1X,I3,' Rain Optics as function of Reff: ',6E11.4)
1384!     ENDDO
1385
1386!     ----------------------------------------------------------------
1387
1388RETURN
1389END SUBROUTINE SUSWN
Note: See TracBrowser for help on using the repository browser.