GCC Code Coverage Report
Directory: ./ Exec Total Coverage
File: phylmd/rrtm/su_aerop.F90 Lines: 298 298 100.0 %
Date: 2023-06-30 12:51:15 Branches: 578 580 99.7 %

Line Branch Exec Source
1
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
1
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
1
ALF_FA(:) =0.0_JPRB
81
1
ASY_FA(:) =0.0_JPRB
82
1
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
1
  &  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
1
  &   .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
1
  &   .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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &  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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
13
  &   .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
4
ALF_DD( :,1) = (/ 2.210_JPRB, 1.011_JPRB,  .398_JPRB /)
426
4
ALF_DD( :,2) = (/ 2.210_JPRB, 1.011_JPRB,  .398_JPRB /)
427
4
ALF_DD( :,3) = (/ 2.210_JPRB, 1.011_JPRB,  .398_JPRB /)
428
4
ALF_DD( :,4) = (/ 2.357_JPRB,  .985_JPRB,  .398_JPRB /)
429
4
ALF_DD( :,5) = (/ 2.453_JPRB, 0.949_JPRB, 0.404_JPRB /)
430
4
ALF_DD( :,6) = (/ 2.549_JPRB,  .901_JPRB,  .414_JPRB /)
431
4
ALF_DD( :,7) = (/ 2.605_JPRB,  .877_JPRB,  .423_JPRB /)
432
4
ALF_DD( :,8) = (/ 2.635_JPRB, 0.869_JPRB, 0.428_JPRB /)
433
4
ALF_DD( :,9) = (/ 2.596_JPRB, 0.978_JPRB, 0.428_JPRB /)
434
4
ALF_DD(:,10) = (/ 2.555_JPRB, 1.039_JPRB,  .424_JPRB /)
435
4
ALF_DD(:,11) = (/ 2.227_JPRB, 1.360_JPRB,  .404_JPRB /)
436
4
ALF_DD(:,12) = (/ 2.058_JPRB, 1.460_JPRB, 0.405_JPRB /)
437
4
ALF_DD(:,13) = (/ 2.041_JPRB, 1.474_JPRB,  .405_JPRB /)
438
4
ALF_DD(:,14) = (/ 1.117_JPRB, 1.523_JPRB, 0.502_JPRB /)
439
4
ALF_DD(:,15) = (/ 0.583_JPRB, 1.117_JPRB, 0.596_JPRB /)
440
4
ALF_DD(:,16) = (/ 0.263_JPRB, 0.705_JPRB, 0.607_JPRB /)
441
442
4
ASY_DD( :,1) = (/  .673_JPRB,  .745_JPRB,  .799_JPRB /)
443
4
ASY_DD( :,2) = (/  .673_JPRB,  .745_JPRB,  .799_JPRB /)
444
4
ASY_DD( :,3) = (/  .673_JPRB,  .745_JPRB,  .799_JPRB /)
445
4
ASY_DD( :,4) = (/  .691_JPRB,  .712_JPRB,  .792_JPRB /)
446
4
ASY_DD( :,5) = (/ 0.704_JPRB, 0.679_JPRB, 0.791_JPRB /)
447
4
ASY_DD( :,6) = (/  .715_JPRB,  .640_JPRB,  .786_JPRB /)
448
4
ASY_DD( :,7) = (/  .726_JPRB,  .605_JPRB,  .782_JPRB /)
449
4
ASY_DD( :,8) = (/ 0.731_JPRB, 0.591_JPRB, 0.782_JPRB /)
450
4
ASY_DD( :,9) = (/ 0.741_JPRB, 0.606_JPRB, 0.756_JPRB /)
451
4
ASY_DD(:,10) = (/  .740_JPRB,  .617_JPRB,  .738_JPRB /)
452
4
ASY_DD(:,11) = (/  .732_JPRB,  .698_JPRB,  .687_JPRB /)
453
4
ASY_DD(:,12) = (/ 0.725_JPRB, 0.722_JPRB, 0.679_JPRB /)
454
4
ASY_DD(:,13) = (/  .723_JPRB,  .722_JPRB,  .669_JPRB /)
455
4
ASY_DD(:,14) = (/ 0.645_JPRB, 0.755_JPRB, 0.682_JPRB /)
456
4
ASY_DD(:,15) = (/ 0.556_JPRB, 0.716_JPRB, 0.718_JPRB /)
457
4
ASY_DD(:,16) = (/ 0.392_JPRB, 0.652_JPRB, 0.730_JPRB /)
458
459
4
OMG_DD( :,1) = (/  .974_JPRB,  .944_JPRB,  .894_JPRB /)
460
4
OMG_DD( :,2) = (/  .974_JPRB,  .944_JPRB,  .894_JPRB /)
461
4
OMG_DD( :,3) = (/  .974_JPRB,  .944_JPRB,  .894_JPRB /)
462
4
OMG_DD( :,4) = (/  .978_JPRB,  .947_JPRB,  .900_JPRB /)
463
4
OMG_DD( :,5) = (/ 0.980_JPRB, 0.947_JPRB, 0.905_JPRB /)
464
4
OMG_DD( :,6) = (/  .985_JPRB,  .957_JPRB,  .926_JPRB /)
465
4
OMG_DD( :,7) = (/  .988_JPRB,  .963_JPRB,  .938_JPRB /)
466
4
OMG_DD( :,8) = (/ 0.987_JPRB, 0.960_JPRB, 0.933_JPRB /)
467
4
OMG_DD( :,9) = (/ 0.990_JPRB, 0.973_JPRB, 0.946_JPRB /)
468
4
OMG_DD(:,10) = (/  .996_JPRB,  .990_JPRB,  .978_JPRB /)
469
4
OMG_DD(:,11) = (/  .997_JPRB,  .994_JPRB,  .982_JPRB /)
470
4
OMG_DD(:,12) = (/ 0.992_JPRB, 0.988_JPRB, 0.958_JPRB /)
471
4
OMG_DD(:,13) = (/  .997_JPRB,  .996_JPRB,  .984_JPRB /)
472
4
OMG_DD(:,14) = (/ 0.996_JPRB, 0.997_JPRB, 0.989_JPRB /)
473
4
OMG_DD(:,15) = (/ 0.995_JPRB, 0.997_JPRB, 0.994_JPRB /)
474
4
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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  & 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
13
  &  .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
13
  &  .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
13
  & 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
13
  &  .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
13
  & 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
13
  & 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
13
  & 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB /)
969
970
!     ----------------------------------------------------------------
971
1
IF (LHOOK) CALL DR_HOOK('SU_AEROP',1,ZHOOK_HANDLE)
972
1
END SUBROUTINE SU_AEROP