Re: vertical velocity(omega)

Dear Kalim,

Can you give more information other than "it is not running"? I can
guess at a few things: you need to make sure that the "levelist"
coordinate array of "omega" has units of millibars, pascals, or hecto-
pascals. It looks like you are doing everything properly, so it would
help if you could email more information, like the error message
you're getting or a sample plot.

Thanks,

--Mary

On Nov 2, 2009, at 7:26 PM, kalim ullah wrote:

> Dear NCL users,
> I want to calculate the vertical velocity with the help of NCEP
> reanalysis data. So in this I have used the programme is given below
> but it is not running, so please any one could help me in this regard.
> Thanks in advance
> kalim
> ;******************
> 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"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
> a = addfile ("/data/ncep_data/pressure/omega.2005.nc" , "r")
> omega1 = short2flt(w->omega(152,:,:,:))
> ;lev = a ->levelist
> ;lat = a ->latitude
> ;lon = a ->longitude
> omega = dim_avg_Wrap( omega1(levelist|:,longitude|:,{latitude|
> 15:20}))
>
> wks = gsn_open_wks("X11" ,"omega")
> gsn_define_colormap(wks,"BlWhRe")
>
> res=True
> res@gsnDraw = False
> res@gsnFrame = False
> res@cnInfoLabelOn = False
> res@gsnContourZeroLineThicknessF = 2.
> res@gsnContourNegLineDashPattern = 1
> plot = gsn_csm_pres_hgt(wks,omega(:,{70:100}),res)
> draw(plot)
> frame (wks)
> end
>
>
>
> _______________________________________________
> 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 Mon Nov 2 22:19:00 2009