/[MITgcm]/MITgcm/verification/flt_example/input/read_flt_traj.m
ViewVC logotype

Diff of /MITgcm/verification/flt_example/input/read_flt_traj.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by jmc, Wed Dec 3 01:50:15 2008 UTC revision 1.3 by jmc, Sun Jan 4 01:03:59 2009 UTC
# Line 8  function [flt,data,header] = read_flt_tr Line 8  function [flt,data,header] = read_flt_tr
8  % >> plot( flts(3).time, flts(3).x/1e3 )  % >> plot( flts(3).time, flts(3).x/1e3 )
9  % >> for k=1:126;plot(flts(k).x/1e3,flts(k).y/1e3);hold on;end;hold off  % >> for k=1:126;plot(flts(k).x/1e3,flts(k).y/1e3);hold on;end;hold off
10    
11    % $Header$
12    % $Name$
13    
14  imax=10;                % record size  imax=10;                % record size
15  ieee='b';               % IEEE big-endian format  ieee='b';               % IEEE big-endian format
16  bytesPerRec=imax*8;     % 8 bytes per real*8  bytesPerRec=imax*8;     % 8 bytes per real*8
17    
18  fls=dir([fName '.*data']);  [I]=strfind(fName,'/');
19    if length(I) == 0,
20     bDr='';
21    else
22     fprintf(' found Dir Sep in file name (');
23     fprintf(' %i',I);
24     bDr=fName(1:I(end));
25     fprintf(' ) ; load files from Dir "%s"\n',bDr);
26    end
27    
28    fls=dir([fName,'.*data']);
29    
30  data=zeros(imax,0);  data=zeros(imax,0);
31  header=zeros(imax,0);  header=zeros(imax,0);
32    
33  % Read everything  % Read everything
34  for k=1:size(fls,1)  for k=1:size(fls,1)
35   fid=fopen(fls(k).name,'r',ieee);   fid=fopen([bDr,fls(k).name],'r',ieee);
36    %fprintf('fid= %i\n',fid);
37   nrecs=fls(k).bytes/bytesPerRec;   nrecs=fls(k).bytes/bytesPerRec;
38   ldata=fread(fid,[imax nrecs],'real*8');   ldata=fread(fid,[imax nrecs],'real*8');
39   fclose(fid);   fclose(fid);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22