[MITgcm-support] Matlab for post-processing!

[Dimitris Menemenlis]

Van Thinh, how big is your domain, Nx * Ny * Nz

Can you fit a single slice Nx * Ny in memory?

There is another routine that I use to read MITgcm output files:
MITgcm/utils/matlab/cs_grid/read_cs/readbin.m
also available here:
ftp://ecco2.jpl.nasa.gov/matlab/mitgcm/readbin.m

For example, if you have Nx=1000; Ny=1000; Nz=50;
you would:

for k=1:Nz
  tmp=readbin(filename,[Nx Ny],1,'real*4',k-1);
  tmp=single(tmp);    % if needed
  % do whatever you need to do with tmp
end