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

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

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

# kappa=Rd/cpp : real (default=0.2857143)
kappa = 0.2857143

# Heat capacity at constant pressure : real (default=1004.70885)
cpp = 1004.5

# Reference pressure : real (default=101325)
preff = 1e5

# Planetary rotation rate : real (default=7.292E-5)
omega = 0


#---------------- Constants ----------------

# Reference temperature T0 for DCMIP test cases : real (default=300)
ncar_T0 = 300

# Reference pressure p0 for DCMIP test cases, used by disvert if disver=ncar : real (default=1e5)
ncar_p0 = 1e5

# Exponent for B(eta), used by disvert if disvert=ncar : integer (default=1)
ncar_disvert_c = 1


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

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

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

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

# Model layer thickness in m, used if disvert=ncar : real (default=300)
ncar_dz = 200

# 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 ----------------

# Equations solved : [gcm|adv] (default=gcm)
caldyn = adv

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


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

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

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


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

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

# Type of guiding terms : [none|dcmip1] (default=none)
guided_type = dcmip1

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

# DCMIP 1.1 advection test, initial tracer field : 
#    [const|slotted_cyl|cos_bell|dbl_cos_bell_q1|dbl_cos_bell_q2|complement|hadley|dcmip11] (default=cos_bell)
dcmip1_shape = dcmip11

# DCMIP 1.1 advection test, wind field : [solid|deform|hadley] (default=deform)
dcmip1_wind = deform