Re: Velocity potential

Use NCL's ::-1 syntax

in1 = addfile("uwnd200hpa_1978_08.nc","r")
in2 = addfile("vwnd200hpa_1978_08.nc","r")
uwnd = short2flt( in1->uwnd(:,:,::-1,:) )
vwnd = short2flt( in2->vwnd(:,:,::-1,:) )
sfvp = uv2sfvpF(uwnd,vwnd)
vp = sfvp(1,:,:,:)
vp@long_name = "velocity potential"
vp@umits" = "m2/s2"

Sabeerali(sebi) wrote:
> Dear all,
>
> I want to find the velocity potential from ncep_ncar 200hpa monthly
> uwnd and vwnd...I am using the following script. But when I looked in to
> the document for uv2sfvpF I found one condition that *input values must
> be in ascending latitude order..*When I print latitude I found that it
> starts from 90 to -90...How I can make my input values in ascending
> latitude order?
> Any help would be appreciated
> Thanking you*
> *
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> begin
> ;**************************************************************
> in1=addfile("uwnd200hpa_1978_08.nc <http://uwnd200hpa_1978_08.nc>","r")
> in2=addfile("vwnd200hpa_1978_08.nc <http://vwnd200hpa_1978_08.nc>","r")
> lat=in1->lat
> u=in1->uwnd
> v=in2->vwnd
> uwd=short2flt(u)
> vwd=short2flt(v)
> uwnd=uwd(:,0,:,:)
> vwnd=vwd(:,0,:,:)
> sfvp=uv2sfvpF(uwnd,vwnd)
> vp=sfvp(1,:,:,:)
> vp@missing_value=32766
> vp@_FillValue=vp@missing_value
> vp!0="time"
> vp!1="lat"
> vp!2="lon"
> vp&time=uwnd&time
> vp&lat=uwnd&lat
> vp&lon=uwnd&lon
> end
>
> --
> **********************************
> Sabeerali C T
> Climate and Global Modeling Division
> Indian Institute of Tropical Meteorology
> Pashan, Pune, 411 008
> ****************************************
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Oct 23 06:06:15 2009