matlab_class/ecco_v4/v4_read_data.m

function [fld]=v4_read_data(fileName,irec);
%usage: fld=v4_read_data(fileName,irec);  'fast' read of 2D fields no irec in all [fileName '*.data'] 

gcmfaces_global;

dir0=strfind(fileName,'/'); if isempty(dir0); dir0='./'; else; dir0=fileName(1:dir0(end)); end;
fileList=dir([fileName '.data']);

nn=length(fileList);
fld=zeros(90,1170,nn);
for ii=1:nn; 
fid_cur=fopen([dir0 fileList(ii).name],'r','b');
recl=90*1170*4; position0=recl*(irec-1);
status=fseek(fid_cur,position0,'bof');
fld(:,:,ii)=fread(fid_cur,[90 1170],'float32');
fclose(fid_cur);
end;

fld=convert2gcmfaces(fld);

1	gforget	1.1	function [fld]=v4_read_data(fileName,irec);
2			%usage: fld=v4_read_data(fileName,irec); 'fast' read of 2D fields no irec in all [fileName '*.data']
3
4			gcmfaces_global;
5
6			dir0=strfind(fileName,'/'); if isempty(dir0); dir0='./'; else; dir0=fileName(1:dir0(end)); end;
7			fileList=dir([fileName '.data']);
8
9			nn=length(fileList);
10			fld=zeros(90,1170,nn);
11			for ii=1:nn;
12			fid_cur=fopen([dir0 fileList(ii).name],'r','b');
13			recl=9011704; position0=recl*(irec-1);
14			status=fseek(fid_cur,position0,'bof');
15			fld(:,:,ii)=fread(fid_cur,[90 1170],'float32');
16			fclose(fid_cur);
17			end;
18
19			fld=convert2gcmfaces(fld);
20