... | ... | @@ -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/).
|
|
|
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").
|
|
|
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:
|
|
|
- numpy
|
... | ... | @@ -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:
|
|
|
```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.
|
|
|
|
... | ... | @@ -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:
|
|
|
```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
|
|
|
```text
|
... | ... | |