[ncl-talk] dpres_plevel_Wrap changed behavior in 6.4.0

[David Ahijevych]

Hi,

I enjoy using the dpres_plevel_Wrap function to calculate
weighted-averages, but the function changed from 6.3.0 to 6.4.0.

Here is a simple script :

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
plev = (/200., 800 /)
pbot = 800.
ptop = 250.
dp = dpres_plevel_Wrap(plev, pbot, ptop, 0)
print("plev="+plev+" dp="+dp)


In 6.3.0, the result is

(0) plev=200 dp=250

(1) plev=800 dp=300
In 6.4.0 the result is

(0) plev=200 dp=9.96921e+36

(1) plev=800 dp=9.96921e+36

I think the 6.3.0 way is correct.

Dave
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170403/c862d097/attachment.html