# Options of the 1D-simulation 
# http://projects.knmi.nl/rico/setup1d_composite.html
#--------------------------------------------------------------------------------------------#
# Large Scale Subsidence w [m/s] Apply the subsidence on the prognostic 
# fields of q_v, q_l, T, u and v of the model.
# Large Scale Temperature Forcing due to advection and radiation, dT/dt [K/s] 
# Large Scale Horizontal Moisture Advection dqv/dt [(g/kg)/s] 
# The zonal u-component of the geostrophic wind is decreasing with 2.0 * 10^-3 s-1
# corresponding with the observed wind. 
# The geostrophic v-component is assumed to be equal to the meridional wind v. 
# The only surface conditions that are prescribed are the SST and the mean surface pressure.
#---------------------------------------------------------------------------------------------#
forcing_type =  104
# Forcage par tsol
ok_flux_surf = n
iflag_radia=0
time_ini = 0    

tadv = 1
tadvh = 1  ! Large Scale Temp Forcing due to advection and radiation
tadvv = 0
qadv = 1 
qadvh = 1  ! Large Scale Horizontal Moisture Advection
qadvv = 0 
thadv = 0
thadvh = 0
thadvv = 0
trad = 0
forc_omega = 0
forc_u = 1   
forc_v = 1    
forc_w = 1     ! Large scale subsidence
forc_geo = 1
forc_ustar = 0
nudging_u = 0
nudging_v = 0
nudging_t = 0
nudging_q = 0
## latitude 
rlat = 18.49
## longitude 
rlon = 0.0   
#rlon = -149.0   
## Grid box area 
#airephy=25e10
airephy=1.
## fixe le cycle annuel au point vernal (degres)
##solarlong0 = 81.
##type of surface: 0=ocean,1=land,2=glacier,3=banquise
nat_surf = 0
## SST (not used if type_ts_forcing=1 in lmdz1d.F)
tsurf = 299.8
##  surface pressure
psurf = 101540.
##  surface altitude
zsurf =  0.
## rogosite
rugos = 0.0001   
##
wtsurf =  0.0
##
wqsurf =  0.0
##
albedo = 0.19
## age de la neige
agesno = 50.0    
##
restart_runoff = 0.0    
