Update Getting started authored by Doriann Blain's avatar Doriann Blain
...@@ -104,6 +104,7 @@ This is nice, but the resolution is quite low. ...@@ -104,6 +104,7 @@ This is nice, but the resolution is quite low.
# More precision ! # More precision !
This time, our goal will be to have more precise results. We will use our calculated temperature profile as input, and a higher resolution power. We will also add a stellar spectrum, and use an advanced mode to calculate the eddy diffusion coefficient. To keep it simple, we will consider only KCl and Na2S clouds. This time, our goal will be to have more precise results. We will use our calculated temperature profile as input, and a higher resolution power. We will also add a stellar spectrum, and use an advanced mode to calculate the eddy diffusion coefficient. To keep it simple, we will consider only KCl and Na2S clouds.
## Setup
1. Download the R500 compressed *k*-tables [here](https://gitlab.obspm.fr/dblain/exorem/-/tree/master/data/k_coefficients_tables). 1. Download the R500 compressed *k*-tables [here](https://gitlab.obspm.fr/dblain/exorem/-/tree/master/data/k_coefficients_tables).
2. Decompress them inside the _data/k_coefficients_tables_ directory executing e.g. `tar xJvf R500.tar.xz R500`. 2. Decompress them inside the _data/k_coefficients_tables_ directory executing e.g. `tar xJvf R500.tar.xz R500`.
3. It is always better to use a stellar spectrum rather than a blackbody spectrum. Here we will download a [BT-Settl](http://svo2.cab.inta-csic.es/theory/newov2/index.php) spectrum model. Take T_eff = 6200 K, Log(g) = 4.5, and a metallicity of 0. Put the file into the _data/stellar_spectra_ directory and rename it e.g. "spectrum_BTSettl_6200K_logg4.5_met0.dat" (mind the .dat extension). Replace the header of the file by the following : 3. It is always better to use a stellar spectrum rather than a blackbody spectrum. Here we will download a [BT-Settl](http://svo2.cab.inta-csic.es/theory/newov2/index.php) spectrum model. Take T_eff = 6200 K, Log(g) = 4.5, and a metallicity of 0. Put the file into the _data/stellar_spectra_ directory and rename it e.g. "spectrum_BTSettl_6200K_logg4.5_met0.dat" (mind the .dat extension). Replace the header of the file by the following :
...@@ -141,10 +142,10 @@ This time, our goal will be to have more precise results. We will use our calcul ...@@ -141,10 +142,10 @@ This time, our goal will be to have more precise results. We will use our calcul
wavenumber_max = 50010 ! (cm-1) last wavenumber wavenumber_max = 50010 ! (cm-1) last wavenumber
wavenumber_step = 20 ! (cm-1) wavenumber step size wavenumber_step = 20 ! (cm-1) wavenumber step size
``` ```
9. Add the required cloud information. We will use the "fixed sedimentation parameter mode", based on a model from [Ackerman and Marley, 2001](https://iopscience.iop.org/article/10.1086/321540/meta), and a cloud coverage of 50%: 9. Add the required cloud information. We will use the "fixed sedimentation parameter mode", based on a model from [Ackerman and Marley, 2001](https://iopscience.iop.org/article/10.1086/321540/meta), and a cloud coverage of 15%:
```text ```text
cloud_mode = 'fixedSedimentation' ! cloud mode ('fixedRadius'|'fixedSedimentation'|'fixedRadiusCondensation'|'fixedRadiusTime') cloud_mode = 'fixedSedimentation' ! cloud mode ('fixedRadius'|'fixedSedimentation'|'fixedRadiusCondensation'|'fixedRadiusTime')
cloud_fraction = 0.5 ! cloud cover fraction cloud_fraction = 0.15 ! cloud cover fraction
cloud_names = 'KCl', 'Na2S' ! condensing species forming the clouds cloud_names = 'KCl', 'Na2S' ! condensing species forming the clouds
cloud_particle_radius = 5e-6, 5e-6 ! (m) mean radius of the cloud particles (fixed radius modes) cloud_particle_radius = 5e-6, 5e-6 ! (m) mean radius of the cloud particles (fixed radius modes)
...@@ -160,11 +161,15 @@ This time, our goal will be to have more precise results. We will use our calcul ...@@ -160,11 +161,15 @@ This time, our goal will be to have more precise results. We will use our calcul
and and
```text ```text
n_iterations = 20 ! number of iterations n_iterations = 20 ! number of iterations
n_non_adiabatic_iterations = 0 ! number of iterations without including the adiabatic correction (necessary for convergence) n_non_adiabatic_iterations = 10 ! number of iterations without including the adiabatic correction (necessary for convergence)
``` ```
11. Finally, update your paths so that you load the right k-tables and the right temperature profile: 11. Finally, update your paths so that you load the right k-tables and the right temperature profile:
```text ```text
path_k_coefficients = '../data/k_coefficients_tables/R500/' ! path to the k coefficients files path_k_coefficients = '../data/k_coefficients_tables/R500/' ! path to the k coefficients files
path_temperature_profile = '../outputs/exorem/' ! path to the a-priori temperature profile file path_temperature_profile = '../outputs/exorem/' ! path to the a-priori temperature profile file
``` ```
12. This should take longer than before, because we added clouds and because of the increase in resolution.
\ No newline at end of file ## Running
As before, go to the _bin_ directory and execute `./exorem.exe ../inputs/corot-4-b.nml`. This should take longer than before, because we added clouds (2 times longer) and because of the increase in resolution (10 times longer !).
## Plotting
\ No newline at end of file