source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/su_aerop.F90 @ 5313

Last change on this file since 5313 was 3331, checked in by acozic, 6 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 48.4 KB
Line 
1SUBROUTINE SU_AEROP
2
3!**** *SU_AEROP*   - INITIALIZE MODULES YOEAEROP
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOEAEROP THE MODULE THAT CONTAINS THE
8!           OPTICAL PROPERTIES FOR THE PROGNOSTIC AEROSOLS
9
10!**   INTERFACE.
11!     ----------
12!        *CALL* *SU_AEROP
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        YOEAEROP
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!        from O.BOUCHER's stand-alone codes (LOA, 2005-03)
37
38!     MODIFICATIONS.
39!     --------------
40!        ORIGINAL : 2005-03-25
41
42!     ------------------------------------------------------------------
43
44USE PARKIND1  ,ONLY : JPRB
45USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
46
47USE YOEAEROP  , ONLY : ALF_SU, ASY_SU, OMG_SU, ALF_OM, ASY_OM, OMG_OM, &
48  & ALF_DD, ASY_DD, OMG_DD, ALF_SS, ASY_SS, OMG_SS, ALF_BC, ASY_BC, OMG_BC, &
49  & ALF_FA, ASY_FA, OMG_FA
50
51IMPLICIT NONE
52
53REAL(KIND=JPRB) :: ZHOOK_HANDLE
54!     ----------------------------------------------------------------
55IF (LHOOK) CALL DR_HOOK('SU_AEROP',0,ZHOOK_HANDLE)
56
57
58!*      1.    PARAMETERS RELATED TO OPTICAL PROPERTIES
59!             ----------------------------------------
60
61! All arrays are given in three blocks:
62! first, for 16 wavelengths relevant to observations
63! 340, 380, 400, 440, 469, 500, 532, 555, 645, 670, 800, 858, 865, 1240, 1640, 2130
64
65!(  not done yet, depends whether SW6 or RRTM_SW is used, to be computed for RRTM_LW)
66! second, for the spectral intervals of the SW radiation scheme
67! third, for the spectral intervals of the LW radiation scheme
68
69! For sulfates (_SU) and organic matters (_OM), values are given as function
70! of RH (12 values) for the 8 wavelengths
71! For sea-salt (_SS), 3 bins are considered      (0.03, 0.50, 5.0, 20.)
72! For desert dust (_DD), 3 bins are considered   (0.03, 0.55, 0.9, 20.)
73!
74! IF BIN LIMITS ARE CHANGED, MAKE SURE THAT THE RELEVANT SEDIMENTATION SPEEDS ARE
75! RECOMPUTED ACCORDINGLY
76
77! RH  0.         10.         20.         30.         40.         50.
78!    60.         70.         80.         85.         90.         95.   
79
80ALF_FA(:) =0.0_JPRB
81ASY_FA(:) =0.0_JPRB
82OMG_FA(:) =0.0_JPRB
83
84! alpha_aer_bc
85ALF_BC( :) = (/ &
86  & 14.237_JPRB, 14.237_JPRB, &
87  & 14.237_JPRB, 12.619_JPRB, 12.000_JPRB, 10.684_JPRB,  9.840_JPRB,  9.412_JPRB,  9.000_JPRB, &
88  &  7.212_JPRB,  5.677_JPRB,  5.500_JPRB,  5.113_JPRB,  5.000_JPRB,  5.000_JPRB,  5.000_JPRB /)
89
90! cg_aer_bc
91ASY_BC( :) = (/ &
92  &   .398_JPRB,   .398_JPRB, &
93  &   .398_JPRB,   .379_JPRB,   .353_JPRB,   .350_JPRB,   .341_JPRB,   .335_JPRB,   .300_JPRB, &
94  &   .296_JPRB,   .263_JPRB,   .250_JPRB,   .248_JPRB,   .200_JPRB,   .200_JPRB,   .200_JPRB /)
95
96! om_aer_bc
97OMG_BC( :) = (/ &
98  &   .270_JPRB,   .270_JPRB, &
99  &   .270_JPRB,   .251_JPRB,   .240_JPRB,   .226_JPRB,   .213_JPRB,   .206_JPRB,   .180_JPRB, &
100  &   .166_JPRB,   .132_JPRB,   .120_JPRB,   .118_JPRB,   .100_JPRB,   .100_JPRB,   .100_JPRB  /)
101
102! alpha_aer_sulfate
103ALF_SU( :,1) = (/ &
104  &  6.823_JPRB, 6.823_JPRB, 6.823_JPRB, 6.823_JPRB, 9.168_JPRB,10.246_JPRB, &
105  & 11.685_JPRB,13.796_JPRB,17.504_JPRB,20.617_JPRB,26.335_JPRB,42.574_JPRB /)
106ALF_SU( :,2) = (/ &
107  &  6.823_JPRB, 6.823_JPRB, 6.823_JPRB, 6.823_JPRB, 9.168_JPRB,10.246_JPRB, &
108  & 11.685_JPRB,13.796_JPRB,17.504_JPRB,20.617_JPRB,26.335_JPRB,42.574_JPRB /)
109ALF_SU( :,3) = (/ &
110  &  6.823_JPRB, 6.823_JPRB, 6.823_JPRB, 6.823_JPRB, 9.168_JPRB,10.246_JPRB, &
111  & 11.685_JPRB,13.796_JPRB,17.504_JPRB,20.617_JPRB,26.335_JPRB,42.574_JPRB /)
112ALF_SU( :,4) = (/ &
113  &  6.023_JPRB, 6.023_JPRB, 6.023_JPRB, 6.023_JPRB, 8.115_JPRB, 9.094_JPRB, &
114  & 10.407_JPRB,12.344_JPRB,15.773_JPRB,18.661_JPRB,24.014_JPRB,39.387_JPRB /)
115ALF_SU( :,5) = (/ &
116  &  6.023_JPRB, 6.023_JPRB, 6.023_JPRB, 6.023_JPRB, 8.115_JPRB, 9.094_JPRB, &
117  & 10.407_JPRB,12.344_JPRB,15.773_JPRB,18.661_JPRB,24.014_JPRB,39.387_JPRB /)
118
119ALF_SU( :,6) = (/ &
120  &  5.009_JPRB, 5.009_JPRB, 5.009_JPRB, 5.009_JPRB, 6.779_JPRB, 7.626_JPRB, &
121  &  8.769_JPRB,10.466_JPRB,13.501_JPRB,16.072_JPRB,20.889_JPRB,34.944_JPRB /)
122ALF_SU( :,7) = (/ &
123  &  4.548_JPRB, 4.548_JPRB, 4.548_JPRB, 4.548_JPRB, 6.171_JPRB, 6.956_JPRB, &
124  &  8.016_JPRB, 9.597_JPRB,12.438_JPRB,14.852_JPRB,19.397_JPRB,32.764_JPRB /)
125
126ALF_SU( :,8) = (/ &
127  &  4.311_JPRB, 4.311_JPRB, 4.311_JPRB, 4.311_JPRB, 5.858_JPRB, 6.609_JPRB, &
128  &  7.626_JPRB, 9.145_JPRB,11.882_JPRB,14.211_JPRB,18.609_JPRB,31.596_JPRB /)
129ALF_SU( :,9) = (/ &
130  &  4.311_JPRB, 4.311_JPRB, 4.311_JPRB, 4.311_JPRB, 5.858_JPRB, 6.609_JPRB, &
131  &  7.626_JPRB, 9.145_JPRB,11.882_JPRB,14.211_JPRB,18.609_JPRB,31.596_JPRB /)
132
133ALF_SU(:,10) = (/ &
134  &  3.053_JPRB, 3.053_JPRB, 3.053_JPRB, 3.053_JPRB, 4.191_JPRB, 4.758_JPRB, &
135  &  5.532_JPRB, 6.700_JPRB, 8.836_JPRB,10.671_JPRB,14.193_JPRB,24.859_JPRB /)
136ALF_SU(:,11) = (/ &
137  &  2.155_JPRB, 2.155_JPRB, 2.155_JPRB, 2.155_JPRB, 2.992_JPRB, 3.416_JPRB, &
138  &  3.999_JPRB, 4.886_JPRB, 6.532_JPRB, 7.958_JPRB,10.735_JPRB,19.336_JPRB /)
139ALF_SU(:,12) = (/ &
140  &  2.155_JPRB, 2.155_JPRB, 2.155_JPRB, 2.155_JPRB, 2.992_JPRB, 3.416_JPRB, &
141  &  3.999_JPRB, 4.886_JPRB, 6.532_JPRB, 7.958_JPRB,10.735_JPRB,19.336_JPRB /)
142
143ALF_SU(:,13) = (/ &
144  &  1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 2.551_JPRB, 2.920_JPRB, &
145  &  3.429_JPRB, 4.207_JPRB, 5.657_JPRB, 6.920_JPRB, 9.392_JPRB,17.128_JPRB /)
146ALF_SU(:,14) = (/ &
147  &  1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 2.551_JPRB, 2.920_JPRB, &
148  &  3.429_JPRB, 4.207_JPRB, 5.657_JPRB, 6.920_JPRB, 9.392_JPRB,17.128_JPRB /)
149ALF_SU(:,15) = (/ &
150  &  1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 2.551_JPRB, 2.920_JPRB, &
151  &  3.429_JPRB, 4.207_JPRB, 5.657_JPRB, 6.920_JPRB, 9.392_JPRB,17.128_JPRB /)
152ALF_SU(:,16) = (/ &
153  &  1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 1.827_JPRB, 2.551_JPRB, 2.920_JPRB, &
154  &  3.429_JPRB, 4.207_JPRB, 5.657_JPRB, 6.920_JPRB, 9.392_JPRB,17.128_JPRB /)
155
156! cg_aer_sulfate
157ASY_SU( :,1) = (/ &
158  &  0.640_JPRB,  .640_JPRB,  .640_JPRB,  .640_JPRB,  .695_JPRB,  .704_JPRB, &
159  &   .713_JPRB,  .723_JPRB,  .734_JPRB,  .744_JPRB,  .754_JPRB,  .775_JPRB /)
160ASY_SU( :,2) = (/ &
161  &  0.640_JPRB,  .640_JPRB,  .640_JPRB,  .640_JPRB,  .695_JPRB,  .704_JPRB, &
162  &   .713_JPRB,  .723_JPRB,  .734_JPRB,  .744_JPRB,  .754_JPRB,  .775_JPRB /)
163ASY_SU( :,3) = (/ &
164  &  0.640_JPRB,  .640_JPRB,  .640_JPRB,  .640_JPRB,  .695_JPRB,  .704_JPRB, &
165  &   .713_JPRB,  .723_JPRB,  .734_JPRB,  .744_JPRB,  .754_JPRB,  .775_JPRB /)
166ASY_SU( :,4) = (/ &
167  &   .632_JPRB,  .632_JPRB,  .632_JPRB,  .632_JPRB,  .687_JPRB,  .696_JPRB, &
168  &   .705_JPRB,  .716_JPRB,  .728_JPRB,  .738_JPRB,  .749_JPRB,  .771_JPRB /)
169ASY_SU( :,5) = (/ &
170  &   .632_JPRB,  .632_JPRB,  .632_JPRB,  .632_JPRB,  .687_JPRB,  .696_JPRB, &
171  &   .705_JPRB,  .716_JPRB,  .728_JPRB,  .738_JPRB,  .749_JPRB,  .771_JPRB /)
172
173ASY_SU( :,6) = (/ &
174  &   .620_JPRB,  .620_JPRB,  .620_JPRB,  .620_JPRB,  .674_JPRB,  .683_JPRB, &
175  &   .693_JPRB,  .704_JPRB,  .717_JPRB,  .728_JPRB,  .740_JPRB,  .764_JPRB /)
176ASY_SU( :,7) = (/ &
177  &   .613_JPRB,  .613_JPRB,  .613_JPRB,  .613_JPRB,  .667_JPRB,  .677_JPRB, &
178  &   .687_JPRB,  .698_JPRB,  .711_JPRB,  .722_JPRB,  .735_JPRB,  .760_JPRB /)
179
180ASY_SU( :,8) = (/ &
181  &   .609_JPRB,  .609_JPRB,  .609_JPRB,  .609_JPRB,  .664_JPRB,  .673_JPRB, &
182  &   .683_JPRB,  .694_JPRB,  .708_JPRB,  .719_JPRB,  .733_JPRB,  .758_JPRB /)
183ASY_SU( :,9) = (/ &
184  &   .609_JPRB,  .609_JPRB,  .609_JPRB,  .609_JPRB,  .664_JPRB,  .673_JPRB, &
185  &   .683_JPRB,  .694_JPRB,  .708_JPRB,  .719_JPRB,  .733_JPRB,  .758_JPRB /)
186
187ASY_SU(:,10) = (/ &
188  &   .584_JPRB,  .584_JPRB,  .584_JPRB,  .584_JPRB,  .637_JPRB,  .647_JPRB, &
189  &   .658_JPRB,  .670_JPRB,  .686_JPRB,  .698_JPRB,  .713_JPRB,  .742_JPRB /)
190ASY_SU(:,11) = (/ &
191  &   .558_JPRB,  .558_JPRB,  .558_JPRB,  .558_JPRB,  .609_JPRB,  .620_JPRB, &
192  &   .631_JPRB,  .644_JPRB,  .661_JPRB,  .674_JPRB,  .691_JPRB,  .723_JPRB /)
193ASY_SU(:,12) = (/ &
194  &   .558_JPRB,  .558_JPRB,  .558_JPRB,  .558_JPRB,  .609_JPRB,  .620_JPRB, &
195  &   .631_JPRB,  .644_JPRB,  .661_JPRB,  .674_JPRB,  .691_JPRB,  .723_JPRB /)
196
197ASY_SU(:,13) = (/ &
198  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .596_JPRB,  .606_JPRB, &
199  &   .618_JPRB,  .632_JPRB,  .649_JPRB,  .663_JPRB,  .680_JPRB,  .714_JPRB /)
200ASY_SU(:,14) = (/ &
201  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .596_JPRB,  .606_JPRB, &
202  &   .618_JPRB,  .632_JPRB,  .649_JPRB,  .663_JPRB,  .680_JPRB,  .714_JPRB /)
203ASY_SU(:,15) = (/ &
204  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .596_JPRB,  .606_JPRB, &
205  &   .618_JPRB,  .632_JPRB,  .649_JPRB,  .663_JPRB,  .680_JPRB,  .714_JPRB /)
206ASY_SU(:,16) = (/ &
207  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .596_JPRB,  .606_JPRB, &
208  &   .618_JPRB,  .632_JPRB,  .649_JPRB,  .663_JPRB,  .680_JPRB,  .714_JPRB /)
209
210! om_aer_sulfate
211OMG_SU( :,1) = (/ &
212  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
213  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
214OMG_SU( :,2) = (/ &
215  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
216  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
217OMG_SU( :,3) = (/ &
218  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
219  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
220OMG_SU( :,4) = (/ &
221  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
222  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
223OMG_SU( :,5) = (/ &
224  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
225  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
226OMG_SU( :,6) = (/ &
227  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
228  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
229OMG_SU( :,7) = (/ &
230  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
231  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /) 
232
233OMG_SU( :,8) = (/ &
234  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
235  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
236OMG_SU( :,9) = (/ &
237  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
238  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
239OMG_SU(:,10) = (/ &
240  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, & 
241  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
242OMG_SU(:,11) = (/ &
243  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
244  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
245OMG_SU(:,12) = (/ &
246  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
247  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
248OMG_SU(:,13) = (/ &
249  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
250  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
251OMG_SU(:,14) = (/ &
252  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
253  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
254OMG_SU(:,15) = (/ &
255  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
256  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
257OMG_SU(:,16) = (/ &
258  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, &
259  &  1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
260
261
262! alpha_aer_om (organic matter)
263ALF_OM( :,1) = (/ &
264  &  4.973_JPRB, 4.973_JPRB, 4.973_JPRB, 4.973_JPRB, 7.408_JPRB, 8.294_JPRB, &
265  &  9.468_JPRB,11.168_JPRB,14.087_JPRB,16.688_JPRB,21.321_JPRB,24.076_JPRB /)
266ALF_OM( :,2) = (/ &
267  &  4.973_JPRB, 4.973_JPRB, 4.973_JPRB, 4.973_JPRB, 7.408_JPRB, 8.294_JPRB, &
268  &  9.468_JPRB,11.168_JPRB,14.087_JPRB,16.688_JPRB,21.321_JPRB,24.076_JPRB /)
269ALF_OM( :,3) = (/ &
270  &  4.973_JPRB, 4.973_JPRB, 4.973_JPRB, 4.973_JPRB, 7.408_JPRB, 8.294_JPRB, &
271  &  9.468_JPRB,11.168_JPRB,14.087_JPRB,16.688_JPRB,21.321_JPRB,24.076_JPRB /)
272ALF_OM( :,4) = (/ &
273  &  4.396_JPRB, 4.396_JPRB, 4.396_JPRB, 4.396_JPRB, 6.612_JPRB, 7.425_JPRB, &
274  &  8.508_JPRB,10.084_JPRB,12.810_JPRB,15.257_JPRB,19.652_JPRB,22.764_JPRB /)
275ALF_OM( :,5) = (/ &
276  &  4.396_JPRB, 4.396_JPRB, 4.396_JPRB, 4.396_JPRB, 6.612_JPRB, 7.425_JPRB, &
277  &  8.508_JPRB,10.084_JPRB,12.810_JPRB,15.257_JPRB,19.652_JPRB,22.764_JPRB /)
278
279ALF_OM( :,6) = (/ &
280  &  3.664_JPRB, 3.664_JPRB, 3.664_JPRB, 3.664_JPRB, 5.585_JPRB, 6.298_JPRB, &
281  &  7.253_JPRB, 8.652_JPRB,11.097_JPRB,13.314_JPRB,17.340_JPRB,20.901_JPRB /)
282ALF_OM( :,7) = (/ &
283  &  3.331_JPRB, 3.331_JPRB, 3.331_JPRB, 3.331_JPRB, 5.111_JPRB, 5.775_JPRB, &
284  &  6.667_JPRB, 7.979_JPRB,10.281_JPRB,12.380_JPRB,16.211_JPRB,19.964_JPRB /)
285
286ALF_OM( :,8) = (/ &
287  &  3.159_JPRB, 3.159_JPRB, 3.159_JPRB, 3.159_JPRB, 4.865_JPRB, 5.502_JPRB, &
288  &  6.361_JPRB, 7.625_JPRB, 9.851_JPRB,11.885_JPRB,15.608_JPRB,19.454_JPRB /)
289ALF_OM( :,9) = (/ &
290  &  3.159_JPRB, 3.159_JPRB, 3.159_JPRB, 3.159_JPRB, 4.865_JPRB, 5.502_JPRB, &
291  &  6.361_JPRB, 7.625_JPRB, 9.851_JPRB,11.885_JPRB,15.608_JPRB,19.454_JPRB /)
292
293ALF_OM(:,10) = (/ &
294  &  2.248_JPRB, 2.248_JPRB, 2.248_JPRB, 2.248_JPRB, 3.537_JPRB, 4.027_JPRB, &
295  &  4.692_JPRB, 5.682_JPRB, 7.451_JPRB, 9.094_JPRB,12.150_JPRB,16.386_JPRB /)
296ALF_OM(:,11) = (/ &
297  &  1.595_JPRB, 1.595_JPRB, 1.595_JPRB, 1.595_JPRB, 2.559_JPRB, 2.931_JPRB, &
298  &  3.439_JPRB, 4.203_JPRB, 5.587_JPRB, 6.888_JPRB, 9.348_JPRB,13.658_JPRB /)
299ALF_OM(:,12) = (/ &
300  &  1.595_JPRB, 1.595_JPRB, 1.595_JPRB, 1.595_JPRB, 2.559_JPRB, 2.931_JPRB, &
301  &  3.439_JPRB, 4.203_JPRB, 5.587_JPRB, 6.888_JPRB, 9.348_JPRB,13.658_JPRB /)
302
303ALF_OM(:,13) = (/ &
304  &  1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 2.195_JPRB, 2.520_JPRB, &
305  &  2.967_JPRB, 3.641_JPRB, 4.867_JPRB, 6.028_JPRB, 8.236_JPRB,12.495_JPRB /)
306ALF_OM(:,14) = (/ &
307  &  1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 2.195_JPRB, 2.520_JPRB, &
308  &  2.967_JPRB, 3.641_JPRB, 4.867_JPRB, 6.028_JPRB, 8.236_JPRB,12.495_JPRB /)
309ALF_OM(:,15) = (/ &
310  &  1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 2.195_JPRB, 2.520_JPRB, &
311  &  2.967_JPRB, 3.641_JPRB, 4.867_JPRB, 6.028_JPRB, 8.236_JPRB,12.495_JPRB /)
312ALF_OM(:,16) = (/ &
313  &  1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 1.356_JPRB, 2.195_JPRB, 2.520_JPRB, &
314  &  2.967_JPRB, 3.641_JPRB, 4.867_JPRB, 6.028_JPRB, 8.236_JPRB,12.495_JPRB /)
315
316ASY_OM( :,1) = (/ &
317  &   .644_JPRB,  .644_JPRB,  .644_JPRB,  .644_JPRB,  .676_JPRB,  .684_JPRB, &
318  &   .692_JPRB,  .700_JPRB,  .711_JPRB,  .718_JPRB,  .726_JPRB,  .543_JPRB /)
319ASY_OM( :,2) = (/ &
320  &   .644_JPRB,  .644_JPRB,  .644_JPRB,  .644_JPRB,  .676_JPRB,  .684_JPRB, &
321  &   .692_JPRB,  .700_JPRB,  .711_JPRB,  .718_JPRB,  .726_JPRB,  .543_JPRB /)
322ASY_OM( :,3) = (/ &
323  &   .644_JPRB,  .644_JPRB,  .644_JPRB,  .644_JPRB,  .676_JPRB,  .684_JPRB, &
324  &   .692_JPRB,  .700_JPRB,  .711_JPRB,  .718_JPRB,  .726_JPRB,  .543_JPRB /)
325ASY_OM( :,4) = (/ &
326  &   .635_JPRB,  .635_JPRB,  .635_JPRB,  .635_JPRB,  .669_JPRB,  .676_JPRB, &
327  &   .685_JPRB,  .694_JPRB,  .705_JPRB,  .713_JPRB,  .721_JPRB,  .538_JPRB /)
328ASY_OM( :,5) = (/ &
329  &   .635_JPRB,  .635_JPRB,  .635_JPRB,  .635_JPRB,  .669_JPRB,  .676_JPRB, &
330  &   .685_JPRB,  .694_JPRB,  .705_JPRB,  .713_JPRB,  .721_JPRB,  .538_JPRB /)
331
332ASY_OM( :,6) = (/ &
333  &   .623_JPRB,  .623_JPRB,  .623_JPRB,  .623_JPRB,  .657_JPRB,  .665_JPRB, &
334  &   .674_JPRB,  .684_JPRB,  .696_JPRB,  .704_JPRB,  .714_JPRB,  .529_JPRB /)
335ASY_OM( :,7) = (/ &
336  &   .616_JPRB,  .616_JPRB,  .616_JPRB,  .616_JPRB,  .651_JPRB,  .659_JPRB, &
337  &   .668_JPRB,  .678_JPRB,  .691_JPRB,  .699_JPRB,  .710_JPRB,  .524_JPRB /)
338
339ASY_OM( :,8) = (/ &
340  &   .612_JPRB,  .612_JPRB,  .612_JPRB,  .612_JPRB,  .647_JPRB,  .655_JPRB, &
341  &   .664_JPRB,  .675_JPRB,  .688_JPRB,  .696_JPRB,  .707_JPRB,  .521_JPRB /)
342ASY_OM( :,9) = (/ &
343  &   .612_JPRB,  .612_JPRB,  .612_JPRB,  .612_JPRB,  .647_JPRB,  .655_JPRB, &
344  &   .664_JPRB,  .675_JPRB,  .688_JPRB,  .696_JPRB,  .707_JPRB,  .521_JPRB /)
345
346ASY_OM(:,10) = (/ &
347  &   .586_JPRB,  .586_JPRB,  .586_JPRB,  .586_JPRB,  .623_JPRB,  .632_JPRB, &
348  &   .641_JPRB,  .653_JPRB,  .667_JPRB,  .677_JPRB,  .690_JPRB,  .502_JPRB /)
349ASY_OM(:,11) = (/ &
350  &   .559_JPRB,  .559_JPRB,  .559_JPRB,  .559_JPRB,  .597_JPRB,  .606_JPRB, &
351  &   .617_JPRB,  .629_JPRB,  .645_JPRB,  .656_JPRB,  .670_JPRB,  .482_JPRB /)
352ASY_OM(:,12) = (/ &
353  &   .559_JPRB,  .559_JPRB,  .559_JPRB,  .559_JPRB,  .597_JPRB,  .606_JPRB, &
354  &   .617_JPRB,  .629_JPRB,  .645_JPRB,  .656_JPRB,  .670_JPRB,  .482_JPRB /)
355
356ASY_OM(:,13) = (/ &
357  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .584_JPRB,  .594_JPRB, &
358  &   .605_JPRB,  .617_JPRB,  .634_JPRB,  .645_JPRB,  .660_JPRB,  .472_JPRB /)
359ASY_OM(:,14) = (/ &
360  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .584_JPRB,  .594_JPRB, &
361  &   .605_JPRB,  .617_JPRB,  .634_JPRB,  .645_JPRB,  .660_JPRB,  .472_JPRB /)
362ASY_OM(:,15) = (/ &
363  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .584_JPRB,  .594_JPRB, &
364  &   .605_JPRB,  .617_JPRB,  .634_JPRB,  .645_JPRB,  .660_JPRB,  .472_JPRB /)
365ASY_OM(:,16) = (/ &
366  &   .546_JPRB,  .546_JPRB,  .546_JPRB,  .546_JPRB,  .584_JPRB,  .594_JPRB, &
367  &   .605_JPRB,  .617_JPRB,  .634_JPRB,  .645_JPRB,  .660_JPRB,  .472_JPRB /)
368
369OMG_OM( :,1) = (/ &
370  &   .971_JPRB,  .971_JPRB,  .971_JPRB,  .971_JPRB,  .980_JPRB,  .982_JPRB, &
371  &   .984_JPRB,  .987_JPRB,  .989_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
372OMG_OM( :,2) = (/ &
373  &   .971_JPRB,  .971_JPRB,  .971_JPRB,  .971_JPRB,  .980_JPRB,  .982_JPRB, &
374  &   .984_JPRB,  .987_JPRB,  .989_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
375OMG_OM( :,3) = (/ &
376  &   .971_JPRB,  .971_JPRB,  .971_JPRB,  .971_JPRB,  .980_JPRB,  .982_JPRB, &
377  &   .984_JPRB,  .987_JPRB,  .989_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
378OMG_OM( :,4) = (/ &
379  &   .971_JPRB,  .971_JPRB,  .971_JPRB,  .971_JPRB,  .980_JPRB,  .982_JPRB, &
380  &   .985_JPRB,  .987_JPRB,  .990_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
381OMG_OM( :,5) = (/ &
382  &   .971_JPRB,  .971_JPRB,  .971_JPRB,  .971_JPRB,  .980_JPRB,  .982_JPRB, &
383  &   .985_JPRB,  .987_JPRB,  .990_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
384
385OMG_OM( :,6) = (/ &
386  &   .970_JPRB,  .970_JPRB,  .970_JPRB,  .970_JPRB,  .980_JPRB,  .982_JPRB, &
387  &   .984_JPRB,  .987_JPRB,  .990_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
388OMG_OM( :,7) = (/ &
389  &   .970_JPRB,  .970_JPRB,  .970_JPRB,  .970_JPRB,  .980_JPRB,  .982_JPRB, &
390  &   .984_JPRB,  .987_JPRB,  .990_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
391
392OMG_OM( :,8) = (/ &
393  &   .969_JPRB,  .969_JPRB,  .969_JPRB,  .969_JPRB,  .980_JPRB,  .982_JPRB, &
394  &   .984_JPRB,  .987_JPRB,  .990_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
395OMG_OM( :,9) = (/ &
396  &   .969_JPRB,  .969_JPRB,  .969_JPRB,  .969_JPRB,  .980_JPRB,  .982_JPRB, &
397  &   .984_JPRB,  .987_JPRB,  .990_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
398
399OMG_OM(:,10) = (/ &
400  &   .967_JPRB,  .967_JPRB,  .967_JPRB,  .967_JPRB,  .979_JPRB,  .981_JPRB, &
401  &   .984_JPRB,  .986_JPRB,  .989_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
402OMG_OM(:,11) = (/ &
403  &   .964_JPRB,  .964_JPRB,  .964_JPRB,  .964_JPRB,  .977_JPRB,  .979_JPRB, &
404  &   .982_JPRB,  .985_JPRB,  .989_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
405OMG_OM(:,12) = (/ &
406  &   .964_JPRB,  .964_JPRB,  .964_JPRB,  .964_JPRB,  .977_JPRB,  .979_JPRB, &
407  &   .982_JPRB,  .985_JPRB,  .989_JPRB,  .991_JPRB,  .993_JPRB, 1.000_JPRB /)
408
409OMG_OM(:,13) = (/ &
410  &   .962_JPRB,  .962_JPRB,  .962_JPRB,  .962_JPRB,  .976_JPRB,  .979_JPRB, &
411  &   .982_JPRB,  .985_JPRB,  .988_JPRB,  .990_JPRB,  .993_JPRB,  .999_JPRB /)
412OMG_OM(:,14) = (/ &
413  &   .962_JPRB,  .962_JPRB,  .962_JPRB,  .962_JPRB,  .976_JPRB,  .979_JPRB, &
414  &   .982_JPRB,  .985_JPRB,  .988_JPRB,  .990_JPRB,  .993_JPRB,  .999_JPRB /)
415OMG_OM(:,15) = (/ &
416  &   .962_JPRB,  .962_JPRB,  .962_JPRB,  .962_JPRB,  .976_JPRB,  .979_JPRB, &
417  &   .982_JPRB,  .985_JPRB,  .988_JPRB,  .990_JPRB,  .993_JPRB,  .999_JPRB /)
418OMG_OM(:,16) = (/ &
419  &   .962_JPRB,  .962_JPRB,  .962_JPRB,  .962_JPRB,  .976_JPRB,  .979_JPRB, &
420  &   .982_JPRB,  .985_JPRB,  .988_JPRB,  .990_JPRB,  .993_JPRB,  .999_JPRB /)
421
422
423!alpha_aer_dd (desert dust)
424!                    bin1,       bin2,       bin3
425ALF_DD( :,1) = (/ 2.210_JPRB, 1.011_JPRB,  .398_JPRB /)
426ALF_DD( :,2) = (/ 2.210_JPRB, 1.011_JPRB,  .398_JPRB /)
427ALF_DD( :,3) = (/ 2.210_JPRB, 1.011_JPRB,  .398_JPRB /)
428ALF_DD( :,4) = (/ 2.357_JPRB,  .985_JPRB,  .398_JPRB /)
429ALF_DD( :,5) = (/ 2.453_JPRB, 0.949_JPRB, 0.404_JPRB /)
430ALF_DD( :,6) = (/ 2.549_JPRB,  .901_JPRB,  .414_JPRB /)
431ALF_DD( :,7) = (/ 2.605_JPRB,  .877_JPRB,  .423_JPRB /)
432ALF_DD( :,8) = (/ 2.635_JPRB, 0.869_JPRB, 0.428_JPRB /)
433ALF_DD( :,9) = (/ 2.596_JPRB, 0.978_JPRB, 0.428_JPRB /)
434ALF_DD(:,10) = (/ 2.555_JPRB, 1.039_JPRB,  .424_JPRB /)
435ALF_DD(:,11) = (/ 2.227_JPRB, 1.360_JPRB,  .404_JPRB /)
436ALF_DD(:,12) = (/ 2.058_JPRB, 1.460_JPRB, 0.405_JPRB /)
437ALF_DD(:,13) = (/ 2.041_JPRB, 1.474_JPRB,  .405_JPRB /)
438ALF_DD(:,14) = (/ 1.117_JPRB, 1.523_JPRB, 0.502_JPRB /)
439ALF_DD(:,15) = (/ 0.583_JPRB, 1.117_JPRB, 0.596_JPRB /)
440ALF_DD(:,16) = (/ 0.263_JPRB, 0.705_JPRB, 0.607_JPRB /)
441
442ASY_DD( :,1) = (/  .673_JPRB,  .745_JPRB,  .799_JPRB /)
443ASY_DD( :,2) = (/  .673_JPRB,  .745_JPRB,  .799_JPRB /)
444ASY_DD( :,3) = (/  .673_JPRB,  .745_JPRB,  .799_JPRB /)
445ASY_DD( :,4) = (/  .691_JPRB,  .712_JPRB,  .792_JPRB /)
446ASY_DD( :,5) = (/ 0.704_JPRB, 0.679_JPRB, 0.791_JPRB /)
447ASY_DD( :,6) = (/  .715_JPRB,  .640_JPRB,  .786_JPRB /)
448ASY_DD( :,7) = (/  .726_JPRB,  .605_JPRB,  .782_JPRB /)
449ASY_DD( :,8) = (/ 0.731_JPRB, 0.591_JPRB, 0.782_JPRB /)
450ASY_DD( :,9) = (/ 0.741_JPRB, 0.606_JPRB, 0.756_JPRB /)
451ASY_DD(:,10) = (/  .740_JPRB,  .617_JPRB,  .738_JPRB /)
452ASY_DD(:,11) = (/  .732_JPRB,  .698_JPRB,  .687_JPRB /)
453ASY_DD(:,12) = (/ 0.725_JPRB, 0.722_JPRB, 0.679_JPRB /)
454ASY_DD(:,13) = (/  .723_JPRB,  .722_JPRB,  .669_JPRB /)
455ASY_DD(:,14) = (/ 0.645_JPRB, 0.755_JPRB, 0.682_JPRB /)
456ASY_DD(:,15) = (/ 0.556_JPRB, 0.716_JPRB, 0.718_JPRB /)
457ASY_DD(:,16) = (/ 0.392_JPRB, 0.652_JPRB, 0.730_JPRB /)
458
459OMG_DD( :,1) = (/  .974_JPRB,  .944_JPRB,  .894_JPRB /)
460OMG_DD( :,2) = (/  .974_JPRB,  .944_JPRB,  .894_JPRB /)
461OMG_DD( :,3) = (/  .974_JPRB,  .944_JPRB,  .894_JPRB /)
462OMG_DD( :,4) = (/  .978_JPRB,  .947_JPRB,  .900_JPRB /)
463OMG_DD( :,5) = (/ 0.980_JPRB, 0.947_JPRB, 0.905_JPRB /)
464OMG_DD( :,6) = (/  .985_JPRB,  .957_JPRB,  .926_JPRB /)
465OMG_DD( :,7) = (/  .988_JPRB,  .963_JPRB,  .938_JPRB /)
466OMG_DD( :,8) = (/ 0.987_JPRB, 0.960_JPRB, 0.933_JPRB /)
467OMG_DD( :,9) = (/ 0.990_JPRB, 0.973_JPRB, 0.946_JPRB /)
468OMG_DD(:,10) = (/  .996_JPRB,  .990_JPRB,  .978_JPRB /)
469OMG_DD(:,11) = (/  .997_JPRB,  .994_JPRB,  .982_JPRB /)
470OMG_DD(:,12) = (/ 0.992_JPRB, 0.988_JPRB, 0.958_JPRB /)
471OMG_DD(:,13) = (/  .997_JPRB,  .996_JPRB,  .984_JPRB /)
472OMG_DD(:,14) = (/ 0.996_JPRB, 0.997_JPRB, 0.989_JPRB /)
473OMG_DD(:,15) = (/ 0.995_JPRB, 0.997_JPRB, 0.994_JPRB /)
474OMG_DD(:,16) = (/ 0.994_JPRB, 0.997_JPRB, 0.996_JPRB /)
475
476
477!-- For sea-salt, XX( RH, wavelength, binNo)
478
479! BIN No. 1
480! optical thickness
481ALF_SS( :, 1, 1) = (/ &
482  & 1.299_JPRB, 1.299_JPRB, 1.299_JPRB, 1.299_JPRB, 3.482_JPRB, 3.940_JPRB, &
483  & 4.572_JPRB, 5.358_JPRB, 6.455_JPRB, 7.261_JPRB, 8.565_JPRB,11.796_JPRB /)
484ALF_SS( :, 2, 1) = (/ &
485  & 1.299_JPRB, 1.299_JPRB, 1.299_JPRB, 1.299_JPRB, 3.482_JPRB, 3.940_JPRB, &
486  & 4.572_JPRB, 5.358_JPRB, 6.455_JPRB, 7.261_JPRB, 8.565_JPRB,11.796_JPRB /)
487ALF_SS( :, 3, 1) = (/ &
488  & 1.299_JPRB, 1.299_JPRB, 1.299_JPRB, 1.299_JPRB, 3.482_JPRB, 3.940_JPRB, &
489  & 4.572_JPRB, 5.358_JPRB, 6.455_JPRB, 7.261_JPRB, 8.565_JPRB,11.796_JPRB /)
490ALF_SS( :, 4, 1) = (/ &
491  & 1.143_JPRB, 1.143_JPRB, 1.143_JPRB, 1.143_JPRB, 3.393_JPRB, 3.761_JPRB, &
492  & 4.436_JPRB, 5.291_JPRB, 6.530_JPRB, 7.452_JPRB, 8.917_JPRB,12.189_JPRB /)
493ALF_SS( :, 5, 1) = (/ &
494  & 1.031_JPRB, 1.031_JPRB, 1.031_JPRB, 1.031_JPRB, 3.266_JPRB, 3.589_JPRB, &
495  & 4.286_JPRB, 5.165_JPRB, 6.472_JPRB, 7.496_JPRB, 9.102_JPRB,12.541_JPRB /)
496
497ALF_SS( :, 6, 1) = (/ &
498  &  .927_JPRB,  .927_JPRB,  .927_JPRB,  .927_JPRB, 3.124_JPRB, 3.387_JPRB, &
499  & 4.083_JPRB, 4.985_JPRB, 6.361_JPRB, 7.426_JPRB, 9.175_JPRB,12.904_JPRB /)
500ALF_SS( :, 7, 1) = (/ &
501  &  .830_JPRB,  .830_JPRB,  .830_JPRB,  .830_JPRB, 2.943_JPRB, 3.174_JPRB, &
502  & 3.855_JPRB, 4.761_JPRB, 6.156_JPRB, 7.300_JPRB, 9.170_JPRB,13.236_JPRB /)
503
504ALF_SS( :, 8, 1) = (/ &
505  & 0.763_JPRB, 0.763_JPRB, 0.763_JPRB, 0.763_JPRB, 2.806_JPRB, 3.012_JPRB, &
506  & 3.691_JPRB, 4.578_JPRB, 5.999_JPRB, 7.146_JPRB, 9.085_JPRB,13.416_JPRB /)
507ALF_SS( :, 9, 1) = (/ &
508  & 0.545_JPRB, 0.545_JPRB, 0.545_JPRB, 0.545_JPRB, 2.301_JPRB, 2.434_JPRB, &
509  & 3.035_JPRB, 3.857_JPRB, 5.227_JPRB, 6.386_JPRB, 8.466_JPRB,13.570_JPRB /)
510
511ALF_SS( :,10, 1) = (/ &
512  &  .500_JPRB,  .500_JPRB,  .500_JPRB,  .500_JPRB, 2.162_JPRB, 2.283_JPRB, &
513  & 2.869_JPRB, 3.662_JPRB, 5.004_JPRB, 6.147_JPRB, 8.239_JPRB,13.447_JPRB /)
514ALF_SS( :,11, 1) = (/ &
515  &  .314_JPRB,  .314_JPRB,  .314_JPRB,  .314_JPRB, 1.575_JPRB, 1.648_JPRB, &
516  & 2.106_JPRB, 2.754_JPRB, 3.898_JPRB, 4.910_JPRB, 6.855_JPRB,12.275_JPRB /)
517ALF_SS( :,12, 1) = (/ &
518  & 0.252_JPRB, 0.252_JPRB, 0.252_JPRB, 0.252_JPRB, 1.350_JPRB, 1.417_JPRB, &
519  & 1.836_JPRB, 2.425_JPRB, 3.458_JPRB, 4.416_JPRB, 6.229_JPRB,11.532_JPRB /)
520
521ALF_SS( :,13, 1) = (/ &
522  &  .246_JPRB,  .246_JPRB,  .246_JPRB,  .246_JPRB, 1.325_JPRB, 1.391_JPRB, &
523  & 1.805_JPRB, 2.388_JPRB, 3.408_JPRB, 4.357_JPRB, 6.157_JPRB,11.438_JPRB /)
524ALF_SS( :,14, 1) = (/ &
525  & 0.078_JPRB, 0.078_JPRB, 0.078_JPRB, 0.078_JPRB, 0.516_JPRB, 0.557_JPRB, &
526  & 0.759_JPRB, 1.059_JPRB, 1.601_JPRB, 2.125_JPRB, 3.232_JPRB, 6.884_JPRB /)
527ALF_SS( :,15, 1) = (/ &
528  & 0.028_JPRB, 0.028_JPRB, 0.028_JPRB, 0.028_JPRB, 0.213_JPRB, 0.237_JPRB, &
529  & 0.327_JPRB, 0.466_JPRB, 0.746_JPRB, 1.037_JPRB, 1.666_JPRB, 3.881_JPRB /)
530ALF_SS( :,16, 1) = (/ &
531  & 0.010_JPRB, 0.010_JPRB, 0.010_JPRB, 0.010_JPRB, 0.086_JPRB, 0.099_JPRB, &
532  & 0.141_JPRB, 0.207_JPRB, 0.338_JPRB, 0.472_JPRB, 0.773_JPRB, 2.013_JPRB /)
533
534
535
536
537
538
539! assymetry factor
540ASY_SS( :, 1, 1) = (/ &
541  &  .698_JPRB,  .698_JPRB,  .698_JPRB,  .698_JPRB,  .707_JPRB,  .754_JPRB, &
542  &  .756_JPRB,  .755_JPRB,  .750_JPRB,  .744_JPRB,  .732_JPRB,  .717_JPRB /)
543ASY_SS( :, 2, 1) = (/ &
544  &  .698_JPRB,  .698_JPRB,  .698_JPRB,  .698_JPRB,  .707_JPRB,  .754_JPRB, &
545  &  .756_JPRB,  .755_JPRB,  .750_JPRB,  .744_JPRB,  .732_JPRB,  .717_JPRB /)
546ASY_SS( :, 3, 1) = (/ &
547  &  .698_JPRB,  .698_JPRB,  .698_JPRB,  .698_JPRB,  .707_JPRB,  .754_JPRB, &
548  &  .756_JPRB,  .755_JPRB,  .750_JPRB,  .744_JPRB,  .732_JPRB,  .717_JPRB /)
549ASY_SS( :, 4, 1) = (/ &
550  &  .682_JPRB,  .682_JPRB,  .682_JPRB,  .682_JPRB,  .711_JPRB,  .755_JPRB, &
551  &  .759_JPRB,  .762_JPRB,  .760_JPRB,  .757_JPRB,  .747_JPRB,  .724_JPRB /)
552ASY_SS( :, 5, 1) = (/ &
553  & 0.670_JPRB, 0.670_JPRB, 0.670_JPRB, 0.670_JPRB, 0.711_JPRB, 0.753_JPRB, &
554  & 0.759_JPRB, 0.763_JPRB, 0.765_JPRB, 0.763_JPRB, 0.756_JPRB, 0.733_JPRB /)
555
556ASY_SS( :, 6, 1) = (/ &
557  &  .659_JPRB,  .659_JPRB,  .659_JPRB,  .659_JPRB,  .709_JPRB,  .749_JPRB, &
558  &  .756_JPRB,  .763_JPRB,  .767_JPRB,  .767_JPRB,  .764_JPRB,  .743_JPRB /)
559ASY_SS( :, 7, 1) = (/ &
560  &  .645_JPRB,  .645_JPRB,  .645_JPRB,  .645_JPRB,  .705_JPRB,  .743_JPRB, &
561  &  .752_JPRB,  .761_JPRB,  .767_JPRB,  .770_JPRB,  .768_JPRB,  .751_JPRB /)
562
563ASY_SS( :, 8, 1) = (/ &
564  & 0.633_JPRB, 0.633_JPRB, 0.633_JPRB, 0.633_JPRB, 0.701_JPRB, 0.738_JPRB, &
565  & 0.749_JPRB, 0.758_JPRB, 0.767_JPRB, 0.769_JPRB, 0.771_JPRB, 0.758_JPRB /)
566ASY_SS( :, 9, 1) = (/ &
567  & 0.581_JPRB, 0.581_JPRB, 0.581_JPRB, 0.581_JPRB, 0.679_JPRB, 0.715_JPRB, &
568  & 0.729_JPRB, 0.743_JPRB, 0.757_JPRB, 0.764_JPRB, 0.772_JPRB, 0.772_JPRB /)
569
570ASY_SS( :,10, 1) = (/ &
571  &  .569_JPRB,  .569_JPRB,  .569_JPRB,  .569_JPRB,  .672_JPRB,  .707_JPRB, &
572  &  .723_JPRB,  .738_JPRB,  .753_JPRB,  .761_JPRB,  .771_JPRB,  .774_JPRB /)
573ASY_SS( :,11, 1) = (/ &
574  &  .513_JPRB,  .513_JPRB,  .513_JPRB,  .513_JPRB,  .633_JPRB,  .667_JPRB, &
575  &  .685_JPRB,  .704_JPRB,  .726_JPRB,  .740_JPRB,  .756_JPRB,  .774_JPRB /)
576ASY_SS( :,12, 1) = (/ &
577  & 0.475_JPRB, 0.475_JPRB, 0.475_JPRB, 0.475_JPRB, 0.609_JPRB, 0.644_JPRB, &
578  & 0.667_JPRB, 0.688_JPRB, 0.712_JPRB, 0.727_JPRB, 0.746_JPRB, 0.770_JPRB /)
579
580ASY_SS( :,13, 1) = (/ &
581  &  .470_JPRB,  .470_JPRB,  .470_JPRB,  .470_JPRB,  .606_JPRB,  .641_JPRB, &
582  &  .665_JPRB,  .686_JPRB,  .710_JPRB,  .726_JPRB,  .745_JPRB,  .769_JPRB /)
583ASY_SS( :,14, 1) = (/ &
584  & 0.230_JPRB, 0.230_JPRB, 0.230_JPRB, 0.230_JPRB, 0.473_JPRB, 0.506_JPRB, &
585  & 0.539_JPRB, 0.568_JPRB, 0.601_JPRB, 0.627_JPRB, 0.664_JPRB, 0.716_JPRB /)
586ASY_SS( :,15, 1) = (/ &
587  & 0.128_JPRB, 0.128_JPRB, 0.128_JPRB, 0.128_JPRB, 0.276_JPRB, 0.311_JPRB, &
588  & 0.357_JPRB, 0.413_JPRB, 0.482_JPRB, 0.522_JPRB, 0.566_JPRB, 0.635_JPRB /)
589ASY_SS( :,16, 1) = (/ &
590  & 0.076_JPRB, 0.076_JPRB, 0.076_JPRB, 0.076_JPRB, 0.158_JPRB, 0.179_JPRB, &
591  & 0.206_JPRB, 0.242_JPRB, 0.297_JPRB, 0.343_JPRB, 0.417_JPRB, 0.540_JPRB /)
592
593! single scattering albedo
594OMG_SS( :, 1, 1) = (/ &
595  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
596  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
597OMG_SS( :, 2, 1) = (/ &
598  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
599  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
600OMG_SS( :, 3, 1) = (/ &
601  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
602  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
603OMG_SS( :, 4, 1) = (/ &
604  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
605  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
606OMG_SS( :, 5, 1) = (/ &
607  & 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 1.000_JPRB, &
608  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
609
610OMG_SS( :, 6, 1) = (/ &
611  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
612  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
613OMG_SS( :, 7, 1) = (/ &
614  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
615  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
616
617OMG_SS( :, 8, 1) = (/ &
618  & 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 1.000_JPRB, &
619  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
620OMG_SS( :, 9, 1) = (/ &
621  & 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 1.000_JPRB, &
622  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
623
624OMG_SS( :,10, 1) = (/ &
625  &  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB, 1.000_JPRB, &
626  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
627OMG_SS( :,11, 1) = (/ &
628  &  .998_JPRB,  .998_JPRB,  .998_JPRB,  .998_JPRB,  .999_JPRB, 1.000_JPRB, &
629  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
630OMG_SS( :,12, 1) = (/ &
631  & 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB, 1.000_JPRB, &
632  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
633
634OMG_SS( :,13, 1) = (/ &
635  &  .998_JPRB,  .998_JPRB,  .998_JPRB,  .998_JPRB,  .999_JPRB, 1.000_JPRB, &
636  & 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
637OMG_SS( :,14, 1) = (/ &
638  & 0.997_JPRB, 0.997_JPRB, 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.999_JPRB, &
639  & 0.999_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
640OMG_SS( :,15, 1) = (/ &
641  & 0.994_JPRB, 0.994_JPRB, 0.994_JPRB, 0.994_JPRB, 0.997_JPRB, 0.999_JPRB, &
642  & 0.999_JPRB, 0.999_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
643OMG_SS( :,16, 1) = (/ &
644  & 0.988_JPRB, 0.988_JPRB, 0.988_JPRB, 0.988_JPRB, 0.995_JPRB, 0.998_JPRB, &
645  & 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 1.000_JPRB, 1.000_JPRB, 1.000_JPRB /)
646
647! BIN No. 2
648ALF_SS( :, 1, 2) = (/ &
649  &  .137_JPRB,  .137_JPRB,  .137_JPRB,  .137_JPRB,  .264_JPRB,  .309_JPRB, &
650  &  .351_JPRB,  .406_JPRB,  .494_JPRB,  .567_JPRB,  .697_JPRB, 1.031_JPRB /)
651ALF_SS( :, 2, 2) = (/ &
652  &  .137_JPRB,  .137_JPRB,  .137_JPRB,  .137_JPRB,  .264_JPRB,  .309_JPRB, &
653  &  .351_JPRB,  .406_JPRB,  .494_JPRB,  .567_JPRB,  .697_JPRB, 1.031_JPRB /)
654ALF_SS( :, 3, 2) = (/ &
655  &  .137_JPRB,  .137_JPRB,  .137_JPRB,  .137_JPRB,  .264_JPRB,  .309_JPRB, &
656  &  .351_JPRB,  .406_JPRB,  .494_JPRB,  .567_JPRB,  .697_JPRB, 1.031_JPRB /)
657ALF_SS( :, 4, 2) = (/ &
658  &  .140_JPRB,  .140_JPRB,  .140_JPRB,  .140_JPRB,  .268_JPRB,  .313_JPRB, &
659  &  .356_JPRB,  .411_JPRB,  .496_JPRB,  .568_JPRB,  .697_JPRB, 1.035_JPRB /)
660ALF_SS( :, 5, 2) = (/ &
661  & 0.141_JPRB, 0.141_JPRB, 0.141_JPRB, 0.141_JPRB, 0.272_JPRB, 0.317_JPRB, &
662  & 0.360_JPRB, 0.415_JPRB, 0.500_JPRB, 0.571_JPRB, 0.697_JPRB, 1.033_JPRB /)
663
664ALF_SS( :, 6, 2) = (/ &
665  &  .143_JPRB,  .143_JPRB,  .143_JPRB,  .143_JPRB,  .274_JPRB,  .322_JPRB, &
666  &  .365_JPRB,  .420_JPRB,  .505_JPRB,  .574_JPRB,  .699_JPRB, 1.033_JPRB /)
667ALF_SS( :, 7, 2) = (/ &
668  &  .144_JPRB,  .144_JPRB,  .144_JPRB,  .144_JPRB,  .279_JPRB,  .326_JPRB, &
669  &  .370_JPRB,  .425_JPRB,  .510_JPRB,  .578_JPRB,  .701_JPRB, 1.032_JPRB /)
670
671ALF_SS( :, 8, 2) = (/ &
672  & 0.144_JPRB, 0.144_JPRB, 0.144_JPRB, 0.144_JPRB, 0.281_JPRB, 0.329_JPRB, &
673  & 0.373_JPRB, 0.429_JPRB, 0.514_JPRB, 0.584_JPRB, 0.707_JPRB, 1.032_JPRB /)
674ALF_SS( :, 9, 2) = (/ &
675  & 0.146_JPRB, 0.146_JPRB, 0.146_JPRB, 0.146_JPRB, 0.290_JPRB, 0.337_JPRB, &
676  & 0.384_JPRB, 0.443_JPRB, 0.530_JPRB, 0.602_JPRB, 0.724_JPRB, 1.041_JPRB /)
677
678ALF_SS( :,10, 2) = (/ &
679  &  .147_JPRB,  .147_JPRB,  .147_JPRB,  .147_JPRB,  .292_JPRB,  .340_JPRB, &
680  &  .385_JPRB,  .446_JPRB,  .533_JPRB,  .606_JPRB,  .728_JPRB, 1.047_JPRB /)
681ALF_SS( :,11, 2) = (/ &
682  &  .146_JPRB,  .146_JPRB,  .146_JPRB,  .146_JPRB,  .299_JPRB,  .345_JPRB, &
683  &  .395_JPRB,  .459_JPRB,  .551_JPRB,  .628_JPRB,  .758_JPRB, 1.078_JPRB /)
684ALF_SS( :,12, 2) = (/ &
685  & 0.147_JPRB, 0.147_JPRB, 0.147_JPRB, 0.147_JPRB, 0.300_JPRB, 0.348_JPRB, &
686  & 0.397_JPRB, 0.459_JPRB, 0.559_JPRB, 0.634_JPRB, 0.768_JPRB, 1.095_JPRB /)
687
688ALF_SS( :,13, 2) = (/ &
689  &  .147_JPRB,  .147_JPRB,  .147_JPRB,  .147_JPRB,  .301_JPRB,  .348_JPRB, &
690  &  .398_JPRB,  .459_JPRB,  .559_JPRB,  .635_JPRB,  .769_JPRB, 1.098_JPRB /)
691ALF_SS( :,14, 2) = (/ &
692  & 0.148_JPRB, 0.148_JPRB, 0.148_JPRB, 0.148_JPRB, 0.305_JPRB, 0.349_JPRB, &
693  & 0.403_JPRB, 0.466_JPRB, 0.566_JPRB, 0.651_JPRB, 0.802_JPRB, 1.167_JPRB /)
694ALF_SS( :,15, 2) = (/ &
695  & 0.146_JPRB, 0.146_JPRB, 0.146_JPRB, 0.146_JPRB, 0.302_JPRB, 0.351_JPRB, &
696  & 0.399_JPRB, 0.461_JPRB, 0.563_JPRB, 0.651_JPRB, 0.805_JPRB, 1.183_JPRB /)
697ALF_SS( :,16, 2) = (/ &
698  & 0.155_JPRB, 0.155_JPRB, 0.155_JPRB, 0.155_JPRB, 0.305_JPRB, 0.348_JPRB, &
699  & 0.402_JPRB, 0.472_JPRB, 0.576_JPRB, 0.657_JPRB, 0.793_JPRB, 1.187_JPRB /)
700
701! assymetry factor
702ASY_SS( :, 1, 2) = (/ &
703  &  .741_JPRB,  .741_JPRB,  .741_JPRB,  .741_JPRB,  .755_JPRB,  .778_JPRB, &
704  &  .782_JPRB,  .789_JPRB,  .796_JPRB,  .801_JPRB,  .809_JPRB,  .819_JPRB /)
705ASY_SS( :, 2, 2) = (/ &
706  &  .741_JPRB,  .741_JPRB,  .741_JPRB,  .741_JPRB,  .755_JPRB,  .778_JPRB, &
707  &  .782_JPRB,  .789_JPRB,  .796_JPRB,  .801_JPRB,  .809_JPRB,  .819_JPRB /)
708ASY_SS( :, 3, 2) = (/ &
709  &  .741_JPRB,  .741_JPRB,  .741_JPRB,  .741_JPRB,  .755_JPRB,  .778_JPRB, &
710  &  .782_JPRB,  .789_JPRB,  .796_JPRB,  .801_JPRB,  .809_JPRB,  .819_JPRB /)
711ASY_SS( :, 4, 2) = (/ &
712  &  .739_JPRB,  .739_JPRB,  .739_JPRB,  .739_JPRB,  .750_JPRB,  .776_JPRB, &
713  &  .779_JPRB,  .784_JPRB,  .792_JPRB,  .795_JPRB,  .802_JPRB,  .816_JPRB /)
714ASY_SS( :, 5, 2) = (/ &
715  & 0.737_JPRB, 0.737_JPRB, 0.737_JPRB, 0.737_JPRB, 0.747_JPRB, 0.773_JPRB, &
716  & 0.778_JPRB, 0.782_JPRB, 0.788_JPRB, 0.791_JPRB, 0.798_JPRB, 0.812_JPRB /)
717
718ASY_SS( :, 6, 2) = (/ &
719  &  .734_JPRB,  .734_JPRB,  .734_JPRB,  .734_JPRB,  .745_JPRB,  .772_JPRB, &
720  &  .775_JPRB,  .780_JPRB,  .785_JPRB,  .789_JPRB,  .795_JPRB,  .807_JPRB /)
721ASY_SS( :, 7, 2) = (/ &
722  &  .732_JPRB,  .732_JPRB,  .732_JPRB,  .732_JPRB,  .744_JPRB,  .770_JPRB, &
723  &  .775_JPRB,  .778_JPRB,  .784_JPRB,  .787_JPRB,  .793_JPRB,  .804_JPRB /)
724
725ASY_SS( :, 8, 2) = (/ &
726  & 0.730_JPRB, 0.730_JPRB, 0.730_JPRB, 0.730_JPRB, 0.742_JPRB, 0.770_JPRB, &
727  & 0.773_JPRB, 0.779_JPRB, 0.782_JPRB, 0.785_JPRB, 0.790_JPRB, 0.802_JPRB /)
728ASY_SS( :, 9, 2) = (/ &
729  & 0.723_JPRB, 0.723_JPRB, 0.723_JPRB, 0.723_JPRB, 0.739_JPRB, 0.766_JPRB, &
730  & 0.771_JPRB, 0.776_JPRB, 0.780_JPRB, 0.783_JPRB, 0.786_JPRB, 0.794_JPRB /)
731
732ASY_SS( :,10, 2) = (/ &
733  &  .720_JPRB,  .720_JPRB,  .720_JPRB,  .720_JPRB,  .737_JPRB,  .765_JPRB, &
734  &  .770_JPRB,  .774_JPRB,  .780_JPRB,  .782_JPRB,  .785_JPRB,  .793_JPRB /)
735ASY_SS( :,11, 2) = (/ &
736  &  .705_JPRB,  .705_JPRB,  .705_JPRB,  .705_JPRB,  .730_JPRB,  .758_JPRB, &
737  &  .763_JPRB,  .770_JPRB,  .775_JPRB,  .777_JPRB,  .781_JPRB,  .787_JPRB /)
738ASY_SS( :,12, 2) = (/ &
739  & 0.701_JPRB, 0.701_JPRB, 0.701_JPRB, 0.701_JPRB, 0.727_JPRB, 0.757_JPRB, &
740  & 0.760_JPRB, 0.767_JPRB, 0.773_JPRB, 0.776_JPRB, 0.780_JPRB, 0.788_JPRB /)
741
742ASY_SS( :,13, 2) = (/ &
743  &  .701_JPRB,  .701_JPRB,  .701_JPRB,  .701_JPRB,  .726_JPRB,  .756_JPRB, &
744  &  .760_JPRB,  .766_JPRB,  .773_JPRB,  .776_JPRB,  .780_JPRB,  .787_JPRB /)
745ASY_SS( :,14, 2) = (/ &
746  & 0.671_JPRB, 0.671_JPRB, 0.671_JPRB, 0.671_JPRB, 0.701_JPRB, 0.731_JPRB, &
747  & 0.738_JPRB, 0.747_JPRB, 0.753_JPRB, 0.759_JPRB, 0.767_JPRB, 0.779_JPRB /)
748ASY_SS( :,15, 2) = (/ &
749  & 0.642_JPRB, 0.642_JPRB, 0.642_JPRB, 0.642_JPRB, 0.675_JPRB, 0.711_JPRB, &
750  & 0.716_JPRB, 0.723_JPRB, 0.733_JPRB, 0.742_JPRB, 0.751_JPRB, 0.764_JPRB /)
751ASY_SS( :,16, 2) = (/ &
752  & 0.662_JPRB, 0.662_JPRB, 0.662_JPRB, 0.662_JPRB, 0.653_JPRB, 0.690_JPRB, &
753  & 0.699_JPRB, 0.709_JPRB, 0.720_JPRB, 0.724_JPRB, 0.730_JPRB, 0.748_JPRB /)
754
755! single scattering albedo
756OMG_SS( :, 1, 2) = (/ &
757  &  .990_JPRB,  .990_JPRB,  .990_JPRB,  .990_JPRB,  .986_JPRB,  .994_JPRB, &
758  &  .995_JPRB,  .996_JPRB,  .997_JPRB,  .997_JPRB,  .998_JPRB,  .999_JPRB /)
759OMG_SS( :, 2, 2) = (/ &
760  &  .990_JPRB,  .990_JPRB,  .990_JPRB,  .990_JPRB,  .986_JPRB,  .994_JPRB, &
761  &  .995_JPRB,  .996_JPRB,  .997_JPRB,  .997_JPRB,  .998_JPRB,  .999_JPRB /)
762OMG_SS( :, 3, 2) = (/ &
763  &  .990_JPRB,  .990_JPRB,  .990_JPRB,  .990_JPRB,  .986_JPRB,  .994_JPRB, &
764  &  .995_JPRB,  .996_JPRB,  .997_JPRB,  .997_JPRB,  .998_JPRB,  .999_JPRB /)
765OMG_SS( :, 4, 2) = (/ &
766  &  .991_JPRB,  .991_JPRB,  .991_JPRB,  .991_JPRB,  .987_JPRB,  .995_JPRB, &
767  &  .995_JPRB,  .996_JPRB,  .997_JPRB,  .997_JPRB,  .998_JPRB,  .999_JPRB /)
768OMG_SS( :, 5, 2) = (/ &
769  & 0.992_JPRB, 0.992_JPRB, 0.992_JPRB, 0.992_JPRB, 0.988_JPRB, 0.995_JPRB, &
770  & 0.996_JPRB, 0.996_JPRB, 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.999_JPRB /)
771
772OMG_SS( :, 6, 2) = (/ &
773  &  .992_JPRB,  .992_JPRB,  .992_JPRB,  .992_JPRB,  .989_JPRB,  .996_JPRB, &
774  &  .996_JPRB,  .997_JPRB,  .997_JPRB,  .997_JPRB,  .998_JPRB,  .999_JPRB /)
775OMG_SS( :, 7, 2) = (/ &
776  &  .993_JPRB,  .993_JPRB,  .993_JPRB,  .993_JPRB,  .989_JPRB,  .996_JPRB, &
777  &  .996_JPRB,  .997_JPRB,  .997_JPRB,  .998_JPRB,  .998_JPRB,  .999_JPRB /)
778
779OMG_SS( :, 8, 2) = (/ &
780  & 0.993_JPRB, 0.993_JPRB, 0.993_JPRB, 0.993_JPRB, 0.990_JPRB, 0.996_JPRB, &
781  & 0.996_JPRB, 0.997_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB /)
782OMG_SS( :, 9, 2) = (/ &
783  & 0.994_JPRB, 0.994_JPRB, 0.994_JPRB, 0.994_JPRB, 0.991_JPRB, 0.997_JPRB, &
784  & 0.997_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB /)
785
786OMG_SS( :,10, 2) = (/ &
787  &  .994_JPRB,  .994_JPRB,  .994_JPRB,  .994_JPRB,  .992_JPRB,  .997_JPRB, &
788  &  .997_JPRB,  .997_JPRB,  .998_JPRB,  .998_JPRB,  .999_JPRB,  .999_JPRB /)
789OMG_SS( :,11, 2) = (/ &
790  &  .995_JPRB,  .995_JPRB,  .995_JPRB,  .995_JPRB,  .993_JPRB,  .997_JPRB, &
791  &  .998_JPRB,  .998_JPRB,  .998_JPRB,  .998_JPRB,  .999_JPRB,  .999_JPRB /)
792OMG_SS( :,12, 2) = (/ &
793  & 0.995_JPRB, 0.995_JPRB, 0.995_JPRB, 0.995_JPRB, 0.993_JPRB, 0.998_JPRB, &
794  & 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB /)
795
796OMG_SS( :,13, 2) = (/ &
797  &  .995_JPRB,  .995_JPRB,  .995_JPRB,  .995_JPRB,  .993_JPRB,  .998_JPRB, &
798  &  .998_JPRB,  .998_JPRB,  .998_JPRB,  .999_JPRB,  .999_JPRB,  .999_JPRB /)
799OMG_SS( :,14, 2) = (/ &
800  & 0.997_JPRB, 0.997_JPRB, 0.997_JPRB, 0.997_JPRB, 0.995_JPRB, 0.998_JPRB, &
801  & 0.998_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB /)
802OMG_SS( :,15, 2) = (/ &
803  & 0.997_JPRB, 0.997_JPRB, 0.997_JPRB, 0.997_JPRB, 0.996_JPRB, 0.999_JPRB, &
804  & 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 1.000_JPRB /)
805OMG_SS( :,16, 2) = (/ &
806  & 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.997_JPRB, 0.999_JPRB, &
807  & 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 0.999_JPRB, 1.000_JPRB, 1.000_JPRB /)
808
809! BIN No. 3
810ALF_SS( :, 1, 3) = (/ &
811  &  .038_JPRB,  .038_JPRB,  .038_JPRB,  .038_JPRB,  .078_JPRB,  .091_JPRB, &
812  &  .105_JPRB,  .122_JPRB,  .148_JPRB,  .170_JPRB,  .208_JPRB,  .308_JPRB /)
813ALF_SS( :, 2, 3) = (/ &
814  &  .038_JPRB,  .038_JPRB,  .038_JPRB,  .038_JPRB,  .078_JPRB,  .091_JPRB, &
815  &  .105_JPRB,  .122_JPRB,  .148_JPRB,  .170_JPRB,  .208_JPRB,  .308_JPRB /)
816ALF_SS( :, 3, 3) = (/ &
817  &  .038_JPRB,  .038_JPRB,  .038_JPRB,  .038_JPRB,  .078_JPRB,  .091_JPRB, &
818  &  .105_JPRB,  .122_JPRB,  .148_JPRB,  .170_JPRB,  .208_JPRB,  .308_JPRB /)
819ALF_SS( :, 4, 3) = (/ &
820  &  .038_JPRB,  .038_JPRB,  .038_JPRB,  .038_JPRB,  .079_JPRB,  .091_JPRB, &
821  &  .105_JPRB,  .122_JPRB,  .149_JPRB,  .171_JPRB,  .209_JPRB,  .308_JPRB /)
822ALF_SS( :, 5, 3) = (/ &
823  & 0.039_JPRB, 0.039_JPRB, 0.039_JPRB, 0.039_JPRB, 0.079_JPRB, 0.092_JPRB, &
824  & 0.105_JPRB, 0.122_JPRB, 0.149_JPRB, 0.171_JPRB, 0.210_JPRB, 0.309_JPRB /)
825
826ALF_SS( :, 6, 3) = (/ &
827  &  .039_JPRB,  .039_JPRB,  .039_JPRB,  .039_JPRB,  .079_JPRB,  .092_JPRB, &
828  &  .105_JPRB,  .123_JPRB,  .149_JPRB,  .171_JPRB,  .209_JPRB,  .309_JPRB /)
829ALF_SS( :, 7, 3) = (/ &
830  &  .039_JPRB,  .039_JPRB,  .039_JPRB,  .039_JPRB,  .079_JPRB,  .092_JPRB, &
831  &  .105_JPRB,  .123_JPRB,  .149_JPRB,  .172_JPRB,  .210_JPRB,  .310_JPRB /)
832
833ALF_SS( :, 8, 3) = (/ &
834  & 0.039_JPRB, 0.039_JPRB, 0.039_JPRB, 0.039_JPRB, 0.079_JPRB, 0.092_JPRB, &
835  & 0.106_JPRB, 0.123_JPRB, 0.150_JPRB, 0.172_JPRB, 0.210_JPRB, 0.310_JPRB /)
836ALF_SS( :, 9, 3) = (/ &
837  & 0.039_JPRB, 0.039_JPRB, 0.039_JPRB, 0.039_JPRB, 0.080_JPRB, 0.093_JPRB, &
838  & 0.106_JPRB, 0.123_JPRB, 0.151_JPRB, 0.172_JPRB, 0.211_JPRB, 0.311_JPRB /)
839
840ALF_SS( :,10, 3) = (/ &
841  &  .039_JPRB,  .039_JPRB,  .039_JPRB,  .039_JPRB,  .080_JPRB,  .093_JPRB, &
842  &  .107_JPRB,  .124_JPRB,  .151_JPRB,  .173_JPRB,  .212_JPRB,  .311_JPRB /)
843ALF_SS( :,11, 3) = (/ &
844  &  .040_JPRB,  .040_JPRB,  .040_JPRB,  .040_JPRB,  .081_JPRB,  .094_JPRB, &
845  &  .108_JPRB,  .124_JPRB,  .152_JPRB,  .174_JPRB,  .212_JPRB,  .314_JPRB /)
846ALF_SS( :,12, 3) = (/ &
847  & 0.040_JPRB, 0.040_JPRB, 0.040_JPRB, 0.040_JPRB, 0.081_JPRB, 0.094_JPRB, &
848  & 0.108_JPRB, 0.126_JPRB, 0.152_JPRB, 0.174_JPRB, 0.213_JPRB, 0.314_JPRB /)
849
850ALF_SS( :,13, 3) = (/ &
851  &  .040_JPRB,  .040_JPRB,  .040_JPRB,  .040_JPRB,  .081_JPRB,  .094_JPRB, &
852  &  .108_JPRB,  .126_JPRB,  .152_JPRB,  .174_JPRB,  .213_JPRB,  .314_JPRB /)
853ALF_SS( :,14, 3) = (/ &
854  & 0.040_JPRB, 0.040_JPRB, 0.040_JPRB, 0.040_JPRB, 0.082_JPRB, 0.096_JPRB, &
855  & 0.109_JPRB, 0.127_JPRB, 0.156_JPRB, 0.178_JPRB, 0.215_JPRB, 0.318_JPRB /)
856ALF_SS( :,15, 3) = (/ &
857  & 0.043_JPRB, 0.043_JPRB, 0.043_JPRB, 0.043_JPRB, 0.084_JPRB, 0.097_JPRB, &
858  & 0.112_JPRB, 0.131_JPRB, 0.159_JPRB, 0.180_JPRB, 0.218_JPRB, 0.324_JPRB /)
859ALF_SS( :,16, 3) = (/ &
860  & 0.042_JPRB, 0.042_JPRB, 0.042_JPRB, 0.042_JPRB, 0.086_JPRB, 0.102_JPRB, &
861  & 0.115_JPRB, 0.131_JPRB, 0.157_JPRB, 0.181_JPRB, 0.225_JPRB, 0.326_JPRB /)
862
863! assymetry factor
864ASY_SS( :, 1, 3) = (/ &
865  &  .803_JPRB,  .803_JPRB,  .803_JPRB,  .803_JPRB,  .810_JPRB,  .828_JPRB, &
866  &  .831_JPRB,  .834_JPRB,  .837_JPRB,  .839_JPRB,  .842_JPRB,  .845_JPRB /)
867ASY_SS( :, 2, 3) = (/ &
868  &  .803_JPRB,  .803_JPRB,  .803_JPRB,  .803_JPRB,  .810_JPRB,  .828_JPRB, &
869  &  .831_JPRB,  .834_JPRB,  .837_JPRB,  .839_JPRB,  .842_JPRB,  .845_JPRB /)
870ASY_SS( :, 3, 3) = (/ &
871  &  .803_JPRB,  .803_JPRB,  .803_JPRB,  .803_JPRB,  .810_JPRB,  .828_JPRB, &
872  &  .831_JPRB,  .834_JPRB,  .837_JPRB,  .839_JPRB,  .842_JPRB,  .845_JPRB /)
873ASY_SS( :, 4, 3) = (/ &
874  &  .800_JPRB,  .800_JPRB,  .800_JPRB,  .800_JPRB,  .808_JPRB,  .827_JPRB, &
875  &  .830_JPRB,  .832_JPRB,  .836_JPRB,  .838_JPRB,  .840_JPRB,  .844_JPRB /)
876ASY_SS( :, 5, 3) = (/ &
877  & 0.799_JPRB, 0.799_JPRB, 0.799_JPRB, 0.799_JPRB, 0.807_JPRB, 0.826_JPRB, &
878  & 0.829_JPRB, 0.832_JPRB, 0.834_JPRB, 0.837_JPRB, 0.840_JPRB, 0.845_JPRB /)
879ASY_SS( :, 6, 3) = (/ &
880  &  .796_JPRB,  .796_JPRB,  .796_JPRB,  .796_JPRB,  .806_JPRB,  .824_JPRB, &
881  &  .828_JPRB,  .830_JPRB,  .834_JPRB,  .836_JPRB,  .839_JPRB,  .843_JPRB /)
882ASY_SS( :, 7, 3) = (/ &
883  &  .795_JPRB,  .795_JPRB,  .795_JPRB,  .795_JPRB,  .804_JPRB,  .823_JPRB, &
884  &  .826_JPRB,  .829_JPRB,  .833_JPRB,  .835_JPRB,  .837_JPRB,  .842_JPRB /)
885
886ASY_SS( :, 8, 3) = (/ &
887  & 0.794_JPRB, 0.794_JPRB, 0.794_JPRB, 0.794_JPRB, 0.803_JPRB, 0.822_JPRB, &
888  & 0.826_JPRB, 0.829_JPRB, 0.832_JPRB, 0.834_JPRB, 0.837_JPRB, 0.841_JPRB /)
889ASY_SS( :, 9, 3) = (/ &
890  & 0.788_JPRB, 0.788_JPRB, 0.788_JPRB, 0.788_JPRB, 0.800_JPRB, 0.819_JPRB, &
891  & 0.823_JPRB, 0.825_JPRB, 0.829_JPRB, 0.832_JPRB, 0.835_JPRB, 0.840_JPRB /)
892
893ASY_SS( :,10, 3) = (/ &
894  &  .786_JPRB,  .786_JPRB,  .786_JPRB,  .786_JPRB,  .800_JPRB,  .818_JPRB, &
895  &  .822_JPRB,  .824_JPRB,  .829_JPRB,  .831_JPRB,  .834_JPRB,  .838_JPRB /)
896ASY_SS( :,11, 3) = (/ &
897  &  .779_JPRB,  .779_JPRB,  .779_JPRB,  .779_JPRB,  .794_JPRB,  .813_JPRB, &
898  &  .817_JPRB,  .820_JPRB,  .824_JPRB,  .827_JPRB,  .830_JPRB,  .836_JPRB /)
899ASY_SS( :,12, 3) = (/ &
900  & 0.774_JPRB, 0.774_JPRB, 0.774_JPRB, 0.774_JPRB, 0.791_JPRB, 0.810_JPRB, &
901  & 0.815_JPRB, 0.819_JPRB, 0.822_JPRB, 0.826_JPRB, 0.829_JPRB, 0.833_JPRB /)
902
903ASY_SS( :,13, 3) = (/ &
904  &  .774_JPRB,  .774_JPRB,  .774_JPRB,  .774_JPRB,  .791_JPRB,  .809_JPRB, &
905  &  .814_JPRB,  .819_JPRB,  .821_JPRB,  .825_JPRB,  .829_JPRB,  .834_JPRB /)
906ASY_SS( :,14, 3) = (/ &
907  & 0.754_JPRB, 0.754_JPRB, 0.754_JPRB, 0.754_JPRB, 0.774_JPRB, 0.795_JPRB, &
908  & 0.799_JPRB, 0.805_JPRB, 0.813_JPRB, 0.816_JPRB, 0.819_JPRB, 0.826_JPRB /)
909ASY_SS( :,15, 3) = (/ &
910  & 0.751_JPRB, 0.751_JPRB, 0.751_JPRB, 0.751_JPRB, 0.762_JPRB, 0.784_JPRB, &
911  & 0.790_JPRB, 0.797_JPRB, 0.802_JPRB, 0.804_JPRB, 0.808_JPRB, 0.819_JPRB /)
912ASY_SS( :,16, 3) = (/ &
913  & 0.708_JPRB, 0.708_JPRB, 0.708_JPRB, 0.708_JPRB, 0.752_JPRB, 0.778_JPRB, &
914  & 0.781_JPRB, 0.783_JPRB, 0.788_JPRB, 0.792_JPRB, 0.801_JPRB, 0.808_JPRB /)
915
916! single scattering albedo
917OMG_SS( :, 1, 3) = (/ &
918  &  .980_JPRB,  .980_JPRB,  .980_JPRB,  .980_JPRB,  .982_JPRB,  .992_JPRB, &
919  &  .993_JPRB,  .994_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
920OMG_SS( :, 2, 3) = (/ &
921  &  .980_JPRB,  .980_JPRB,  .980_JPRB,  .980_JPRB,  .982_JPRB,  .992_JPRB, &
922  &  .993_JPRB,  .994_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
923OMG_SS( :, 3, 3) = (/ &
924  &  .980_JPRB,  .980_JPRB,  .980_JPRB,  .980_JPRB,  .982_JPRB,  .992_JPRB, &
925  &  .993_JPRB,  .994_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
926OMG_SS( :, 4, 3) = (/ &
927  &  .981_JPRB,  .981_JPRB,  .981_JPRB,  .981_JPRB,  .981_JPRB,  .991_JPRB, &
928  &  .993_JPRB,  .994_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
929OMG_SS( :, 5, 3) = (/ &
930  & 0.981_JPRB, 0.981_JPRB, 0.981_JPRB, 0.981_JPRB, 0.981_JPRB, 0.991_JPRB, &
931  & 0.993_JPRB, 0.994_JPRB, 0.995_JPRB, 0.996_JPRB, 0.997_JPRB, 0.998_JPRB /)
932
933OMG_SS( :, 6, 3) = (/ &
934  &  .981_JPRB,  .981_JPRB,  .981_JPRB,  .981_JPRB,  .981_JPRB,  .991_JPRB, &
935  &  .993_JPRB,  .994_JPRB,  .995_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
936OMG_SS( :, 7, 3) = (/ &
937  &  .982_JPRB,  .982_JPRB,  .982_JPRB,  .982_JPRB,  .990_JPRB,  .992_JPRB, &
938  &  .993_JPRB,  .994_JPRB,  .995_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
939
940OMG_SS( :, 8, 3) = (/ &
941  & 0.982_JPRB, 0.982_JPRB, 0.982_JPRB, 0.982_JPRB, 0.990_JPRB, 0.992_JPRB, &
942  & 0.993_JPRB, 0.994_JPRB, 0.995_JPRB, 0.996_JPRB, 0.997_JPRB, 0.998_JPRB /)
943OMG_SS( :, 9, 3) = (/ &
944  & 0.982_JPRB, 0.982_JPRB, 0.982_JPRB, 0.982_JPRB, 0.990_JPRB, 0.992_JPRB, &
945  & 0.993_JPRB, 0.994_JPRB, 0.995_JPRB, 0.996_JPRB, 0.997_JPRB, 0.998_JPRB /)
946
947OMG_SS( :,10, 3) = (/ &
948  &  .983_JPRB,  .983_JPRB,  .983_JPRB,  .983_JPRB,  .987_JPRB,  .992_JPRB, &
949  &  .993_JPRB,  .994_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
950OMG_SS( :,11, 3) = (/ &
951  &  .985_JPRB,  .985_JPRB,  .985_JPRB,  .985_JPRB,  .990_JPRB,  .993_JPRB, &
952  &  .994_JPRB,  .995_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
953OMG_SS( :,12, 3) = (/ &
954  & 0.985_JPRB, 0.985_JPRB, 0.985_JPRB, 0.985_JPRB, 0.990_JPRB, 0.993_JPRB, &
955  & 0.994_JPRB, 0.995_JPRB, 0.996_JPRB, 0.997_JPRB, 0.997_JPRB, 0.998_JPRB /)
956
957OMG_SS( :,13, 3) = (/ &
958  &  .986_JPRB,  .986_JPRB,  .986_JPRB,  .986_JPRB,  .990_JPRB,  .993_JPRB, &
959  &  .994_JPRB,  .995_JPRB,  .996_JPRB,  .996_JPRB,  .997_JPRB,  .998_JPRB /)
960OMG_SS( :,14, 3) = (/ &
961  & 0.989_JPRB, 0.989_JPRB, 0.989_JPRB, 0.989_JPRB, 0.991_JPRB, 0.994_JPRB, &
962  & 0.995_JPRB, 0.996_JPRB, 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.999_JPRB /)
963OMG_SS( :,15, 3) = (/ &
964  & 0.992_JPRB, 0.992_JPRB, 0.992_JPRB, 0.992_JPRB, 0.994_JPRB, 0.996_JPRB, &
965  & 0.996_JPRB, 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB /)
966OMG_SS( :,16, 3) = (/ &
967  & 0.993_JPRB, 0.993_JPRB, 0.993_JPRB, 0.993_JPRB, 0.995_JPRB, 0.997_JPRB, &
968  & 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB /)
969
970!     ----------------------------------------------------------------
971IF (LHOOK) CALL DR_HOOK('SU_AEROP',1,ZHOOK_HANDLE)
972END SUBROUTINE SU_AEROP
Note: See TracBrowser for help on using the repository browser.