changeset
Error fixes in depolarization.rst
Tue, 25 Oct 2016 12:51:13 +0300
- author
- Michael Kottas <mike.kottas@gmail.com>
- date
- Tue, 25 Oct 2016 12:51:13 +0300
- changeset 70
- 31295554bf20
- parent 69
- 380d41368260
- child 71
- 26d9dac079e9
Error fixes in depolarization.rst
| docs/depolarization/depolarization.rst | file | annotate | diff | comparison | revisions |
--- a/docs/depolarization/depolarization.rst Tue Oct 25 12:50:55 2016 +0300 +++ b/docs/depolarization/depolarization.rst Tue Oct 25 12:51:13 2016 +0300 @@ -1,3 +1,10 @@ +.. role:: underline + :class: underline + +.. role:: red + + + 1. Particle Linear Depolarization Ratio Implementation ====================================================== @@ -27,7 +34,7 @@ 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: -.. code-block:: python +:: double VolumeDepol(Length) ; double ErrorVolumeDepol(Length) ; @@ -52,7 +59,7 @@ .. 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 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 :underline:`ONLY` the method b) is considered. 1.3 SCC procedure to calculate the PLDRP ---------------------------------------- @@ -61,11 +68,11 @@ #. 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; +#. 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 :underline:`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; -#. 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. +#. In SCC v4.0 the polarization channels are :underline:`NOT` labeled on the base of their polarization state (as it was done in the SCC v3.11) but :underline:`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. -: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: +:WARNING: In switching from the SCC v3.11 to SCC v4.0 the following modifications have been made on :underline:`ALL` channels of :underline:`ALL` registered configurations: *elPP→elPR* *elCP→elPT* @@ -167,9 +174,7 @@ | 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: - -testing the inline code :code:`test` +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: - :code:`Polarization_Channel_Gain_Factor` (*apparent calibration factor* - :math:`\eta^*` ) - :code:`Polarization_Channel_Gain_Factor_Correction` (*calib. factor corr.* – *K*) @@ -185,8 +190,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 :underline:`REMOVED`; +#. The :underline:`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` @@ -256,26 +261,26 @@ :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: It this variable is found in the SCC input file the corresponding settings in the SCC database will be :underline:`OVERWRITTEN`. Unless you don't have any valid reason to overwrite the database value this variable should not be used. 3. The variables: - .. code-block:: python + :: 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; + have been added. Both these variable are :underline:`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; -4. The variable :code:`Depolarization_Factor` has been *REMOVED*. +4. The variable :code:`Depolarization_Factor` has been :underline:`REMOVED`. 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. - 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). + 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 :underline:`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). - :WARNING: After this transition period *ONLY* automatic calibration will be allowed! + :WARNING: After this transition period :underline:`ONLY` automatic calibration will be allowed! -5. The new *OPTIONAL* variable: +5. The new :underline:`OPTIONAL` variable: :code:`string channel\_string\_ID(channels)` @@ -287,84 +292,55 @@ :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! -Real Example -============ - -This section describes all the practical steps the users need to follow -to switch from SCC v3.11 to new SCC v4.0. - -**IMPORTANT:** +3. Real Example +=============== -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. +This section describes all the practical steps the users need to follow to switch from SCC v3.11 to new SCC v4.0. + +:IMPORTANT: + If your lidar system is not equipped with any polarization channels :underline:`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. -Modification of polarization channel parameters ------------------------------------------------ +3.1 Modification of polarization channel parameters +--------------------------------------------------- -In what it follows it is assumed you already have registered one or more -lidar configurations in the SCC database and that such configurations -have been already used to produce optical products (aerosol extinction -and/or backscatter coefficients) by means of the SCC v3.11. +In what it follows it is assumed you already have registered one or more lidar configurations in the SCC database and that such configurations have been already used to produce optical products (aerosol extinction and/or backscatter coefficients) by means of the SCC v3.11. -Let's assume your 3+2 system is registered in the SCC database and the -settings used by the SCC v3.11 are the ones summarized in table 3.1. +Let's assume your 3+2 system is registered in the SCC database and the settings used by the SCC v3.11 are the ones summarized in table 3.1. -**Table 3.1:** Example of configuration in SCC v3.11 +:Table 3.1: Example of configuration in SCC v3.11 +----------------+--------------+----------------+-------------+-----------+ | Channel Name | Channel ID | Channel Type | nighttime | daytime | +----------------+--------------+----------------+-------------+-----------+ -| 355 | 1 | elT | | | +| 355 | 1 | elT | x | x | +----------------+--------------+----------------+-------------+-----------+ -| 387 | 2 | vrRN2 | | | +| 387 | 2 | vrRN2 | x | | +----------------+--------------+----------------+-------------+-----------+ -| 532 cross | 3 | elCP | | | +| 532 cross | 3 | elCP | x | x | +----------------+--------------+----------------+-------------+-----------+ -| 532 parallel | 4 | elPP | | | +| 532 parallel | 4 | elPP | x | x | +----------------+--------------+----------------+-------------+-----------+ -| 607 | 5 | vrRN2 | | | +| 607 | 5 | vrRN2 | x | | +----------------+--------------+----------------+-------------+-----------+ -| 1064 | 6 | elT | | | +| 1064 | 6 | elT | x | x | +----------------+--------------+----------------+-------------+-----------+ -We assume there are 2 system configurations called “nighttime” and -“daytime”. The nighttime configuration contains all the available lidar -channels (in order to calculate, for example, the aerosol extinction at -355 and 532nm and the aerosol backscatter at 355, 532 and 1064nm) while -in daytime conditions only elastic channels are used (only elastic -backscatter coefficients are generated). +We assume there are 2 system configurations called “nighttime” and “daytime”. The nighttime configuration contains all the available lidar channels (in order to calculate, for example, the aerosol extinction at 355 and 532nm and the aerosol backscatter at 355, 532 and 1064nm) while in daytime conditions only elastic channels are used (only elastic backscatter coefficients are generated). + +To make these settings working with SCC v4.0 it is needed to modify :underline:ONLY` the products properties involving the polarization channels (532 cross and parallel). All the products not involving the polarization channels :underline:`DO NOT` need any modification and should work in the SCC v4.0 exactly as they did in SCC v3.11. In the example above the aerosol extinction and backscatter coefficient at 355nm, the extinction at 532nm as well as the backscatter coefficient at 1064nm do not required any +modification. Let's focus on the modifications needed for the calculation of backscatter at 532nm. -To make these settings working with SCC v4.0 it is needed to modify -*ONLY* the products properties involving the polarization channels (532 -cross and parallel). All the products not involving the polarization -channels *DO NOT* need any modification and should work in the SCC v4.0 -exactly as they did in SCC v3.11. In the example above the aerosol -extinction and backscatter coefficient at 355nm, the extinction at 532nm -as well as the backscatter coefficient at 1064nm do not required any -modification. Let's focus on the modifications needed for the -calculation of backscatter at 532nm. +.. figure:: figure3.1.png + :height: 369 + :width: 1037 + :scale: 100 % + :align: center -|image0| How to select signal types + **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 elCP to elPT and and in the same way the 532 -parallel channel should be changed from elPP to 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 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). **Table 3.2:** The same of table 3.1 but with new channel types introduced in SCC v4.0 @@ -372,487 +348,347 @@ +----------------+--------------+----------------+-------------+-----------+ | Channel Name | Channel ID | Channel Type | nighttime | daytime | +----------------+--------------+----------------+-------------+-----------+ -| 355 | 1 | elT | | | +| 355 | 1 | elT | x | x | +----------------+--------------+----------------+-------------+-----------+ -| 387 | 2 | vrRN2 | | | +| 387 | 2 | vrRN2 | x | | +----------------+--------------+----------------+-------------+-----------+ -| 532 cross | 3 | **elPT** | | | +| 532 cross | 3 | :red:`elPT` | x | x | +----------------+--------------+----------------+-------------+-----------+ -| 532 parallel | 4 | **elPR** | | | +| 532 parallel | 4 | :red:`elPR` | x | x | +----------------+--------------+----------------+-------------+-----------+ -| 607 | 5 | vrRN2 | | | +| 607 | 5 | vrRN2 | x | | +----------------+--------------+----------------+-------------+-----------+ -| 1064 | 6 | elT | | | +| 1064 | 6 | elT | x | x | +----------------+--------------+----------------+-------------+-----------+ -The other change about the polarization channels required to run the SCC -v4.0 is the definition of the polarization crosstalk parameters for all -the polarization channels available. Such parameters can be defined for -each polarization channel using the SCC web interface (see figure 3.2). -In particular among the channel parameters there is a new tab called -*Polarization crosstalk parameters* where it is possible to insert the -values from for the parameters *G* and *H* and the corresponding -statistical and systematic errors if available. In case you have -measured *G* and *H* for your polarization channels please insert the -corresponding values there. Otherwise you can insert the ideal values as -reported in table 1.1. +The other change about the polarization channels required to run the SCC v4.0 is the definition of the polarization crosstalk parameters for all the polarization channels available. Such parameters can be defined for each polarization channel using the SCC web interface (see figure 3.2). In particular among the channel parameters there is a new tab called *Polarization crosstalk parameters* where it is possible to insert the values from for the parameters *G* and *H* and the corresponding statistical and systematic errors if available. In case you have measured *G* and *H* for your polarization channels please insert the corresponding values there. Otherwise you can insert the ideal values as reported in table 1.1. -|image1| *Polarization crosstalk parameters* tab in channel properties -(SCC v4.0). - -Definition of new calibration configuration and product -------------------------------------------------------- +.. figure:: figure3.2.png + :height: 479 + :width: 1890 + :scale: 100 % + :align: center -In this section we will see how to set the polarization calibration -parameters: the calibration constant (called h\ :sup:`\*` 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 D90 calibration method. + **Figure 3.2:** Polarization crosstalk parameters tab in channel properties (SCC v4.0). -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.2 Definition of new calibration configuration and product +----------------------------------------------------------- + +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 :underline:`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 3.1). **Table 3.3:** Polarization calibration configurations assuming D90 calibration method +----------------------------+--------------+----------------+----------------------+ -| Channel Name | Channel ID | Channel Type | depol\_calibration | +| Channel Name | Channel ID | Channel Type | depol_calibration | +----------------------------+--------------+----------------+----------------------+ -| 532 cross +45 degrees | 10 | +45elPT | | +| 532 cross +45 degrees | 10 | +45elPT | x | +----------------------------+--------------+----------------+----------------------+ -| 532 parallel +45 degrees | 11 | +45elPR | | +| 532 parallel +45 degrees | 11 | +45elPR | x | +----------------------------+--------------+----------------+----------------------+ -| 532 cross -45 degrees | 12 | -45elPT | | +| 532 cross -45 degrees | 12 | -45elPT | x | +----------------------------+--------------+----------------+----------------------+ -| 532 parallel -45 degrees | 13 | -45elPR | | +| 532 parallel -45 degrees | 13 | -45elPR | x | +----------------------------+--------------+----------------+----------------------+ -Finally we should add to the configuration “depol\_calibration” a -product “\ *Linear polarization calibration”* to be used for the -calibration. According to the example given above and to the usecase -document attached we should use an usecase=4 for this example. +Finally we should add to the configuration “depol_calibration” a product “*Linear polarization calibration”* to be used for the calibration. According to the example given above and to the usecase document attached we should use an usecase=4 for this example. -Other “\ *Linear polarization calibration”* options to be specified are -reported in figure 3.3. The most important factor you should insert here -is the *Pol calibration correction factor* (K). The ideal value for this -parameter is 1. Anyway if you have measured the parameter K please fill -in the measured value and the corresponding measurement errors. - -|image2| Options for *Linear polarization calibration product*. +Other “*Linear polarization calibration”* options to be specified are reported in figure 3.3. The most important factor you should insert here is the *Pol calibration correction factor* (K). The ideal value for this parameter is 1. Anyway if you have measured the parameter K please fill in the measured value and the corresponding measurement errors. -As you can see it is possible to fill in only the K correction factor -and not the calibration constant h\ :sup:`\*`. +.. figure:: figure3.3.png + :height: 495 + :width: 1887 + :scale: 100 % + :align: center -Actually for a *LIMITED* period of time it will be possible to fill in -also the constant h\ :sup:`\*` 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 h\ :sup:`\*` -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. + **Figure 3.3:** Options for *Linear polarization calibration product*. + +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 :underline:`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 :underline:`NOT` possible to provide calibration constant using this procedure and the parameter :math:`\eta^*` can be calculated :underline:`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_Min(channels) + double Pol_Calib\_Range_Max(channels) -where to specify the altitude ranges in meters in which the polarization -calibration should be done. +where to specify the altitude ranges in meters in which the polarization calibration should be done. According to the table 3.3 this file should be something similar to: - -dimensions: - -channels = 4 ; - -nb\_of\_time\_scales = 1 ; - -points = 16380 ; - -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) ; + dimensions: + channels = 4 ; + nb\_of\_time\_scales = 1 ; + points = 16380 ; + 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) ; -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: + // 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 ; -:System = "mysystem" ; - -:Longitude\_degrees\_east = 15.723771 ; - -:RawData\_Start\_Time\_UT = "220000" ; + data: + channel\_ID = 10, 11, 12, 13 ; -:RawData\_Start\_Date = "20130620" ; + Background\_Low = 30000, 30000, 30000, 30000 ; -:Measurement\_ID = "20130620po00" ; + Background\_High = 50000, 50000, 50000, 50000 ; -:Altitude\_meter\_asl = 760. ; + id\_timescale = 0, 0, 0, 0 ; -:RawData\_Stop\_Time\_UT = "230333" ; + Laser\_Pointing\_Angle = 0 ; -:Latitude\_degrees\_north = 40.601039 ; - -data: + Molecular\_Calc = 0 ; -channel\_ID = 10, 11, 12, 13 ; - -Background\_Low = 30000, 30000, 30000, 30000 ; - -Background\_High = 50000, 50000, 50000, 50000 ; + Laser\_Pointing\_Angle\_of\_Profiles = + 0, + 0, + 0 ; -id\_timescale = 0, 0, 0, 0 ; - -Laser\_Pointing\_Angle = 0 ; - -Molecular\_Calc = 0 ; - -Laser\_Pointing\_Angle\_of\_Profiles = - -0, - -0, + Raw\_Data\_Start\_Time = + 0, + 300, + 600 ; -0 ; - -Raw\_Data\_Start\_Time = - -0, + Raw\_Data\_Stop\_Time = + 210, + 510, + 810 ; -300, - -600 ; + Laser\_Shots = + 1200, 1200, 1200, 1200, + 1200, 1200, 1200, 1200, + 1200, 1200, 1200, 1200 ; -Raw\_Data\_Stop\_Time = + Pressure\_at\_Lidar\_Station = 1010 ; -210, + Temperature\_at\_Lidar\_Station = 14 ; -510, + Pol\_Calib\_Range\_Min = 1000, 1000, 1000, 1000 ; -810 ; + Pol\_Calib\_Range\_Min = 2000, 2000, 2000, 2000 ; -Laser\_Shots = + Raw\_Lidar\_Data = …...; -1200, 1200, 1200, 1200, - -1200, 1200, 1200, 1200, +The file above assume the following calibration measurements have been done: -1200, 1200, 1200, 1200 ; +1. First +45 degrees acquisition followed by a corresponding -45 degrees acquisition -Pressure\_at\_Lidar\_Station = 1010 ; - -Temperature\_at\_Lidar\_Station = 14 ; + a. Measurement at +45 degrees -Pol\_Calib\_Range\_Min = 1000, 1000, 1000, 1000 ; + Start Time: 20130620 22:00:00 -Pol\_Calib\_Range\_Min = 2000, 2000, 2000, 2000 ; + Stop Time: 20130620 22:01:00 -Raw\_Lidar\_Data = …...; + Shots: 1200 -The file above assume the following calibration measurements have been -done: + b. Measurement at -45 degrees -1. First +45 degrees acquisition followed by a corresponding -45 degrees - acquisition + Start Time: 20130620 22:02:30 - a. Measurement at +45 degrees + Stop Time: 20130620 22:03:30 -Start Time: 20130620 22:00:00 - -Stop Time: 20130620 22:01:00 + Shots: 1200 -Shots: 1200 +2. Second +45 degrees acquisition followed by a corresponding -45 degrees acquisition -a. Measurement at -45 degrees - -Start Time: 20130620 22:02:30 + a. Measurement at +45 degrees -Stop Time: 20130620 22:03:30 + Start Time: 20130620 22:05:00 -Shots: 1200 + Stop Time: 20130620 22:06:00 -1. Second +45 degrees acquisition followed by a corresponding -45 - degrees acquisition + Shots: 1200 - a. Measurement at +45 degrees + b. Measurement at -45 degrees -Start Time: 20130620 22:05:00 - -Stop Time: 20130620 22:06:00 + Start Time: 20130620 22:07:30 -Shots: 1200 + Stop Time: 20130620 22:08:30 -a. Measurement at -45 degrees - -Start Time: 20130620 22:07:30 + Shots: 1200 -Stop Time: 20130620 22:08:30 +3. Third +45 degrees acquisition followed by a corresponding -45 degrees acquisition -Shots: 1200 + a. Measurement at +45 degrees -1. Third +45 degrees acquisition followed by a corresponding -45 degrees - acquisition + Start Time: 20130620 22:10:00 - a. Measurement at +45 degrees + Stop Time: 20130620 22:11:00 -Start Time: 20130620 22:10:00 - -Stop Time: 20130620 22:11:00 + Shots: 1200 -Shots: 1200 + b. Measurement at -45 degrees -a. Measurement at -45 degrees + Start Time: 20130620 22:12:30 -Start Time: 20130620 22:12:30 + Stop Time: 20130620 22:13:30 -Stop Time: 20130620 22:13:30 - -Shots: 1200 + Shots: 1200 -As you can see there are 3 cycles of consecutive measurements at +45 and --45 degrees. That's way the dimension time is set to 3. +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 RawData\_Start\_Date and -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`) -Raw\_Data\_Start\_Time[0]=0 (start of the first +45 measurement in -seconds since RawData\_Start\_Time\_UT) - -Raw\_Data\_Stop\_Time[0]=210 (stop of the first -45 measurement in -seconds since 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`) Following a similar procedure for the other 2 cycles we have: -Raw\_Data\_Start\_Time[1]=300 (start of the second +45 measurement in -seconds since RawData\_Start\_Time\_UT) +:code:`Raw_Data_Start_Time[1]=300` (start of the second +45 measurement in seconds since :code:`RawData_Start_Time_UT`) -Raw\_Data\_Stop\_Time[1]=510 (stop of the second -45 measurement in -seconds since 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`) -Raw\_Data\_Start\_Time[2]=600 (start of the third +45 measurement in -seconds since 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`) -Raw\_Data\_Stop\_Time[2]=810 (stop of the third -45 measurement in -seconds since RawData\_Start\_Time\_UT) - -Moreover, according to the order of the channels in the channel\_ID -variable, the Raw\_Lidar\_Data array should be filled as it follows: +:code:`Raw_Data_Stop_Time[2]=810` (stop of the third -45 measurement in seconds since :code:`RawData_Start_Time_UT`) -Raw\_Lidar\_Data[0][0][points] → 1\ :sup:`st` measured transmitted -signal at +45 degrees +Moreover, according to the order of the channels in the :code:`channel_ID` variable, the :code:`Raw_Lidar_Data` array should be filled as it follows: -Raw\_Lidar\_Data[0][1][points] → 1\ :sup:`st` measured reflected signal -at +45 degrees +:code:`Raw_Lidar_Data[0][0][points]` :math:`\rightarrow` 1\ :sup:`st` measured transmitted signal at +45 degrees -Raw\_Lidar\_Data[0][2][points] → 1\ :sup:`st` measured transmitted -signal at -45 degrees +:code:`Raw_Lidar_Data[0][1][points]` :math:`\rightarrow` 1\ :sup:`st` measured transmitted signal at +45 degrees -Raw\_Lidar\_Data[0][3][points] → 1\ :sup:`st` measured reflected signal -at -45 degrees +:code:`Raw_Lidar_Data[0][2][points]` :math:`\rightarrow` 1\ :sup:`st` measured transmitted signal at -45 degrees -Raw\_Lidar\_Data[1][0][points] → 2\ :sup:`nd` measured transmitted -signal at +45 degrees +:code:`Raw_Lidar_Data[0][3][points]` :math:`\rightarrow` 1\ :sup:`st` measured transmitted signal at -45 degrees -Raw\_Lidar\_Data[1][1][points] → 2\ :sup:`nd` measured reflected signal -at +45 degrees +:code:`Raw_Lidar_Data[1][0][points]` :math:`\rightarrow` 2\ :sup:`nd` measured transmitted signal at +45 degrees -Raw\_Lidar\_Data[1][2][points] → 2\ :sup:`nd` measured transmitted -signal at -45 degrees +:code:`Raw_Lidar_Data[1][1][points]` :math:`\rightarrow` 2\ :sup:`nd` measured transmitted signal at +45 degrees -Raw\_Lidar\_Data[1][3][points] → 2\ :sup:`nd` measured reflected signal -at -45 degrees +:code:`Raw_Lidar_Data[1][2][points]` :math:`\rightarrow` 2\ :sup:`nd` measured transmitted signal at -45 degrees -Raw\_Lidar\_Data[2][0][points] → 3\ :sup:`rd` measured transmitted -signal at +45 degrees +:code:`Raw_Lidar_Data[1][3][points]` :math:`\rightarrow` 2\ :sup:`nd` measured transmitted signal at -45 degrees -Raw\_Lidar\_Data[2][1][points] → 3\ :sup:`rd` measured reflected signal -at +45 degrees +:code:`Raw_Lidar_Data[2][0][points]` :math:`\rightarrow` 3\ :sup:`rd` measured transmitted signal at +45 degrees -Raw\_Lidar\_Data[2][2][points] → 3\ :sup:`rd` measured transmitted -signal at -45 degrees +:code:`Raw_Lidar_Data[2][1][points]` :math:`\rightarrow` 3\ :sup:`rd` measured transmitted signal at +45 degrees -Raw\_Lidar\_Data[2][3][points] → 3\ :sup:`rd` measured reflected signal -at -45 degrees +:code:`Raw_Lidar_Data[2][2][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). +:code:`Raw_Lidar_Data[2][3][points]` :math:`\rightarrow` 3\ :sup:`rd` measured transmitted signal at -45 degrees -**Definition of “Raman/Elastic backscatter and linear depolarization ratio”** +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). + +3.3 Definition of “Raman/Elastic backscatter and linear depolarization ratio” ----------------------------------------------------------------------------- -In order to calculate the *PLDR* we need to modify the polarization -related products linked to the “standard” measurement configurations -(the configuration called “nighttime” and/or “daytime” in table 3.2). +In order to calculate the *PLDR* we need to modify the polarization related products linked to the “standard” measurement configurations (the configuration called “nighttime” and/or “daytime” in table 3.2). -Let's suppose we have defined the following products (defined already in -SCC v3.11): +Let's suppose we have defined the following products (defined already in SCC v3.11): **Table 3.4:** Example of products configuration in SCC v3.11 +-----------------------+--------------+-----------------------+-------------+-----------+ | Product Name | Product ID | Product Type | nighttime | daytime | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Raman backscatter | 1 | Raman backscatter | | | +| Raman backscatter | 1 | Raman backscatter | x | | | | | | | | | 355nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Extinction | 2 | Extinction | | | +| Extinction | 2 | Extinction | x | | | | | | | | | 387nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Raman backscatter | 3 | Raman backscatter | | | +| Raman backscatter | 3 | Raman backscatter | x | | | | | | | | | 532nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Extinction | 4 | Extinction | | | +| Extinction | 4 | Extinction | x | | | | | | | | | 532nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Elastic backscatter | 5 | Elastic backscatter | | | +| Elastic backscatter | 5 | Elastic backscatter | | x | | | | | | | | 355nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Elastic backscatter | 6 | Elastic backscatter | | | +| Elastic backscatter | 6 | Elastic backscatter | | x | | | | | | | | 532nm | | | | | +-----------------------+--------------+-----------------------+-------------+-----------+ -| Elastic backscatter | 7 | Elastic backscatter | | | +| Elastic backscatter | 7 | Elastic backscatter | x | x | | | | | | | | 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 bakscatter 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 +**Table 3.5:** The same of table 3.4 but with new product types introduced in SCC v4.0 +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ | Product Name | Product ID | Product Type | nighttime | daytime | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Raman backscatter | 1 | Raman backscatter | | | +| Raman backscatter | 1 | Raman backscatter | x | | | | | | | | | 355nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Extinction | 2 | Extinction | | | +| Extinction | 2 | Extinction | x | | | | | | | | | 387nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Raman backscatter | 10 | **Raman backscatter and linear depolarization ratio** | | | +| Raman backscatter | 10 | :red:`Raman backscatter and linear depolarization ratio` | x | | | | | | | | | 532nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Extinction | 4 | Extinction | | | +| Extinction | 4 | Extinction | x | | | | | | | | | 532nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Elastic backscatter | 5 | Elastic backscatter | | | +| Elastic backscatter | 5 | Elastic backscatter | | x | | | | | | | | 355nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Elastic backscatter | 11 | **Elastic backscatter and linear depolarization ratio** | | | +| Elastic backscatter | 11 | :red:`Elastic backscatter and linear depolarization ratio`| | x | | | | | | | | 532nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -| Elastic backscatter | 7 | Elastic backscatter | | | +| Elastic backscatter | 7 | Elastic backscatter | x | x | | | | | | | | 1064nm | | | | | +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ -As you can see in table 3.5, the old product IDs=3 and 6 (present in -table 3.4) have been replaced with the new product ID=10 and 11 to -guarantee the calculation of *PLDR*. +As you can see in table 3.5, the old product IDs=3 and 6 (present in table 3.4) have been replaced with the new product ID=10 and 11 to guarantee the calculation of *PLDR*. -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. - -|image3| How to link a product to calibrate with a calibration product. +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. -**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. +.. figure:: figure3.4.png + :height: 102 + :width: 1895 + :scale: 100 % + :align: center -.. |image0| image:: ./media/image1.png - :width: 6.69514in - :height: 2.40764in -.. |image1| image:: ./media/image2.png - :width: 6.69306in - :height: 1.71458in -.. |image2| image:: ./media/image3.png - :width: 6.69306in - :height: 1.77431in -.. |image3| image:: ./media/image4.png - :width: 6.69306in - :height: 0.36389in + **Figure 3.4:** 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
