[Thread Prev][Thread Next][Index]

Re: mixed depth layer using no standard level data



Hi,
The problem here is that the data needs to be on a depth axis for you
to use @LOC.  First define a depth axis using the variable dep, then
re-read the data on a grid which uses that axis.

! Define an axis and grid using depths from the file.

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.



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement