----------------------------------------------------
Grid input
nx = 205 
ny  = 59 
xori = 493324.42 
yori = 5676182.02 
alfa = -51.15 
depfile = bathy_2d.dep 
vardx = 1 
xfile = x_2d.dep 
yfile = y_2d.dep 
posdwn = -1 
----------------------------------------------------
General constants
rho=1025
g=9.81
----------------------------------------------------
Wave discretisation
thetanaut=0
thetamin = -57.5
thetamax = 57.5
dtheta = 15
----------------------------------------------------
Numerics input
CFL = 0.7
eps = 0.01
hmin = 0.2
gammax = 5
umin = 0.0
Hwci = 1.0
----------------------------------------------------
Boundary condition options
front = 1
back = 1
left = 0
right = 0
----------------------------------------------------
Tidal boundary condition 
tideloc = 1
zs0file = tideinput.txt
tidelen = 99
paulrevere = 0
----------------------------------------------------
Wave boundary conditions
instat = 41
Required for instat 3,4,5,7, 41
bcfile = waveinput_2d.txt
taper = 100
----------------------------------------------------
Wave calculation options
form = 2
wci = 1
turb = 2
swtable = RF_table.txt
Tsmin = 1
tsfac = 0.1
facua = 0.1
dzsdx_cr = 0.1
facthr = 7.5
 ----------------------------------------------------
Flow calculation options
C=60.
nuh=0.1
----------------------------------------------------
Sediment transport calculation options
dico = 1
D50 = 0.000214
D90 = 0.000300
struct=1
ne_layer=hardlayer_2d.dep
----------------------------------------------------
Morphological calculation options
morfac = 1.0
---------------------------------------------------- 
Time input
tstart = 0
#tsglobal = timeseries.txt
tint  = 172800
tstop = 172800
----------------------------------------------------
nglobalvar = 1
zb