#---------------------------------#
# Run parameters for the 1D model #                                         
#---------------------------------#

## Planetary constants
## ~~~~~~~~~~~~~~~~~~~~~~~
## NB: those are mandatory
## ~~~~~~~~~~~~~~~~~~~~~~~
## LENGTH OF A DAY in s
daysec         = 35733. 
## GRAVITY in m s-2
g              = 24.79 
##PLANETARY RADIUS in m
rad            = 7.1492E+07
## LENGTH OF A YEAR in days
year_day       = 10475.5 
## MIN DIST STAR-PLANET in AU [periastron]
periastr       = 4.95
## MAX DIST STAR-PLANET in AU [apoastron]
apoastr        = 5.46
## DATE OF PERIASTRON in days
peri_day       = 1520
## OBLIQUITY in deg
obliquit       = 3.13
## SURFACE PRESSURE in Pa
psurf          = 3.e5
#psurf          = 66000

## Time integration parameters
## ~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Initial date (in solar days,=0 at Ls=0)
day0           = 0
# Initial local time (in hours, between 0 and 24)
time           = 12 
# Number of time steps per sol
day_step       = 1 
# Number of sols to run 
ndt            = 63000
# ecriture physiq
ecritphy = 200

## Vertical levels
## ~~~~~~~~~~~~~~~
# hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels)
hybrid         = .false.
# autocompute vertical discretisation? (useful for exoplanet runs)
autozlevs      = .false.
# Ceiling pressure (Pa) ? (only used if autozlevs = .true.)
pceil          = 0.2

## Thermal properties
## ~~~~~~~~~~~~~~~~~~
# Simulate global averaged conditions ? 
global1d       = .false.
# Latitude (deg) [only used if global1d = F]
latitude       = 13.
# Solar Zenith angle (deg) [only used if global1d = T]
szangle        = 60
# Albedo of bare ground
albedo         = 0.
# Emissivity of bare ground
emis           = 1.0
# Soil thermal inertia (SI)
inertia        = 0.001
# Initial CO2 ice on the surface (kg.m-2)
co2ice         = 0

## Wind profile
## ~~~~~~~~~~~~
## zonal eastward component of the geostrophic wind (m/s)
u              = 0.
# 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
#         [Mars] ichoice=2   Savidjari profile (as Seiff but with dT/dz=cte)
#         [Mars] ichoice=3   Lindner (polar profile)
#         [Mars] ichoice=4   inversion
#         [Mars] 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           = 170
# 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.