mike@68: 1. Particle Linear Depolarization Ratio Implementation mike@68: ====================================================== ulalume3@67: mike@68: 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. ulalume3@67: mike@68: 1.1 Background ulalume3@67: -------------- ulalume3@67: mike@68: The calculation of the volume linear depolarization ratio profile (*VLDR*) and particle linear depolarization ratio profile (*PLDR*) needs two different steps: ulalume3@67: mike@68: #. the calibration of the polarization sensitive lidar channels; mike@68: #. the calculation of the *VLDR* or *PLDR* itself. mike@68: mike@68: 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. mike@68: mike@68: 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). ulalume3@67: mike@68: New signal types have been introduced to take into account special channel configurations used for calibration purposes. ulalume3@67: mike@68: 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: ulalume3@67: mike@68: #. *Linear polarization calibration (factor* h) *(product\_type\_id=6);* mike@68: #. *Raman backscatter and linear depolarization ratio ulalume3@67: (product\_type\_id=7);* mike@68: #. *Elastic backscatter and linear depolarization ratio ulalume3@67: (product\_type\_id=8).* ulalume3@67: mike@68: 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: ulalume3@67: mike@68: .. code-block:: python ulalume3@67: mike@68: double VolumeDepol(Length) ; mike@68: double ErrorVolumeDepol(Length) ; mike@68: ErrorVolumeDepol:long\_name = "absolute error of VolumeDepol" ; mike@68: double ParticleDepol(Length) ; mike@68: double ErrorParticleDepol(Length) ; mike@68: ErrorParticleDepol:long\_name = "absolute error of ParticleDepol" ; ulalume3@67: mike@68: 1.2 Polarization calibration ulalume3@67: ---------------------------- ulalume3@67: mike@68: An important point is the definition of reliable *PLDR* calibration procedures. Within EARLINET the following calibration procedures are currently used: ulalume3@67: mike@68: a) Rayleigh calibration; mike@68: b) +45 calibration method, or D90 calibration method (made by +45 and -45 measurements); mike@68: c) 3 signals (total, cross and parallel). ulalume3@67: mike@68: 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. ulalume3@67: mike@68: For what it concerns the method c) it, basically, requires to solve the equation: ulalume3@67: mike@68: .. math:: mike@68: \alpha_s P_s + \alpha_p P_p = P ulalume3@67: mike@68: 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. mike@68: mike@68: 1.3 SCC procedure to calculate the PLDRP mike@68: ---------------------------------------- mike@68: mike@68: According to what mentioned before the SCC calculates the *PLDR* through the following steps: ulalume3@67: mike@68: #. 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); mike@68: mike@68: #. 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; ulalume3@67: mike@68: #. 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. ulalume3@67: mike@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: mike@68: *elPP→elPR* ulalume3@67: mike@68: *elCP→elPT* ulalume3@67: mike@68: *elPPnr→elPRnr* ulalume3@67: mike@68: *elPPfr→ elPRfr* ulalume3@67: mike@68: *elCPnr→ elPTnr* ulalume3@67: mike@68: *elCPfr→ elPTfr* mike@68: mike@68: Please be sure these modifications reflect to your actual lidar setup(cross channels are transmitted and parallel channels are reflected); ulalume3@67: mike@68: 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); mike@68: #. The file at point 2 is pre-processed by **ELPP** module which applies the standard pre-processing procedures applied to “standard” lidar data; mike@68: #. 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; mike@68: #. 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; mike@68: #. 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); mike@68: #. The user needs to submit another SCC raw data file containing the “standard” measurements; mike@68: #. 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; ulalume3@67: mike@68: #. 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. mike@68: mike@68: 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. ulalume3@67: mike@68: In case of +45 calibration method :math:`\eta^*` is calculated by: ulalume3@67: mike@68: .. math:: mike@68: \eta^* = \frac{I_R}{I_T}(+45) ulalume3@67: ulalume3@67: While in case of D90 calibration method: ulalume3@67: mike@68: .. math:: mike@68: \eta^* = \sqrt{\frac{I_R}{I_T}(+45) \frac{I_R}{I_T}(-45)} ulalume3@67: ulalume3@67: **ELDA** module calculates the “apparent” *VLDR*: ulalume3@67: mike@68: .. math:: mike@68: \delta^* = \frac{K}{\eta^*} \cdot \frac{I_R}{I_T} ulalume3@67: ulalume3@67: the *VLDR* ulalume3@67: mike@68: .. math:: mike@68: \delta = \frac{\delta^*(G_T + H_T)-(G_R + H_R)}{(G_R - H_R) - \delta^*(G_T - H_T)} ulalume3@67: ulalume3@67: and the *PLDR* ulalume3@67: mike@68: .. math:: mike@68: \delta_{\alpha} = \frac{(1 + \delta_m)\delta R - (1 + \delta)\delta_m}{(1 + \delta_m)R - (1 + \delta)} ulalume3@67: ulalume3@67: where: ulalume3@67: mike@68: - :math:`\eta^*` is the *apparent calibration factor* calculated by **scc\_calibrator** mike@68: mike@68: - *K* is the *calibration factor correction* defined as polarization product option ulalume3@67: mike@68: - :math:`I_T` and :math:`I_R` are the transmitted and the reflected signals in the polarization detection set-up ulalume3@67: mike@68: - :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. ulalume3@67: mike@68: - :math:`\delta_m` is the molecular linear depolarization ratio calculated by ELPP mike@68: mike@68: - *R* is the backscatter ratio ulalume3@67: mike@68: 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. ulalume3@67: mike@68: 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: ulalume3@67: mike@68: .. math:: mike@68: I_{total} \propto \frac{\eta^*}{K}H_R I_T - H_T I_R ulalume3@67: mike@68: 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). ulalume3@67: mike@68: 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: ulalume3@67: mike@68: .. math:: mike@68: G_T=1 , \qquad H_T=-1, \qquad G_R=1, \qquad H_R=1 mike@68: mike@68: 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: ulalume3@67: mike@68: .. math:: mike@68: G_T=1 , \qquad H_T=0, \qquad G_R=1, \qquad H_R=-1 ulalume3@67: ulalume3@67: mike@68: **Table 1.1:** Polarization cross-talk correction parameters for ideal systems ulalume3@67: ulalume3@67: +----------------------+-----------------------------+-----------------+-----------------+-----------------+ ulalume3@67: | Laser polarization | Detected in lidar channel | mike@68: + +-----------------------------+-----------------+-----------------+-----------------+ ulalume3@67: | | Transmitted | Reflected | mike@68: + +-----------------------------+-----------------+-----------------+-----------------+ mike@68: | | :math:`G_T` | :math:`H_T` | :math:`G_R` | :math:`H_R` | ulalume3@67: +----------------------+-----------------------------+-----------------+-----------------+-----------------+ ulalume3@67: | total | 1 | 0 | 1 | 0 | ulalume3@67: +----------------------+-----------------------------+-----------------+-----------------+-----------------+ ulalume3@67: | parallel | 1 | 1 | 1 | 1 | ulalume3@67: +----------------------+-----------------------------+-----------------+-----------------+-----------------+ ulalume3@67: | cross | 1 | -1 | 1 | -1 | ulalume3@67: +----------------------+-----------------------------+-----------------+-----------------+-----------------+ ulalume3@67: mike@68: 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: ulalume3@67: mike@68: testing the inline code :code:`test` ulalume3@67: mike@68: - :code:`Polarization_Channel_Gain_Factor` (*apparent calibration factor* - :math:`\eta^*` ) mike@68: - :code:`Polarization_Channel_Gain_Factor_Correction` (*calib. factor corr.* – *K*) mike@68: - :code:`G_T` mike@68: - :code:`H_T` mike@68: - :code:`G_R` mike@68: - :code:`H_R` ulalume3@67: mike@68: 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. ulalume3@67: mike@68: 2. Changes of the SCC input format mike@68: ================================== ulalume3@67: ulalume3@67: The following minor changes have been applied to raw SCC data format: ulalume3@67: mike@68: #. The optional variable *ID\_Range* has been *REMOVED*; mike@68: #. The *OPTIONAL* variable :code:`int Signal\_Type(channels)` has been added. The possible values are the same available in the SCC\_DB: ulalume3@67: mike@68: :code:`0` :math:`\rightarrow` :code:`elT` ulalume3@67: mike@68: :code:`1` :math:`\rightarrow` :code:`elTnr` ulalume3@67: mike@68: :code:`2` :math:`\rightarrow` :code:`elTfr` ulalume3@67: mike@68: :code:`3` :math:`\rightarrow` :code:`vrRN2` ulalume3@67: mike@68: :code:`4` :math:`\rightarrow` :code:`vrRN2nr` ulalume3@67: mike@68: :code:`5` :math:`\rightarrow` :code:`vrRN2fr` ulalume3@67: mike@68: :code:`6` :math:`\rightarrow` :code:`elPR` ulalume3@67: mike@68: :code:`7` :math:`\rightarrow` :code:`elPT` ulalume3@67: mike@68: :code:`8` :math:`\rightarrow` :code:`pRRlow` ulalume3@67: mike@68: :code:`9` :math:`\rightarrow` :code:`pRRhigh` ulalume3@67: mike@68: :code:`10` :math:`\rightarrow` :code:`elPRnr` ulalume3@67: mike@68: :code:`11` :math:`\rightarrow` :code:`elPRfr` ulalume3@67: mike@68: :code:`12` :math:`\rightarrow` :code:`elPTnr` ulalume3@67: mike@68: :code:`13` :math:`\rightarrow` :code:`elPTfr` ulalume3@67: mike@68: :code:`14` :math:`\rightarrow` :code:`vrRH2O` ulalume3@67: mike@68: :code:`15` :math:`\rightarrow` :code:`pRRhighnr` ulalume3@67: mike@68: :code:`16` :math:`\rightarrow` :code:`pRRhighfr` ulalume3@67: mike@68: :code:`17` :math:`\rightarrow` :code:`pRRlownr` ulalume3@67: mike@68: :code:`18` :math:`\rightarrow` :code:`pRRlowfr` ulalume3@67: mike@68: :code:`19` :math:`\rightarrow` :code:`vrRH2Onr` ulalume3@67: mike@68: :code:`20` :math:`\rightarrow` :code:`vrRH2Ofr` ulalume3@67: mike@68: :code:`21` :math:`\rightarrow` :code:`elTunr` ulalume3@67: mike@68: :code:`22` :math:`\rightarrow` :code:`+45elPT` ulalume3@67: mike@68: :code:`23` :math:`\rightarrow` :code:`+45elPR` ulalume3@67: mike@68: :code:`24` :math:`\rightarrow` :code:`-45elPT` ulalume3@67: mike@68: :code:`25` :math:`\rightarrow` :code:`-45elPR` ulalume3@67: mike@68: :code:`26` :math:`\rightarrow` :code:`+45elPTnr` ulalume3@67: mike@68: :code:`27` :math:`\rightarrow` :code:`+45elPTfr` ulalume3@67: mike@68: :code:`28` :math:`\rightarrow` :code:`+45elPRnr` mike@68: mike@68: :code:`29` :math:`\rightarrow` :code:`+45elPRfr` ulalume3@67: mike@68: :code:`30` :math:`\rightarrow` :code:`-45elPTnr` ulalume3@67: mike@68: :code:`31` :math:`\rightarrow` :code:`-45elPTfr` ulalume3@67: mike@68: :code:`32` :math:`\rightarrow` :code:`-45elPRnr` ulalume3@67: mike@68: :code:`33` :math:`\rightarrow` :code:`-45elPRfr` ulalume3@67: mike@68: :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. ulalume3@67: mike@68: 3. The variables: ulalume3@67: mike@68: .. code-block:: python mike@68: mike@68: double Pol\_Calib\_Range\_Min(channels) mike@68: double Pol\_Calib\_Range\_Max(channels) ulalume3@67: mike@68: 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; mike@68: mike@68: 4. The variable :code:`Depolarization_Factor` has been *REMOVED*. ulalume3@67: mike@68: 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. mike@68: mike@68: 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). ulalume3@67: mike@68: :WARNING: After this transition period *ONLY* automatic calibration will be allowed! ulalume3@67: mike@68: 5. The new *OPTIONAL* variable: ulalume3@67: mike@68: :code:`string channel\_string\_ID(channels)` ulalume3@67: mike@68: has been introduced. mike@68: mike@68: 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. ulalume3@67: mike@68: 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). ulalume3@67: mike@68: :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! ulalume3@67: ulalume3@67: Real Example ulalume3@67: ============ ulalume3@67: ulalume3@67: This section describes all the practical steps the users need to follow ulalume3@67: to switch from SCC v3.11 to new SCC v4.0. ulalume3@67: ulalume3@67: **IMPORTANT:** ulalume3@67: ulalume3@67: If your lidar system is not equipped with any polarization channels *NO* ulalume3@67: changes are required. In this case, the SCC v4.0 should work using the ulalume3@67: same input files and the same database configurations you have used with ulalume3@67: the SCC v3.11. Anyway as in the SCC v4.0 several bugs have been fixed, ulalume3@67: it is recommended to re-run all the measurement IDs you have submitted. ulalume3@67: For doing that you just need to reprocess all your data without the need ulalume3@67: to submit raw data files already uploaded on the server. ulalume3@67: ulalume3@67: The practical example reported below describes the modifications ulalume3@67: required to use the SCC v4.0 for lidar systems equipped with ulalume3@67: polarization channels. ulalume3@67: ulalume3@67: Modification of polarization channel parameters ulalume3@67: ----------------------------------------------- ulalume3@67: ulalume3@67: In what it follows it is assumed you already have registered one or more ulalume3@67: lidar configurations in the SCC database and that such configurations ulalume3@67: have been already used to produce optical products (aerosol extinction ulalume3@67: and/or backscatter coefficients) by means of the SCC v3.11. ulalume3@67: ulalume3@67: Let's assume your 3+2 system is registered in the SCC database and the ulalume3@67: settings used by the SCC v3.11 are the ones summarized in table 3.1. ulalume3@67: ulalume3@67: **Table 3.1:** Example of configuration in SCC v3.11 ulalume3@67: ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | Channel Name | Channel ID | Channel Type | nighttime | daytime | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 355 | 1 | elT |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 387 | 2 | vrRN2 |  | | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 532 cross | 3 | elCP |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 532 parallel | 4 | elPP |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 607 | 5 | vrRN2 |  | | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 1064 | 6 | elT |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: ulalume3@67: We assume there are 2 system configurations called “nighttime” and ulalume3@67: “daytime”. The nighttime configuration contains all the available lidar ulalume3@67: channels (in order to calculate, for example, the aerosol extinction at ulalume3@67: 355 and 532nm and the aerosol backscatter at 355, 532 and 1064nm) while ulalume3@67: in daytime conditions only elastic channels are used (only elastic ulalume3@67: backscatter coefficients are generated). ulalume3@67: ulalume3@67: To make these settings working with SCC v4.0 it is needed to modify ulalume3@67: *ONLY* the products properties involving the polarization channels (532 ulalume3@67: cross and parallel). All the products not involving the polarization ulalume3@67: channels *DO NOT* need any modification and should work in the SCC v4.0 ulalume3@67: exactly as they did in SCC v3.11. In the example above the aerosol ulalume3@67: extinction and backscatter coefficient at 355nm, the extinction at 532nm ulalume3@67: as well as the backscatter coefficient at 1064nm do not required any ulalume3@67: modification. Let's focus on the modifications needed for the ulalume3@67: calculation of backscatter at 532nm. ulalume3@67: ulalume3@67: |image0| How to select signal types ulalume3@67: ulalume3@67: The first modification concerns the settings of the channel type for the ulalume3@67: 532 cross and 532 parallel polarization channels. Starting from SCC v4.0 ulalume3@67: polarization channels are identified as transmitted and reflected ulalume3@67: polarization channels and not on the base of their polarization state. ulalume3@67: So suppose if we suppose the cross polarized channel is transmitted by a ulalume3@67: polarizer beam splitter cube and the parallel is reflected the value ulalume3@67: reported in table 3.1 should be modified as they appear in table 3.2. So ulalume3@67: using the SCC web interface, the signal type of the 532 cross channel ulalume3@67: should be changed from elCP to elPT and and in the same way the 532 ulalume3@67: parallel channel should be changed from elPP to elPR (see figure 3.1). ulalume3@67: ulalume3@67: **Table 3.2:** The same of table 3.1 but with new channel types ulalume3@67: introduced in SCC v4.0 ulalume3@67: ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | Channel Name | Channel ID | Channel Type | nighttime | daytime | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 355 | 1 | elT |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 387 | 2 | vrRN2 |  | | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 532 cross | 3 | **elPT** |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 532 parallel | 4 | **elPR** |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 607 | 5 | vrRN2 |  | | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: | 1064 | 6 | elT |  |  | ulalume3@67: +----------------+--------------+----------------+-------------+-----------+ ulalume3@67: ulalume3@67: The other change about the polarization channels required to run the SCC ulalume3@67: v4.0 is the definition of the polarization crosstalk parameters for all ulalume3@67: the polarization channels available. Such parameters can be defined for ulalume3@67: each polarization channel using the SCC web interface (see figure 3.2). ulalume3@67: In particular among the channel parameters there is a new tab called ulalume3@67: *Polarization crosstalk parameters* where it is possible to insert the ulalume3@67: values from for the parameters *G* and *H* and the corresponding ulalume3@67: statistical and systematic errors if available. In case you have ulalume3@67: measured *G* and *H* for your polarization channels please insert the ulalume3@67: corresponding values there. Otherwise you can insert the ideal values as ulalume3@67: reported in table 1.1. ulalume3@67: ulalume3@67: |image1| *Polarization crosstalk parameters* tab in channel properties ulalume3@67: (SCC v4.0). ulalume3@67: ulalume3@67: Definition of new calibration configuration and product ulalume3@67: ------------------------------------------------------- ulalume3@67: ulalume3@67: In this section we will see how to set the polarization calibration ulalume3@67: parameters: the calibration constant (called h\ :sup:`\*` in section ulalume3@67: 1.3) and the correction to calibration constant (called K in section ulalume3@67: 1.3). ulalume3@67: ulalume3@67: In order to provide such parameters you need to define a new system ulalume3@67: configuration to be used *only* for calibration purposes. Such new ulalume3@67: configuration should include the polarization channels in the ulalume3@67: measurement configuration used for the calibration. Let's suppose we ulalume3@67: want to use the D90 calibration method. ulalume3@67: ulalume3@67: In this case we need to define a new configuration (called for example ulalume3@67: “depol\_calibration”) as reported in the table 3.3. As you can see the ulalume3@67: configuration “depol\_calibration” includes 4 “new” channels. Actually ulalume3@67: the channels “532 cross +45 degrees” (channel ID=10) and “532 cross -45 ulalume3@67: degrees” (channel ID=12) refer to the same physical channel “532 cross” ulalume3@67: reported with channel ID=3 in table 3.2. Anyway we need to define two ulalume3@67: new channel IDs to identify the “532 cross” channel in the two ulalume3@67: polarization rotated configurations (+45 and -45 degrees) needed to ulalume3@67: apply the D90 calibration method. The same is true for the “532 ulalume3@67: parallel” channel. The polarization rotated channels should be labeled ulalume3@67: with the corresponding signal type as reported in table 3.3 (see figure ulalume3@67: 3.1). ulalume3@67: ulalume3@67: **Table 3.3:** Polarization calibration configurations assuming D90 ulalume3@67: calibration method ulalume3@67: ulalume3@67: +----------------------------+--------------+----------------+----------------------+ ulalume3@67: | Channel Name | Channel ID | Channel Type | depol\_calibration | ulalume3@67: +----------------------------+--------------+----------------+----------------------+ ulalume3@67: | 532 cross +45 degrees | 10 | +45elPT |  | ulalume3@67: +----------------------------+--------------+----------------+----------------------+ ulalume3@67: | 532 parallel +45 degrees | 11 | +45elPR |  | ulalume3@67: +----------------------------+--------------+----------------+----------------------+ ulalume3@67: | 532 cross -45 degrees | 12 | -45elPT |  | ulalume3@67: +----------------------------+--------------+----------------+----------------------+ ulalume3@67: | 532 parallel -45 degrees | 13 | -45elPR |  | ulalume3@67: +----------------------------+--------------+----------------+----------------------+ ulalume3@67: ulalume3@67: Finally we should add to the configuration “depol\_calibration” a ulalume3@67: product “\ *Linear polarization calibration”* to be used for the ulalume3@67: calibration. According to the example given above and to the usecase ulalume3@67: document attached we should use an usecase=4 for this example. ulalume3@67: ulalume3@67: Other “\ *Linear polarization calibration”* options to be specified are ulalume3@67: reported in figure 3.3. The most important factor you should insert here ulalume3@67: is the *Pol calibration correction factor* (K). The ideal value for this ulalume3@67: parameter is 1. Anyway if you have measured the parameter K please fill ulalume3@67: in the measured value and the corresponding measurement errors. ulalume3@67: ulalume3@67: |image2| Options for *Linear polarization calibration product*. ulalume3@67: ulalume3@67: As you can see it is possible to fill in only the K correction factor ulalume3@67: and not the calibration constant h\ :sup:`\*`. ulalume3@67: ulalume3@67: Actually for a *LIMITED* period of time it will be possible to fill in ulalume3@67: also the constant h\ :sup:`\*` using a temporary tab called ulalume3@67: *Polarization calibration constant*. This has been done to provide the ulalume3@67: users with the possibility to continue to use the SCC even if an ulalume3@67: automatic calibration made by the SCC was not submitted yet. Anyway ulalume3@67: after a transition period it will be *not* possible to provide ulalume3@67: calibration constant using this procedure and the parameter h\ :sup:`\*` ulalume3@67: can be calculated *ONLY* by the SCC as result of the submission of a ulalume3@67: proper calibration raw input dataset. The format of this input file is ulalume3@67: the same as the standard SCC input file. The only difference is that is ulalume3@67: should contain calibration measurements instead of standard ulalume3@67: measurements. Following our example, such file should contain the ulalume3@67: measurement performed at +45 and -45 degrees at 532nm. Also the channel ulalume3@67: IDs in the file should reflect the ones reported in table 3.3. ulalume3@67: ulalume3@67: Moreover this raw input file has to contain the variables: ulalume3@67: ulalume3@67: *double Pol\_Calib\_Range\_Min(channels)* ulalume3@67: ulalume3@67: *double Pol\_Calib\_Range\_Max(channels) * ulalume3@67: ulalume3@67: where to specify the altitude ranges in meters in which the polarization ulalume3@67: calibration should be done. ulalume3@67: ulalume3@67: According to the table 3.3 this file should be something similar to: ulalume3@67: ulalume3@67: dimensions: ulalume3@67: ulalume3@67: channels = 4 ; ulalume3@67: ulalume3@67: nb\_of\_time\_scales = 1 ; ulalume3@67: ulalume3@67: points = 16380 ; ulalume3@67: ulalume3@67: scan\_angles = 1 ; ulalume3@67: ulalume3@67: time = UNLIMITED ; // (3 currently) ulalume3@67: ulalume3@67: variables: ulalume3@67: ulalume3@67: int channel\_ID(channels) ; ulalume3@67: ulalume3@67: double Background\_Low(channels) ; ulalume3@67: ulalume3@67: double Background\_High(channels) ; ulalume3@67: ulalume3@67: int id\_timescale(channels) ; ulalume3@67: ulalume3@67: double Laser\_Pointing\_Angle(scan\_angles) ; ulalume3@67: ulalume3@67: int Molecular\_Calc ; ulalume3@67: ulalume3@67: int Laser\_Pointing\_Angle\_of\_Profiles(time, nb\_of\_time\_scales) ; ulalume3@67: ulalume3@67: int Raw\_Data\_Start\_Time(time, nb\_of\_time\_scales) ; ulalume3@67: ulalume3@67: int Raw\_Data\_Stop\_Time(time, nb\_of\_time\_scales) ; ulalume3@67: ulalume3@67: int Laser\_Shots(time, channels) ; ulalume3@67: ulalume3@67: double Raw\_Lidar\_Data(time, channels, points) ; ulalume3@67: ulalume3@67: double Pressure\_at\_Lidar\_Station ; ulalume3@67: ulalume3@67: double Temperature\_at\_Lidar\_Station ; ulalume3@67: ulalume3@67: double Pol\_Calib\_Range\_Min(channels) ; ulalume3@67: ulalume3@67: double Pol\_Calib\_Range\_Max(channels) ; ulalume3@67: ulalume3@67: // global attributes: ulalume3@67: ulalume3@67: :System = "mysystem" ; ulalume3@67: ulalume3@67: :Longitude\_degrees\_east = 15.723771 ; ulalume3@67: ulalume3@67: :RawData\_Start\_Time\_UT = "220000" ; ulalume3@67: ulalume3@67: :RawData\_Start\_Date = "20130620" ; ulalume3@67: ulalume3@67: :Measurement\_ID = "20130620po00" ; ulalume3@67: ulalume3@67: :Altitude\_meter\_asl = 760. ; ulalume3@67: ulalume3@67: :RawData\_Stop\_Time\_UT = "230333" ; ulalume3@67: ulalume3@67: :Latitude\_degrees\_north = 40.601039 ; ulalume3@67: ulalume3@67: data: ulalume3@67: ulalume3@67: channel\_ID = 10, 11, 12, 13 ; ulalume3@67: ulalume3@67: Background\_Low = 30000, 30000, 30000, 30000 ; ulalume3@67: ulalume3@67: Background\_High = 50000, 50000, 50000, 50000 ; ulalume3@67: ulalume3@67: id\_timescale = 0, 0, 0, 0 ; ulalume3@67: ulalume3@67: Laser\_Pointing\_Angle = 0 ; ulalume3@67: ulalume3@67: Molecular\_Calc = 0 ; ulalume3@67: ulalume3@67: Laser\_Pointing\_Angle\_of\_Profiles = ulalume3@67: ulalume3@67: 0, ulalume3@67: ulalume3@67: 0, ulalume3@67: ulalume3@67: 0 ; ulalume3@67: ulalume3@67: Raw\_Data\_Start\_Time = ulalume3@67: ulalume3@67: 0, ulalume3@67: ulalume3@67: 300, ulalume3@67: ulalume3@67: 600 ; ulalume3@67: ulalume3@67: Raw\_Data\_Stop\_Time = ulalume3@67: ulalume3@67: 210, ulalume3@67: ulalume3@67: 510, ulalume3@67: ulalume3@67: 810 ; ulalume3@67: ulalume3@67: Laser\_Shots = ulalume3@67: ulalume3@67: 1200, 1200, 1200, 1200, ulalume3@67: ulalume3@67: 1200, 1200, 1200, 1200, ulalume3@67: ulalume3@67: 1200, 1200, 1200, 1200 ; ulalume3@67: ulalume3@67: Pressure\_at\_Lidar\_Station = 1010 ; ulalume3@67: ulalume3@67: Temperature\_at\_Lidar\_Station = 14 ; ulalume3@67: ulalume3@67: Pol\_Calib\_Range\_Min = 1000, 1000, 1000, 1000 ; ulalume3@67: ulalume3@67: Pol\_Calib\_Range\_Min = 2000, 2000, 2000, 2000 ; ulalume3@67: ulalume3@67: Raw\_Lidar\_Data = …...; ulalume3@67: ulalume3@67: The file above assume the following calibration measurements have been ulalume3@67: done: ulalume3@67: ulalume3@67: 1. First +45 degrees acquisition followed by a corresponding -45 degrees ulalume3@67: acquisition ulalume3@67: ulalume3@67: a. Measurement at +45 degrees ulalume3@67: ulalume3@67: Start Time: 20130620 22:00:00 ulalume3@67: ulalume3@67: Stop Time: 20130620 22:01:00 ulalume3@67: ulalume3@67: Shots: 1200 ulalume3@67: ulalume3@67: a. Measurement at -45 degrees ulalume3@67: ulalume3@67: Start Time: 20130620 22:02:30 ulalume3@67: ulalume3@67: Stop Time: 20130620 22:03:30 ulalume3@67: ulalume3@67: Shots: 1200 ulalume3@67: ulalume3@67: 1. Second +45 degrees acquisition followed by a corresponding -45 ulalume3@67: degrees acquisition ulalume3@67: ulalume3@67: a. Measurement at +45 degrees ulalume3@67: ulalume3@67: Start Time: 20130620 22:05:00 ulalume3@67: ulalume3@67: Stop Time: 20130620 22:06:00 ulalume3@67: ulalume3@67: Shots: 1200 ulalume3@67: ulalume3@67: a. Measurement at -45 degrees ulalume3@67: ulalume3@67: Start Time: 20130620 22:07:30 ulalume3@67: ulalume3@67: Stop Time: 20130620 22:08:30 ulalume3@67: ulalume3@67: Shots: 1200 ulalume3@67: ulalume3@67: 1. Third +45 degrees acquisition followed by a corresponding -45 degrees ulalume3@67: acquisition ulalume3@67: ulalume3@67: a. Measurement at +45 degrees ulalume3@67: ulalume3@67: Start Time: 20130620 22:10:00 ulalume3@67: ulalume3@67: Stop Time: 20130620 22:11:00 ulalume3@67: ulalume3@67: Shots: 1200 ulalume3@67: ulalume3@67: a. Measurement at -45 degrees ulalume3@67: ulalume3@67: Start Time: 20130620 22:12:30 ulalume3@67: ulalume3@67: Stop Time: 20130620 22:13:30 ulalume3@67: ulalume3@67: Shots: 1200 ulalume3@67: ulalume3@67: As you can see there are 3 cycles of consecutive measurements at +45 and ulalume3@67: -45 degrees. That's way the dimension time is set to 3. ulalume3@67: ulalume3@67: The first +/-45 degrees measurement starts at “20130620 22:00:00” (start ulalume3@67: time of the first +45 measurement) and stops at “20130620 22:03:30” ulalume3@67: (stop time of the fist -45 measurement). As a consequence, according to ulalume3@67: the values of the global attributes RawData\_Start\_Date and ulalume3@67: RawData\_Start\_Time\_UT we have to set: ulalume3@67: ulalume3@67: Raw\_Data\_Start\_Time[0]=0 (start of the first +45 measurement in ulalume3@67: seconds since RawData\_Start\_Time\_UT) ulalume3@67: ulalume3@67: Raw\_Data\_Stop\_Time[0]=210 (stop of the first -45 measurement in ulalume3@67: seconds since RawData\_Start\_Time\_UT) ulalume3@67: ulalume3@67: Following a similar procedure for the other 2 cycles we have: ulalume3@67: ulalume3@67: Raw\_Data\_Start\_Time[1]=300 (start of the second +45 measurement in ulalume3@67: seconds since RawData\_Start\_Time\_UT) ulalume3@67: ulalume3@67: Raw\_Data\_Stop\_Time[1]=510 (stop of the second -45 measurement in ulalume3@67: seconds since RawData\_Start\_Time\_UT) ulalume3@67: ulalume3@67: Raw\_Data\_Start\_Time[2]=600 (start of the third +45 measurement in ulalume3@67: seconds since RawData\_Start\_Time\_UT) ulalume3@67: ulalume3@67: Raw\_Data\_Stop\_Time[2]=810 (stop of the third -45 measurement in ulalume3@67: seconds since RawData\_Start\_Time\_UT) ulalume3@67: ulalume3@67: Moreover, according to the order of the channels in the channel\_ID ulalume3@67: variable, the Raw\_Lidar\_Data array should be filled as it follows: ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[0][0][points] → 1\ :sup:`st` measured transmitted ulalume3@67: signal at +45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[0][1][points] → 1\ :sup:`st` measured reflected signal ulalume3@67: at +45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[0][2][points] → 1\ :sup:`st` measured transmitted ulalume3@67: signal at -45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[0][3][points] → 1\ :sup:`st` measured reflected signal ulalume3@67: at -45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[1][0][points] → 2\ :sup:`nd` measured transmitted ulalume3@67: signal at +45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[1][1][points] → 2\ :sup:`nd` measured reflected signal ulalume3@67: at +45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[1][2][points] → 2\ :sup:`nd` measured transmitted ulalume3@67: signal at -45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[1][3][points] → 2\ :sup:`nd` measured reflected signal ulalume3@67: at -45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[2][0][points] → 3\ :sup:`rd` measured transmitted ulalume3@67: signal at +45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[2][1][points] → 3\ :sup:`rd` measured reflected signal ulalume3@67: at +45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[2][2][points] → 3\ :sup:`rd` measured transmitted ulalume3@67: signal at -45 degrees ulalume3@67: ulalume3@67: Raw\_Lidar\_Data[2][3][points] → 3\ :sup:`rd` measured reflected signal ulalume3@67: at -45 degrees ulalume3@67: ulalume3@67: Once this file has been created it needs to be submitted to the SCC and ulalume3@67: linked to the configuration “depol\_calibration”. The result of the SCC ulalume3@67: analysis on this file will be the calculation of the calibration ulalume3@67: constant h\ :sup:`\*` that will be logged into the SCC database and can ulalume3@67: be used to calibrate Raman/Elastic backscat ter products (see section ulalume3@67: 3.3). ulalume3@67: ulalume3@67: **Definition of “Raman/Elastic backscatter and linear depolarization ratio”** ulalume3@67: ----------------------------------------------------------------------------- ulalume3@67: ulalume3@67: In order to calculate the *PLDR* we need to modify the polarization ulalume3@67: related products linked to the “standard” measurement configurations ulalume3@67: (the configuration called “nighttime” and/or “daytime” in table 3.2). ulalume3@67: ulalume3@67: Let's suppose we have defined the following products (defined already in ulalume3@67: SCC v3.11): ulalume3@67: ulalume3@67: **Table 3.4:** Example of products configuration in SCC v3.11 ulalume3@67: ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Product Name | Product ID | Product Type | nighttime | daytime | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Raman backscatter | 1 | Raman backscatter |  | | ulalume3@67: | | | | | | ulalume3@67: | 355nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Extinction | 2 | Extinction |  | | ulalume3@67: | | | | | | ulalume3@67: | 387nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Raman backscatter | 3 | Raman backscatter |  | | ulalume3@67: | | | | | | ulalume3@67: | 532nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Extinction | 4 | Extinction |  | | ulalume3@67: | | | | | | ulalume3@67: | 532nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Elastic backscatter | 5 | Elastic backscatter | |  | ulalume3@67: | | | | | | ulalume3@67: | 355nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Elastic backscatter | 6 | Elastic backscatter | |  | ulalume3@67: | | | | | | ulalume3@67: | 532nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: | Elastic backscatter | 7 | Elastic backscatter |  |  | ulalume3@67: | | | | | | ulalume3@67: | 1064nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------+-------------+-----------+ ulalume3@67: ulalume3@67: Product ID=1, 2, 4, 5, 7 do not need any modification as they do not ulalume3@67: involve polarization channels. The only product that need to be modified ulalume3@67: are the Product ID=3 and 6. To produce b532 files containing also *PLDR* ulalume3@67: we need to modify the “nighttime” and “daytime” configurations to ulalume3@67: include a product of type “Raman bakscatter and linear depolarization ulalume3@67: ratio” or “Elastic bakscatter and linear depolarization ratio” ulalume3@67: respectively. So the configuration reported in table 3.4 should be ulalume3@67: changed to match what is included in table 3.5. ulalume3@67: ulalume3@67: **Table 3.5:** The same of table 3.4 but with new product types ulalume3@67: introduced in SCC v4.0 ulalume3@67: ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Product Name | Product ID | Product Type | nighttime | daytime | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Raman backscatter | 1 | Raman backscatter |  | | ulalume3@67: | | | | | | ulalume3@67: | 355nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Extinction | 2 | Extinction |  | | ulalume3@67: | | | | | | ulalume3@67: | 387nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Raman backscatter | 10 | **Raman backscatter and linear depolarization ratio** |  | | ulalume3@67: | | | | | | ulalume3@67: | 532nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Extinction | 4 | Extinction |  | | ulalume3@67: | | | | | | ulalume3@67: | 532nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Elastic backscatter | 5 | Elastic backscatter | |  | ulalume3@67: | | | | | | ulalume3@67: | 355nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Elastic backscatter | 11 | **Elastic backscatter and linear depolarization ratio** | |  | ulalume3@67: | | | | | | ulalume3@67: | 532nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: | Elastic backscatter | 7 | Elastic backscatter |  |  | ulalume3@67: | | | | | | ulalume3@67: | 1064nm | | | | | ulalume3@67: +-----------------------+--------------+-----------------------------------------------------------+-------------+-----------+ ulalume3@67: ulalume3@67: As you can see in table 3.5, the old product IDs=3 and 6 (present in ulalume3@67: table 3.4) have been replaced with the new product ID=10 and 11 to ulalume3@67: guarantee the calculation of *PLDR*. ulalume3@67: ulalume3@67: It is important to set among the product options of the product ID=10 ulalume3@67: and 11 which calibration product we want to use for calibration (see ulalume3@67: section 3.2). This can be done using the SCC web interface setting the ulalume3@67: appropriate setting in the tab *Polarization calibration products* (see ulalume3@67: figure 3.4). According to the current example you should set here the ulalume3@67: calibration product defined in section 3.2. ulalume3@67: ulalume3@67: |image3| How to link a product to calibrate with a calibration product. ulalume3@67: ulalume3@67: **WARNING:** Please not that also *Raman/Elastic backscatter products* ulalume3@67: need to be linked to a calibration product because the calibration ulalume3@67: constant and the corresponding correction factor is needed to calculate ulalume3@67: the total signal out of the two polarization components even if the ulalume3@67: *PLDR* is not involved in the product calculation. ulalume3@67: ulalume3@67: .. |image0| image:: ./media/image1.png ulalume3@67: :width: 6.69514in ulalume3@67: :height: 2.40764in ulalume3@67: .. |image1| image:: ./media/image2.png ulalume3@67: :width: 6.69306in ulalume3@67: :height: 1.71458in ulalume3@67: .. |image2| image:: ./media/image3.png ulalume3@67: :width: 6.69306in ulalume3@67: :height: 1.77431in ulalume3@67: .. |image3| image:: ./media/image4.png ulalume3@67: :width: 6.69306in ulalume3@67: :height: 0.36389in