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');
  delete(tfile)
 catch
  delete(tfile)
  error(sprintf('An error occured while scanning for: %s',varargin{k}));
 end
end
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			eval(['!grep ' varargin{k} ' ' file ' \| sed s/.\*=// \| sed s/NAN/1.23456789/ >! ' tfile])
23			vals(:,k)=textread(tfile,'%f');
24			delete(tfile)
25			catch
26			delete(tfile)
27			error(sprintf('An error occured while scanning for: %s',varargin{k}));
28			end
29			end