blob: a08129254a6c4000cd169c7814931d13d183c4fa (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
begin thorn FunwaveSources
# parameters section?
inherit Funwave
begin variables
end variables
begin temporaries
end temporaries
begin derivatives
Dp u_i = (u_(i+1) - u_(i))/h
Dz v_i = (v_(i+1) - v_(i-1))/(2*h)
end derivatives
begin calculation sourceterms scheduled in funw_source
sourcex = gamma1*mask9*(-u1p + u4)*(Dp_1 p + Dp_2 q) +
gamma1*mask9*(-(gamma2*mask9*(u1pp + u2 + u3)) + u4*Dz_1 u + u*Dz_1 u4 +
v4*Dz_2 u + v*Dz_2 u4)*max(h,mindepthfrc) -
cd*u*sqrt(u**2 + v**2)
end calculation
end thorn
|