function ep_lookOcr100(id, NumOfDays)
% EP_lookOcr100 - The OCR-100's are the Satlantic radiometers that rest in the
% buoy bridle at about 1.5 m depth and record upwelling radiance. 
%
% Requires: READEQPAC

% Peter Strutton
% 29 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
switch lower(id)
case 'ep1'
   loc = '0 155W';
   FileName = 'ep1_specprr_full.nc';
   lon = -155;
case 'ep2'
   loc = '2S 170W';
   FileName = 'ep2_specprr_full.nc';
   lon = -170;
otherwise
   loc = [];
   return
end

FileLoc = fullfile(getOasisPath, 'eqpac', 'netcdf', FileName);

ep     = netcdf(FileLoc);
t      = nctime(ep);
lambda = [412 443 490 510 555];
satA   = [ep{'Lu412_1_1'}(:) ep{'Lu443_1_2'}(:) ep{'Lu490_1_3'}(:) ep{'Lu510_1_4'}(:) ep{'Lu555_1_5'}(:)];
satB   = [ep{'Lu412_1_9'}(:) ep{'Lu443_1_10'}(:) ep{'Lu490_1_11'}(:) ep{'Lu510_1_12'}(:) ep{'Lu555_1_13'}(:)];
satA(find(satA == -999)) = NaN;
satB(find(satB == -999)) = NaN;

ncquit(ep);

% Trim to NumOfDays
DEnd    = max(t);							% Get the last date
i       = find(t >= (DEnd - NumOfDays));	% Find all dates within specified period
t       = t(i,:);
satA    = satA(i,:);
satB    = satB(i,:);

[Year Month Day Hour Minute Second] = datevec(t + gmtoffset(lon));
%N = find(Hour >= 11 & Hour <= 13);
N = find(Hour);
M = find(Hour >= 23 | Hour <= 1);
M = N;

%==================================
% Open and read the file spec.noon
%==================================
fig1    = figure;
set(fig1,'Units','normalized','Position',[0,0.4,0.5,0.5],'Name','Satlantic A ')
subplot('position', [0.13 0.47 0.7750 0.42])
%good = find(~isnan(satA(N)));
plot(t(N),satA(N,:))
legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',0);
ylabel('Radiance [\muW cm^-^2 nm^-^1 str^-1]');
title(['Lu[1.5m] from a OCR-100 (A) at 12:00[local] ' loc])
if NumOfDays < 50
   dtick = 6;
else
   dtick = 3;
end

datetick('x',dtick)

%=====================================
% Open and read the file spec.0m.dark
%=====================================
subplot('position', [0.13 0.11 0.775 0.25])
x = readeqpac(fullfile(getOasisPath, 'eqpac', id, 'spec.0m.dark'));
plot(t(M),satA(M,:))
ylabel('Radiance [\muW cm^-^2 nm^-^1 sr^-^1]');
xlabel('Date')
title(['Lu[1.5m] from a OCR-100 (A) at 00:00[local] ' loc])
if NumOfDays < 50
   dtick = 6;
else
   dtick = 3;
end
datetick('x',dtick)
setfontsize(7)

%=====================================
% Open and read the file spec.0m.1230
%=====================================
fig4 = figure;
set(fig4,'Units','normalized','Position',[0.5,0.1,0.5,0.5],'Name','Satlantic B')
subplot('position', [0.13 0.47 0.7750 0.42])
plot(t(N),satB(N,:))
legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',0);
ylabel('Radiance [v]');
title(['Lu[1.5m] from a OCR-100 (B) at 12:30[local] ' loc])
if NumOfDays < 50
   dtick = 6;
else
   dtick = 3;
end

datetick('x',dtick)

%=====================================
% Open and read the file spec.0m.dark
%=====================================
subplot('position', [0.13 0.11 0.775 0.25])
plot(t(M), satB(M,:))
ylabel('Radiance [\muW cm^-^2 nm^-^1 sr^-^1]');
xlabel('Date')
title(['Lu[1.5m] from a OCR-100 (B) at 00:00[local] ' loc])
if NumOfDays < 50
   dtick = 6;
else
   dtick = 3;
end
datetick('x',dtick)
setfontsize(7)
