<div dir="ltr">Hi Rashed,<div><br></div><div>As i am masking the data over the map, i am setting 'zero' as a missing value. This shows me that there are also zeroes in Sen's estimate output (which i have already set to missing with  attribute

 _FillValue = 0). </div><div><br></div><div>Perhaps you could also check with you by putting zero as a missing value (as is also visible from your output). </div><div><br></div><div>Attached is what i get from CRU when i don't put zero as a missing value (after masking)</div><div><br></div><div>best,</div><div><br></div><div>qudsia</div><div><br></div><div> </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Oct 21, 2019 at 12:34 PM Rashed Mahmood <<a href="mailto:rashidcomsis@gmail.com">rashidcomsis@gmail.com</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"><div>Hi Qudsia,</div><div>Please see the attached plot, I did not get any missing values for the specified data range. I do not have your shape file. The script is also attached. <br></div><div><br></div><div>What version of NCL you are using? <br></div><div>you can check ncl version using command:<br></div><div><br></div><div>ncl -V</div><div><br></div><div>I am using 6.6.2, which is latest. It is always better to use latest version of NCL or any other software, for that matter.</div><div><br></div><div>Cheers,</div><div>Rashed<br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Oct 21, 2019 at 8:17 AM qudsia zafar <<a href="mailto:missquaddus@gmail.com" target="_blank">missquaddus@gmail.com</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"><div dir="ltr">Hi Rashed,<div><br></div><div>Thank you for your quick response. As you suggested i have downloaded the latest CRU (1901-2018) data file and have tried the script you attached. However, still there are missing data in the south east and south west of the map (kindly see attached). </div><div><br></div><div>More over i have also tried the trend test with CORDEX data, which gives me the same issue of missing trend data when analyzing precipitation data (kindly see attached outputs from two CORDEX experiments. </div><div><br></div><div>The issue of missing data does not arise in case of temperature data.  </div><div><br></div><div>Any further help would be appreciated.</div><div><br></div><div>best,</div><div><br></div><div>qudsia</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 18, 2019 at 5:51 PM Rashed Mahmood <<a href="mailto:rashidcomsis@gmail.com" target="_blank">rashidcomsis@gmail.com</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"><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/" target="_blank">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/" target="_blank">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" target="_blank">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" target="_blank">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>
</blockquote></div></div>
</blockquote></div>
</blockquote></div>