function ep1ocr100look(NumOfDays) % EP1OCR100LOOK - The OCR-100's are the Satlantic radiometers that rest in the % buoy bridle at about 1.5 m depth and record upwelling radiance. This script reads % the files \oasis\eqpac\ep1\spec.0m.dark, sat1.dark, spec.noon, spec.0m.1230 and % splots the radiance data (7 wavelengths in each figure). % % Requires: READEQPAC % Peter Strutton % 29 Dec 1998; Brian Schlining % Prompt for input if no input arguments are specified if ~nargin NumOfDays = str2num(input('How many days of data would you like to view? ','s')); end % Set the path to the data file depending on what type of machine this function is run from switch computer case 'PCWIN' % Assumes the pwd is /oasis/epqac/brian/ %PathS = ['..\ep1']; PathS = ['//tsunami.shore.mbari.org/oasis/eqpac/ep1/']; otherwise % Assigns an absolute path in UNIX PathS = ['/hosts/tsunami/oasis/eqpac/ep1']; end loc = '0 155W'; %================================== % Open and read the file spec.noon %================================== FileLoc = [PathS filesep 'spec.noon']; X = readeqpac(FileLoc, NumOfDays); 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]) plot(X(:,1),X(:,9:15)) legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',0); ylabel('Radiance [\muW cm^-^2 nm^-^1 str^-1]'); %xlabel('Date') 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 %===================================== FileLoc = [PathS filesep 'spec.0m.dark']; X = readeqpac(FileLoc, NumOfDays); subplot('position', [0.13 0.11 0.775 0.25]) plot(X(:,1),X(:,9:15)) %legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',0); 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 %===================================== FileLoc = [PathS filesep 'spec.0m.1230']; X = readeqpac(FileLoc, NumOfDays); 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(X(:,1),X(:,9:15)) legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',0); ylabel('Radiance [v]'); %xlabel('Date') 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 %===================================== FileLoc = [PathS filesep 'sat1.dark']; X = readeqpac(FileLoc, NumOfDays); subplot('position', [0.13 0.11 0.775 0.25]) plot(X(:,1),X(:,2:8)) 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)