Latest News: CDIP buoy station 111 Anacapa Passage, CA was decommissioned May 14, 2019. If you relied on this data, we would like to hear from you.

User Preferences

Timezone
Units
Lat/Lon format
Defaults

Geo and time bounds

073 - SCRIPPS PIER, LA JOLLA CA   Notice:04/16/2019 Wave sensors cleaned 1815 to 1900 UTC Data set:

ENERGY/DIRECTION MONTHLY VECTOR PLOT

Energy and direction spectrum data - derived from a fast fourier transform of approximately 30 minutes of wave data - wave data) is plotted with frequency on the x axis and time on the y axis. The length of a vector corresponds to the energy density (with scale to the right of the plot) and the wavefront direction is in the direction of the arrow. True North is at the top of the page. Every fourth data run is plotted.


Matlab code for generating the plot (no longer used)

%-------------------------------------------------------------------------------
% Name:
%   .dbase/web_programs/plotpro/matlab_scripts/feather.m
% Description:
%   Creates the monthly vector plot (feathers).
%   A few lines are pre-pended to this script.
% Called by:
%   .f90/plot_utils.make_plot
%-------------------------------------------------------------------------------


path(path,'/project/dbase/web_programs/plotpro/matlab_scripts');

wkg = wkg_path;

% give short names for loaded variables.
% NOTE: x = frequency data, y = energy density data, z = direction

eval(['freq=load(''',wkg,'freq.',pid,'01'');']);
eval(['ener=load(''',wkg,'ener.',pid,'01'');']);
eval(['dire=load(''',wkg,'dire.',pid,'01'');']);

x = freq(1,:);
y = ener;
z = dire;

% check to see if we have data
if ~ (any(x>0) & any(y>0) & any(z>0)) 
  unix(['\rm ',wkg,'plot_',pid,'*']); % for pp when run mountain with archive.
  clear all;
  quit;
end

dt = read_date_file([pid,'01'],wkg);
tim = dt.day + (dt.hour + dt.min/60)/24; % in days

yearchar = num2str(dt.year(1));
monthnum = dt.month(1);

clear freq ener date

% set energy to zero where direction is -1
I = find(z<0);
y(I) = 0;
clear I;

% set figure location on paper 
ppleft = .05;
ppbottom = .05;
ppwidth = .9;
ppheight = .93;

% set figure window
figure('units','normal','paperunits','normal', 'position',[.1,.1,.5,.8], ...
                      'paperposition',[ppleft,ppbottom,ppwidth,ppheight],'visible','off');

% to get arrow directions right get paper dimensions
xpage = 8.5;
ypage = 11;

% setup axes for text outside the plot region
tah = axes('position',[0 0 1 1],'visible','off');

% setup axes for data plot.
axleft = 0.1;
axbottom = 0.1;
axwidth = 0.7;
axheight = 0.76;
dah = axes('position',[axleft, axbottom, axwidth, axheight],'visible','off');

% setup plot scale parameters.
uxlim = 0.25;
lxlim = 0.04;
uylim = 33;
lylim = 1;

sy= size(y);
sz= size(z);

% get date information and index I for first of every
% day and every six hours.
pday = 0;
ptime = 0;
j=1;
k=1;
for i=1:sz(1)
  if ( (dt.day(i) ~= pday) | (tim(i) - ptime > 5/24))
    pday = dt.day(i);
    ptime = tim(i);
    I(k) = i;
    k = k + 1;
  end
end


% relate paper dimensions to axes scales
J = find( x<=uxlim & x >= lxlim);
Imin = min(J);
Imax = max(J);
xscale =  ( Imax - Imin )/(xpage*ppwidth*axwidth);
%xscale =  ( uxlim - lxlim )/(xpage*ppwidth*axwidth);
yscale =  ( uylim - lylim )/(ypage*ppheight*axheight);

% remove records not used
y = y(I,J);
ymax = max(max(y));
z = z(I,J);
tim = tim(I);
clear I;

%-------------------------------------------------------------------------------
%-- PLOT VECTORS ---------------------------------------------------------------
%-------------------------------------------------------------------------------

stretch = 7;
%ymax=30.;  % Use for fixed scale
%stretch = 0;
[u,v]=pol2cart(-1*(z+90)*(2*pi/360),y./ymax);
quiver(J,tim,xscale.*u,yscale.*v,stretch);


% LABEL X AXIS

numticks = 5;
diff = floor((Imax - Imin)/(numticks-1));
xt = Imin + (0:(numticks-1))*diff;
fxtl = num2str(x(xt)',2);
pxtl = num2str(1./x(xt)',2);
set(gca,'xlim',[Imin,Imax],'xtick',xt,'xticklabel',fxtl);
text(xt,-.5*ones(1,numticks),pxtl,'horizontalalignment','center');
text(Imax+.5,0.5,'Frequency (Hz)','horizontalalignment','left',...
     'fontsize',10);
text(Imax+.5,-.37,'Period (s)','horizontalalignment','left',...
     'fontsize',10);
 

% LABEL Y AXIS

ylabel('Day of Month (UTC)','fontsize',12);
set(gca,'ylim',[lylim,uylim],'ytick',[lylim:uylim-2],...
   'yticklabel',num2str([lylim:uylim-2]'));


% GET STATION INFO

info = read_info_file([pid,'01'],wkg);


% WRITE STATION INFO ON PLOT

set(gcf,'currentaxes',tah);

text(.45,.955,[info.stnid,' ',info.stname],'horizontalalignment',...
  'center','fontsize',14);

stream_name = ['(',char(info.stream_name),')'];
if ( length(stream_name) > 5 )
  text(.45, .932, stream_name, 'horizontalalignment',...
  'center','fontsize',12);
end

text(.45,.903,'ENERGY/DIRECTION SPECTRUM','horizontalalignment',...
  'center','fontsize',14);

mname = get_monthname(monthnum,'long');
text(.45,.883,[char(mname),' ',yearchar],'horizontalalignment',...
  'center','fontsize',14);



% PUT SCALE ON FIGURE

plotmax = 0.197; % magic number: length of ymax (inches) on plot per unit of stretch

ymaxlength = plotmax*stretch; % length in inches of ymax on plot
%ymaxlength = plotmax; % use for fixed scale

scalelength = 2; % desired approximate length of scale in inches

ytop = ymax*scalelength/ymaxlength;

if ytop < .1
  rndval = .01;
elseif ytop < 1
  rndval = .1;
elseif ytop < 10
  rndval = 1;
else  
  rndval = 10;
end

scaleytop = rndval*round(ytop/rndval); % top value on scale

% y scale height in normalized units

yscaleheight = (scalelength/ypage)*(scaleytop/ytop);

axes('position',[.870,.60,.05,yscaleheight]);

hold on; 
 

%draw the scale line (4 units long)
ytick = [0,scaleytop/4,scaleytop/2,3*scaleytop/4,scaleytop];
yticklabs = num2str(ytick',4);
feather(0,scaleytop);
hold on;
set(gca,'ylim',[0,scaleytop],'ytick',ytick,'yticklabel',yticklabs,'xticklabel',[]);
xlim=get(gca,'xlim');
ylim=get(gca,'ylim');
xtext = (xlim(2)-xlim(1))/7 + xlim(2);
text(xtext,ytop/2,'Energy Density (m^2/Hz)','rotation',90, ...
  'horizontalalignment','center','fontsize',12);
text((xlim(2)-xlim(1))/2 + xlim(1), ylim(2) + ylim(2)/5,...
  'TRUE','horizontalalignment','center','fontsize',12);
text((xlim(2)-xlim(1))/2 + xlim(1), ylim(2) + ylim(2)/10,...
  'NORTH','horizontalalignment','center','fontsize',12);

blue_box(0.01)

% Print the plot

%saveas(gcf,'/home/corey/feather.fig')

set(gcf, 'visible','off');
png_name = [wkg,'plot_',pid,'.png'];
print('-dpng','-r75',png_name);


quit