Commit 3d025eb7 authored by Rory O'Hara Murray's avatar Rory O'Hara Murray

Remove redundant linse from turbine_area_sigma.m

parent 87c3b7e9
......@@ -31,16 +31,12 @@ else
splot = true;
end
% Turbine and sigma layer parameters
elev = 0; % water elevation above/below MSL - change this to see how the sigma layer occupation fraction changes with the tide
depth = H + elev;
dT = depth - Ht; % turbine hub depth
dT = H - Ht; % turbine hub depth
assert(dT>r, 'Turbine will stick out of water');
dLay = depth./sigLay;
zLev = [0:-dLay:-depth]';
dLay = H./sigLay;
zLev = [0:-dLay:-H]';
% what sigma layer is the hub in?
drsl = zLev+dT; % depth of hub relative to each sigma level
......@@ -56,7 +52,7 @@ if plot_fig
plot([-r r], zLev*[1 1])
xlabel('Distance (m)')
ylabel('Depth (m)')
title([num2str(depth, '%2.0f') ' m water depth'])
title([num2str(H, '%2.0f') ' m water depth'])
% draw rotor area
a=0;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment