Update Getting started authored by Doriann Blain's avatar Doriann Blain
......@@ -104,6 +104,7 @@ This is nice, but the resolution is quite low.
# 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.
## Setup
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`.
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
wavenumber_max = 50010 ! (cm-1) last wavenumber
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
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_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
and
```text
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:
```text
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
```
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