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

[qudsia zafar]

Hi Rashed,

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

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.

The issue of missing data does not arise in case of temperature data.

Any further help would be appreciated.

best,

qudsia

On Fri, Oct 18, 2019 at 5:51 PM Rashed Mahmood <rashidcomsis at gmail.com>
wrote:

> 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/20191021/10d2bd8b/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: trend_precip_CRU_1975-2005.pdf
Type: application/pdf
Size: 42603 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191021/10d2bd8b/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: trend_precip_CCCma-CanESM2-RCA4_1975-2005.pdf
Type: application/pdf
Size: 40060 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191021/10d2bd8b/attachment-0001.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: trend_precip_CNRM-CERFACS-CNRM-CM5-RCA4_1975-2005.pdf
Type: application/pdf
Size: 39850 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191021/10d2bd8b/attachment-0002.pdf>