Velocity potential

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","r")
in2=addfile("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
Received on Fri Oct 23 01:04:06 2009