function com_plot(nfiles)
%A function to overplot data  
%to compare vAriables from EPIC files.
% Input argument is how many files.
cnames = ['rgbymckrgbymckrgbymck'];
for i=1:nfiles
   nc =netcdf
	v = uigetvar(nc)
	eptime1 = nc{'time'}(:);
	eptime2 = nc{'time2'}(:);
   eval(['time' num2str(i) ' = ep_datenum([eptime1 eptime2]);']);
   eval(['first(i) = time' num2str(i) '(1);']);
   eval(['last(i) = time' num2str(i) '(length(eptime1));']);
   eval(['data' num2str(i) ' = v(:);']);
   eval(['vmin(i) = min(data' num2str(i) ');']);
   eval(['vmax(i) = max(data' num2str(i) ');']);
   vlabel{i} = justname(name(nc));
   which_var = name(v);
	ncclose   
end

start = min(first);
stop = max(last);
lowest = min(vmin);
highest = max(vmax);
startday = floor(start);
axis ([start stop lowest highest])
hold off
for i = 1:nfiles
   eval(['plot (time' num2str(i) ', data' num2str(i) ', ''color'', cnames(i))' ])
	eval(['text (start,data' num2str(i) '(1),vlabel{i})'])
	hold on
end
xlabel ('time')
ylabel (which_var)
%dateaxis('x',0,startday)
dateaxis('x',0,0)
ncclose