all[]
2
1
1
Hs = 0.2; #wave height
g = 9.81; #gravity
lambda = 2.0; #wavelength
h = 0.4; #channel depth
freq = 1; #frequency
#(Hs/2)*sigma*(cosh(k*(h+z))/sinh(k*h))*sin(k*x);
#-(Hs/2)*sigma*(sinh(k*(h+z))/sinh(k*h))*cos(k*x);
#(Hs/2)*sigma*(cosh(k*(h+zcoord))/sinh(k*h))*sin(k*x+(pi/2))
# -(Hs/2)*sigma*(sinh(k*(h+zcoord))/sinh(k*h))*cos(k*x+(pi/2))
k = 2*pi/lambda;
sigma = 2*pi*freq;
zcoord = z-h;
velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;
101325
1
1273.15
1
1
0.5
0
2
0
mesh_viscosity_1 = 0.001;
20
1e-08
1
1
0
1000
0.001
1017.24
0.02495
0
0
-9.81
0
0
0
bk
bot
fr
in
out
top
0.01
200
0