... | ... | @@ -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 |