itgplot

D-ITG2.4源代码

#! /usr/bin/octave -qf


graphicformat="postscript eps color";
extension="eps";

if nargin>0
   inputf=argv{1};
else
   printf("%s\n","Input file name missing, exiting...");
   return;
endif

if (f=fopen(inputf,"r"))<0
   printf("%s%s%s\n","File ",inputf," does not exist, exiting...");
   return;
endif

if size(s=getenv("DITG_PLOT_TITLE"))
   title(s);
endif

# read labels
labels=blanks(0);
fscanf(f,"%s",1);    # we do not need the first string (Time)
while !strcmp(s=fscanf(f,"%s",1),"Aggregate-Flow")
labels=[labels; strcat("-@;",s,";") ];
endwhile
labels=[labels; strcat("-@;",s,";") ];    # add the last string (Aggregate-Flow)

if nargin>1
   mask=intersection(str2num(argv{2}),1:size(labels,1));
else
   mask=1:size(labels,1);
endif

if !size(mask)
   printf("%s\n","No flow selected to be plotted");
   return;
endif

# read data 
data=fscanf(f,"%f",[size(labels,1)+1,Inf])';

temp=split(strrep(inputf,"\\","/"),"/");
outputf=split(deblank(temp(size(temp,1),:)),".");
ylbl=outputf(1,:);

eval(strcat("gset term ",graphicformat));
gset key below box
eval(strcat("gset output '",ylbl,".",extension,"'"));   # set output file

ylbl(1,1)=toupper(ylbl(1,1));
ylabel(ylbl);
xlabel("Time (s)");

plot(data(:,1), data(:,mask+1), labels(mask,:));

gset output 

printf("%s\n","Done!");