mpi_threading_mode=funneled
INCLUDEDEF = const.def

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

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

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

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


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

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

omp_level_size=1

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

# Time step in s : real (default=480)
dt = 300

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

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

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

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

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

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


#---------------- Time ----------------

# Run length in s : real (default=??)
#            30 days
#run_length = 2592000
#run_length = 604800
run_length = 10800

# Interval in s between two outputs : integer (default=??)
write_period = 10800


#---------------- Physical parameters ----------------

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

# Initial state : 
#   [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06)
#etat0 = dcmip2016_baroclinic_wave
etat0 = start_file

restart_file_name=test_restart

itau_check_conserv=1

itau_physics=3
physics=dcmip2016
#physics_dcmip2016=moist_baroclinic
physics_dcmip2016=dry_baroclinic
physics_dcmip2016_PBL=no
#caldyn_thermo=theta_fake_moist
thermo=theta