global_ocean.128x60x15/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	%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);