#
#  $Id: default.par,v 1.11 2004/09/22 15:13:47 schatz Exp $
#
#   Example parameter file for our dune simulation software
#


################################
# grid
#
NX = 160 #200 #500	#  number of grid columns
NY = 4 #8		#  number of grid lines
dx = 1

Nt = 20000 #200000 
Nt0 = 0000

save.every  = 1000
save.x-line = 1
save.dir    = DAT

# suppress saving of some variables:
# dontsave.<truncated file name> = 1
# for example:
dontsave.u = 1
dontsave.flux= 1
dontsave.shear= 1
dontsave.shear_pert= 1
dontsave.stall= 1
dontsave.rho = 1
dontsave.h_deposit= 1
dontsave.h_nonerod= 1
dontsave.h_sep= 1
dontsave.dhdt= 1

#salt.dt_h_max = 2000 #2300               # 10000 * dx
#salt.dh_max   = 0.         # 0.05 * dx

update.xmin = 0 #5
update.xmin_veget = 0 #300

calc.shift_back = 1
shift_back_target.cm0	= 0
shift_back_centrex.cmT	= 1

calc.x_periodic       = 0
calc.y_periodic       = 1
calc.volume.correction = 1
calc.analyze = 1
################################
# influx:  const or outflux

influx = const
calc.flux_in_eq_out   = 0

q_in = 0.0 # 0.0 #0.001

qinout.avg = 1
qinout.fact = 1
qinout.yoff = -10

q_in_min = 0.00

###############################################################################
# Wind model: const, flatrand (flat probability distribution), sine
# (sinusoidal variation) or bi (bidirectional)

wind = const

constwind.direction = 0
constwind.u* = 0.38

flatwind.du* = 0
flatwind.avgdir = 0
flatwind.ddir = 30

sinewind.du* = 0
sinewind.avgdir = 0
sinewind.ddir = 0.1

biwind.dir1= 60
biwind.dir2= -60
biwind.u*1 = 0.38
biwind.u*2 = 0.38
biwind.interval1 = 125000
biwind.interval2 = 125000

Bagnold-shields_parameter= 0.11 #Bagnold-Shields parameter to calcuNas praias e dunas, sob o efeito contínuo da água e dos ventos marinhos, a vegetação tem aspecto peculiar, destacando-se entre as espécies o campim-da-areia (remiria maritima), carrapicho-da-praia ou espinho-de-roseta (Acicarpha spathulata) e pimenteira (Cordia curassavica).late u*ft (this parameter should be a function of the grain Re.)

#################################################################################
# Wind field after Hunt et al. (HLR) (1988), Weng et al. (1991)

hlr.z0    = 1e-3 #0.5e-3  # after Rassmussen (1996) z0 = 0.045 * u*^2 / (2g);
                    # z0(u*=0.5) = 6e-4

hlr.tau_y = 1       # calculate y-components of the shear

# Dumps high frequencies by multiplying the
                    # shear (in Fourier space) with a Gaussian.
hlr.cut_k  = 2      # sigma = k_cut
###############################################################################
# paramters for separation bubble

sepbubble = parabolic

sep.angle = 20

sepbub.parabolic= 0 #1

bubble.length = 0 #5.7 #8.5 #4.5

sepbub.smooth = 4   # smooth the sep. bubble in order to get the polynoms
                       # laterally smooth. (0 = off, increasing number smooths
                       # stronger, negative numbers include ky-filtering)

sepbub.tau   = 0.05     # Shear stress pertubation below the separation bubble:
                       # tau_hlr * (1-(h_sepbub-h) / (sepbub.tau*crit.slope))

# parameters for traditional separation bubble: sepbubble = P3

sepbub.slope = 0.2 #0.25   # maximum slope at the inflection point of the 3rd
                       # order polynom. (-0.25 = 14 degree)

#sepbub.k_cut = 2.0     # smoothing for slope calculation dh0
#sepbub.filter_k= 1
#sepbub.k_cut_y = 2.0   # smoothing in y-direction

###############################################################################
# Saltation model:
# For details see Sauermann et al. (2001) or PhD, G. Sauermann (2001)

salt.g  = 9.81         # gravity
fluid_viscosity= 1.8e-5
salt.d  = 250e-6       # grain diameter 250um
salt.D  = 0.0          # diffusion constant
salt.Z/zm = 0	       # focal height Z and zm ratio

beta  = 5.7e-4         # direct entrainment according to Anderson (1991): beta = 5.7e-4
gamma = 0.2            # defines saturation transients

###############################################################################
# Vegetation model:

calc.vegetation= 0

vegetation.root= 0	#root size proportional to veget height? Answer: 0/1
vegetation.m = 0.16 	# accounts for the non-uniformity of surface stress
			# [value for creosote, Wyatt & Nickling (1997)]
vegetation.sigma = 1.5	# ratio of vegetation width to height (sigma=wv/hv)
			# [value for creosote, Wyatt & Nickling (1997)]
vegetation.beta = 200	# ratio of vegetation to surface drag coefficients
			# [value for creosote, Wyatt & Nickling (1997)]
vegetation.concent = 0.44   # concentration of plants 
		           # number of plants per m^2... dx_min=1.5
vegetation.Max.height = 1.0   # maximum vegetation height (Hv)

vegetation.init.value = 0 #1

vegetation.omega = 0	# frequency of sinusoidal variations of vegetation.t

vegetation.t = 2.1024e7	# time to reach the max height without sand erosion or
			# deposition.
vegetation.erosion = 1  	# exposition to sand erosion
vegetation.deposition= 1 	# exposition to sand deposition
vegetation.erosion.limit= 1 
vegetation.grow.eros= 1 # perturbation of the veg. growth rate due to erosion
			# 1 (no perturbation) - 0 (no growth)

###############################################################################
# Watertable model

calc.watertable = 0		# 0 (1) turns the watertable off (on)
watertable.meanheight = 2.0	# average height of the watertable level Hw
watertable.amplitude = 0	# amplitude of the sinusoidal variations of Hw
watertable.time.cycle = 1e7	# frequency of the sinusoidal variations of Hw
calc.capillary.fringe = 0	# calculation with/without capillary fringe
capillary.fringe.height = 0	# height of the capillary fringe

###############################################################################
# avalanches

avalanche = flow
aval.new.maxiter = 20 #50
aval.new.relax = 0.9 # this is the only parameter of the algorithm and controls the slope relaxing rate, for a high enough value the slope relax independently of it

# static and dynamic angle of repose:
aval.angle_repose_stat = 34
aval.angle_repose_dyn = 33

##########################################
# initial sand surface
#

Init-Surf = beach #plain
init.peaks = 0.0

# ---- flat surface, initialise with constant: Init-Surf = plain

plain.Height = 5. #2.0

# ---- aleatory surface: Init-Surf = alea

alea.Height = 3.0
alea.Fluctuations = 0.02

# ----  initialisation from file: Init-Surf = init_h ----

init_h.file = h.00000.dat

#Barchan_init/init_barchan_dx1.5.dat

#DataFiles/veget_test_1/h.20000.dat

#veget_no_sep/h.00001.dat

#veget_test_2/h.20000.dat

#Barchan_init/init_barchan_1.dat

init_h.sym  = 0
init_h.x-line = 0
init_h.plain.height= 0.0

# ---- gaussian surface: Init-Surf = gauss ----
gauss.volume = 0
#gauss.x_0   = 30
#gauss.y_0   = 400
gauss.h_0 = 4 
gauss.sigma = 0 #66.66
gauss.max.slope = 10 #15 #20
gauss.scale_y = 1.0
gauss.plain_height= 0 #1.0
gauss.plain_xslope= 0.0
gauss.plain_yslope= 0.0

# ---- beach surface: Init-Surf = beach ----

beach.h_0 = 4. #1. #2.0
#beach.x_0 = 30
#beach.y_0 = 500
#beach.Xsigma = 0
beach.sigma = 0
beach.max.slope = 10 #5
beach.scale_y = 0.0
beach.scale_x = 1.0

# ---- conical surface: Init-Surf = cone ----

cone.volume = 0
cone.h_0 = 20
cone.x_0   = 100
cone.y_0   = 80

# ---- Paraboloid: Init-Surf = paraboloid ----

parabola.theta = 34
parabola.h = 6
parabola.sym = 1

# parabola.l = 5.4*h + 10.7
# parabola.w = 5.8*h + 1.9
# parabola.b0 = -2.4*h + 7.9
# parabola.b2 = 1.46 * l / PTsqr(w) 

# parabola.x_0 = (NX-1)dx/2
# parabola.y_0 = (NY-1)dx/2 


# ---- initial surface defined by a matlab file: Init-Surf: matlab ----

matlab.file = ./zi.dat

##########################################
# initialisation of rock surface under the sand
# (names of parameters are as above but prepended by "nonerod.")

nonerod.Init-Surf = plain #init_h

nonerod.plain.Height = 1e-10

nonerod.gauss.volume = 0
nonerod.gauss.x_0   = 600
nonerod.gauss.y_0   = 90
nonerod.gauss.h_0   = 10
nonerod.gauss.sigma = 0
nonerod.gauss.Xsigma = 0.25
nonerod.gauss.max.slope = 3
nonerod.gauss.scale_y = 100
nonerod.gauss.plain_height= 0
nonerod.gauss.plain_xslope= 0.0
nonerod.gauss.plain_yslope= 0.0

nonerod.init_h.file = Barchans/Barchan_96/h.100000.dat
nonerod.init_h.sym  = 0
nonerod.init_h.x-line = 1

##########################################
# initialisation of the vegetation
# (names of parameters are as above but prepended by "veget.")

veget.Init-Surf = plain #init_h
veget.Pos.Init-Surf = plain #init_h

veget.plain.Height = 0

veget.init_h.file = veget.55000.dat
veget.init_h.sym  = 0
veget.init_h.x-line = 1

veget.Pos.plain.Height = 0

veget.Pos.init_h.file = veget_pos.55000.dat
veget.Pos.init_h.sym  = 0
veget.Pos.init_h.x-line = 1

##########################################
# fixed constants / materials, etc

rho_air    = 1.225
rho_sand   = 1650
rho_quartz = 2650

##########################################
#  material constants

fluid_density    = 1.225
grain_density    = 2650
packing		 = 0.6226	# sand density / grain density = 1 - void ratio