--- a/docs/depolarization/depolarization.rst Mon Dec 05 13:25:30 2016 +0200 +++ b/docs/depolarization/depolarization.rst Mon Dec 05 13:26:38 2016 +0200 @@ -3,6 +3,9 @@ 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. +.. important:: + If your lidar system is not equipped with any polarization channels **NO** changes are required. In this case, the SCC v4.0 should work using the same input files and the same database configurations you have used with the SCC v3.11. Anyway as in the SCC v4.0 several bugs have been fixed, it is recommended to re-run all the measurement IDs you have submitted. For doing that you just need to reprocess all your data without the need to submit raw data files already uploaded on the server. + 1.1 Background -------------- @@ -11,7 +14,7 @@ #. the calibration of the polarization sensitive lidar channels; #. the calculation of the *VLDR* or *PLDR* itself. -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. +The SCC allows the user to make both the above points. In particular the calibration step is made by a completely new module called **ELDEC** (Earlinet Lidar Depolarization Calibrator) which computes the *apparent calibration factor* :math:`\eta^*` 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. 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). @@ -19,11 +22,11 @@ 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: -#. *Linear polarization calibration (factor* h) *(product\_type\_id=6);* +#. *Linear polarization calibration (factor* :math:`\eta`) *(product_type_id=6);* #. *Raman backscatter and linear depolarization ratio - (product\_type\_id=7);* + (product_type_id=7);* #. *Elastic backscatter and linear depolarization ratio - (product\_type\_id=8).* + (product_type_id=8).* 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: @@ -31,10 +34,10 @@ double VolumeDepol(Length) ; double ErrorVolumeDepol(Length) ; - ErrorVolumeDepol:long\_name = "absolute error of VolumeDepol" ; + ErrorVolumeDepol:long_name = "absolute error of VolumeDepol" ; double ParticleDepol(Length) ; double ErrorParticleDepol(Length) ; - ErrorParticleDepol:long\_name = "absolute error of ParticleDepol" ; + ErrorParticleDepol:long_name = "absolute error of ParticleDepol" ; 1.2 Polarization calibration ---------------------------- @@ -42,7 +45,7 @@ An important point is the definition of reliable *PLDR* calibration procedures. Within EARLINET the following calibration procedures are currently used: a) Rayleigh calibration; - b) +45 calibration method, or D90 calibration method (made by +45 and -45 measurements); + b) +45 calibration method, or :math:`\Delta90` calibration method (made by +45 and -45 measurements); c) 3 signals (total, cross and parallel). 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. @@ -52,14 +55,14 @@ .. math:: \alpha_s P_s + \alpha_p P_p = P -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. +in two different 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. 1.3 SCC procedure to calculate the PLDRP ---------------------------------------- According to what mentioned before the SCC calculates the *PLDR* through the following steps: -#. 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); +#. 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); #. 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; @@ -81,23 +84,23 @@ Please be sure these modifications reflect to your actual lidar setup(cross channels are transmitted and parallel channels are reflected); 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); -#. The file at point 2 is pre-processed by **ELPP** module which applies the standard pre-processing procedures applied to “standard” lidar data; -#. 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; -#. 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; -#. 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); +#. The file at point 4 is pre-processed by **ELPP** module which applies the standard pre-processing procedures applied to “standard” lidar data; +#. The pre-processed files are then processed by the new modules **ELDEC** which calculates :math:`\eta^*` *the apparent calibration factor* and logs it into the SCC_DB; +#. 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; +#. The product defined at point 7 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); #. The user needs to submit another SCC raw data file containing the “standard” measurements; -#. 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; +#. 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; -#. 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. +#. 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. -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. +The *apparent calibration factor* :math:`\eta^*` is calculated by the **ELDEC** 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. In case of +45 calibration method :math:`\eta^*` is calculated by: .. math:: \eta^* = \frac{I_R}{I_T}(+45) -While in case of D90 calibration method: +While in case of :math:`\Delta90` calibration method: .. math:: \eta^* = \sqrt{\frac{I_R}{I_T}(+45) \frac{I_R}{I_T}(-45)} @@ -119,13 +122,13 @@ where: - - :math:`\eta^*` is the *apparent calibration factor* calculated by **scc\_calibrator** + - :math:`\eta^*` is the *apparent calibration factor* calculated by **ELDEC** - *K* is the *calibration factor correction* defined as polarization product option - :math:`I_T` and :math:`I_R` are the transmitted and the reflected signals in the polarization detection set-up - - :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. + - :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. - :math:`\delta_m` is the molecular linear depolarization ratio calculated by ELPP @@ -145,10 +148,10 @@ .. math:: G_T=1 , \qquad H_T=-1, \qquad G_R=1, \qquad H_R=1 -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: +If, on the other hand, we have the perpendicular polarized lidar signal on reflected channel and the total polarized on the transmitted for and ideal system we have: .. math:: - G_T=1 , \qquad H_T=0, \qquad G_R=1, \qquad H_R=-1 + G_T=1 , \qquad H_T=1, \qquad G_R=1, \qquad H_R=-1 **Table 1.1:** Polarization cross-talk correction parameters for ideal systems @@ -162,9 +165,9 @@ +----------------------+-----------------------------+-----------------+-----------------+-----------------+ | total | 1 | 0 | 1 | 0 | +----------------------+-----------------------------+-----------------+-----------------+-----------------+ -| parallel | 1 | 1 | 1 | 1 | +| parallel | 1 | 1 | 1 | -1 | +----------------------+-----------------------------+-----------------+-----------------+-----------------+ -| cross | 1 | -1 | 1 | -1 | +| cross | 1 | -1 | 1 | 1 | +----------------------+-----------------------------+-----------------+-----------------+-----------------+ 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: @@ -183,8 +186,8 @@ The following minor changes have been applied to raw SCC data format: -#. The optional variable *ID\_Range* has been **REMOVED**; -#. The **OPTIONAL** variable :code:`int Signal\_Type(channels)` has been added. The possible values are the same available in the SCC\_DB: +#. The optional variable *ID_Range* has been **REMOVED**; +#. The **OPTIONAL** variable :code:`int Signal_Type(channels)` has been added. The possible values are the same available in the SCC_DB: :code:`0` :math:`\rightarrow` :code:`elT` @@ -254,14 +257,14 @@ :code:`33` :math:`\rightarrow` :code:`-45elPRfr` - .. 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. + .. warning:: 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. 3. The variables: :: - double Pol\_Calib\_Range\_Min(channels) - double Pol\_Calib\_Range\_Max(channels) + double Pol_Calib_Range_Min(channels) + double Pol_Calib_Range_Max(channels) 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; @@ -275,7 +278,7 @@ 5. The new **OPTIONAL** variable: - :code:`string channel\_string\_ID(channels)` + :code:`string channel_string_ID(channels)` has been introduced. @@ -293,7 +296,7 @@ :IMPORTANT: If your lidar system is not equipped with any polarization channels **NO** changes are required. In this case, the SCC v4.0 should work using the same input files and the same database configurations you have used with the SCC v3.11. Anyway as in the SCC v4.0 several bugs have been fixed,it is recommended to re-run all the measurement IDs you have submitted. For doing that you just need to reprocess all your data without the need to submit raw data files already uploaded on the server. -The practical example reported below describes the modifications required to use the SCC v4.0 for lidar systems equipped with polarization channels. +The practical example reported below describes the modifications required to use the SCC v4.0 for lidar systems equipped with polarization channels. Lidar systems not equipped with polarization channels do not require any modification to switch to SCC v4.0. 3.1 Modification of polarization channel parameters --------------------------------------------------- @@ -333,7 +336,7 @@ **Figure 3.1**: How to select signal types -The first modification concerns the settings of the channel type for the 532 cross and 532 parallel polarization channels. Starting from SCC v4.0 polarization channels are identified as transmitted and reflected polarization channels and not on the base of their polarization state. So suppose if we suppose the cross polarized channel is transmitted by a polarizer beam splitter cube and the parallel is reflected the value reported in table 3.1 should be modified as they appear in table 3.2. So using the SCC web interface, the signal type of the 532 cross channel should be changed from :code:`elCP` to :code:`elPT` and in the same way the 532 parallel channel should be changed from :code:`elPP` to :code:`elPR` (see figure 3.1). +The first modification concerns the settings of the channel type for the 532 cross and 532 parallel polarization channels. Starting from SCC v4.0 polarization channels are identified as transmitted and reflected polarization channels and not on the base of their polarization state. So if we suppose that the cross polarized channel is transmitted by a polarizer beam splitter cube, and the parallel is reflected, the value reported in table 3.1 should be modified as they appear in table 3.2. So using the SCC web interface, the signal type of the 532 cross channel should be changed from :code:`elCP` to :code:`elPT` and in the same way the 532 parallel channel should be changed from :code:`elPP` to :code:`elPR` (see figure 3.1). **Table 3.2:** The same of table 3.1 but with new channel types introduced in SCC v4.0 @@ -369,11 +372,10 @@ In this section we will see how to set the polarization calibration parameters: the calibration constant (called :math:`\eta^*` in section 1.3) and the correction to calibration constant (called K in section 1.3). In order to provide such parameters you need to define a new system configuration to be used **ONLY** for calibration purposes. Such new configuration should include the polarization channels in the measurement configuration used for the calibration. Let's suppose we want to use the :math:`\Delta90` calibration method. -In this case we need to define a new configuration (called for example “depol_calibration”) as reported in the table 3.3. As you can see the configuration “depol\_calibration” includes 4 “new” channels. Actually the channels “532 cross +45 degrees” (channel ID=10) and “532 cross -45 degrees” (channel ID=12) refer to the same physical channel “532 cross” reported with channel ID=3 in table 3.2. Anyway we need to define two new channel IDs to identify the “532 cross” channel in the two polarization rotated configurations (+45 and -45 degrees) needed to apply the D90 calibration method. The same is true for the “532 parallel” channel. The polarization rotated channels should be labeled with the corresponding signal type as reported in table 3.3 (see figure +In this case we need to define a new configuration (called for example “depol_calibration”) as reported in the table 3.3. As you can see the configuration “depol_calibration” includes 4 “new” channels. Actually the channels “532 cross +45 degrees” (channel ID=10) and “532 cross -45 degrees” (channel ID=12) refer to the same physical channel “532 cross” reported with channel ID=3 in table 3.2. Anyway we need to define two new channel IDs to identify the “532 cross” channel in the two polarization rotated configurations (+45 and -45 degrees) needed to apply the D90 calibration method. The same is true for the “532 parallel” channel. The polarization rotated channels should be labeled with the corresponding signal type as reported in table 3.3 (see figure 3.1). -**Table 3.3:** Polarization calibration configurations assuming D90 -calibration method +**Table 3.3:** Polarization calibration configurations assuming :math:`\Delta90` calibration method +----------------------------+--------------+----------------+----------------------+ | Channel Name | Channel ID | Channel Type | depol_calibration | @@ -401,94 +403,102 @@ As you can see it is possible to fill in only the K correction factor and not the calibration constant :math:`\eta^*`. -Actually for a **LIMITED** period of time it will be possible to fill in also the constant :math:`\eta^*` using a temporary tab called *Polarization calibration constant*. This has been done to provide the users with the possibility to continue to use the SCC even if an automatic calibration made by the SCC was not submitted yet. Anyway after a transition period it will be **NOT** possible to provide calibration constant using this procedure and the parameter :math:`\eta^*` can be calculated **ONLY** by the SCC as result of the submission of a proper calibration raw input dataset. The format of this input file is the same as the standard SCC input file. The only difference is that is should contain calibration measurements instead of standard measurements. Following our example, such file should contain the measurement performed at +45 and -45 degrees at 532nm. Also the channel IDs in the file should reflect the ones reported in table 3.3. +Actually for a **LIMITED** period of time it will be possible to fill in also the constant :math:`\eta^*` using a temporary option shown in figure 3.4. This has been done to provide the users with the possibility to continue to use the SCC even if an automatic calibration made by the SCC was not submitted yet. Anyway after a transition period it will be **NOT** possible to provide calibration constant using this procedure and the parameter :math:`\eta^*` can be calculated **ONLY** by the SCC as result of the submission of a proper calibration raw input dataset. The format of this input file is the same as the standard SCC input file. The only difference is that is should contain calibration measurements instead of standard measurements. Following our example, such file should contain the measurement performed at +45 and -45 degrees at 532nm. Also the channel IDs in the file should reflect the ones reported in table 3.3. Moreover this raw input file has to contain the variables: :: double Pol_Calib_Range_Min(channels) - double Pol_Calib\_Range_Max(channels) + double Pol_Calib_Range_Max(channels) where to specify the altitude ranges in meters in which the polarization calibration should be done. +.. figure:: media/figure3.4.png + :height: 806 + :width: 1896 + :scale: 100 % + :align: center + + **Figure 3.4:** To provide polarization calibration (:math:`\eta^*`) values manually just use the button “Add polarization calibration” in the upper-right corner. This option will be available only for a limited period of time. After that only SCC calculated calibration constants will be accepted. + According to the table 3.3 this file should be something similar to: :: dimensions: channels = 4 ; - nb\_of\_time\_scales = 1 ; + nb_of_time_scales = 1 ; points = 16380 ; - scan\_angles = 1 ; + scan_angles = 1 ; time = UNLIMITED ; // (3 currently) variables: - int channel\_ID(channels) ; - double Background\_Low(channels) ; - double Background\_High(channels) ; - int id\_timescale(channels) ; - double Laser\_Pointing\_Angle(scan\_angles) ; - int Molecular\_Calc ; - int Laser\_Pointing\_Angle\_of\_Profiles(time, nb\_of\_time\_scales) ; - int Raw\_Data\_Start\_Time(time, nb\_of\_time\_scales) ; - int Raw\_Data\_Stop\_Time(time, nb\_of\_time\_scales) ; - int Laser\_Shots(time, channels) ; - double Raw\_Lidar\_Data(time, channels, points) ; - double Pressure\_at\_Lidar\_Station ; - double Temperature\_at\_Lidar\_Station ; - double Pol\_Calib\_Range\_Min(channels) ; - double Pol\_Calib\_Range\_Max(channels) ; + int channel_ID(channels) ; + double Background_Low(channels) ; + double Background_High(channels) ; + int id_timescale(channels) ; + double Laser_Pointing_Angle(scan_angles) ; + int Molecular_Calc ; + int Laser_Pointing_Angle_of_Profiles(time, nb_of_time_scales) ; + int Raw_Data_Start_Time(time, nb_of_time_scales) ; + int Raw_Data_Stop_Time(time, nb_of_time_scales) ; + int Laser_Shots(time, channels) ; + double Raw_Lidar_Data(time, channels, points) ; + double Pressure_at_Lidar_Station ; + double Temperature_at_Lidar_Station ; + double Pol_Calib_Range_Min(channels) ; + double Pol_Calib_Range_Max(channels) ; // global attributes: :System = "mysystem" ; - :Longitude\_degrees\_east = 15.723771 ; - :RawData\_Start\_Time\_UT = "220000" ; - :RawData\_Start\_Date = "20130620" ; - :Measurement\_ID = "20130620po00" ; - :Altitude\_meter\_asl = 760. ; - :RawData\_Stop\_Time\_UT = "230333" ; - :Latitude\_degrees\_north = 40.601039 ; + :Longitude_degrees_east = 15.723771 ; + :RawData_Start_Time_UT = "220000" ; + :RawData_Start_Date = "20130620" ; + :Measurement_ID = "20130620po00" ; + :Altitude_meter_asl = 760. ; + :RawData_Stop_Time_UT = "230333" ; + :Latitude_degrees_north = 40.601039 ; data: - channel\_ID = 10, 11, 12, 13 ; + channel_ID = 10, 11, 12, 13 ; - Background\_Low = 30000, 30000, 30000, 30000 ; + Background_Low = 30000, 30000, 30000, 30000 ; - Background\_High = 50000, 50000, 50000, 50000 ; + Background_High = 50000, 50000, 50000, 50000 ; - id\_timescale = 0, 0, 0, 0 ; + id_timescale = 0, 0, 0, 0 ; - Laser\_Pointing\_Angle = 0 ; + Laser_Pointing_Angle = 0 ; - Molecular\_Calc = 0 ; + Molecular_Calc = 0 ; - Laser\_Pointing\_Angle\_of\_Profiles = + Laser_Pointing_Angle_of_Profiles = 0, 0, 0 ; - Raw\_Data\_Start\_Time = + Raw_Data_Start_Time = 0, 300, 600 ; - Raw\_Data\_Stop\_Time = + Raw_Data_Stop_Time = 210, 510, 810 ; - Laser\_Shots = + Laser_Shots = 1200, 1200, 1200, 1200, 1200, 1200, 1200, 1200, 1200, 1200, 1200, 1200 ; - Pressure\_at\_Lidar\_Station = 1010 ; + Pressure_at_Lidar_Station = 1010 ; - Temperature\_at\_Lidar\_Station = 14 ; + Temperature_at_Lidar_Station = 14 ; - Pol\_Calib\_Range\_Min = 1000, 1000, 1000, 1000 ; + Pol_Calib_Range_Min = 1000, 1000, 1000, 1000 ; - Pol\_Calib\_Range\_Min = 2000, 2000, 2000, 2000 ; + Pol_Calib_Range_Min = 2000, 2000, 2000, 2000 ; - Raw\_Lidar\_Data = …...; + Raw_Lidar_Data = …...; The file above assume the following calibration measurements have been done: @@ -548,10 +558,10 @@ As you can see there are 3 cycles of consecutive measurements at +45 and -45 degrees. That way the dimension :code:`time` is set to 3. -The first +/-45 degrees measurement starts at “20130620 22:00:00” (start time of the first +45 measurement) and stops at “20130620 22:03:30” (stop time of the fist -45 measurement). As a consequence, according to the values of the global attributes :code:`RawData\_Start\_Date` and :code:`RawData_Start_Time_UT` we have to set: +The first +/-45 degrees measurement starts at “20130620 22:00:00” (start time of the first +45 measurement) and stops at “20130620 22:03:30” (stop time of the fist -45 measurement). As a consequence, according to the values of the global attributes :code:`RawData_Start_Date` and :code:`RawData_Start_Time_UT` we have to set: :code:`Raw_Data_Start_Time[0]=0` (start of the first +45 measurement in -seconds since :code:`RawData_Start_Time\_UT`) +seconds since :code:`RawData_Start_Time_UT`) :code:`Raw_Data_Stop_Time[0]=210` (stop of the first -45 measurement in seconds since :code:`RawData_Start_Time_UT`) @@ -560,7 +570,7 @@ :code:`Raw_Data_Start_Time[1]=300` (start of the second +45 measurement in seconds since :code:`RawData_Start_Time_UT`) -:code:`Raw_Data_Stop_Time[1]=510` (stop of the second -45 measurement in seconds since :code:`RawData_Start_Time\_UT`) +:code:`Raw_Data_Stop_Time[1]=510` (stop of the second -45 measurement in seconds since :code:`RawData_Start_Time_UT`) :code:`Raw_Data_Start_Time[2]=600` (start of the third +45 measurement in seconds since :code:`RawData_Start_Time_UT`) @@ -592,7 +602,7 @@ :code:`Raw_Lidar_Data[2][3][points]` :math:`\rightarrow` 3\ :sup:`rd` measured transmitted signal at -45 degrees -Once this file has been created it needs to be submitted to the SCC and linked to the configuration “depol\_calibration”. The result of the SCC analysis on this file will be the calculation of the calibration constant h\ :sup:`\*` that will be logged into the SCC database and can be used to calibrate Raman/Elastic backscat ter products (see section 3.3). +Once this file has been created it needs to be submitted to the SCC and linked to the configuration “depol_calibration”. The result of the SCC analysis on this file will be the calculation of the calibration constant h\ :sup:`\*` that will be logged into the SCC database and can be used to calibrate Raman/Elastic backscatter products (see section 3.3). 3.3 Definition of “Raman/Elastic backscatter and linear depolarization ratio” ----------------------------------------------------------------------------- @@ -635,7 +645,7 @@ | 1064nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -Product ID=1, 2, 4, 5, 7 do not need any modification as they do not involve polarization channels. The only product that need to be modified are the Product ID=3 and 6. To produce b532 files containing also *PLDR* we need to modify the “nighttime” and “daytime” configurations to include a product of type “Raman bakscatter and linear depolarization ratio” or “Elastic bakscatter and linear depolarization ratio” respectively. So the configuration reported in table 3.4 should be +Product ID=1, 2, 4, 5, 7 do not need any modification as they do not involve polarization channels. The only product that need to be modified are the Product ID=3 and 6. To produce b532 files containing also *PLDR* we need to modify the “nighttime” and “daytime” configurations to include a product of type “Raman backscatter and linear depolarization ratio” or “Elastic bakscatter and linear depolarization ratio” respectively. So the configuration reported in table 3.4 should be changed to match what is included in table 3.5. **Table 3.5:** The same of table 3.4 but with new product types introduced in SCC v4.0 @@ -676,12 +686,12 @@ It is important to set among the product options of the product ID=10 and 11 which calibration product we want to use for calibration (see section 3.2). This can be done using the SCC web interface setting the appropriate setting in the tab *Polarization calibration products* (see figure 3.4). According to the current example you should set here the calibration product defined in section 3.2. -.. figure:: media/figure3.4.png +.. figure:: media/figure3.5.png :height: 102 :width: 1895 :scale: 100 % :align: center - **Figure 3.4:** How to link a product to calibrate with a calibration product. + **Figure 3.5:** How to link a product to calibrate with a calibration product. -.. warning:: Please not that also *Raman/Elastic backscatter products* need to be linked to a calibration product because the calibration constant and the corresponding correction factor is needed to calculate the total signal out of the two polarization components even if the *PLDR* is not involved in the product calculation. \ No newline at end of file +.. warning:: Please note that also *Raman/Elastic backscatter products* need to be linked to a calibration product because the calibration constant and the corresponding correction factor is needed to calculate the total signal out of the two polarization components even if the *PLDR* is not involved in the product calculation. \ No newline at end of file