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

#### Time integration parameters
#
# Initial date (in martian sols ; =0 at Ls=0)
day0=1500
# Initial local time (in hours, between 0 and 24)
time=0
# Number of time steps per sol
day_step=10
# Number of sols to run
ndt = 2000
ecritphy=20

cpp_mugaz_mode = 2

#### Physical parameters
#
# Day length, in s
daysec = 551854.0625
# length of a year, in days
year_day = 14178.3037109375
# Radius of the planet, in m
rad = 1187000
# Surface pressure (Pa)
psurf=1.1
# Gravity (ms^-2)
g=0.617188990116119
# Molar mass of atmosphere (g)
mugaz=28
# Specific heat capacity of atmosphere (J kg-1 K-1) # ignored if cpp_mugaz_mode = 1
cpp=1040.0
# latitude (in degrees)
latitude=0.0
# orbital distance at perihelion (AU)
periastr=29.65831
# orbital distance at aphelion (AU)
apoastr=49.3050467
# date of periastron, in days
peri_day = 87.3600006103516
# obliquity (degrees)
obliquit=119.591

# Albedo of bare ground
albedo=0.5
# Emissivity of bare ground
emis=0.85
# Soil thermal inertia (SI)
inertia=20
# zonal eastward component of the geostrophic wind (m/s)
u=0
# meridional northward component of the geostrophic wind (m/s)
v=0.
# Initial N2 ice on the surface (kg.m-2)
n2ice=0
# hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels)
hybrid=.false.

# CH4 VMR (%)
ch4ref=0.5
# CO VMR (%)
coref=0.046
# haze ref (kg/kg)
hazeref=0.000000385
# prec haze ref (kg/kg)
prechazeref=0.0000004

# read haze from profhaze
lecthaze = 1

###### 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=1
ichoice=8
# Reference temperature tref (K)
tref=38
# 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