function ep1prr600look(NumOfDays) % EP1PRR600LOOK - Opens and graphs the files \oasis\eqpac\ep1\spec.20m, % spec.20m.1230, spec.20m.dark. % % 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 if NumOfDays < 45 Placement = 0; else Placement = 2; 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'; % Read spec.20m FileLoc = [PathS 'spec.20m']; X = readeqpac(FileLoc, NumOfDays); fig1 = figure; set(fig1,'Units','normalized','Position',[0,0.4,0.5,0.5],'Name','Ed 20m at noon local') subplot(2,1,1) plot(X(:,1),X(:,2:8)) legend('Ed412', 'Ed443', 'Ed490', 'Ed510', 'Ed555', 'Ed670', 'PAR',Placement); ylabel('Irradiance [\muW cm^-^2 nm^-^1]'); xlabel('Date'); title(['Ed[20m] from a PRR-600 at 12:00[local] ' loc]) if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) %fig2 = figure; %set(fig2,'Units','normalized','Position',[0.5,0.4,0.5,0.5],'Name','Lu 20m at noon local') subplot(2,1,2) plot(X(:,1),X(:,9:15)) legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',Placement); ylabel('Radiance [\muW cm^-^2 nm^-^1 sr^-^1]'); xlabel('Date'); title(['Lu[20m] from a PRR-600 at 12:00[local] ' loc]) if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) setfontsize(7) % Read spec.20m.1230 FileLoc = [PathS 'spec.20m.1230']; X = readeqpac(FileLoc, NumOfDays); fig3 = figure; set(fig3,'Units','normalized','Position',[0,0.25,0.5,0.5],'Name','Ed 20m at 1230 local') subplot(2,1,1) plot(X(:,1),X(:,2:8)) legend('Ed412', 'Ed443', 'Ed490', 'Ed510', 'Ed555', 'Ed670', 'PAR',Placement); ylabel('Irradiance [\muW cm^-^2 nm^-^1]'); xlabel('Date'); title(['Ed[20m] from a PRR-600 at 12:30[local] ' loc]) if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) %fig4 = figure; %set(fig4,'Units','normalized','Position',[0.5,0.25,0.5,0.5],'Name','Lu 20m at 1230 local') subplot(2,1,2) plot(X(:,1),X(:,9:15)) title(['Lu[20m] from a PRR-600 at 12:00[local] ' loc]) legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',Placement); ylabel('Radiance [\muW cm^-^2 nm^-^1 sr^-^1]'); xlabel('Date'); title('Lu[20m] from a PRR-600 at 12:30[local] 0 150W') if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) setfontsize(7) % Read spec.20m.dark FileLoc = [PathS 'spec.20m.dark']; X = readeqpac(FileLoc, NumOfDays); fig5 = figure; set(fig5,'Units','normalized','Position',[0,0.1,0.5,0.5],'Name','Ed 20m at 0000 local') subplot(2,1,1) plot(X(:,1),X(:,2:8)) legend('Ed412', 'Ed443', 'Ed490', 'Ed510', 'Ed555', 'Ed670', 'PAR',Placement); ylabel('Irradiance [\muW cm^-^2 nm^-^1]'); xlabel('Date'); title(['Ed[20m] from a PRR-600 at 00:00[local] ' loc]) if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) %fig6 = figure; %set(fig6,'Units','normalized','Position',[0.5,0.1,0.5,0.5],'Name','Lu 20m at 0000 local') subplot(2,1,2) plot(X(:,1),X(:,9:15)) legend('Lu412', 'Lu443', 'Lu490', 'Lu510', 'Lu555', 'Lu670', 'Lu683',Placement); ylabel('Radiance [\muW cm^-^2 nm^-^1 sr^-^1]'); xlabel('Date'); title(['Lu[20m] from a PRR-600 at 00:00[local] ' loc]) if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) setfontsize(7)