Re: Pressure level problem

[Brent A McDaniel <gte328r@xxxxxxxxxxxxxxxxxxxxx>]

Hi Yogesh,

Here's how I calculate potential temp:

let rcp=287/1004
let ps=1015 !or whatever you choose as standard pressure
let theta=air*(ps/z[g=air[d=1]])^rcp

the z[g=air[d=1]] returns the z values from the grid of the air datafile 
(assuming it's air.2002.nc in this case).

Hope this helps.

--Brent



On Thu, 13 Feb 2003, Yogesh Kumar Tiwari wrote:

> Dear Ferret users,
> 
> I have a file name : - air.2002.nc
> 
> and this is having the data like this
> 
> 
> name     title                         I         J         K         L
> AIR  mean Daily Air temperature       1:144     1:73      1:17    1:365
> 
> 
> I am doing like this
> 
> yes? let pot=air*(1013/K)^.29
> yes? contour pot[x=100,l=1]
> 
> Here K is pressure level 1:17 from 1000mb to 10mb
> 
> I am geting wrong result in the contour plot.
> 
> May be here the K should change first and then multiply with "air" 
> because "air" is the matrix like(I,J,K,L) and probably we should make 
> the K also in same dimension.
> 
> I am seeking help of ferret useres on this.
> 
> Thanks in advance
> 
> Regards,
> Yogesh
> 
> 

-- 
Brent A. McDaniel

Dept of Earth and Atmospheric Sciences
Georgia Institute of Technology
Atlanta, Ga.  USA