[MITgcm-support] Matlab for post-processing!

Dimitris Menemenlis menemenlis at sbcglobal.net
Thu Apr 5 17:36:47 EDT 2007


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



More information about the MITgcm-support mailing list