J/A+A/562/A133 Analytical model for irradiated atmospheres (Parmentier+, 2014)
A non-grey analytical model for irradiated atmospheres. I: Derivation.
Parmentier V., Guillot G.
<Astron. Astrophys. 562, A133 (2014)>
=2014A&A...562A.133P 2014A&A...562A.133P
ADC_Keywords: Models, atmosphere ; Planets
Keywords: radiative transfer - planets and satellites: atmospheres -
stars: atmospheres - planetary systems
Abstract:
We present a new analytical model to calculate the thermal structure
of an irradiated atmosphere. Conversely to previous semi-grey model
this model takes into account two thermal opacities and two visible
opacities. The resulting temperature profiles are more realistic than
the ones computed from semi-grey models.
Description:
The model have six parameters to describe the opacities:
- Kappa(N) is the Rosseland mean opacity at each levels of the atmosphere
it does not have to be constant with depth
- Gp is the ratio of the thermal Plank mean opacity to the thermal
Rosseland mean opacity
- Beta is the width ratio of the two thermal bands in the frequency space
- Gv1 is the ratio of the visible opacity in the first visible band to
the thermal Rosseland mean opacity
- Gv2 is the ratio of the visible opacity in the second visible band to
the thermal Rosseland mean opacity
- Betav is the width ratio of the two visible band in the frequency space
Additional parameters describe the physical setting:
- Tirr is the irradiation temperature, given by the stellar flux
- mu is the angle between the vertical direction and the stellar direction
- Tint is the internal temperature, given by the internal luminosity
- P(i) are the pressure levels where the temperature is computed
- grav is the gravity of the planet
- N is the number of atmospheric levels
The code and all the outputs uses SI units.
Installation and use : to install the code use the command "make".
The input parameters must be changed inside the file PaperI.f90.
It is necessary to compile the code again each time.
The subroutine Tprofile.f90 can be directly implemented into one's code.
To launch the code, launch the executable file NonGrey.
The output is in the file PTprofile.csv
File Summary:
--------------------------------------------------------------------------------
FileName Lrecl Records Explanations
--------------------------------------------------------------------------------
ReadMe 80 . This file
paper1.f90 120 61 Main program
tprofile.f90 512 13 Subroutine calculating the profile
test.csv 129 101 Example of output
Makefile 219 52 Makefile
--------------------------------------------------------------------------------
Acknowledgements:
Vivien Parmentier, vivien.parmentier(at)oca.eu
(End) Patricia Vannier [CDS] 26-Nov-2013