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,filesep); 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	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,filesep); 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