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	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