RT Dataset -- Distribution and intensity of water ice signature in South Xanadu and Tui Regio

This dataset contains all Radiative Transfer (RT) results made for the paper:

    > Coutelier et al. 2021 - Icarus.
    > Distribution and intensity of water ice signature in South Xanadu and Tui Regio
    > doi: 10.1016/j.icarus.2021.114464

The data are stored in 4 zipped-folders named with the VIMS cube id
and the tested model (when available):

    - C1481624349_1_FM_diff_FH.zip
    - C1481624349_1_FM_eq_FH.zip
    - C1590648776_1.zip
    - C1809727868_1.zip

In these folders, each VIMS pixel is stored in 2 csv-files, with the following patterns:

    1. <CUBE_ID>_<PIXEL_SAMPLE>-<PIXEL_LINE>[_<TESTED_MODEL>]_spectra.csv
    2. <CUBE_ID>_<PIXEL_SAMPLE>-<PIXEL_LINE>[_<TESTED_MODEL>]_albedo.csv

The first one, starts with a commented header describing the radiative transfer model inputs
and outputs:

    - CUBE_ID: the VIMS cube id (`C1234567890_1` format)
    - SAMPLE: the pixel sample number.
    - LINE: the pixel line number.
    - NKS: the number of correlated-k values used in the RT model.
    - MAX_VAL_NLEGP: the maximum number Legendre polynomial coefficient in the phase function.
    - SURFACE: the surface reflection type (Lambert in our case).
    - SOLARFLUX: the incident flux (adjusted from the incident angle).
    - ALT_MAX_M: the max altitude of the model grid (in meter).
    - RADIUS_MONO: The aerosols monomer radius (in meter).
    - N_MONO: the number of monomers.
    - DF: the aerosols fractal dimension.
    - FH: the haze scaling factor.
    - FM: the mist scaling factor.
    - HAZE_ALT_TRANS: the altitude of transition between mist and haze (in meter).
    - X2: the quadratic error between the data and the output.
    - PHASE: the surface phase angle (in degree).
    - EMI: the surface emergent angle (in degree).
    - INC: the surface incident angle (in degree).
    - SOLAR_MU: Opposite of the cosine of the incident angle (SHDOMPP). Cosine of the incident angle (SPSDISORT)
    - MU_OUT: the cosine of emergent angle.
    - PHI_OUT: the azimuth of the emergent angle (in degrees).

Then, it contains the original and the retrieved spectra:

    - Column 1: the VIMS band channel number.
    - Column 2: the VIMS channel central wavelength (in meter).
    - Column 3: the original VIMS pixel I/F.
    - Column 4: the RT model retrieved I/F.

The second file, contains only the surface albedo retrieved with the RT model with its associated error:

    - Column 1: the VIMS channel central wavelength (in meter).
    - Column 2: the original surface albedo retrieved (negative values were set at 0).
    - Column 3: the adjusted surface albedo retrieved with wider atmospheric windows. [*]
    - Column 4: the upper error value of the surface albedo retrieved in the model.
    - Column 5: the lower error value of the surface albedo retrieved in the model.

[*]: In order to keep a record of edge effects. All null values were replaced with NaN but
the negative values are kept.

Notes
-----

For the VIMS cube `C1481624349_1`, we performed two different tests:

    1. with different haze and mist parameters FH == FM (default).
    2. with the same haze and mist parameters FH != FM.

The radiative transfer retrievals were performed with `SHDOMPP` solver except
on the VIMS cube `C1809727868_1` where we used `SPSDISORT` solver (suffixed with `_SPS` in the filename).

Finally, we also provide a Jupyter Notebook (`Display_cubes.ipynb`) to help visualize the data.
This file can be open locally or with Jupyter Notebook Viewer at: https://nbviewer.jupyter.org

We provide as well the three calibrated VIMS-IR cube files ('C<CUBE_ID>_ir.cub') we worked with.