function ep2prr620look(NumOfDays) % EP2PRR600LOOK - Opens and graphs the files \oasis\eqpac\ep1\spec.noon, % spec.0m.1230, spec.0m.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 % 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 = ['//tsunami.shore.mbari.org/oasis/eqpac/ep2/']; otherwise % Assigns an absolute path in UNIX PathS = ['/hosts/tsunami/oasis/eqpac/ep2/']; end if NumOfDays < 45 Placement = 2; else Placement = 2; end % Read spec.noon FileLoc = [PathS 'spec.noon']; X = readeqpac(FileLoc, NumOfDays); fig1 = figure; set(fig1,'Units','normalized','Position',[0,0.4,0.5,0.5],'Name','PRR-620 at noon local') 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[3m+] from a PRR-620 at 12:00[local] 2S 170W') if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) setfontsize(7) % Read spec.0m.1230 FileLoc = [PathS 'spec.0m.1230']; X = readeqpac(FileLoc, NumOfDays); fig2 = figure; set(fig2,'Units','normalized','Position',[0.25,0.25,0.5,0.5],'Name','PRR-620 at 1230 local') 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[3m+] from a PRR-620 at 12:30[local] 2S 170W') if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) setfontsize(7) % Read spec.0m.dark FileLoc = [PathS 'spec.0m.dark']; X = readeqpac(FileLoc, NumOfDays); fig3 = figure; set(fig3,'Units','normalized','Position',[0.4,0.1,0.5,0.5],'Name','PRR-620 at 0000 local') 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[3m+] from a PRR-620 at 00:00[local] 2S 170W') if NumOfDays < 50 dtick = 6; else dtick = 3; end datetick('x',dtick) setfontsize(7)