Update Getting started authored by Doriann Blain's avatar Doriann Blain
...@@ -23,7 +23,7 @@ If you downloaded the archive in the [_dist_](https://gitlab.obspm.fr/dblain/exo ...@@ -23,7 +23,7 @@ If you downloaded the archive in the [_dist_](https://gitlab.obspm.fr/dblain/exo
3. Put k-coefficients of desired absorbers inside the _data/k_coefficients_tables_ directory. You can use directly *CKC* outputs converted into HDF5 using *txt2h5*, or e.g. *[petitRADTRANS](https://petitradtrans.readthedocs.io/en/latest/content/available_opacities.html)* k-tables. If your absorber name is "absorberName", the lines file needs to be named "absorberName.ktable.exorem.h5" (e.g. H2O.ktable.exorem.h5). You can download *Exo-REM* k-tables [here](https://lesia.obspm.fr/exorem/). 3. Put k-coefficients of desired absorbers inside the _data/k_coefficients_tables_ directory. You can use directly *CKC* outputs converted into HDF5 using *txt2h5*, or e.g. *[petitRADTRANS](https://petitradtrans.readthedocs.io/en/latest/content/available_opacities.html)* k-tables. If your absorber name is "absorberName", the lines file needs to be named "absorberName.ktable.exorem.h5" (e.g. H2O.ktable.exorem.h5). You can download *Exo-REM* k-tables [here](https://lesia.obspm.fr/exorem/).
4. If you intend to use stellar spectra, you can download them e.g. [here](http://svo2.cab.inta-csic.es/theory/newov2/index.php?models=bt-settl). You **must** modify them in order to respect the *Exo-REM* data format (see step 3 of section [More precision: setup](#setup-1)). 4. If you intend to use stellar spectra, you can download them e.g. [here](http://svo2.cab.inta-csic.es/theory/newov2/index.php?models=bt-settl). You **must** modify them in order to respect the *Exo-REM* data format (see step 3 of section [More precision: setup](#setup-1)).
5. Verify if all the condensates and gases thermochemical tables are in the _data/thermochemical_tables_ directory. If you want to add more species to the chemical model, respect the same format and use "speciesName.tct.dat" as file name (e.g. "H2O.tct.dat"). 5. Verify if all the condensates and gases thermochemical tables are in the _data/thermochemical_tables_ directory. If you want to add more species to the chemical model, respect the same format and use "speciesName.tct.dat" as file name (e.g. "H2O.tct.dat").
6. Put a temperature profile as a priori inside the _inputs/atmospheres/temperature_profiles_ directory. The *Exo-REM* data format must be respected. You should have received an example of such a file with your *Exo-REM* distribution. 6. Put a temperature profile as a priori inside the _inputs/atmospheres/temperature_profiles_ directory. The *Exo-REM* data format must be respected. You can also use an *Exo-REM* HDF5 output file. You should have received an example of such a file with your *Exo-REM* distribution.
To plot the figures using the provided plot functions, you will need Python3, and the following Python packages: To plot the figures using the provided plot functions, you will need Python3, and the following Python packages:
- numpy - numpy
...@@ -109,7 +109,7 @@ By default, at the end of the calculations, the results are stored in the _outpu ...@@ -109,7 +109,7 @@ By default, at the end of the calculations, the results are stored in the _outpu
Let's take a look at our results. Go back to the _exorem_ directory and execute: Let's take a look at our results. Go back to the _exorem_ directory and execute:
```bash ```bash
python exorem_plot.py corot-4b python exorem_plot.py ./outputs/exorem/corot-4b.h5
``` ```
The figures should be in the _outputs/figures_ directory. The figures should be in the _outputs/figures_ directory.
...@@ -180,7 +180,7 @@ This time, our goal will be to have more precise results. We will use our calcul ...@@ -180,7 +180,7 @@ This time, our goal will be to have more precise results. We will use our calcul
``` ```
10. Update the retrieval parameters, we will use our previously retrieved temperature profile, and a lower number of iterations since we should be close to a solution: 10. Update the retrieval parameters, we will use our previously retrieved temperature profile, and a lower number of iterations since we should be close to a solution:
```text ```text
temperature_profile_file = 'temperature_profile_corot-4b.dat' ! a-priori temperature profile file temperature_profile_file = 'corot-4b.h5' ! a-priori temperature profile file
``` ```
and and
```text ```text
... ...
......