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

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

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

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

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


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

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

# Discretization of Coriolis term : [energy|enstophy] (default=energy)
caldyn_conserv = enstrophy


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

# Pressure level to which interpolate T : real (default=0)
out_pression_level = 85000

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

# Run length in s : real (default=??)
run_length = 864000


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

# Planetary radius : real (default=6371220)
radius = 6.37122e6

# Gravity : real (default=9.80616)
g = 9.80616

# Vertical coordinate : [eta_mass|eta_lag] (default=mass)
caldyn_eta = eta_lag

# Boussinesq approximation, requires eta_lag : [.FALSE.|.TRUE.] (default=.FALSE.)
boussinesq = .TRUE.

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

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

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

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

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

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

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