Changes
Page history
Update Getting started
authored
Dec 04, 2020
by
Doriann Blain
Show whitespace changes
Inline
Side-by-side
Getting-started.md
View page @
416d3244
...
...
@@ -115,6 +115,11 @@ python exorem_plot.py ./outputs/exorem/corot-4b.h5
```
The figures should be in the _outputs/figures_ directory.
Alternatively, if you set
`output_hdf5 = False`
to have ASCII files, you need to execute:
```
bash
python exorem_plot.py corot-4b
--inputs_extension
dat
```
The temperature profile figure (temperature_profile_corot-4b) should look like this:

...
...
@@ -201,7 +206,7 @@ As before, go to the _bin_ directory and execute `./exorem.exe ../inputs/corot-4
## Plotting
Again, let's take a look at our results. Go back to the _exorem_ directory and execute:
```
bash
python exorem_plot.py corot-4b_R500
python exorem_plot.py
./outputs/exorem/
corot-4b_R500
.h5
```
The temperature profile figure (temperature_profile_corot-4b) should look like this:
...
...
@@ -226,7 +231,7 @@ Not happy with the figures you get ? What if for example you wanted to see the c
3.
Then, simply do:
```
python
plot_contribution_transmission_spectra
(
'
./outputs/exorem/
spectra_
corot-4b_R500.
dat
'
,
'
./outputs/exorem/corot-4b_R500.
h5
'
,
wvn2wvl
=
True
,
xmin
=
0.5e-6
,
xmax
=
1.5e-6
,
...
...
@@ -235,9 +240,12 @@ Not happy with the figures you get ? What if for example you wanted to see the c
```
4.
If you want to see the effect of the clouds on the transmission spectrum, you can try:
```
python
plot_transmission_spectrum
(
'
./outputs/exorem/spectra_corot-4b_R500.dat
'
,
wvn2wvl
=
True
,
star_radius
=
814e6
,
cloud_coverage
=
0
,
color
=
'
r
'
,
label
=
'
no cloud
'
)
plot_transmission_spectrum
(
'
./outputs/exorem/spectra_corot-4b_R500.dat
'
,
wvn2wvl
=
True
,
star_radius
=
814e6
,
cloud_coverage
=
0.5
,
color
=
'
g
'
,
label
=
'
50% cover
'
)
plot_transmission_spectrum
(
'
./outputs/exorem/spectra_corot-4b_R500.dat
'
,
wvn2wvl
=
True
,
star_radius
=
814e6
,
cloud_coverage
=
1
,
color
=
'
b
'
,
label
=
'
full cover
'
)
plot_transmission_spectrum
(
'
./outputs/exorem/example_R50_beta8.h5
'
,
wvn2wvl
=
True
,
cloud_coverage
=
0
,
color
=
'
r
'
,
label
=
'
no cloud
'
)
plot_transmission_spectrum
(
'
./outputs/exorem/example_R50_beta8.h5
'
,
wvn2wvl
=
True
,
cloud_coverage
=
0.5
,
color
=
'
g
'
,
label
=
'
50% cover
'
)
plot_transmission_spectrum
(
'
./outputs/exorem/example_R50_beta8.h5
'
,
wvn2wvl
=
True
,
cloud_coverage
=
1
,
color
=
'
b
'
,
label
=
'
full cover
'
)
plt
.
xlim
([
0.3e-6
,
2e-6
])
plt
.
gca
().
set_xscale
(
'
log
'
)
plt
.
legend
()
...
...
...
...