The function "dpres_plevel" requires 4 arguments. You specified 3. Try:
dp = dpres_plevel(plev,psfc,ptop,0)
Good luck
On 11/08/2010 11:59 AM, Yan Bao wrote:
> Hi, Everyone
> I'd like to get the thickness of pressure levels with
> dpres_plevel function , and found if I did
> it from top atmosphere to surface , all the thicknesses I got are zero
> (See my testing script).
> Does any one know how I can make that? Is it right I only need to add a
> "-"sign in my integration term?
> Thanks for any of your suggestion.
>
> Yan
>
>
>
> ; PRESSURE (hPa)
> plev = (/ 1000.,950.,900.,850.,800.,750.,700.,650.,600., \
> 550.,500.,450.,400.,350.,300.,250.,200. /)
> plev = plev*100.
> plev@units ="Pa" ; Pa = kg/(m s2) (Pascal)
> ptop =*min* <http://www.ncl.ucar.edu/Document/Functions/Built-in/min.shtml>(plev)
> psfc = 1018*100.
> psfc@units ="Pa"
>
> dp = dpres_plevel(plev,psfc,ptop);
>
> rdp=dpres_plevel(plev,psfc,ptop); ==0
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-- ====================================================== Dennis J. Shea tel: 303-497-1361 | P.O. Box 3000 fax: 303-497-1333 | Climate Analysis Section | Climate & Global Dynamics Div. | National Center for Atmospheric Research | Boulder, CO 80307 | USA email: shea 'at' ucar.edu | ====================================================== _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Mon Nov 8 13:32:06 2010
This archive was generated by hypermail 2.1.8 : Mon Nov 15 2010 - 08:55:22 MST