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');
   if abs(length(tmp)-lngt)>1
    % try to read one line at a time in order to deal with special case
    % of values like, e.g.: "  -2.9248686233802-321"
    fid=fopen(tfile);
    n=0;
    while(~feof(fid))
     n=n+1;
     tmp2=fgetl(fid);
     val=sscanf(tmp2,'%f');
     if length(val)>1
      tmp(n)=0;
     else
      tmp(n)=val;
     end
    end
    fid=fclose(fid);
    tmp=tmp(1:n);
   end
   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	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
24	% vals(:,k)=textread(tfile,'%f');
25	% When output file is from an ongoing integration, one or more of
26	% the diagnostics may be missing at the last available time step.
27	% The code below accomodates this difference in length.
28
29	if k==1
30	vals(:,k)=textread(tfile,'%f');
31	lngt=length(vals(:,k));
32	else
33	tmp=textread(tfile,'%f');
34	if abs(length(tmp)-lngt)>1
35	% try to read one line at a time in order to deal with special case
36	% of values like, e.g.: " -2.9248686233802-321"
37	fid=fopen(tfile);
38	n=0;
39	while(~feof(fid))
40	n=n+1;
41	tmp2=fgetl(fid);
42	val=sscanf(tmp2,'%f');
43	if length(val)>1
44	tmp(n)=0;
45	else
46	tmp(n)=val;
47	end
48	end
49	fid=fclose(fid);
50	tmp=tmp(1:n);
51	end
52	if abs(length(tmp)-lngt)>1
53	error(sprintf('An error occured while scanning for: %s',varargin{k}));
54	else
55	lngt=min(lngt,length(tmp));
56	vals(1:lngt,k)=tmp(1:lngt);
57	end
58
59	end
60	delete(tfile)
61	catch
62	delete(tfile)
63	error(sprintf('An error occured while scanning for: %s',varargin{k}));
64	end
65	end
66	vals=vals(1:lngt,:);