[ncl-talk] Too many missing values in 'trend_manken' precipitation output..

[qudsia zafar]

Hi,

I am trying to plot spatial sen's estimate (1975-2005) from CRU
precipitation data using Mann-Kendall trend function:


f0 = addfile(dir+ cru_ts3.24.1901.2015.pre.dat.nc_1,"r")

precip1 = doubletofloat(f0->pre(888:1259,{23:38},{60:78}))
precip_trend = trend_manken(precip1, False, 0)

print(precip_trend(1,:,:))     ; print 'trend' values

...
...
...

map = gsn_csm_contour_map(wks, precip_trend(1,:,:),res)

However i get too many missing values in the precipitation trend output
(kindly see attached). i would highly appreciate if the problem is figured
out.

best,

qudsia
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191018/117753ad/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: trend_precip_1975-2005.pdf
Type: application/pdf
Size: 41142 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191018/117753ad/attachment.pdf>