mpi_threading_mode=funneled
# ------------------------------- Mesh ---------------------------------

metric_type = icosa_area
area_center_lon = 3
area_center_lat = 47
area_radius=800

#---------------- guided ----------------
#guided_type = none

guided_type = nudging
itau_nudging = 120
nudging_zone = area


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

# nbp                 20  40  80 160
# T-edge length (km) 500 250 120  60

# sub splitting of main rhombus : integer (default=1)
nsplit_i=4
nsplit_j=4
omp_level_size=1

# read_metric : read grid description from start.nc file 
read_metric= n

# optim_it : mesh optimisation : number of iteration : integer (default=0)
optim_it= 0

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

# timestep : real (default=480s)
dt = 30.

# advection called every itau_adv time steps : integer (default 2)
# standard : umax=100m/s vs c=340m/s (ratio 1:3)
# in JW06 umax=35m/s vs c=340m/s (ratio 1:10)
itau_adv=1

# run length : real (default=dt*itaumax)
#run_length=864000
run_length=2592000

# etat0_start_iteration_reset : set y to reset iteration from start.nc
etat0_start_iteration_reset= y

# output field period : integer (default none)
write_period=10800

itau_check_conserv=100

# ---------------------------------- Misc --------------------------------
# number of tracer (nqtot) : integer (default 1)
nqtot=4


# pression value where output is interpolated : real (default=0, no output)
out_pression_level=85000

start_file_name=start
restart_file_name=restart

# Use of oasis or not (needed to initialize and finalize coupling mode)
using_oasis= n

# ---------------------------------- Dissipation ---------------------------------

# dissipation time graddiv : real (default=5000)
tau_graddiv = 360

# number of iteration for graddiv : integer (default=1)
nitergdiv = 1

# dissipation time nxgradrot (default=5000)
tau_gradrot = 2160

# number of iteration for nxgradrot : integer (default=1)
nitergrot=2

# dissipation time divgrad (theta) (default=5000)
tau_divgrad=1080

# number of iteration for divgrad : integer (default=1)
niterdivgrad=2

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

# kind of physics : string : none, dcmip (default=none)
#physics=phys_external
physics=none
# itau_physics: dt_phys=itau_physics*dt here dt=450s => dt_phys=900s
itau_physics=15



# Following parameters are changed by options set in dynamico.card
#create_etat0_limit= y
#etat0= database
etat0 = start_file
iflag_phys= 1

# heldsz_p0 is only read if etat0=held_suarez (aquaplanet simulation)
heldsz_p0=101080
