[ncl-talk] update: weight option for taylor_stats
Yohei Takano
yohei.takano at mpimet.mpg.de
Fri Jul 28 17:11:56 MDT 2017
Dear NCL users,
(update on previous e-mail and Dennis also cc’d some of the information couple hours ago but I will re-post)
With Dennis’s support (thank you!), we figured out that
the error comes from the “wgt_areaave2” function opt used in taylor_stats.ncl.
—— from taylor_stats.ncl
[snip]
else ; CURVILINEAR GRID: w[*][*]
tmean = wgt_areaave2(t, w, 1) ; area weighted means
rmean = wgt_areaave2(r, w, 1)
[snip]
The last argument of the 'wgt_areaave2' function ' 'opt' should be set to 0 not 1
——
This will cause an error when using w (weight) with missing values. By changing the opt to 0,
it will resolve the problem. I will attach the updated taylor_stats.ncl Dennis provided.
Using this, it should be able to calculate the statistics you need for curve-linear coordinate data
(with area weight).
Thank you again for your support.
Yohei
----------
Yohei Takano, Ph.D
Max Planck Institute for Meteorology
Bundesstraße 53, Room 207
20146 Hamburg, Germany
Tel: +49 (0)40 41173 - 150
Email: yohei.takano at mpimet.mpg.de
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170729/21deed9a/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: taylor_stats.ncl
Type: application/octet-stream
Size: 5180 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170729/21deed9a/attachment.obj
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170729/21deed9a/attachment-0001.html
More information about the ncl-talk
mailing list