MITgcm_contrib/mitgcm_tools/mitgcmhistory.m

function [vals] = mitgcmhistory(file,varargin)
%vals = mitgcmhistory(FILE,EXPR1,...);
%
%Extracts the expressions "expr1","expr2",... from the file "file".
%This assumes output in the standard form defined by the MITgcm
%monitor package and is not a replacement for TEXTREAD.
% 
%e.g.
%>> vals=mitgcmhistory('output.txt','time_secondsf','ke_mean','ke_max');
%>> plot(vals(:,1)/86400,vals(:,2:3));
%
% Written by adcroft@mit.edu, 2001
%$Header:

if nargin<2
 error('You must supply a filename and at least one search expression!')
end

tfile=sprintf('/tmp/grepexpr%15.15f',rand);
for k=1:nargin-1;
 try
  eval(['!grep ' varargin{k} ' ' file ' | sed s/.\*=// | sed s/NAN/1.23456789/ > ' tfile])

  % vals(:,k)=textread(tfile,'%f');

  % When output file is from an ongoing integration, one or more of
  % the diagnostics may be missing at the last available time step.
  % The code below accomodates this difference in length.

  if k==1
   vals(:,k)=textread(tfile,'%f');
   lngt=length(vals(:,k));
  else
   tmp=textread(tfile,'%f');
   % allow 
   if abs(length(tmp)-lngt)>1
    error(sprintf('An error occured while scanning for: %s',varargin{k}));
   else
    lngt=min(lngt,length(tmp));
    vals(1:lngt,k)=tmp(1:lngt);
   end

  end
  delete(tfile)
 catch
  delete(tfile)
  error(sprintf('An error occured while scanning for: %s',varargin{k}));
 end
end
vals=vals(1:lngt,:);
1	adcroft	1.1	function [vals] = mitgcmhistory(file,varargin)
2			%vals = mitgcmhistory(FILE,EXPR1,...);
3			%
4			%Extracts the expressions "expr1","expr2",... from the file "file".
5			%This assumes output in the standard form defined by the MITgcm
6			%monitor package and is not a replacement for TEXTREAD.
7			%
8			%e.g.
9			%>> vals=mitgcmhistory('output.txt','time_secondsf','ke_mean','ke_max');
10			%>> plot(vals(:,1)/86400,vals(:,2:3));
11			%
12			% Written by adcroft@mit.edu, 2001
13			%$Header:
14
15			if nargin<2
16			error('You must supply a filename and at least one search expression!')
17			end
18
19			tfile=sprintf('/tmp/grepexpr%15.15f',rand);
20			for k=1:nargin-1;
21			try
22	dimitri	1.2	eval(['!grep ' varargin{k} ' ' file ' \| sed s/.\*=// \| sed s/NAN/1.23456789/ > ' tfile])
23	dimitri	1.3
24			% vals(:,k)=textread(tfile,'%f');
25
26			% When output file is from an ongoing integration, one or more of
27			% the diagnostics may be missing at the last available time step.
28			% The code below accomodates this difference in length.
29
30			if k==1
31			vals(:,k)=textread(tfile,'%f');
32			lngt=length(vals(:,k));
33			else
34			tmp=textread(tfile,'%f');
35			% allow
36			if abs(length(tmp)-lngt)>1
37			error(sprintf('An error occured while scanning for: %s',varargin{k}));
38			else
39			lngt=min(lngt,length(tmp));
40			vals(1:lngt,k)=tmp(1:lngt);
41			end
42
43			end
44	adcroft	1.1	delete(tfile)
45			catch
46			delete(tfile)
47			error(sprintf('An error occured while scanning for: %s',varargin{k}));
48			end
49			end
50	dimitri	1.3	vals=vals(1:lngt,:);