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
================================================================================
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
