From iusia79 at hotmail.com  Wed Jan 12 10:39:38 2005
From: iusia79 at hotmail.com (nadia lobue)
Date: Wed Jan 12 11:01:27 2005
Subject: [Sea-mat] matlab stick plot 
Message-ID: <BAY102-F3557B2C9ECCA8032F231B5A6890@phx.gbl>

An HTML attachment was scrubbed...
URL: http://mailman.whoi.edu/pipermail/sea-mat/attachments/20050112/a92afea2/attachment.htm
From jhyatt at whoi.edu  Wed Jan 12 16:29:47 2005
From: jhyatt at whoi.edu (Jason Hyatt)
Date: Wed Jan 12 16:29:59 2005
Subject: [Sea-mat] matlab stick plot
In-Reply-To: <BAY102-F3557B2C9ECCA8032F231B5A6890@phx.gbl>
References: <BAY102-F3557B2C9ECCA8032F231B5A6890@phx.gbl>
Message-ID: <41E596CB.20003@whoi.edu>

timeplt, just watch to make sure you are using the correct julian day.
http://woodshole.er.usgs.gov/operations/sea-mat/timeplt-html/

nadia lobue wrote:

>
>
> Hello everybody,
>
> I'm a new user to matlab. I have many oceanographic data in particular 
> I have currentometer data and I trying to make a simple stick plot 
> with matlab program.Someone knows a specific function to make stick 
> plot or a  method in order to make this type of diagrams?
>
> Thank's a lot
>
> Nadia
>
>  
>
>
>
> ------------------------------------------------------------------------
> Scarica gratuitamente MSN Toolbar 
> <http://g.msn.com/8HMAITIT/2728??PS=47575> !
>
>------------------------------------------------------------------------
>
>_______________________________________________
>Sea-mat maillist  -  Sea-mat@whoi.edu
>http://mailman.whoi.edu/mailman/listinfo/sea-mat
>  
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.whoi.edu/pipermail/sea-mat/attachments/20050112/fc2e4118/attachment.htm
From iusia79 at hotmail.com  Thu Jan 13 09:29:22 2005
From: iusia79 at hotmail.com (nadia lobue)
Date: Thu Jan 13 09:35:21 2005
Subject: [Sea-mat] matlab stick plot
In-Reply-To: <41E596CB.20003@whoi.edu>
Message-ID: <BAY102-F21455CED859C1852E95FACA68A0@phx.gbl>

An HTML attachment was scrubbed...
URL: http://mailman.whoi.edu/pipermail/sea-mat/attachments/20050113/6e7b16d7/attachment.htm
From fmachin at becarios.ulpgc.es  Thu Jan 13 11:17:58 2005
From: fmachin at becarios.ulpgc.es (Francisco =?iso-8859-1?Q?Jos=E9=20Mach=EDn=20Jim=E9nez?=)
Date: Thu Jan 13 11:20:03 2005
Subject: [Sea-mat] matlab stick plot
References: <BAY102-F21455CED859C1852E95FACA68A0@phx.gbl>
Message-ID: <41E69F36.F6CF6957@becarios.ulpgc.es>

    Hi, Nadia:
    Forget the example given in the Matlab help. theta, r and so on are
variables necessary to built a vector time serie, whose components are
finally u and v.

    As you described your case, you have the module vector (u), its
angle with respect to the east and the time when data were recorded. You
need again the components, say u_x and u_y, which are computed as:

u_x=u*cos(angle in radians)
u_y=u*sin(angle in radians)

And finally, you "feather" them:

feather(u_x,u_y)

    Let me know of your progress.

    Regards,

    Francis.

**********************************
Francisco Jos? Mach?n Jim?nez
Dpto. F?sica-U.L.P.G.C.
Campus Universitario de Tafira
35017 Las Palmas de Gran Canaria
Espa?a
Tfn:34-928-454521
FAX:34-928-452922
**********************************


From zhaocf at mail.ouc.edu.cn  Wed Feb  2 21:14:29 2005
From: zhaocf at mail.ouc.edu.cn (zhaocf-ouc)
Date: Thu Feb  3 07:28:33 2005
Subject: [Sea-mat] SEA-MAT toolbox
Message-ID: <307396818.23541@ouc.edu.cn>

Dear Sir,

I would like to use SEA-MAT toolbox in ocean science researches.
Could you please add my email addresss in your email list? Where 
can I download the toolbox?

Best regards,

Chaofang ZHAO
Ocean Remote Sensing Institute 
Ocean University of China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.whoi.edu/pipermail/sea-mat/attachments/20050203/0de116ca/attachment.htm
From klinck at stockmann.ccpo.odu.edu  Thu May 19 09:17:19 2005
From: klinck at stockmann.ccpo.odu.edu (John Klinck)
Date: Thu May 19 10:21:40 2005
Subject: [Sea-mat] update for matlab function
Message-ID: <200505191317.j4JDHJkf015579@stockmann.ccpo.odu.edu>


Attached is  the matlab function for calculating dynamic normal modes 
(dynmodes.m) which I submitted a few years ago. The header now contains a more 
extensive explanation of what the routine does. I get a few requests per year to 
explain the routine, which was poorly documented in its original form. Please 
replace the klinck.tar.Z and klinck.zip with this file (unless you use tar or 
zip wrappers to avoid differences in end of file for different computers).

Let me know if you have any questions. 

John

John M. Klinck
Center for Coastal Physical Oceanography
Department of Ocean, Earth and Atmospheric Sciences
Crittenton Hall
Old Dominion University
Norfolk, Virginia 23529
phone: 757 683-6005, fax: 757 683-5550 email: klinck@ccpo.odu.edu
-------------- next part --------------
function [wmodes,pmodes,ce]=dynmodes(Nsq,p,nmodes)
% DYNMODES calculates ocean dynamic vertical modes
%  taking a column vector of Brunt-Vaisala values (Nsq) at
%  different pressures (p) and calculating some number of 
%  dynamic modes (nmodes). 
%  Note: The input pressures need not be uniformly spaced, 
%    and the deepest pressure is assumed to be the bottom.
%
%  The program is based on the developement in section 6.13 
%  (p. 167ff) of Atmosphere-Ocean Dynamics by Adrian Gill, 
%  Academic Press, 1982.
%
%  The variables describing flow in the ocean or atmosphere are separated
%  into vertical (depending only on height or depth) and horizontal and
%  time. For example, the vertical velocity (w) is represented as a sum of
%  different vertical structures as 
%  $$ 
%  w = \sum_{n=0}^\infty \hat h_n(z) \hat w_n (x,y,t).  
%  $$
%
%  There is a related set of vertical functions ($\hat p_n(z)$) which
%  represent the vertical structure of horizontal velocity and pressure
%  anomaly. In the program below, wmodes is a matrix containing the $\hat
%  h_n(z)$ functions and pmodes are the $\hat p_n(z)$ functions.  wmodes(:,n)
%  or pmodes(:,n) is the nth vertical structure at the depths (p) provided as
%  input.
%
%  Each set of functions is mutually orthogonal with appropriate weighting. 
%
%  The governing equation for the vertical structures (6.11.18 in Gill) is 
%  $$ 
%    d^2/dz^2 \hat h_n + (N^2(z) /c_e) \hat h_n = 0, n = 1, 2, ... ,
%  $$
%  where $N^2(z)$ is the Brunt-Vaiasala freqency squared and $c_e$ is the
%  eigenvalue, which is the phase speed of the internal gravity wave with the
%  indicated vertical structure. The boundary conditions are that $\hat h$
%  vanishes at the surface and bottom (the rigid lid assumption). The
%  solution for n=0 is a special case described in Gill which is very weakly
%  affected by the density structure.
%
%  The vertical velocity structure (wmodes) is calculated at the depths (p)
%  for which $N^2$ is provided on input. The differential equation above is
%  converted to a set of algebraic equations using second order central
%  finite differences. The linear algebraic equation is 
%  $A * wmodes + \lambda B * wmodes = 0$, where B contains $N^2$ along its
%  diagonal. The matlab routine eig(A,B) calculate the eigenvalues $\lambda$
%  which are the reciprocals of the wave speed $c_e$.
%
%  USAGE: [wmodes,pmodes,ce]=dynmodes(Nsq,p,nmodes);
%                               or
%                            dynmodes;  % to show demo 
%
%   Inputs:   Nsq = column vector of Brunt-Vaisala buoyancy frequency (s^-2)
%               p = column vector of pressure (decibars)
%          nmodes = number of vertical modes to calculate 
%  
%   Outputs:   wmodes = vertical velocity structure
%              pmodes = horizontal velocity structure
%                  ce = modal speed (m/s)
%  developed by John Klinck. July, 1999
%  send comments and corrections to klinck@ccpo.odu.edu

if nargin<1
   help(mfilename);
   nplot=3;
%      Calling the function with no input arguements runs
%      one of these two test problems. 
%      This also serves as an example driver to calculate 
%      vertical modes.
%
%      problem 1
%    solution is h = ho sin(z /ce) where ce = 1 / n pi
%     ce = 0.3183 / n = [ 0.3183 0.1591 0.1061]
%  p=0:.05:1;
%  z=-p;
%  n=length(p);
%  Nsq(1:n)=1;
%  nmodes=3;
%
%      problem 2
%    solution is h = ho sin(No z /ce) where ce = No H / n pi
%    for No=1.e-3 and H = 400, the test values are 
%     ce = 0.127324 / n = [ 0.127324, 0.063662, 0.042441]
%
  p=0:10:400;
  z=-p;
  n=length(p);
  Nsq(1:n)=1.e-6;
  nmodes=3;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
%
  [wmodes,pmodes,ce]=dynmodes(Nsq,p,nmodes);

  figure(1)
  plot(Nsq,z);
  title('Buoyancy Frequency Squared (s^{-2})')

  figure(2)
  plot(ce(1:nplot),'r:o');
  title(' Modal Speed (m/s)')

  figure(3)
  plot(wmodes(:,1:nplot),z);
  title('Vertical Velocity Structure')

  figure(4)
  plot(pmodes(:,1:nplot),z);
  title('Horizontal Velocity Structure')

  figure(gcf)
  return
end
  
rho0=1028;
%    convert to column vector if necessary
[m,n] = size(p);
if n == 1
   p=p';
end
[m,n] = size(Nsq);
if n == 1
   Nsq=Nsq';
   n=m;
end

%                 check for surface value
if p(1) > 0
%             add surface pressure with top Nsq value
    z(1)=0;
    z(2:n+1)=-p(1:n);
    N2(1)=Nsq(1);
    N2(2:n+1)=Nsq(1:n);
    nz=n+1;
else
    z=-p;
    N2=Nsq;
    nz=n;
end

%          calculate depths and spacing
%        spacing
dz(1:nz-1)=z(1:nz-1)-z(2:nz);
%        midpoint depth
zm(1:nz-1)=z(1:nz-1)-.5*dz(1:nz-1)'';
%        midpoint spacing
dzm=zeros(1,nz);
dzm(2:nz-1)=zm(1:nz-2)-zm(2:nz-1);
dzm(1)=dzm(2);
dzm(nz)=dzm(nz-1);

%        get dynamic modes
A = zeros(nz,nz);
B = zeros(nz,nz);
%             create matrices   
for i=2:nz-1
  A(i,i) = 1/(dz(i-1)*dzm(i))  + 1/(dz(i)*dzm(i));
  A(i,i-1) = -1/(dz(i-1)*dzm(i));
  A(i,i+1) = -1/(dz(i)*dzm(i));
end
for i=1:nz
  B(i,i)=N2(i);
end
%             set boundary conditions
A(1,1)=-1.;
A(nz,1)=-1.;

[wmodes,e] = eig(A,B);

%          extract eigenvalues
e=diag(e);
%
ind=find(imag(e)==0);
e=e(ind);
wmodes=wmodes(:,ind);
%
ind=find(e>=1.e-10);
e=e(ind);
wmodes=wmodes(:,ind);
%
[e,ind]=sort(e);
wmodes=wmodes(:,ind);

nm=length(e);
ce=1./sqrt(e);
%                   create pressure structure
pmodes=zeros(size(wmodes));

for i=1:nm
%           calculate first deriv of vertical modes
  pr=diff(wmodes(:,i));   
  pr(1:nz-1)= pr(1:nz-1)./dz(1:nz-1)';
  pr=pr*rho0*ce(i)*ce(i);
%       linearly interpolate back to original depths
  pmodes(2:nz-1,i)=.5*(pr(2:nz-1)+pr(1:nz-2));
  pmodes(1,i)=pr(1);
  pmodes(nz,i)=pr(nz-1);
end
From AAllen at rdc.uscg.mil  Tue Jun 21 11:32:49 2005
From: AAllen at rdc.uscg.mil (Allen, Art)
Date: Tue Jun 21 11:34:00 2005
Subject: [Sea-mat] Wave height calculator
Message-ID: <0B81C49C79E8D611AA330002B3B79A0D2947DF@rdcex.rdc.uscg.mil>

I am looking for a simple  wave height calculator that would generate
significant wave height from wind speed,  duration and fetch.


Art Allen
Arthur A. Allen
U.S. Coast Guard
Office of Search and Rescue (G-OPR-1)
located at:
U.S. Coast Guard International Ice Patrol
1082 Shennecossett Road
Groton, CT, 06340-6096
(860) 441-2747
(860) 441-2773 FAX
aallen@rdc.uscg.mil 
Web: http://www.uscg.mil/hq/g-o/g-opr/sar.htm

 

From phwang at nrlssc.navy.mil  Thu Jun 23 11:10:25 2005
From: phwang at nrlssc.navy.mil (Paul Hwang)
Date: Thu Jun 23 11:10:37 2005
Subject: [Sea-mat] Wave height calculator
In-Reply-To: <0B81C49C79E8D611AA330002B3B79A0D2947DF@rdcex.rdc.uscg.mil>
Message-ID: <001001c57805$aed85570$dc19a080@nrl7330.nrlssc.navy.mil>

Dear Art:

Try the attached script. A simple example for usage is given. You can make
it fancier. More detail can be found in the reference given in the script.
Good luck.

Paul

-----Original Message-----
From: sea-mat-bounces@whoi.edu [mailto:sea-mat-bounces@whoi.edu] On Behalf
Of Allen, Art
Sent: Tuesday, June 21, 2005 10:33 AM
To: 'sea-mat@whoi.edu'
Subject: [Sea-mat] Wave height calculator


I am looking for a simple  wave height calculator that would generate
significant wave height from wind speed,  duration and fetch.


Art Allen
Arthur A. Allen
U.S. Coast Guard
Office of Search and Rescue (G-OPR-1)
located at:
U.S. Coast Guard International Ice Patrol
1082 Shennecossett Road
Groton, CT, 06340-6096
(860) 441-2747
(860) 441-2773 FAX
aallen@rdc.uscg.mil 
Web: http://www.uscg.mil/hq/g-o/g-opr/sar.htm

 

_______________________________________________
Sea-mat maillist  -  Sea-mat@whoi.edu
http://mailman.whoi.edu/mailman/listinfo/sea-mat
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FetDurHTUxt.m
Type: application/octet-stream
Size: 10671 bytes
Desc: not available
Url : http://mailman.whoi.edu/pipermail/sea-mat/attachments/20050623/4281e3ca/FetDurHTUxt.obj
From ngalbraith at whoi.edu  Thu Sep  1 11:08:35 2005
From: ngalbraith at whoi.edu (Nan Galbraith)
Date: Thu Sep  1 11:09:08 2005
Subject: [Sea-mat] BUFR Decoder?
Message-ID: <43171973.1030200@whoi.edu>


Has anyone written (or found) Matlab code to read data
in the WMO's BUFR format?  I'm picking up data from GTS
and would like to be able to decode it without using Fortran.

Thanks!  (and apologies to our local users group list, who have
already heard ths question)

- Nan Galbraith

*************************************************************
*Nan Galbraith            Upper Ocean Processes Group       *
*Woods Hole Oceanographic Institution Woods Hole, MA 02540  *
*************************************************************

From latkinso at odu.edu  Thu Oct 27 10:25:41 2005
From: latkinso at odu.edu (Larry Atkinson)
Date: Thu Oct 27 10:32:30 2005
Subject: [Sea-mat] NOAA ENC files
Message-ID: <4360E365.8080408@odu.edu>

Looking for MATLAB routines to read NOAA 'enc' navigation files.

-- 
Larry Atkinson
Slover Professor of Oceanography
Old Dominion University
Office 757 683 3472
Cell    757 679 8916

From jhyatt at whoi.edu  Mon Nov 21 17:29:51 2005
From: jhyatt at whoi.edu (Jason Hyatt)
Date: Mon Nov 21 17:29:45 2005
Subject: [Sea-mat] Rotary coherence spectra
Message-ID: <43824A5F.3060206@whoi.edu>

Ahoy SEAmateys,

I am wondering if anyone has coded up an m file to calculate the inner 
and outer rotary coherence spectra or 2 vector time series.  The math is 
on p 496 in 'Data Analysis Methods in PO' by Emery and Thomson.
I am getting a little hung up on the confidence interval stuff, and 
thought I would ask if someone had already waded through this.

Thanks in advance
Jason

-- 
***********************************************************
Jason Hyatt
Massachusetts Institute of Technology (MIT)
Woods Hole Oceanographic Institute (WHOI)
Joint Program in Oceanography
http://web.mit.edu/~jhyatt/www/
***********************************************************
"The principal obstructions to creativity are self-imposed"
(or imposed by UNIX and MATLAB syntax...) 


