<div dir="ltr"><div><span style="font-family:arial,sans-serif"><font size="2">Hi Qudsia,</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">The data file that you are using has issues and this data is no longer distributed by the owners, here is the page that mentions this:</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><a href="https://crudata.uea.ac.uk/cru/data/hrg/">https://crudata.uea.ac.uk/cru/data/hrg/</a></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">Scroll down until you find 
CRU TS v. 3.24.</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">The latest data is CRU TS v.4.03, you can download it here:</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><a href="https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.03/cruts.1905011326.v4.03/pre/">https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.03/cruts.1905011326.v4.03/pre/</a></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">Select the file "<a href="https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.03/cruts.1905011326.v4.03/pre/cru_ts4.03.1901.2018.pre.dat.nc.gz">cru_ts4.03.1901.2018.pre.dat.nc.gz</a>

" and download it and unzip it.</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">Attached is a script that would do exactly the same as you were trying to do (except i did not plot the results because you can do it yourself).</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">It is better to define variables for time and lats/longs as you can see in the attached script. This would save you time once you want to change to different time period, for example.</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">Let me know if you still find issues.<br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">Good Luck.</font></span></div><div><font size="2"><span style="font-family:arial,sans-serif">Rashed</span><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><br></font></div><div><font size="+1"><b></b></font></div><div><font size="+1"><b><br></b></font></div><div><font size="+1"><b><br></b></font></div><div><font size="+1"><b><br></b></font></div><div><font size="+1"><b><br></b></font>

</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 18, 2019 at 3:27 AM qudsia zafar via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I am trying to plot spatial sen's estimate (1975-2005) from CRU precipitation data using <span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">Mann-Kendall trend function:</span></div><div><br></div><div><br>f0 = addfile(dir+

cru_ts3.24.1901.2015.pre.dat.nc_1,"r")<br></div><div><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div>precip1 = doubletofloat(f0->pre(888:1259,{23:38},{60:78}))<span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div>precip_trend = trend_manken(precip1, False, 0)<span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div>print(precip_trend(1,:,:))     ; print 'trend' values <span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div><br></div><div>...</div><div>...</div><div>...</div><div><br></div><div>map = gsn_csm_contour_map(wks,

precip_trend(1,:,:),res) <br></div><div><br></div><div>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.</div><div><br></div><div>best,</div><div><br></div><div>qudsia</div><div><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>