#
#----------------------------------------------------------------------- 
#run control parameters:                                          
#------------------------------                                          

# planet type
planet_type = generic

# Number of days to run model for 
#     nday=370
nday=1

# ndynstep (overrides nday): number of dynamical steps to run for
#     ndynstep=60

# Number of dynamical steps per day (must be a multiple of iperiod) (originally, 900 or 960)
 day_step = 960

# Apply a Matsuno step every iperiod dynamical step
  iperiod=5

# Control output information in the dynamics every iconser dynamical steps
  iconser=120

# Apply dissipation every idissip dynamical steps (originally, idissip=5)
  idissip=5
#  dissip_period=1

# dissipation operator to use (star or non-star)
 lstardis=.true.

# use hybrid vertical coordinate (else will use sigma levels)
 hybrid=.true.

# iterate lateral dissipation operator gradiv nitergdiv times
nitergdiv=1

# iterate lateral dissipation operator nxgradrot nitergrot times
nitergrot=2

# iterate lateral dissipation operator divgrad niterh times
   niterh=2

# time scale (s) for shortest wavelengths for u,v (gradiv)
 tetagdiv=4000.

# time scale (s) for shortest wavelengths for u,v (nxgradrot)
 tetagrot=5000.

# time scale (s) for shortest wavelengths for h (divgrad)
 tetatemp=5000.

# coefficient for gamdissip                                            
  coefdis=0.

# time marching scheme (Matsuno if purmats is true, else Matsuno-Leapfrog)
  purmats=.false.

# run with (true) or without (false) physics
   physic=.true.

# call physics every iphysiq dynamical steps (originally, iphysiq=5)
  iphysiq=10

# Use a regular grid
  grireg=.true.

# longitude (degrees) of zoom center
   clon=0.

# latitude (degrees) of zoom center
   clat=0.

# enhancement factor of zoom, along longitudes
  grossismx=1.

# enhancement factor of zoom, along latitudes
 grossismy=1.

#  Use an hyperbolic function f(y) if .true., else use a sine
  fxyhypb=.false.

# extention along longitudes of zoom region (fraction of global domain)
   dzoomx= 0.

# extention along latitudes of zoom region (fraction of global domain)
   dzoomy=0.

# zoom stiffness along longitudes
    taux=2.

# zoom stiffness along latitudes
    tauy=2.

#  Function  f(y) as y = Sin(latitude) if = .true. ,  else  y = latitude
  ysinus= .false.

# Use a sponge layer
  callsponge  = .true.
 
# Sponge layer extends over topmost nsponge layers
  nsponge = 3

# Sponge:  relaxed towards teta=teta_ave
#        and mode0(u=v=0), mode1(u=u_ave,v=0), mode2(u=u_ave,v=v_ave)
  mode_sponge= 2

# Sponge layer time scale (s):  tetasponge
  tetasponge = 5000

# some definitions for the physics, in file 'callphys.def'
INCLUDEDEF=callphys.def
