cpl_gray%2Bocn/input_ocn/gendata.m


nc=32;

%-- make flat bathy (4.km deep, but real depth will be smaller if sum(drF) < 4.km)
%yy=yc*pi/90;
depth=4000.;
var=-depth*ones(6*nc,nc);
fname='flat_4km.bin';
%fid=fopen(fname,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
%fprintf(['write file: ',fname,'\n']);

%var=rdda(fname,[6*nc nc],1,'real*8','b');
%mnV=min(var(:)); MxV=max(var(:));

rDir='../cpl_idl_05.ocn_c300/'; it=2592000;
namf='Ttave';
var=rdmds([rDir,namf],it);
mnV=min(var(:)); MxV=max(var(:));
fprintf(' file: "%s" loaded: min,Max= %8.3f , %8.3f\n',namf,mnV,MxV);
fname='tIni_cpl.bin';
fid=fopen(fname,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
fprintf(['write file: ',fname,'\n']);

namf='Stave';
var=rdmds([rDir,namf],it);
mnV=min(var(:)); MxV=max(var(:));
fprintf(' file: "%s" loaded: min,Max= %8.3f , %8.3f\n',namf,mnV,MxV);
fname='sIni_cpl.bin';
fid=fopen(fname,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
fprintf(['write file: ',fname,'\n']);

return
1	jmc	1.1
2			nc=32;
3
4			%-- make flat bathy (4.km deep, but real depth will be smaller if sum(drF) < 4.km)
5			%yy=yc*pi/90;
6			depth=4000.;
7			var=-depthones(6nc,nc);
8			fname='flat_4km.bin';
9			%fid=fopen(fname,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
10			%fprintf(['write file: ',fname,'\n']);
11
12			%var=rdda(fname,[6nc nc],1,'real8','b');
13			%mnV=min(var(:)); MxV=max(var(:));
14
15			rDir='../cpl_idl_05.ocn_c300/'; it=2592000;
16			namf='Ttave';
17			var=rdmds([rDir,namf],it);
18			mnV=min(var(:)); MxV=max(var(:));
19			fprintf(' file: "%s" loaded: min,Max= %8.3f , %8.3f\n',namf,mnV,MxV);
20			fname='tIni_cpl.bin';
21			fid=fopen(fname,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
22			fprintf(['write file: ',fname,'\n']);
23
24			namf='Stave';
25			var=rdmds([rDir,namf],it);
26			mnV=min(var(:)); MxV=max(var(:));
27			fprintf(' file: "%s" loaded: min,Max= %8.3f , %8.3f\n',namf,mnV,MxV);
28			fname='sIni_cpl.bin';
29			fid=fopen(fname,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
30			fprintf(['write file: ',fname,'\n']);
31
32			return