LMDZ
su_aerop.F90
Go to the documentation of this file.
1 SUBROUTINE 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 
44 USE parkind1 ,ONLY : jprb
45 USE yomhook ,ONLY : lhook, dr_hook
46 
47 USE 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 
51 IMPLICIT NONE
52 
53 REAL(KIND=JPRB) :: ZHOOK_HANDLE
54 ! ----------------------------------------------------------------
55 IF (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 
80 alf_fa(:) =0.0_jprb
81 asy_fa(:) =0.0_jprb
82 omg_fa(:) =0.0_jprb
83 
84 ! alpha_aer_bc
85 alf_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
91 asy_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
97 omg_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
103 alf_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 /)
106 alf_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 /)
109 alf_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 /)
112 alf_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 /)
115 alf_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 
119 alf_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 /)
122 alf_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 
126 alf_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 /)
129 alf_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 
133 alf_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 /)
136 alf_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 /)
139 alf_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 
143 alf_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 /)
146 alf_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 /)
149 alf_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 /)
152 alf_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
157 asy_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 /)
160 asy_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 /)
163 asy_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 /)
166 asy_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 /)
169 asy_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 
173 asy_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 /)
176 asy_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 
180 asy_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 /)
183 asy_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 
187 asy_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 /)
190 asy_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 /)
193 asy_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 
197 asy_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 /)
200 asy_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 /)
203 asy_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 /)
206 asy_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
211 omg_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 /)
214 omg_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 /)
217 omg_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 /)
220 omg_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 /)
223 omg_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 /)
226 omg_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 /)
229 omg_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 
233 omg_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 /)
236 omg_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 /)
239 omg_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 /)
242 omg_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 /)
245 omg_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 /)
248 omg_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 /)
251 omg_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 /)
254 omg_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 /)
257 omg_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)
263 alf_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 /)
266 alf_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 /)
269 alf_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 /)
272 alf_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 /)
275 alf_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 
279 alf_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 /)
282 alf_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 
286 alf_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 /)
289 alf_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 
293 alf_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 /)
296 alf_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 /)
299 alf_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 
303 alf_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 /)
306 alf_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 /)
309 alf_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 /)
312 alf_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 
316 asy_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 /)
319 asy_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 /)
322 asy_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 /)
325 asy_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 /)
328 asy_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 
332 asy_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 /)
335 asy_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 
339 asy_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 /)
342 asy_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 
346 asy_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 /)
349 asy_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 /)
352 asy_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 
356 asy_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 /)
359 asy_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 /)
362 asy_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 /)
365 asy_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 
369 omg_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 /)
372 omg_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 /)
375 omg_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 /)
378 omg_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 /)
381 omg_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 
385 omg_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 /)
388 omg_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 
392 omg_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 /)
395 omg_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 
399 omg_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 /)
402 omg_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 /)
405 omg_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 
409 omg_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 /)
412 omg_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 /)
415 omg_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 /)
418 omg_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
425 alf_dd( :,1) = (/ 2.210_jprb, 1.011_jprb, .398_jprb /)
426 alf_dd( :,2) = (/ 2.210_jprb, 1.011_jprb, .398_jprb /)
427 alf_dd( :,3) = (/ 2.210_jprb, 1.011_jprb, .398_jprb /)
428 alf_dd( :,4) = (/ 2.357_jprb, .985_jprb, .398_jprb /)
429 alf_dd( :,5) = (/ 2.453_jprb, 0.949_jprb, 0.404_jprb /)
430 alf_dd( :,6) = (/ 2.549_jprb, .901_jprb, .414_jprb /)
431 alf_dd( :,7) = (/ 2.605_jprb, .877_jprb, .423_jprb /)
432 alf_dd( :,8) = (/ 2.635_jprb, 0.869_jprb, 0.428_jprb /)
433 alf_dd( :,9) = (/ 2.596_jprb, 0.978_jprb, 0.428_jprb /)
434 alf_dd(:,10) = (/ 2.555_jprb, 1.039_jprb, .424_jprb /)
435 alf_dd(:,11) = (/ 2.227_jprb, 1.360_jprb, .404_jprb /)
436 alf_dd(:,12) = (/ 2.058_jprb, 1.460_jprb, 0.405_jprb /)
437 alf_dd(:,13) = (/ 2.041_jprb, 1.474_jprb, .405_jprb /)
438 alf_dd(:,14) = (/ 1.117_jprb, 1.523_jprb, 0.502_jprb /)
439 alf_dd(:,15) = (/ 0.583_jprb, 1.117_jprb, 0.596_jprb /)
440 alf_dd(:,16) = (/ 0.263_jprb, 0.705_jprb, 0.607_jprb /)
441 
442 asy_dd( :,1) = (/ .673_jprb, .745_jprb, .799_jprb /)
443 asy_dd( :,2) = (/ .673_jprb, .745_jprb, .799_jprb /)
444 asy_dd( :,3) = (/ .673_jprb, .745_jprb, .799_jprb /)
445 asy_dd( :,4) = (/ .691_jprb, .712_jprb, .792_jprb /)
446 asy_dd( :,5) = (/ 0.704_jprb, 0.679_jprb, 0.791_jprb /)
447 asy_dd( :,6) = (/ .715_jprb, .640_jprb, .786_jprb /)
448 asy_dd( :,7) = (/ .726_jprb, .605_jprb, .782_jprb /)
449 asy_dd( :,8) = (/ 0.731_jprb, 0.591_jprb, 0.782_jprb /)
450 asy_dd( :,9) = (/ 0.741_jprb, 0.606_jprb, 0.756_jprb /)
451 asy_dd(:,10) = (/ .740_jprb, .617_jprb, .738_jprb /)
452 asy_dd(:,11) = (/ .732_jprb, .698_jprb, .687_jprb /)
453 asy_dd(:,12) = (/ 0.725_jprb, 0.722_jprb, 0.679_jprb /)
454 asy_dd(:,13) = (/ .723_jprb, .722_jprb, .669_jprb /)
455 asy_dd(:,14) = (/ 0.645_jprb, 0.755_jprb, 0.682_jprb /)
456 asy_dd(:,15) = (/ 0.556_jprb, 0.716_jprb, 0.718_jprb /)
457 asy_dd(:,16) = (/ 0.392_jprb, 0.652_jprb, 0.730_jprb /)
458 
459 omg_dd( :,1) = (/ .974_jprb, .944_jprb, .894_jprb /)
460 omg_dd( :,2) = (/ .974_jprb, .944_jprb, .894_jprb /)
461 omg_dd( :,3) = (/ .974_jprb, .944_jprb, .894_jprb /)
462 omg_dd( :,4) = (/ .978_jprb, .947_jprb, .900_jprb /)
463 omg_dd( :,5) = (/ 0.980_jprb, 0.947_jprb, 0.905_jprb /)
464 omg_dd( :,6) = (/ .985_jprb, .957_jprb, .926_jprb /)
465 omg_dd( :,7) = (/ .988_jprb, .963_jprb, .938_jprb /)
466 omg_dd( :,8) = (/ 0.987_jprb, 0.960_jprb, 0.933_jprb /)
467 omg_dd( :,9) = (/ 0.990_jprb, 0.973_jprb, 0.946_jprb /)
468 omg_dd(:,10) = (/ .996_jprb, .990_jprb, .978_jprb /)
469 omg_dd(:,11) = (/ .997_jprb, .994_jprb, .982_jprb /)
470 omg_dd(:,12) = (/ 0.992_jprb, 0.988_jprb, 0.958_jprb /)
471 omg_dd(:,13) = (/ .997_jprb, .996_jprb, .984_jprb /)
472 omg_dd(:,14) = (/ 0.996_jprb, 0.997_jprb, 0.989_jprb /)
473 omg_dd(:,15) = (/ 0.995_jprb, 0.997_jprb, 0.994_jprb /)
474 omg_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
481 alf_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 /)
484 alf_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 /)
487 alf_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 /)
490 alf_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 /)
493 alf_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 
497 alf_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 /)
500 alf_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 
504 alf_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 /)
507 alf_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 
511 alf_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 /)
514 alf_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 /)
517 alf_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 
521 alf_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 /)
524 alf_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 /)
527 alf_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 /)
530 alf_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
540 asy_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 /)
543 asy_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 /)
546 asy_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 /)
549 asy_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 /)
552 asy_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 
556 asy_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 /)
559 asy_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 
563 asy_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 /)
566 asy_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 
570 asy_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 /)
573 asy_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 /)
576 asy_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 
580 asy_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 /)
583 asy_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 /)
586 asy_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 /)
589 asy_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
594 omg_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 /)
597 omg_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 /)
600 omg_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 /)
603 omg_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 /)
606 omg_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 
610 omg_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 /)
613 omg_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 
617 omg_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 /)
620 omg_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 
624 omg_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 /)
627 omg_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 /)
630 omg_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 
634 omg_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 /)
637 omg_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 /)
640 omg_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 /)
643 omg_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
648 alf_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 /)
651 alf_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 /)
654 alf_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 /)
657 alf_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 /)
660 alf_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 
664 alf_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 /)
667 alf_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 
671 alf_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 /)
674 alf_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 
678 alf_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 /)
681 alf_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 /)
684 alf_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 
688 alf_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 /)
691 alf_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 /)
694 alf_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 /)
697 alf_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
702 asy_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 /)
705 asy_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 /)
708 asy_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 /)
711 asy_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 /)
714 asy_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 
718 asy_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 /)
721 asy_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 
725 asy_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 /)
728 asy_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 
732 asy_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 /)
735 asy_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 /)
738 asy_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 
742 asy_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 /)
745 asy_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 /)
748 asy_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 /)
751 asy_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
756 omg_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 /)
759 omg_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 /)
762 omg_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 /)
765 omg_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 /)
768 omg_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 
772 omg_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 /)
775 omg_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 
779 omg_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 /)
782 omg_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 
786 omg_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 /)
789 omg_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 /)
792 omg_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 
796 omg_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 /)
799 omg_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 /)
802 omg_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 /)
805 omg_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
810 alf_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 /)
813 alf_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 /)
816 alf_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 /)
819 alf_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 /)
822 alf_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 
826 alf_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 /)
829 alf_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 
833 alf_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 /)
836 alf_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 
840 alf_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 /)
843 alf_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 /)
846 alf_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 
850 alf_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 /)
853 alf_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 /)
856 alf_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 /)
859 alf_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
864 asy_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 /)
867 asy_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 /)
870 asy_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 /)
873 asy_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 /)
876 asy_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 /)
879 asy_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 /)
882 asy_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 
886 asy_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 /)
889 asy_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 
893 asy_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 /)
896 asy_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 /)
899 asy_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 
903 asy_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 /)
906 asy_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 /)
909 asy_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 /)
912 asy_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
917 omg_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 /)
920 omg_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 /)
923 omg_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 /)
926 omg_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 /)
929 omg_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 
933 omg_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 /)
936 omg_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 
940 omg_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 /)
943 omg_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 
947 omg_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 /)
950 omg_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 /)
953 omg_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 
957 omg_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 /)
960 omg_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 /)
963 omg_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 /)
966 omg_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 ! ----------------------------------------------------------------
971 IF (lhook) CALL dr_hook('SU_AEROP',1,zhook_handle)
972 END SUBROUTINE su_aerop
yoeaerop::alf_om
real(kind=jprb), dimension(12, 16) alf_om
Definition: yoeaerop.F90:16
yoeaerop::asy_su
real(kind=jprb), dimension(12, 16) asy_su
Definition: yoeaerop.F90:18
yoeaerop::asy_dd
real(kind=jprb), dimension(3, 16) asy_dd
Definition: yoeaerop.F90:14
yoeaerop::alf_fa
real(kind=jprb), dimension(16) alf_fa
Definition: yoeaerop.F90:15
yoeaerop::asy_fa
real(kind=jprb), dimension(16) asy_fa
Definition: yoeaerop.F90:15
yoeaerop::alf_ss
real(kind=jprb), dimension(12, 16, 3) alf_ss
Definition: yoeaerop.F90:17
yoeaerop::omg_om
real(kind=jprb), dimension(12, 16) omg_om
Definition: yoeaerop.F90:16
parkind1::jprb
integer, parameter jprb
Definition: parkind1.F90:31
yoeaerop::omg_dd
real(kind=jprb), dimension(3, 16) omg_dd
Definition: yoeaerop.F90:14
yoeaerop::omg_fa
real(kind=jprb), dimension(16) omg_fa
Definition: yoeaerop.F90:15
yoeaerop::asy_ss
real(kind=jprb), dimension(12, 16, 3) asy_ss
Definition: yoeaerop.F90:17
yoeaerop::alf_su
real(kind=jprb), dimension(12, 16) alf_su
Definition: yoeaerop.F90:18
yomhook::lhook
logical lhook
Definition: yomhook.F90:12
yoeaerop::asy_om
real(kind=jprb), dimension(12, 16) asy_om
Definition: yoeaerop.F90:16
yoeaerop::alf_bc
real(kind=jprb), dimension(16) alf_bc
Definition: yoeaerop.F90:13
yomhook::dr_hook
subroutine dr_hook(CDNAME, KSWITCH, PKEY)
Definition: yomhook.F90:17
yoeaerop::omg_bc
real(kind=jprb), dimension(16) omg_bc
Definition: yoeaerop.F90:13
yoeaerop::asy_bc
real(kind=jprb), dimension(16) asy_bc
Definition: yoeaerop.F90:13
yoeaerop::omg_su
real(kind=jprb), dimension(12, 16) omg_su
Definition: yoeaerop.F90:18
su_aerop
subroutine su_aerop
Definition: su_aerop.F90:2
yoeaerop::alf_dd
real(kind=jprb), dimension(3, 16) alf_dd
Definition: yoeaerop.F90:14
yoeaerop::omg_ss
real(kind=jprb), dimension(12, 16, 3) omg_ss
Definition: yoeaerop.F90:17