cs_adj_advdiff/input_ad/mk_bathy4gcm.m

setMinDepth=input('setMinDepth (>0) = ? ');

h0=rdda('bathy_cs32.bin',[6*32 32],1,'real*8','b');
h1=min(0,h0);
[I J]=find(-setMinDepth<h1 & h1<0 );
fprintf('make deeper %i grid points\n',length(I));
h1(find(-setMinDepth<h1 & h1<0 ))=-setMinDepth;

namfil=['bathy_Hmin',int2str(setMinDepth),'.bin'];
fid=fopen(namfil,'w','b'); fwrite(fid,h1,'real*8'); fclose(fid);
fprintf([' write bathy on file:',namfil,' \n']);

return

rac='grid_cs32/';
set_axis 
load_cs
var=min(h1,0); var(find(var==0))=NaN;
shift=-1; cbV=0; AxBx=[-180 180 -90 90];
ccB=[-5000 100];
grph_CS(var,xcs,ycs,xcg,ycg,ccB(1),ccB(2),shift,cbV,AxBx);

1	heimbach	1.1	setMinDepth=input('setMinDepth (>0) = ? ');
2
3			h0=rdda('bathy_cs32.bin',[632 32],1,'real8','b');
4			h1=min(0,h0);
5			[I J]=find(-setMinDepth<h1 & h1<0 );
6			fprintf('make deeper %i grid points\n',length(I));
7			h1(find(-setMinDepth<h1 & h1<0 ))=-setMinDepth;
8
9			namfil=['bathy_Hmin',int2str(setMinDepth),'.bin'];
10			fid=fopen(namfil,'w','b'); fwrite(fid,h1,'real*8'); fclose(fid);
11			fprintf([' write bathy on file:',namfil,' \n']);
12
13			return
14
15			rac='grid_cs32/';
16			set_axis
17			load_cs
18			var=min(h1,0); var(find(var==0))=NaN;
19			shift=-1; cbV=0; AxBx=[-180 180 -90 90];
20			ccB=[-5000 100];
21			grph_CS(var,xcs,ycs,xcg,ycg,ccB(1),ccB(2),shift,cbV,AxBx);
22