function [bathy]=bumpexp(meshnod,dep,hh,hl,hs,hw)
% TWO EXPONENTIAL BUMPS
% [bathy]=bumpexp(meshnod,dep,hh,hl,hs,hw)
% For bathymetry 
%      dep is depth of outer plane, 
%      hh is height of bump above the plane
%     dep = 10, hh = 1 will give a bump which is 9 deep on its top.
% [bathy]=bumpexp(meshnod,depth,humpheight,humplength,humpspace,humpwidth);
% For drag coefficient 
%     cdd is drag coefficient on outer plane
%     cddbump is drag coefficient of mid point of bump
%       So insert  cdd-cddbump for  hh. 
% [meshcd]=bumpexp(meshnod,cdd,cdd-cddbump,humplength,humpspace,humpwidth);

humplength = 12000;
humpheight = 10;
depth = 20;
humpspace =6000;
humpwidth = 3000;


depth = dep;
humpheight = hh;
humplength = hl;
humpspace = hs;
humpwidth = hw;

x = meshnod(:,2)-mean(meshnod(:,2));
y = meshnod(:,3)-mean(meshnod(:,3)); 
bathy = meshnod(:,1:2);

bathy(:,2) = depth - humpheight*(...
  exp(-((x-humpspace/2)/humpwidth ).^2).*exp(-(y/humplength).^2) +...
  exp(-((x+humpspace/2)/humpwidth ).^2).*exp(-(y/humplength).^2) );

