docs/depolarization/depolarization.rst

changeset 68
f697817dad5f
parent 67
ae93ea881019
child 70
31295554bf20
equal deleted inserted replaced
67:ae93ea881019 68:f697817dad5f
1 **Single Calculus Chain ** 1 1. Particle Linear Depolarization Ratio Implementation
2 2 ======================================================
3 **version: 4.0** 3
4 4 The most important improvement included in the SCC v4.0 is the implementation of a new optical product which is the particle linear depolarization ratio.
5 **date: Date (fixed)** 5
6 6 1.1 Background
7 **DRAFT**
8
9 This document describes the main changes implemented in the SCC v4.0
10 with respect to what already provided in the SCC v3.11. It will be also
11 reported the modifications the users need to perform to run the new
12 version of SCC.
13
14 Table of Contents
15
16 1. Particle Linear Depolarization Ratio Implementation 3
17
18 1.1 Background 3
19
20 1.2 Polarization calibration 4
21
22 1.3 SCC procedure to calculate the PLDRP 4
23
24 2.Changes of the SCC input format 8
25
26 3.Real Example 10
27
28 3.1 Modification of polarization channel parameters 10
29
30 3.2 Definition of new calibration configuration and product 12
31
32 3.3 Definition of “Raman/Elastic backscatter and linear depolarization
33 ratio” 16
34
35 Particle Linear Depolarization Ratio Implementation
36 ===================================================
37
38 The most important improvement included in the SCC v4.0 is the
39 implementation of a new optical product which is the particle linear
40 depolarization ratio.
41
42 **Background**
43 -------------- 7 --------------
44 8
45 The calculation of the volume linear depolarization ratio profile 9 The calculation of the volume linear depolarization ratio profile (*VLDR*) and particle linear depolarization ratio profile (*PLDR*) needs two different steps:
46 (*VLDR*) and particle linear depolarization ratio profile (*PLDR*) needs 10
47 two different steps: 11 #. the calibration of the polarization sensitive lidar channels;
48 12 #. the calculation of the *VLDR* or *PLDR* itself.
49 1. the calibration of the polarization sensitive lidar channels; 13
50 14 The SCC allows the user to make both the above points. In particular the calibration step is made by a completely new module called **scc\_calibrator** which computes the *apparent calibration factor* :math:`\beta^*` out of the pre-processed data provided by the standard **ELPP** (Earlinet Lidar Pre-Processor) module and it records it in the SCC database (SCC\_DB). Once logged into the SCC\_DB this factor can be used whenever it is necessary.
51 2. the calculation of the *VLDR* or *PLDR* itself. 15
52 16 The raw lidar calibration measurements should be put in a NetCDF file which has the same structure as the “standard” raw SCC NetCDF input file (for more details see sections 2 and 3.2).
53 The SCC allows the user to make both the above points. In particular the 17
54 calibration step is made by a completely new module called 18 New signal types have been introduced to take into account special channel configurations used for calibration purposes.
55 **scc\_calibrator** which computes the *apparent calibration factor* 19
56 h\ :sup:`\*` out of the pre-processed data provided by the standard 20 Moreover new product types for both calibration and *PLDR* calculation have been defined. As, in principle, it is possible to calculate the *PLDR* only when the aerosol backscatter coefficient profile is available the following new products have been defined:
57 **ELPP** (Earlinet Lidar Pre-Processor) module and it records it in the 21
58 SCC database (SCC\_DB). Once logged into the SCC\_DB this factor can be 22 #. *Linear polarization calibration (factor* h) *(product\_type\_id=6);*
59 used whenever it is necessary. 23 #. *Raman backscatter and linear depolarization ratio
60
61 The raw lidar calibration measurements should be put in a NetCDF file
62 which has the same structure as the “standard” raw SCC NetCDF input file
63 (for more details see sections 2 and 3.2).
64
65 New signal types have been introduced to take into account special
66 channel configurations used for calibration purposes.
67
68 Moreover new product types for both calibration and *PLDR* calculation
69 have been defined. As, in principle, it is possible to calculate the
70 *PLDR* only when the aerosol backscatter coefficient profile is
71 available the following new products have been defined:
72
73 1. *Linear polarization calibration (factor* h) *(product\_type\_id=6);*
74
75 2. *Raman backscatter and linear depolarization ratio
76 (product\_type\_id=7);* 24 (product\_type\_id=7);*
77 25 #. *Elastic backscatter and linear depolarization ratio
78 3. *Elastic backscatter and linear depolarization ratio
79 (product\_type\_id=8).* 26 (product\_type\_id=8).*
80 27
81 The first product in the above list is used only for calibration while 28 The first product in the above list is used only for calibration while the other two are used for the calculation of *PLDR*. Basically, in most of the cases, the products 2 and 3 are equivalent to the corresponding backscatter product types with the exception that also the following new variables are available:
82 the other two are used for the calculation of *PLDR*. Basically, in most 29
83 of the cases, the products 2 and 3 are equivalent to the corresponding 30 .. code-block:: python
84 backscatter product types with the exception that also the following new 31
85 variables are available: 32 double VolumeDepol(Length) ;
86 33 double ErrorVolumeDepol(Length) ;
87 double VolumeDepol(Length) ; 34 ErrorVolumeDepol:long\_name = "absolute error of VolumeDepol" ;
88 35 double ParticleDepol(Length) ;
89 double ErrorVolumeDepol(Length) ; 36 double ErrorParticleDepol(Length) ;
90 37 ErrorParticleDepol:long\_name = "absolute error of ParticleDepol" ;
91 ErrorVolumeDepol:long\_name = "absolute error of VolumeDepol" ; 38
92 39 1.2 Polarization calibration
93 double ParticleDepol(Length) ;
94
95 double ErrorParticleDepol(Length) ;
96
97 ErrorParticleDepol:long\_name = "absolute error of ParticleDepol" ;
98
99 **Polarization calibration**
100 ---------------------------- 40 ----------------------------
101 41
102 An important point is the definition of reliable *PLDR* calibration 42 An important point is the definition of reliable *PLDR* calibration procedures. Within EARLINET the following calibration procedures are currently used:
103 procedures. Within EARLINET the following calibration procedures are 43
104 currently used: 44 a) Rayleigh calibration;
105 45 b) +45 calibration method, or D90 calibration method (made by +45 and -45 measurements);
106 a) Rayleigh calibration; 46 c) 3 signals (total, cross and parallel).
107 47
108 b) +45 calibration method, or D90 calibration method (made by +45 and 48 It is well known that method a) could produce easily large errors on *PLDR* which cannot be controlled. For this reason only the methods b) and c) can be used to provide reliable polarization calibrations and so only those methods will be implemented in the SCC.
109 -45 measurements); 49
110 50 For what it concerns the method c) it, basically, requires to solve the equation:
111 c) 3 signals (total, cross and parallel). 51
112 52 .. math::
113 It is well known that method a) could produce easily large errors on 53 \alpha_s P_s + \alpha_p P_p = P
114 *PLDR* which cannot be controlled. For this reason only the methods b) 54
115 and c) can be used to provide reliable polarization calibrations and so 55 in two different of atmospheric layers with considerably different *VLDR*. So to calibrate in this way the implementation of automatic layer identification in the SCC is required. As at moment this feature is not yet available within the SCC *ONLY* the method b) is considered.
116 only those methods will be implemented in the SCC. 56
117 57 1.3 SCC procedure to calculate the PLDRP
118 For what it concerns the method c) it, basically, requires to solve the 58 ----------------------------------------
119 equation: 59
120 60 According to what mentioned before the SCC calculates the *PLDR* through the following steps:
121 in two different of atmospheric layers with considerably different 61
122 *VLDR*. So to calibrate in this way the implementation of automatic 62 #. The user needs to create a new system configuration in the SCC\_DB including only lidar channels used for the calibration. One (or more) *Linear polarization calibration (product\_type\_id=6)* product should be associated to this new configuration (see section 3.2 for more details);
123 layer identification in the SCC is required. As at moment this feature 63
124 is not yet available within the SCC *ONLY* the method b) is considered. 64 #. This new system configuration should contain only the polarization channels in the configuration used for the calibration (for example rotated in the polarization plane of +45 degrees). A channel in calibration measurement configuration should have a *DIFFERENT* channel ID from the channel ID corresponding to the same channel in standard measurement configuration. For example, if a system has two polarization channels which in standard measurement configuration correspond to the channel ID=1 and 2 respectively, the same physical channels under calibration measurement configuration should correspond to different channel IDs (let's say ID=3 and 4 for the +45 degrees polarization rotated channels and ID=5 and 6 for the -45 degrees polarization rotated ones in case D90 calibration method is used). Moreover, the polarization channels should be labeled correctly using the new signal types available (*+45elPT, +45elPR, -45elPT, -45elPR, +45elPTnr, +45elPTfr, +45elPRnr, +45elPRfr, -45elPTnr, -45elPTfr, -45elPRnr, -45elPRfr).* For more details see section 3.2;
125 65
126 SCC procedure to calculate the PLDRP 66 #. In SCC v4.0 the polarization channels are *NOT* labeled on the base of their polarization state (as it was done in the SCC v3.11) but *ALWAYS* as transmitted and reflected channels. So the channels that in SCC v3.11 were labeled as *elCP, elCPnr, elCPfr, elPP, elPPnr elPPfr* will be labeled in SCC v4.0 as *elPR, elPRnr elPRfr elPT, elPTnr elPTfr* where the letter *T* stands from transmitted and the letter *R* for reflected.
127 ------------------------------------ 67
128 68 :WARNING: In switching from the SCC v3.11 to SCC v4.0 the following modifications have been made on *ALL* channels of *ALL* registered configurations:
129 According to what mentioned before the SCC calculates the *PLDR* through 69 *elPP→elPR*
130 the following steps: 70
131 71 *elCP→elPT*
132 1. The user needs to create a new system configuration in the SCC\_DB 72
133 including only lidar channels used for the calibration. One (or more) 73 *elPPnr→elPRnr*
134 *Linear polarization calibration (product\_type\_id=6)* product 74
135 should be associated to this new configuration (see section 3.2 for 75 *elPPfr→ elPRfr*
136 more details); 76
137 77 *elCPnr→ elPTnr*
138 2. This new system configuration should contain only the polarization 78
139 channels in the configuration used for the calibration (for example 79 *elCPfr→ elPTfr*
140 rotated in the polarization plane of +45 degrees). A channel in 80
141 calibration measurement configuration should have a *DIFFERENT* 81 Please be sure these modifications reflect to your actual lidar setup(cross channels are transmitted and parallel channels are reflected);
142 channel ID from the channel ID corresponding to the same channel in 82
143 standard measurement configuration. For example, if a system has two 83 4. The user needs to submit a file (same format as raw SCC input file) containing the raw data for the lidar channels defined at the point 1 (see section 3.2 for more details);
144 polarization channels which in standard measurement configuration 84 #. The file at point 2 is pre-processed by **ELPP** module which applies the standard pre-processing procedures applied to “standard” lidar data;
145 correspond to the channel ID=1 and 2 respectively, the same physical 85 #. The pre-processed files are then processed by the new modules **scc\_calibrator** which calculates :math:`\eta^*` *the apparent calibration factor* and logs it into the SCC\_DB;
146 channels under calibration measurement configuration should 86 #. The user needs to create a new system configuration in the SCC\_DB (which should be different from the one used for the calibration) and associate it the new product *Raman backscatter and linear depolarization ratio product\_type\_id=7)* or *Elastic backscatter and linear depolarization ratio (product\_type\_id=8).* Alternatively the calculation of those products can be added to an already existing lidar configuration as long as it is different from the calibration one;
147 correspond to different channel IDs (let's say ID=3 and 4 for the +45 87 #. The product defined at point 5 should be linked to the product containing the polarization calibration (defined at point 1) in a way that the *apparent calibration factor* can be selected from the SCC\_DB (see section 3.3 and in particular figure 3.4);
148 degrees polarization rotated channels and ID=5 and 6 for the -45 88 #. The user needs to submit another SCC raw data file containing the “standard” measurements;
149 degrees polarization rotated ones in case D90 calibration method is 89 #. Finally **ELPP** and **ELDA** will produce a b-file containing backscatter coefficient profile and *PLDR*. In particular this calculation is made in two different steps: from the pre-processed lidar polarization signals, and taking into account the *apparent calibration factor* and the *calibration factor correction K* (defined as option of *Linear polarization calibration* product\ *)* written into the SCC\_DB, an “apparent” *VLDR* :math:`\delta^*` is calculated. Even if :math:`\delta^*` is a calibrated quantity it can be still affected by possible systematic errors due to not perfect optics or alignment of the system;
150 used). Moreover, the polarization channels should be labeled 90
151 correctly using the new signal types available (*+45elPT, +45elPR, 91 #. To take into account these errors a corrected *VLDR* (:math:`\delta`) is calculated using the *polarization cross-talk correction parameters* *G* and *H* calculated on the base of Müller matrix formalism. These cross-talk correction parameters (*G* and *H*) are stored in the SCC\_DB for each lidar channels (see section 3.1 in particular figure 3.2). Finally the *PLDR* is calculated using the backscatter coefficient profile and the molecular LDRP calculated by ELPP considering the center wavelength and bandwidth of the channels interference filter.
152 -45elPT, -45elPR, +45elPTnr, +45elPTfr, +45elPRnr, +45elPRfr, 92
153 -45elPTnr, -45elPTfr, -45elPRnr, -45elPRfr).* For more details see 93 The *apparent calibration factor* :math:`\eta^*` is calculated by the **scc\_calibrator** module as the geometrical mean of the ratio of the +/-45 degrees reflected to the +/- 45 degrees transmitted signals within an altitude calibration range defined by the users in the raw data input files.
154 section 3.2; 94
155 95 In case of +45 calibration method :math:`\eta^*` is calculated by:
156 3. In SCC v4.0 the polarization channels are *NOT* labeled on the base 96
157 of their polarization state (as it was done in the SCC v3.11) but 97 .. math::
158 *ALWAYS* as transmitted and reflected channels. So the channels that 98 \eta^* = \frac{I_R}{I_T}(+45)
159 in SCC v3.11 were labeled as *elCP, elCPnr, elCPfr, elPP, elPPnr
160 elPPfr* will be labeled in SCC v4.0 as *elPR, elPRnr elPRfr elPT,
161 elPTnr elPTfr* where the letter *T* stands from transmitted and the
162 letter *R* for reflected.
163
164 **WARNING:** In switching from the SCC v3.11 to SCC v4.0 the following
165 modifications have been made on *ALL* channels of *ALL* registered
166 configurations:
167
168 *elPP→elPR*
169
170 *elCP→elPT*
171
172 *elPPnr→elPRnr*
173
174 *elPPfr→ elPRfr*
175
176 *elCPnr→ elPTnr*
177
178 *elCPfr→ elPTfr*
179
180 Please be sure these modifications reflect to your actual lidar setup
181 (cross channels are transmitted and parallel channels are reflected);
182
183 1. The user needs to submit a file (same format as raw SCC input file)
184 containing the raw data for the lidar channels defined at the point 1
185 (see section 3.2 for more details);
186
187 2. The file at point 2 is pre-processed by **ELPP** module which applies
188 the standard pre-processing procedures applied to “standard” lidar
189 data;
190
191 3. The pre-processed files are then processed by the new modules
192 **scc\_calibrator** which calculates h\ :sup:`\*` *the apparent
193 calibration factor* and logs it into the SCC\_DB;
194
195 4. The user needs to create a new system configuration in the SCC\_DB
196 (which should be different from the one used for the calibration) and
197 associate it the new product *Raman backscatter and linear
198 depolarization ratio (product\_type\_id=7)* or *Elastic backscatter
199 and linear depolarization ratio (product\_type\_id=8).* Alternatively
200 the calculation of those products can be added to an already existing
201 lidar configuration as long as it is different from the calibration
202 one;
203
204 5. The product defined at point 5 should be linked to the product
205 containing the polarization calibration (defined at point 1) in a way
206 that the *apparent calibration factor* can be selected from the
207 SCC\_DB (see section 3.3 and in particular figure 3.4);
208
209 6. The user needs to submit another SCC raw data file containing the
210 “standard” measurements;
211
212 7. Finally **ELPP** and **ELDA** will produce a b-file containing
213 backscatter coefficient profile and *PLDR*. In particular this
214 calculation is made in two different steps: from the pre-processed
215 lidar polarization signals, and taking into account the *apparent
216 calibration factor* and the *calibration factor correction K*
217 (defined as option of *Linear polarization calibration* product\ *)*
218 written into the SCC\_DB, an “apparent” *VLDR* *d\ :sup:`\*`* is
219 calculated. Even if *d\ :sup:`\*`* is a calibrated quantity it can be
220 still affected by possible systematic errors due to not perfect
221 optics or alignment of the system;
222
223 8. To take into account these errors a corrected *VLDR* (*d)* is
224 calculated using the *polarization cross-talk correction parameters*
225 *G* and *H* calculated on the base of Müller matrix formalism. These
226 cross-talk correction parameters (*G* and *H*) are stored in the
227 SCC\_DB for each lidar channels (see section 3.1 in particular figure
228 3.2). Finally the *PLDR* is calculated using the backscatter
229 coefficient profile and the molecular LDRP calculated by ELPP
230 considering the center wavelength and bandwidth of the channels
231 interference filter.
232
233 The *apparent calibration factor* h\ :sup:`\*` is calculated by the
234 **scc\_calibrator** module as the geometrical mean of the ratio of the
235 +/-45 degrees reflected to the +/- 45 degrees transmitted signals within
236 an altitude calibration range defined by the users in the raw data input
237 files.
238
239 In case of +45 calibration method h\ :sup:`\*` is calculated by:
240
241 (1.1)
242 99
243 While in case of D90 calibration method: 100 While in case of D90 calibration method:
244 101
245 (1.2) 102 .. math::
103 \eta^* = \sqrt{\frac{I_R}{I_T}(+45) \frac{I_R}{I_T}(-45)}
246 104
247 **ELDA** module calculates the “apparent” *VLDR*: 105 **ELDA** module calculates the “apparent” *VLDR*:
248 106
249 (1.3) 107 .. math::
108 \delta^* = \frac{K}{\eta^*} \cdot \frac{I_R}{I_T}
250 109
251 the *VLDR* 110 the *VLDR*
252 111
253 (1.4) 112 .. math::
113 \delta = \frac{\delta^*(G_T + H_T)-(G_R + H_R)}{(G_R - H_R) - \delta^*(G_T - H_T)}
254 114
255 and the *PLDR* 115 and the *PLDR*
256 116
257 (1.5) 117 .. math::
118 \delta_{\alpha} = \frac{(1 + \delta_m)\delta R - (1 + \delta)\delta_m}{(1 + \delta_m)R - (1 + \delta)}
258 119
259 where: 120 where:
260 121
261 - h\ :sup:`\*` is the *apparent calibration factor* calculated by 122 - :math:`\eta^*` is the *apparent calibration factor* calculated by **scc\_calibrator**
262 **scc\_calibrator** 123
263 124 - *K* is the *calibration factor correction* defined as polarization product option
264 - *K* is the *calibration factor correction* defined as polarization 125
265 product option 126 - :math:`I_T` and :math:`I_R` are the transmitted and the reflected signals in the polarization detection set-up
266 127
267 - *I\ :sub:`T`* and I\ *:sub:`R`* are the transmitted and the reflected 128 - :math:`G_{T,R}` and :math:`H_{T,R}` are *polarization cross-talk correction parameters* for the transmitted and reflected signals used to correct for systematic errors. Both these factors are defined in the SCC\_DB for each lidar channel.
268 signals in the polarization detection set-up 129
269 130 - :math:`\delta_m` is the molecular linear depolarization ratio calculated by ELPP
270 - *G\ :sub:`T,R`* and *H\ :sub:`T,R`* are *polarization cross-talk 131
271 correction parameters* for the transmitted and reflected signals used 132 - *R* is the backscatter ratio
272 to correct for systematic errors. Both these factors are defined in 133
273 the SCC\_DB for each lidar channel. 134 Please note once again that the polarization channels are described in terms of transmitted and reflected signals. This means that according to different lidar instrumental configurations, the transmitted or the reflected channel can contain total, perpendicular or parallel polarized signals.
274 135
275 - *d\ :sub:`m`* is the molecular linear depolarization ratio calculated 136 In order to retrieve the backscatter profile the total signal must be obtained combining the transmitted and reflected polarized signals. The following formula is used:
276 by ELPP 137
277 138 .. math::
278 - *R* is the backscatter ratio 139 I_{total} \propto \frac{\eta^*}{K}H_R I_T - H_T I_R
279 140
280 Please note once again that the polarization channels are described in 141 The formulas above are general and can be adapted to all possible polarization lidar configurations selecting the right polarization cross-talk correction parameters (see Table 1.1).
281 terms of transmitted and reflected signals. This means that according to 142
282 different lidar instrumental configurations, the transmitted or the 143 Let's suppose, for example, we have the perpendicular polarized lidar signal on the transmitted channel and the parallel polarized on reflected channel. For an ideal system (no diattenuation and cross-talk) we have:
283 reflected channel can contain total, perpendicular or parallel polarized 144
284 signals. 145 .. math::
285 146 G_T=1 , \qquad H_T=-1, \qquad G_R=1, \qquad H_R=1
286 In order to retrieve the backscatter profile the total signal must be 147
287 obtained combining the transmitted and reflected polarized signals. The 148 If, on the other hands, we have the perpendicular polarized lidar signal on reflected channel and the total polarized on the transmitted for and ideal system we have:
288 following formula is used: 149
289 150 .. math::
290 (1.6) 151 G_T=1 , \qquad H_T=0, \qquad G_R=1, \qquad H_R=-1
291 152
292 The formulas above are general and can be adapted to all possible 153
293 polarization lidar configurations selecting the right polarization 154 **Table 1.1:** Polarization cross-talk correction parameters for ideal systems
294 cross-talk correction parameters (see Table 1.1).
295
296 Let's suppose, for example, we have the perpendicular polarized lidar
297 signal on the transmitted channel and the parallel polarized on
298 reflected channel. For an ideal system (no diattenuation and cross-talk)
299 we have:
300
301 If, on the other hands, we have the perpendicular polarized lidar signal
302 on reflected channel and the total polarized on the transmitted for and
303 ideal system we have:
304
305 **Table 1.1:** Polarization cross-talk correction parameters for ideal
306 systems
307 155
308 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 156 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
309 | Laser polarization | Detected in lidar channel | 157 | Laser polarization | Detected in lidar channel |
310 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 158 + +-----------------------------+-----------------+-----------------+-----------------+
311 | | Transmitted | Reflected | 159 | | Transmitted | Reflected |
312 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 160 + +-----------------------------+-----------------+-----------------+-----------------+
313 | | *G\ :sub:`T`* | *H\ :sub:`T`* | *G\ :sub:`R`* | *H\ :sub:`R`* | 161 | | :math:`G_T` | :math:`H_T` | :math:`G_R` | :math:`H_R` |
314 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 162 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
315 | total | 1 | 0 | 1 | 0 | 163 | total | 1 | 0 | 1 | 0 |
316 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 164 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
317 | parallel | 1 | 1 | 1 | 1 | 165 | parallel | 1 | 1 | 1 | 1 |
318 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 166 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
319 | cross | 1 | -1 | 1 | -1 | 167 | cross | 1 | -1 | 1 | -1 |
320 +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 168 +----------------------+-----------------------------+-----------------+-----------------+-----------------+
321 169
322 The *apparent calibration factor* (h:sup:`\*`), *the calibration factor 170 The *apparent calibration factor (:math:`\eta^*`), *the calibration factor correction* (*K*) and the *polarization cross-talk correction parameters* are stored by **ELPP** module in the intermediate NetCDF files using the following variables:
323 correction* (*K*) and the *polarization cross-talk correction 171
324 parameters* are stored by **ELPP** module in the intermediate NetCDF 172 testing the inline code :code:`test`
325 files using the following variables: 173
326 174 - :code:`Polarization_Channel_Gain_Factor` (*apparent calibration factor* - :math:`\eta^*` )
327 - *Polarization\_Channel\_Gain\_Factor (apparent calibration factor* – 175 - :code:`Polarization_Channel_Gain_Factor_Correction` (*calib. factor corr.* – *K*)
328 h\ :sup:`\*`) 176 - :code:`G_T`
329 177 - :code:`H_T`
330 - *Polarization\_Channel\_Gain\_Factor\_Correction (calib. factor 178 - :code:`G_R`
331 corr.* – *K*) 179 - :code:`H_R`
332 180
333 - *G\_T* 181 Finally new usecases have been defined to take into account all the possible lidar configurations. The details on that are provided as a separate file.
334 182
335 - *H\_T* 183 2. Changes of the SCC input format
336 184 ==================================
337 - *G\_R*
338
339 - *H\_R*
340
341 Finally new usecases have been defined to take into account all the
342 possible lidar configurations. The details on that are provided as a
343 separate file.
344
345 Changes of the SCC input format
346 ===============================
347 185
348 The following minor changes have been applied to raw SCC data format: 186 The following minor changes have been applied to raw SCC data format:
349 187
350 1. The optional variable *ID\_Range* has been *REMOVED*; 188 #. The optional variable *ID\_Range* has been *REMOVED*;
351 189 #. The *OPTIONAL* variable :code:`int Signal\_Type(channels)` has been added. The possible values are the same available in the SCC\_DB:
352 2. The *OPTIONAL* variable *int Signal\_Type(channels)* has been added. 190
353 The possible values are the same available in the SCC\_DB: 191 :code:`0` :math:`\rightarrow` :code:`elT`
354 192
355 0 *→* elT 193 :code:`1` :math:`\rightarrow` :code:`elTnr`
356 194
357 1 *→* elTnr 195 :code:`2` :math:`\rightarrow` :code:`elTfr`
358 196
359 2 *→* elTfr 197 :code:`3` :math:`\rightarrow` :code:`vrRN2`
360 198
361 3 *→* vrRN2 199 :code:`4` :math:`\rightarrow` :code:`vrRN2nr`
362 200
363 4 *→* vrRN2nr 201 :code:`5` :math:`\rightarrow` :code:`vrRN2fr`
364 202
365 5 *→* vrRN2fr 203 :code:`6` :math:`\rightarrow` :code:`elPR`
366 204
367 6 *→* elPR 205 :code:`7` :math:`\rightarrow` :code:`elPT`
368 206
369 7 *→* elPT 207 :code:`8` :math:`\rightarrow` :code:`pRRlow`
370 208
371 8 *→* pRRlow 209 :code:`9` :math:`\rightarrow` :code:`pRRhigh`
372 210
373 9 *→* pRRhigh 211 :code:`10` :math:`\rightarrow` :code:`elPRnr`
374 212
375 10 *→* elPRnr 213 :code:`11` :math:`\rightarrow` :code:`elPRfr`
376 214
377 11 *→* elPRfr 215 :code:`12` :math:`\rightarrow` :code:`elPTnr`
378 216
379 12 *→* elPTnr 217 :code:`13` :math:`\rightarrow` :code:`elPTfr`
380 218
381 13 *→* elPTfr 219 :code:`14` :math:`\rightarrow` :code:`vrRH2O`
382 220
383 14 *→* vrRH2O 221 :code:`15` :math:`\rightarrow` :code:`pRRhighnr`
384 222
385 15 *→* pRRhighnr 223 :code:`16` :math:`\rightarrow` :code:`pRRhighfr`
386 224
387 16 *→* pRRhighfr 225 :code:`17` :math:`\rightarrow` :code:`pRRlownr`
388 226
389 17 *→* pRRlownr 227 :code:`18` :math:`\rightarrow` :code:`pRRlowfr`
390 228
391 18 *→* pRRlowfr 229 :code:`19` :math:`\rightarrow` :code:`vrRH2Onr`
392 230
393 19 *→* vrRH2Onr 231 :code:`20` :math:`\rightarrow` :code:`vrRH2Ofr`
394 232
395 20 *→* vrRH2Ofr 233 :code:`21` :math:`\rightarrow` :code:`elTunr`
396 234
397 21 *→* elTunr 235 :code:`22` :math:`\rightarrow` :code:`+45elPT`
398 236
399 *22 → +45elPT* 237 :code:`23` :math:`\rightarrow` :code:`+45elPR`
400 238
401 *23 → +45elPR* 239 :code:`24` :math:`\rightarrow` :code:`-45elPT`
402 240
403 *24 → -45elPT* 241 :code:`25` :math:`\rightarrow` :code:`-45elPR`
404 242
405 *25 → -45elPR* 243 :code:`26` :math:`\rightarrow` :code:`+45elPTnr`
406 244
407 *26 → +45elPTnr* 245 :code:`27` :math:`\rightarrow` :code:`+45elPTfr`
408 246
409 *27 → +45elPTfr* 247 :code:`28` :math:`\rightarrow` :code:`+45elPRnr`
410 248
411 *28 → +45elPRnr* 249 :code:`29` :math:`\rightarrow` :code:`+45elPRfr`
412 250
413 *29 → +45elPRfr* 251 :code:`30` :math:`\rightarrow` :code:`-45elPTnr`
414 252
415 *30 → -45elPTnr* 253 :code:`31` :math:`\rightarrow` :code:`-45elPTfr`
416 254
417 *31 → -45elPTfr* 255 :code:`32` :math:`\rightarrow` :code:`-45elPRnr`
418 256
419 *32 → -45elPRnr* 257 :code:`33` :math:`\rightarrow` :code:`-45elPRfr`
420 258
421 *33 → -45elPRfr* 259 :WARNING: It this variable is found in the SCC input file the corresponding settings in the SCC database will be *overwritten*. Unless you don't have any valid reason to overwrite the database value this variable should not be used.
422 260
423 **WARNING:** It this variable is found in the SCC input file the 261 3. The variables:
424 corresponding settings in the SCC database will be *overwritten*. Unless 262
425 you don't have any valid reason to overwrite the database value this 263 .. code-block:: python
426 variable should not be used. 264
427 265 double Pol\_Calib\_Range\_Min(channels)
428 1. The variables: 266 double Pol\_Calib\_Range\_Max(channels)
429 267
430 *double Pol\_Calib\_Range\_Min(channels)* 268 have been added. Both these variable are *MANDATORY* for any calibration raw dataset. These variable should be included only the polarization calibration measurements and should specify the altitude range (meters) in which the polarization calibration should be made. For more details see section 3.3;
431 269
432 *double Pol\_Calib\_Range\_Max(channels) * 270 4. The variable :code:`Depolarization_Factor` has been *REMOVED*.
433 271
434 have been added. Both these variable are *mandatory* for any calibration 272 The SCC v3.11 used this variable to get polarization calibration factor for the calculation of the total signal out of cross and parallels ones. As the SCC v4.0 is able to calculate the same parameter by itself, the use of this variable is *NOT* possible anymore. The recommended way to get a valid and quality assured depolarization calibration factor is to submit to the SCC v4.0 a polarization calibration dataset and let the SCC to calculate such factor.
435 raw dataset. 273
436 274 To make this change more smooth and to provide the users with the possibility to continue to analyze their data with the SCC v4.0 even if a calibration dataset has not been submitted yet, it will be possible for a *LIMITED* period of time to submit the calibration constant via the SCC web interface. The SCC will keep track of the used calibration method (automatic or manual).
437 These variable should be included only the polarization calibration 275
438 measurements and should specify the altitude range (meters) in which the 276 :WARNING: After this transition period *ONLY* automatic calibration will be allowed!
439 polarization calibration should be made. For more details see section 277
440 3.3; 278 5. The new *OPTIONAL* variable:
441 279
442 1. The variable *Depolarization\_Factor* has been *removed*. 280 :code:`string channel\_string\_ID(channels)`
443 281
444 The SCC v3.11 used this variable to get polarization calibration factor 282 has been introduced.
445 for the calculation of the total signal out of cross and parallels ones. 283
446 As the SCC v4.0 is able to calculate the same parameter by itself, the 284 Starting from SCC v4.0 the lidar channel can be identified not only by using integers (as it happened until SCC v3.11) but also by using strings.
447 use of this variable is *NOT* possible anymore. The recommended way to 285
448 get a valid and quality assured depolarization calibration factor is to 286 The procedure implemented in the SCC v4.0 to recognize the lidar channel within the raw lidar data is fully backward compatible (old format files are accepted as they are by SCC v4.0).
449 submit to the SCC v4.0 a polarization calibration dataset and let the 287
450 SCC to calculate such factor. 288 :WARNING: Please note that the definition of the new string variable requires netCDF-4 format! The type *string* is not supported in netCDF-3 format!
451
452 To make this change more smooth and to provide the users with the
453 possibility to continue to analyze their data with the SCC v4.0 even if
454 a calibration dataset has not been submitted yet, it will be possible
455 for a *LIMITED* period of time to submit the calibration constant via
456 the SCC web interface. The SCC will keep track of the used calibration
457 method (automatic or manual).
458
459 **WARNING:** After this transition period *only* automatic calibration
460 will be allowed!
461
462 1. The new *optional* variable:
463
464 *string channel\_string\_ID(channels)*
465
466 has been introduced.
467
468 Starting from SCC v4.0 the lidar channel can be identified not only by
469 using integers (as it happened until SCC v3.11) but also by using
470 strings.
471
472 The procedure implemented in the SCC v4.0 to recognize the lidar channel
473 within the raw lidar data is fully backward compatible (old format files
474 are accepted as they are by SCC v4.0).
475
476 **WARNING:** Please note that the definition of the new string variable
477 requires netCDF-4 format! The type *string* is not supported in netCDF-3
478 format!
479 289
480 Real Example 290 Real Example
481 ============ 291 ============
482 292
483 This section describes all the practical steps the users need to follow 293 This section describes all the practical steps the users need to follow

mercurial