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	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