Free surface flows 121
Computer code 5.5
% Example 5.5 Wind Circulation in Thermaikos Gulf
% ho = Water depth [m];
% u = Velocity in the x-direction [m/s];
% v = Velocity in the y-direction [m/s];
% wx = Wind velocity along the x-axis [m/s];
% wy = Wind velocity along the y-axis [m/s];
% dif = Lower limit of the eddy viscosity;
% cf = Coriolis coefficient [1/s];
% fs = Water surface friction coefficient;
% fb = Bed friction coefficient;
% sp = Smagorinsky parameter for horizontal eddy coefficent
[dimensionless];
% Dd = Spatial step in both x and y directions [m];
% Dt = Time step [s];
% nx = Number of computational steps in the x-direction;
% ny = Number of computational steps in the y-direction;
22
20
18
16
14
12
10
8
6
4
2
2
4
6
8
10
12
14
16
18
20
22
Figure 5.18 Wind-induced current circulation pattern.