flt_example/input/read_flt_traj.m

function [flt,data,header] = read_flt_traj(varargin)
% Reads the float_trajectories files.
%
% flts=read_flt_traj(File_Names,[Worldlength]);
% input Worldlength (= 4 or 8) is optional
% returns a structured array with fields 'time','x','y','k','u','v','t','s','p'
%
% eg.
% >> flts=read_flt_traj('float_trajectories',4);
% >> plot( flts(3).time, flts(3).x/1e3 )
% >> for k=1:126;plot(flts(k).x/1e3,flts(k).y/1e3);hold on;end;hold off

% $Header: /u/gcmpack/MITgcm/verification/flt_example/input/read_flt_traj.m,v 1.5 2009/02/25 01:41:00 dfer Exp $
% $Name:  $

fName = varargin{1};
imax=13;                  % record size
ieee='b';                 % IEEE big-endian format
WORDLENGTH = 8;           % 8 bytes per real*8
if length(varargin)==2
   WORDLENGTH = varargin{2};
end
bytesPerRec=imax*WORDLENGTH;
rtype =['real*',num2str(WORDLENGTH)];

[I]=strfind(fName,'/');
if length(I) == 0,
 bDr='';
else
 fprintf(' found Dir Sep in file name (');
 fprintf(' %i',I);
 bDr=fName(1:I(end));
 fprintf(' ) ; load files from Dir "%s"\n',bDr);
end

fls=dir([fName,'.*data']);

data=zeros(imax,0);
header=zeros(imax,0);

% Read everything
for k=1:size(fls,1)
 fid=fopen([bDr,fls(k).name],'r',ieee);
%fprintf('fid= %i\n',fid);
 nrecs=fls(k).bytes/bytesPerRec;
 ldata=fread(fid,[imax nrecs],rtype);
 fclose(fid);
 header=[header ldata(:,1)];
 data=[data ldata(:,2:end)];
 clear ldata;
end

flt=struct('numsteps',[],'time',[],'x',[],'y',[],'z',[]);

% Sort it all out
for k=1:max(max(data(1,:)));
 j=find( data(1,:)==k );
 [t,jj]=sort( data(2,j) ); j=j(jj);
 flt(k).time=data(2,j);
 flt(k).x=data( 3,j);
 flt(k).y=data( 4,j);
 flt(k).z=data( 5,j);
 flt(k).i=data( 6,j);
 flt(k).j=data( 7,j);
 flt(k).k=data( 8,j);
 flt(k).p=data( 9,j);
 flt(k).u=data(10,j);
 flt(k).v=data(11,j);
 flt(k).t=data(12,j);
 flt(k).s=data(13,j);
end

return
1	dfer	1.6	function [flt,data,header] = read_flt_traj(varargin)
2	adcroft	1.1	% Reads the float_trajectories files.
3			%
4	dfer	1.6	% flts=read_flt_traj(File_Names,[Worldlength]);
5			% input Worldlength (= 4 or 8) is optional
6	adcroft	1.1	% returns a structured array with fields 'time','x','y','k','u','v','t','s','p'
7			%
8			% eg.
9	dfer	1.6	% >> flts=read_flt_traj('float_trajectories',4);
10	adcroft	1.1	% >> plot( flts(3).time, flts(3).x/1e3 )
11			% >> for k=1:126;plot(flts(k).x/1e3,flts(k).y/1e3);hold on;end;hold off
12
13	dfer	1.6	% $Header: /u/gcmpack/MITgcm/verification/flt_example/input/read_flt_traj.m,v 1.5 2009/02/25 01:41:00 dfer Exp $
14	jmc	1.3	% $Name: $
15
16	dfer	1.6	fName = varargin{1};
17	dfer	1.5	imax=13; % record size
18			ieee='b'; % IEEE big-endian format
19	dfer	1.6	WORDLENGTH = 8; % 8 bytes per real*8
20			if length(varargin)==2
21			WORDLENGTH = varargin{2};
22			end
23			bytesPerRec=imax*WORDLENGTH;
24			rtype =['real*',num2str(WORDLENGTH)];
25	adcroft	1.1
26	jmc	1.3	[I]=strfind(fName,'/');
27			if length(I) == 0,
28			bDr='';
29			else
30			fprintf(' found Dir Sep in file name (');
31	dfer	1.5	fprintf(' %i',I);
32	jmc	1.3	bDr=fName(1:I(end));
33			fprintf(' ) ; load files from Dir "%s"\n',bDr);
34			end
35
36			fls=dir([fName,'.*data']);
37	adcroft	1.1
38			data=zeros(imax,0);
39			header=zeros(imax,0);
40
41			% Read everything
42			for k=1:size(fls,1)
43	jmc	1.3	fid=fopen([bDr,fls(k).name],'r',ieee);
44			%fprintf('fid= %i\n',fid);
45	adcroft	1.1	nrecs=fls(k).bytes/bytesPerRec;
46	dfer	1.6	ldata=fread(fid,[imax nrecs],rtype);
47	adcroft	1.1	fclose(fid);
48			header=[header ldata(:,1)];
49			data=[data ldata(:,2:end)];
50	dfer	1.6	clear ldata;
51	adcroft	1.1	end
52
53	jmc	1.4	flt=struct('numsteps',[],'time',[],'x',[],'y',[],'z',[]);
54	adcroft	1.1
55			% Sort it all out
56			for k=1:max(max(data(1,:)));
57			j=find( data(1,:)==k );
58			[t,jj]=sort( data(2,j) ); j=j(jj);
59			flt(k).time=data(2,j);
60	jmc	1.4	flt(k).x=data( 3,j);
61			flt(k).y=data( 4,j);
62			flt(k).z=data( 5,j);
63			flt(k).i=data( 6,j);
64			flt(k).j=data( 7,j);
65			flt(k).k=data( 8,j);
66			flt(k).p=data( 9,j);
67			flt(k).u=data(10,j);
68			flt(k).v=data(11,j);
69			flt(k).t=data(12,j);
70			flt(k).s=data(13,j);
71	adcroft	1.1	end
72	dfer	1.6
73			return