#---------------- Mesh ----------------

# Number of subdivisions on a main triangle : integer (default=40)
nbp = 60

# Number of vertical layers : integer (default=19)
llm = 54

# Vertical grid : [std|ncar|ncarl30] (default=std)
#disvert=read_apbp
disvert=plugin

# Mesh optimisation : number of iterations : integer (default=0)
optim_it = 1000

# Sub splitting of main rhombus : integer (default=1)
nsplit_i = 1
nsplit_j = 1

#number of openmp task on vertical level
omp_level_size=1

#---------------- Numerics ----------------

# Advection called every itau_adv time steps : integer (default=2)
itau_adv = 1

# Time step in s : real (default=480)
# dt = 720
# Alternative to specifying "dt", specify number of steps per day : day_step
day_step = 960

# Number of tracers : integer (default=1)
nqtot = 11


#---------------- Time and output ----------------

# Time style : [none|dcmip] (default=dcmip)
time_style = none

# Run length in s : real (default=??)
# run_length = 103680000
run_length = 462.369791666
# Alternative to specifying "run_length", specify number of days to run : ndays
# ndays=1

# Interval in s between two outputs : integer (default=0)
#write_period = 10800
# output every dyn time step => dt = day_length/day_step ; 88775/960=92.4739583333
write_period= 92.4739583333333


# Initial state : 
#   [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06)
etat0=isothermal
etat0_isothermal_temp=200
# pertub initial state by a factor of ...
etat0_ps_white_noise=0.01
etat0_theta_rhodz_white_noise=0.01

#etat0=start_file
#etat0_start_file_colocated=true

# start file name (default is start.nc)
# restart file name (default is restart.nc)
start_file_name=start_icosa
restart_file_name=restart_icosa

etat0_start_iteration_reset=true
etat0_start_iteration=0

# Dissipation time for grad(div) : real (default=5000)
tau_graddiv = 18000

# Exponent of grad(div) disspation : integer (default=1)
nitergdiv = 2

# Dissipation time for curl(curl) : real (default=5000)
tau_gradrot = 18000

# Exponent of curl(curl) disspation : integer (default=1)
nitergrot = 2

# Dissipation time for div(grad) : real (default=5000)
tau_divgrad = 18000

# Exponent of div(grad) disspation : integer (default=1)
niterdivgrad = 2

#---------------- Planet ----------------

INCLUDEDEF=mars_const.def

# -------------------------------- Physics -------------------------------------

# Physics package : [none|held_suarez|dcmip] (default=none)
physics = phys_external

# call physics every itau_physics (dynamical) time step
itau_physics=5
#itau_physics=10

# physics initial state (startfi.nc) if startphy_file=.true.
startphy_file=.false.