

fname = ('../frc2008/swf_grid6.nc');
; %copy of swfl_grid4.nc
 
% h_m = nc_varget(fname,'h');
% h_m=h_m.*0.3048; %convert feet to meters

h= nc_varget(fname,'h');
mask = nc_varget(fname,'mask_rho');
hplot=h;
hplot(~mask) = nan;

figure;
pcolor(hplot)
colorbar


% % %increase depth upstream Calusa
% % 
% x = [-81.7471 -81.7617 -81.8085 -81.8091 -81.8259 -81.8394 -81.7545 ...
%     -81.7088 -81.7564 -81.8412 -81.6888];
% y = [ 26.7185 26.7169 26.6994 26.6995 26.6939 26.6921 26.7226 ...
%     26.7271 26.7268 26.6966 26.7318];
% z = [4 5.75 5.75 5 4 3 5 6.5 5 3 6.5];
% 
% [sortx,ix]=sort(x);
% zsort=z(ix);
% figure;plot(sortx,zsort,'*-')

% 
% %increase depth in pass north of lovers key
% 
% 
% i=[133 133 134 134 135 135 136 136 137 137 138 138 139 139 140 140 141 141];
% j=[198 199 198 199 198 199 198 199 198 199 198 199 198 199 198 199 198 199];
% z=[5*ones(1,14) 3 4.5 2 3];
% 
% %increase depth in pass between captiva and Sanibel
% 
% i=[i 88 89 89 90 90 91 91 92 92 93 93 94 94 95 95 96 96];
% j=[j 129 128 129 128 129 128 129 128 129 128 129 128 129 128 129 128 129];
% z=[z 5*ones(1,13) 4 4 3 3];
% 
% 
% %widen channel between north captiva
% 
% i=[i 90 91 92 93 94 95 96 97 98 99 100 101 90 91 92 93];
% j=[j 108 108 108 108 108 108 108 108 108 108 108 108 109 109 109 109];
% z=[z 6.5*ones(1,10) 6 5 4 4 4 3 ];
% 
% 
% for m=1:length(i),
%     h(i(m),j(m))=z(m);
% end
% 
% % h_m(i,j)=z;
% nc_varput(fname,'h',h);
% 
% 
% hplot=h;
% hplot(~mask) = nan;
% 
% figure;
% pcolor(hplot)
% colorbar
% title('fixed')
% 
