#
#----------------------------------------------------------------------- 
# Run parameters for the 1D 'rcm1d.e' model                                          
#-----------------------------------------------------------------------                                         

## ~~~~~~~~~~~~~~~~~~~~~~~~~~~
## Time integration parameters
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Initial date 
day0=0
# Initial local time (in hours, between 0 and 24)
time=0
# Number of time steps per day
day_step=200
# Number of days to run 
ndt =673

## ~~~~~~~~~~~~~~~~~~~~
## Planetary constants
## ~~~~~~~~~~~~~~~~~~~~
# Planetary radius (m)
rad = 2575000
# Length of a year (days)
year_day=673
# Length of a day (s)
daysec = 1378890
# periastron (AU) : minimum Sun-Planet distance
periastr = 9.05426543135603
# apoastron (AU) : maximum Sun-Planet distance
apoastr = 10.0669795050736
# date of periastron (days)
peri_day = 536
# Obliquity (deg)
obliquit = 26.7
# Surface pressure (Pa)
psurf=146700
# Gravity (ms^-2)
g=1.35

## ~~~~~~~~~~~~~~~~
## Vertical levels
## ~~~~~~~~~~~~~~~~
# Hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels)
hybrid=.true
# Autocompute vertical discretisation? (useful for exoplanet runs)
autozlevs=.false
# Ceiling pressure (Pa) 
pceil = 0.135


## ~~~~~~~~~~~~~~~~~~~
## Thermal properties
## ~~~~~~~~~~~~~~~~~~~
# Simulate global averaged conditions ? 
global1d       = .true.
# Solar Zenith angle (deg) [only used if global1d = T]
szangle        = 60. 
# Latitude (deg) [only used if global1d = F]
latitude       = 0.0
# Molar mass of atmosphere (g)
mugaz=27.83246
# Specific heat capacity of atmosphere?
cpp=1038.73
# Albedo of bare ground
albedo=0.15
# Emissivity of bare ground
emis=0.95
# Soil thermal inertia (SI)
inertia=335
# Initial CO2 ice on the surface (kg.m-2)
co2ice=0

## ~~~~~~~~~~~~~
## Wind profile
## ~~~~~~~~~~~~~
# Zonal eastward component of the geostrophic wind (m/s)
u=10.
# meridional northward component of the geostrophic wind (m/s)
v=0.

## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
## Initial atmospheric temperature profile
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
# Type of initial temperature profile
#         ichoice=1   Constant Temperature:  T=tref
#         ichoice=2   Savidjari profile (as Seiff but with dT/dz=cte)
#         ichoice=3   Lindner (polar profile)
#         ichoice=4   inversion
#         ichoice=5   Seiff  (standard profile, based on Viking entry)
#         ichoice=6   constant T  +  gaussian perturbation (levels)
#         ichoice=7   constant T  + gaussian perturbation (km)
#         ichoice=8   Read in an ascii file "profile" 
ichoice=8
# Reference temperature tref (K)
tref=280
# Add a perturbation to profile if isin=1
isin=0
# peak of gaussian perturbation (for ichoice=6 or 7)
pic=26.522
# width of the gaussian perturbation (for ichoice=6 or 7)
largeur=10
# height of the gaussian perturbation (for ichoice=6 or 7)
hauteur=30.

# Some definitions for the physics, in file 'callphys.def'
INCLUDEDEF=callphys.def
