matlab_class/gcmfaces_IO/read_bin.m

function [fld]=read_bin(fileName,varargin);
%usage: fld=read_bin(fileName);         reads 3D field (N levels in 3rd dim)
%usage: fld=read_bin(fileName,trec);    reads 3D field, at trec (50 levels)  
%usage: fld=read_bin(fileName,trec,k);  reads 2D field, at level k, at trec (50 levels)
%usage: fld=read_bin(fileName,trec,0);  reads 2D field, at trec (1 level)

global mygrid;

n1=mygrid.ioSize(1); n2=mygrid.ioSize(2); n3=length(mygrid.RC); nF=mygrid.nFaces;
recl2D=n1*n2*4; recl3D=n1*n2*n3*4; myprec='float32';
%recl2D=n1*n2*8; recl3D=n1*n2*n3*8; myprec='float64';

if nargin==1; n3=dir(fileName); n3=n3.bytes/recl2D; recl3D=recl2D*n3; end; 
if nargin>1; tt=varargin{1}; else; tt=1; end;
if nargin>2; kk=varargin{2}; else; kk=[]; end;

fid=fopen(fileName,'r','b');
if ~isempty(kk);
  if kk==0;%meaning the field is only two dimensional
      status=fseek(fid,(tt-1)*recl2D,'bof');
      fld=reshape(fread(fid,n1*n2,myprec),[n1 n2]);
  else;
      status=fseek(fid,(tt-1)*recl3D+(kk-1)*recl2D,'bof');
      fld=reshape(fread(fid,n1*n2,myprec),[n1 n2]);
  end;
else;
  status=fseek(fid,(tt-1)*recl3D,'bof');
  fld=reshape(fread(fid,n1*n2*n3,myprec),[n1 n2 n3]);
end;
fclose(fid);

fld=convert2gcmfaces(fld);

1	gforget	1.1	function [fld]=read_bin(fileName,varargin);
2			%usage: fld=read_bin(fileName); reads 3D field (N levels in 3rd dim)
3			%usage: fld=read_bin(fileName,trec); reads 3D field, at trec (50 levels)
4			%usage: fld=read_bin(fileName,trec,k); reads 2D field, at level k, at trec (50 levels)
5			%usage: fld=read_bin(fileName,trec,0); reads 2D field, at trec (1 level)
6
7			global mygrid;
8
9			n1=mygrid.ioSize(1); n2=mygrid.ioSize(2); n3=length(mygrid.RC); nF=mygrid.nFaces;
10			recl2D=n1n24; recl3D=n1n2n3*4; myprec='float32';
11			%recl2D=n1n28; recl3D=n1n2n3*8; myprec='float64';
12
13			if nargin==1; n3=dir(fileName); n3=n3.bytes/recl2D; recl3D=recl2D*n3; end;
14			if nargin>1; tt=varargin{1}; else; tt=1; end;
15			if nargin>2; kk=varargin{2}; else; kk=[]; end;
16
17			fid=fopen(fileName,'r','b');
18			if ~isempty(kk);
19			if kk==0;%meaning the field is only two dimensional
20			status=fseek(fid,(tt-1)*recl2D,'bof');
21			fld=reshape(fread(fid,n1*n2,myprec),[n1 n2]);
22			else;
23			status=fseek(fid,(tt-1)recl3D+(kk-1)recl2D,'bof');
24			fld=reshape(fread(fid,n1*n2,myprec),[n1 n2]);
25			end;
26			else;
27			status=fseek(fid,(tt-1)*recl3D,'bof');
28			fld=reshape(fread(fid,n1n2n3,myprec),[n1 n2 n3]);
29			end;
30			fclose(fid);
31
32			fld=convert2gcmfaces(fld);
33