% grid_edit_step1    
%

oldgrid='model_grid';
newgrid='model_grid.new';

% use readgrid to get 2 d depth array 
[d,x,y,i,j,xx,yy,h,data,str]=readgrid(oldgrid);

% load coastline to aid in grid editing

%load('/usr/users/rps/ecohab/gom_coast2.mat');

% make the gridedit object a GLOBAL 
global grid_obj

%gridedit(x,y,d,coast);
gridedit(x,y,d);

% N.B.  When you finish editing the grid, 
% type 'grid_edit_step2' to write out the new grid

