function ep1mcplook(id, NumOfDays)
% EP1MCPLOOK - Opens and graphs the file /oasis/eqpac/ep1/mcp

% Peter Strutton
% 22 Dec 1998; Brian Schlining

% Prompt for input if no input arguments are specified
if nargin < 2 
   NumOfDays = 30;
end

% Set the path to the data file depending on what type of machine this function is run from
FileName = 'mcp';
FileLoc = fullfile(getOasisPath, 'eqpac', id, FileName);

switch lower(id)
case 'ep1'
   loc = '0 155W';
case 'ep2'
   loc = '2S 170W';
otherwise
   loc = [];
end

% Open and read the file
fid = fopen(FileLoc,'r');
OK = 1;
while OK                         % Trash the header
   S = fgetl(fid);
   if ~strcmp(S(1),'#')
      break
   end
end
i = 1;
while ~feof(fid)                 % Read the data
   XS     = fgetl(fid);          % Get a line from the file
   X(i,:) = sscanf(XS,'%f')';    % Read the data from the line
   i      = i + 1;
end
fclose(fid);

[rX cX] = size(X);                  % Get it's size
D       = dy2date(X(:,1)+1900);           % Convert decimal years to serial date numbers
DEnd    = max(D);                   % Get the last date
i       = find(D >= (DEnd - NumOfDays));  % Find all dates within specified period

% Create the graph
figH = figure;
set(figH,'Name','MCP; 10m & 30m')

j       = find((X(i,2) == 10));     % Skip rows of data with NaN's
k       = find((X(i,2) == 30));     % Skip rows of data with NaN's
aH = plotyy(D(i(j)),X(i(j),3),D(i(k)),X(i(k),3));
%set(aH(1),'XLim',[min(X(:,1)) now],'XTick',[]);
set(aH(1),'XTick',[]);
set(figH,'CurrentAxes',aH(1))
ylabel('10m MCP output [volts]');
set(figH,'CurrentAxes',aH(2))
set(aH(1),'XTick',[])
ylabel('30m MCP output [volts]');
xlabel('Date');
title(['MCP output in volts at 10m & 30m for ' loc]);

%set(aH(2),'XLim',[min(D(i,1)) now])
%hold on
if NumOfDays < 40
   dtick = 6;
else
   dtick = 3
end

datetick('x',dtick)

XLim = get(aH(2),'XLim');
set(aH(1),'XLim',XLim)
setfontsize(7)






