... | @@ -20,7 +20,7 @@ Reminder: *Exo-REM* is not made to simulate highly-irradiated planets (equilibri |
... | @@ -20,7 +20,7 @@ Reminder: *Exo-REM* is not made to simulate highly-irradiated planets (equilibri |
|
## Cookbook when in trouble
|
|
## Cookbook when in trouble
|
|
When you have trouble to reach convergence (for example, if the ratio `J_int / J_th` oscillates or if the code crashes), you can do, in that order:
|
|
When you have trouble to reach convergence (for example, if the ratio `J_int / J_th` oscillates or if the code crashes), you can do, in that order:
|
|
|
|
|
|
1. Set `retrieval_flux_error_bottom = 1e-3` and `retrieval_flux_error_top = 1e-5`. Decrease `weight_apriori` roughly by one order of magnitude. If you crash, decrease by another order of magnitude. Repeat until you can reach the end of your run. Then, Increase `weight_apriori` again to 10. Repeat until you are satisfied.
|
|
1. Set `retrieval_flux_error_bottom = 1e-3` and `retrieval_flux_error_top = 1e-5`. Decrease `weight_apriori` roughly by one order of magnitude. If you crash, decrease by another order of magnitude. Repeat until you can reach the end of your run. Take your result temperature profile as input, then increase `weight_apriori` again to 10. Repeat until you are satisfied.
|
|
2. Increase `retrieval_flux_error_bottom`, but keep it lower than 1e-2. When satisfied, set it back to at must 1e-3 and redo a run.
|
|
2. Increase `retrieval_flux_error_bottom`, but keep it lower than 1e-2. When satisfied, set it back to at must 1e-3 and redo a run.
|
|
3. Increase `retrieval_flux_error_top`, but keep it lower than 1e-2. When satisfied, set it back to at must 1e-5 and redo a run.
|
|
3. Increase `retrieval_flux_error_top`, but keep it lower than 1e-2. When satisfied, set it back to at must 1e-5 and redo a run.
|
|
4. Find, or make-up a temperature profile that is close to what you expect, then use it as input.
|
|
4. Find, or make-up a temperature profile that is close to what you expect, then use it as input.
|
... | | ... | |