Directory:  ./ 

File:  rad/su_aerop.f90 
Date:  20220111 19:19:34 
Exec  Total  Coverage  

Lines:  298  298  100.0% 
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  ! JEANJACQUES MORCRETTE *ECMWF*  
36  ! from O.BOUCHER's standalone codes (LOA, 200503)  
37  
38  ! MODIFICATIONS.  
39  !   
40  ! ORIGINAL : 20050325  
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/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.

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 seasalt (_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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,1) = (/ 2.210_JPRB, 1.011_JPRB, .398_JPRB /) 
426 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,2) = (/ 2.210_JPRB, 1.011_JPRB, .398_JPRB /) 
427 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,3) = (/ 2.210_JPRB, 1.011_JPRB, .398_JPRB /) 
428 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,4) = (/ 2.357_JPRB, .985_JPRB, .398_JPRB /) 
429 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,5) = (/ 2.453_JPRB, 0.949_JPRB, 0.404_JPRB /) 
430 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,6) = (/ 2.549_JPRB, .901_JPRB, .414_JPRB /) 
431 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,7) = (/ 2.605_JPRB, .877_JPRB, .423_JPRB /) 
432 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,8) = (/ 2.635_JPRB, 0.869_JPRB, 0.428_JPRB /) 
433 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD( :,9) = (/ 2.596_JPRB, 0.978_JPRB, 0.428_JPRB /) 
434 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,10) = (/ 2.555_JPRB, 1.039_JPRB, .424_JPRB /) 
435 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,11) = (/ 2.227_JPRB, 1.360_JPRB, .404_JPRB /) 
436 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,12) = (/ 2.058_JPRB, 1.460_JPRB, 0.405_JPRB /) 
437 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,13) = (/ 2.041_JPRB, 1.474_JPRB, .405_JPRB /) 
438 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,14) = (/ 1.117_JPRB, 1.523_JPRB, 0.502_JPRB /) 
439 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,15) = (/ 0.583_JPRB, 1.117_JPRB, 0.596_JPRB /) 
440 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ALF_DD(:,16) = (/ 0.263_JPRB, 0.705_JPRB, 0.607_JPRB /) 
441  
442 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,1) = (/ .673_JPRB, .745_JPRB, .799_JPRB /) 
443 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,2) = (/ .673_JPRB, .745_JPRB, .799_JPRB /) 
444 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,3) = (/ .673_JPRB, .745_JPRB, .799_JPRB /) 
445 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,4) = (/ .691_JPRB, .712_JPRB, .792_JPRB /) 
446 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,5) = (/ 0.704_JPRB, 0.679_JPRB, 0.791_JPRB /) 
447 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,6) = (/ .715_JPRB, .640_JPRB, .786_JPRB /) 
448 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,7) = (/ .726_JPRB, .605_JPRB, .782_JPRB /) 
449 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,8) = (/ 0.731_JPRB, 0.591_JPRB, 0.782_JPRB /) 
450 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD( :,9) = (/ 0.741_JPRB, 0.606_JPRB, 0.756_JPRB /) 
451 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,10) = (/ .740_JPRB, .617_JPRB, .738_JPRB /) 
452 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,11) = (/ .732_JPRB, .698_JPRB, .687_JPRB /) 
453 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,12) = (/ 0.725_JPRB, 0.722_JPRB, 0.679_JPRB /) 
454 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,13) = (/ .723_JPRB, .722_JPRB, .669_JPRB /) 
455 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,14) = (/ 0.645_JPRB, 0.755_JPRB, 0.682_JPRB /) 
456 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,15) = (/ 0.556_JPRB, 0.716_JPRB, 0.718_JPRB /) 
457 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  ASY_DD(:,16) = (/ 0.392_JPRB, 0.652_JPRB, 0.730_JPRB /) 
458  
459 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,1) = (/ .974_JPRB, .944_JPRB, .894_JPRB /) 
460 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,2) = (/ .974_JPRB, .944_JPRB, .894_JPRB /) 
461 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,3) = (/ .974_JPRB, .944_JPRB, .894_JPRB /) 
462 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,4) = (/ .978_JPRB, .947_JPRB, .900_JPRB /) 
463 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,5) = (/ 0.980_JPRB, 0.947_JPRB, 0.905_JPRB /) 
464 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,6) = (/ .985_JPRB, .957_JPRB, .926_JPRB /) 
465 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,7) = (/ .988_JPRB, .963_JPRB, .938_JPRB /) 
466 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,8) = (/ 0.987_JPRB, 0.960_JPRB, 0.933_JPRB /) 
467 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD( :,9) = (/ 0.990_JPRB, 0.973_JPRB, 0.946_JPRB /) 
468 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,10) = (/ .996_JPRB, .990_JPRB, .978_JPRB /) 
469 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,11) = (/ .997_JPRB, .994_JPRB, .982_JPRB /) 
470 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,12) = (/ 0.992_JPRB, 0.988_JPRB, 0.958_JPRB /) 
471 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,13) = (/ .997_JPRB, .996_JPRB, .984_JPRB /) 
472 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,14) = (/ 0.996_JPRB, 0.997_JPRB, 0.989_JPRB /) 
473 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,15) = (/ 0.995_JPRB, 0.997_JPRB, 0.994_JPRB /) 
474 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 3 times.

4  OMG_DD(:,16) = (/ 0.994_JPRB, 0.997_JPRB, 0.996_JPRB /) 
475  
476  
477  ! For seasalt, 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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 12 times.

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 
2/2✓ Branch 0 taken 12 times.
✓ Branch 1 taken 1 times.

13  & 0.997_JPRB, 0.997_JPRB, 0.998_JPRB, 0.998_JPRB, 0.998_JPRB, 0.999_JPRB /) 
969  
970  !   
971 
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.

1  IF (LHOOK) CALL DR_HOOK('SU_AEROP',1,ZHOOK_HANDLE) 
972  1  END SUBROUTINE SU_AEROP  
973 