Diagnostics/DiagUtility/grepread.m

function [vals] = grepread(file,varargin)
% VALS = GREPREAD(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=grepread('output.txt','time_secondsf','ke_mean','ke_max');
% >> plot(vals(:,1)/86400,vals(:,2:3));

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] = grepread(file,varargin)
2	% VALS = GREPREAD(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=grepread('output.txt','time_secondsf','ke_mean','ke_max');
10	% >> plot(vals(:,1)/86400,vals(:,2:3));
11
12	if nargin<2
13	error('You must supply a filename and at least one search expression!')
14	end
15
16	tfile=sprintf('/tmp/grepexpr%15.15f',rand);
17	for k=1:nargin-1;
18	try
19	eval(['!grep ' varargin{k} ' ' file ' \| sed s/.\*=// \| sed s/NAN/1.23456789/ >! ' tfile])
20	vals(:,k)=textread(tfile,'%f');
21	delete(tfile)
22	catch
23	delete(tfile)
24	error(sprintf('An error occured while scanning for: %s',varargin{k}));
25	end
26	end