... | ... | @@ -33,7 +33,7 @@ To plot the figures using the provided plot functions, you will need Python3, an |
|
|
## Setup
|
|
|
In this example, we will simulate the atmosphere of CoRoT-4 b, a well studied planet. A good source of planetary information can be found [here](https://exoplanetarchive.ipac.caltech.edu/). We will use the parameters from Moutou et al. 2008.
|
|
|
|
|
|
Copy and edit the file _inputs/example.nml_, rename it _corot-4b.nml_. An extended description of the input parameters is available [here](Input-file).
|
|
|
Copy and edit the file _inputs/example.nml_, rename it _corot-4b.nml_. An extended description of the input parameters is available [here](Documentation).
|
|
|
|
|
|
1. Open _inputs/_corot-4b.nml_ with any notepad editor.
|
|
|
2. Edit the suffix of your future output files:
|
... | ... | @@ -98,7 +98,7 @@ Copy and edit the file _inputs/example.nml_, rename it _corot-4b.nml_. An extend |
|
|
./exorem.exe ../inputs/corot-4b.nml
|
|
|
```
|
|
|
|
|
|
During the run, keep a look in the terminal at how much the temperature vary (`dT`), at the ratio of the calculated internal radiosity over the theoretical radiosity (`J_int / (sigma * T_int_th**4)`), and at the variation of Chi^2 (`Chi2_var`). You should aim for a ratio close to 1, meaning that the current solution satisfy what you required. A low `dT` mean and `Chi2_var` mean that you are well converged (a negative `Chi2_var` mean that your current solution is better that the solution at the last iteration). You should have well converged by iteration 20, so in this case we overestimated the number of iterations we needed. This just mean that we waited longer than necessary to have the results. However, if the convergence was bad, we would like to re-run the code with a higher `n_iterations`, but this may not be sufficient. More on that [here]().
|
|
|
During the run, keep a look in the terminal at how much the temperature vary (`dT`), at the ratio of the calculated internal radiosity over the theoretical radiosity (`J_int / (sigma * T_int_th**4)`), and at the variation of Chi^2 (`Chi2_var`). You should aim for a ratio close to 1, meaning that the current solution satisfy what you required. A low `dT` mean and `Chi2_var` mean that you are well converged (a negative `Chi2_var` mean that your current solution is better that the solution at the last iteration). You should have well converged by iteration 20, so in this case we overestimated the number of iterations we needed. This just mean that we waited longer than necessary to have the results. However, if the convergence was bad, we would like to re-run the code with a higher `n_iterations`, but this may not be sufficient. More on that [here](Advice-and-troubleshooting).
|
|
|
|
|
|
By default, at the end of the calculations, the results are stored in the _outputs/exorem_ directory.
|
|
|
|
... | ... | |