global_ocean.128x64x15/input/mkx60.m

%fni='bathymetry.bin'; nr=1;
%fni='lev_clim_temp.bin'; nr=15;
%fni='lev_clim_salt.bin'; nr=15;
%fni='ncep_taux.bin';nr=12;
%fni='ncep_tauy.bin';nr=12;
%fni='lev_monthly_temp.bin';nr=12;
%fni='lev_monthly_salt.bin';nr=12;

nx=128;nyi=64;cropy=2;
fid=fopen(fni,'r','ieee-be');
phi=fread(fid,nx*nyi*nr,'float32');
fclose(fid);
phi=reshape(phi,[nx nyi nr]);
phio=phi(:,1+cropy:end-cropy,:);
fno=sprintf('%s_06',fni);
fid=fopen(fno,'w','ieee-be');
fwrite(fid,phio,'float32');
fclose(fid);
1	cnh	1.1	%fni='bathymetry.bin'; nr=1;
2			%fni='lev_clim_temp.bin'; nr=15;
3			%fni='lev_clim_salt.bin'; nr=15;
4			%fni='ncep_taux.bin';nr=12;
5			%fni='ncep_tauy.bin';nr=12;
6			%fni='lev_monthly_temp.bin';nr=12;
7			%fni='lev_monthly_salt.bin';nr=12;
8
9			nx=128;nyi=64;cropy=2;
10			fid=fopen(fni,'r','ieee-be');
11			phi=fread(fid,nxnyinr,'float32');
12			fclose(fid);
13			phi=reshape(phi,[nx nyi nr]);
14			phio=phi(:,1+cropy:end-cropy,:);
15			fno=sprintf('%s_06',fni);
16			fid=fopen(fno,'w','ieee-be');
17			fwrite(fid,phio,'float32');
18			fclose(fid);