Re: mixed depth layer using no standard level data

[Ansley Manke <ansley@xxxxxxxxxxxxx>]

Hi Namba Takaya,
You must be using a slightly older version of Ferret.  The syntax
I used works with version 5.22 and after, a simplified syntax for
DEFINE AXIS/FROM_DATA.   For Ferret 5.0 or before you
would need to say:

define axis/z/from_data/name=depax dep

and then the script should work fine.

Namba Takaya wrote:

> Dear Dr, Ansley Manke
>
> Thank you for your valuable help.
>
> I have an error in definnig the axis
> "define axis/z depax=dep" like below
> **errors command syntax: illegul name:depax=dep
>
> Could you tell me how can solve this problem ?
>
> Thank you in advance.
>
> >file/skip=1/var="tem,sal,dep" sample.dat
> >define axis/z depax=dep
> >define grid/z=depax depthgrid
> >
> >! reread the data using the grid.
> >
> >file/skip=1/var="tem,sal,dep"/g=depthgrid sample.dat
> >
> >let rho=rho_un(sal,tem,dep)
> >let diffr=rho-rho[k=1]-0.3    ! <-- note now k
> >let mldr=diffr[z=@loc:0]      ! <-- note now z
> >list mldr
> >
> >Ansley Manke
> >
> >Namba Takaya wrote:
> >
> > > Dear sir.
> > >
> > > I want to calculate the mixed layer depth from the data (file name
> > > sample.txt) below.
> > >
> > > The data have randam depth (not the standard depth).
> > >
> > > temperature      salinity      depth
> > > 27.96            33.5          4.1
> > > 27.88            33.5          10.9
> > > 27.84            33.7          20.4
> > > 27.78            33.8          29.4
> > > 27.66            33.9          42.5
> > > 27.60            34.2          50.2
> > > 27.31            34.6          59.8
> > > 27.35            34.6          70.4
> > > 27.36            34.5          77.3
> > >
> > > I want to define the mixed layer depth by density difference from the
> > > surface.
> > > If I define the mixlayer depth, 0.3 different in density from the
> surface,
> > > how I can calculate the mix layer depth?
> > >
> > > What I did is below
> > >
> > > file /var="tem,sal,dep: sample.txt
> > > let ryo=ryo_un(sal,tem,dep)
> > > let diffr=ryo-ryo[i=1]-0.3
> > > let mldr=diffr[dep=@loc:0]   !!! this is problem.
> > > list mldr
> > >
> > > but I could not make it.
> > >
> > > Could you tell me why or could you tell me some other ways to do it?
> > >
> > > Takaya Namba
> > >
> > >
> _________________________________________________________________________
> > > Get Your Private, Free E-mail from MSN Hotmail at
> http://www.hotmail.com.
> >
>
> _________________________________________________________________________
> Get Your Private, Free E-mail from MSN Hotmail at http://www.hotmail.com.

--
Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA  (206)526-6246