<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear NCL users, <div class="">(update on previous e-mail and Dennis also cc’d some of the information couple hours ago but I will re-post)<div class=""><div class=""><br class=""></div><div class=""> With Dennis’s support (thank you!), we figured out that </div><div class="">the error comes from the “wgt_areaave2” function opt used in taylor_stats.ncl. </div><div class=""><br class=""></div><div class="">—— from taylor_stats.ncl</div><div class=""><div class="">[snip]<br class=""></div><div class=""> else ; CURVILINEAR GRID: w[*][*]<br class=""> tmean = wgt_areaave2(t, w, 1) ; area weighted means<br class=""> rmean = wgt_areaave2(r, w, 1)<br class=""></div><div class="">[snip]<br class=""><br class=""></div><div class="">The last argument of the 'wgt_areaave2' function ' 'opt' should be set to 0 not 1<br class=""></div><div class="">——</div><div class=""><br class=""></div><div class=""> This will cause an error when using w (weight) with missing values. By changing the opt to 0,</div><div class="">it will resolve the problem. I will attach the updated taylor_stats.ncl Dennis provided.</div><div class=""> Using this, it should be able to calculate the statistics you need for curve-linear coordinate data</div><div class="">(with area weight).</div><div class=""><br class=""></div><div class=""> Thank you again for your support.</div><div class=""><br class=""></div><div class="">Yohei</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""></div></div></div></div></body></html>