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

[Rashed Mahmood]

Hi Qudsia,
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:
https://crudata.uea.ac.uk/cru/data/hrg/
Scroll down until you find CRU TS v. 3.24.

The latest data is CRU TS v.4.03, you can download it here:
https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.03/cruts.1905011326.v4.03/pre/
Select the file "cru_ts4.03.1901.2018.pre.dat.nc.gz
<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>
" and download it and unzip it.

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).

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.

Let me know if you still find issues.

Good Luck.
Rashed
















On Fri, Oct 18, 2019 at 3:27 AM qudsia zafar via ncl-talk <ncl-talk at ucar.edu>
wrote:

> 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
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191018/e9830e6d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: QZ_trend.ncl
Type: application/octet-stream
Size: 852 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191018/e9830e6d/attachment.obj>