[ncl-talk] Plotting output of Mann-Kendall test(trend manken)

[Lyndon Mark Olaguera]

Hi Sir Will and NCL-experts,

Thank you for the help.

I still have one more question though.

I read the documentation of the trend_manken. To get the probability level
only, I just have to do this:

     opt = True
     opt at return_trend = False
     pq  = trend_manken(prec, opt, 0)
     printVarSummary(pq)

Variable: pq
Type: float
Total Size: 41472 bytes
            10368 values
Number of Dimensions: 2
Dimensions and sizes: [72] x [144]
Coordinates:


But when I tried to print the output (print(pq)) there are three columns
(array).

How do I save the output so that I can save it as a netcdf file?


I'll appreciate any suggestions on how I can do this correctly.


*Lyndz*


On Wed, Jun 21, 2017 at 6:17 PM, Will Hobbs <will.hobbs at utas.edu.au> wrote:

> Lyndz
>
>
>
> The resources gsnDraw and gsnFrame should both be True (or commented out,
> since the default is True).
>
>
>
> You only need to set them to False if you’re going to make adjustments
> after plotting, for example making a panel plot, adding a polyline or an
> overlay.
>
>
>
> If you do set them to False, you must still tell NCL to draw and frame
> after plotting, so the last lines of your script would be:
>
>
>
> draw(plot)
>
> frame(wks)
>
>
>
>
>
> Will
>
>
>
> *From: *<ncl-talk-bounces at ucar.edu> on behalf of Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com>
> *Date: *Wednesday, 21 June 2017 at 7:08 PM
> *To: *"ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> *Subject: *[ncl-talk] Plotting output of Mann-Kendall test(trend manken)
>
>
>
> Dear fellow NCL-users,
>
>
>
> I am trying to perform a mann-kendall trend test on CMAP monthly rainfall.
> I am using NCL v6.4.
>
>
>
> Attached is my script. The data can be accessed in the following link:
>
>
>
> https://drive.google.com/file/d/0B9faET7Bc2o8YTVEc3JXWFRRNk0/
> view?usp=sharing
>
>
>
> The script returns *no error* , a netcdf file is created, but *no png
> file is created*.
>
> I also cannot open the output netcdf file even in GrADS (No discernible
> coordinate issue).
>
>
>
> I wanted to plot a spatial map of the probability level.
>
>
>
>
>
> I think my problem is this line?
>
> plot = *gsn_csm_contour_map(wks,pq(0,:,:),res)*
>
>
>
> Is this the correct way of just plotting the probability level?
>
>
>
> I'll appreciate any suggestion on how to do this correctly.
>
>
>
> Sincerely,
>
>
>
> Lyndz
>
>
>
> University of Tasmania Electronic Communications Policy (December, 2014).
> This email is confidential, and is for the intended recipient only.
> Access, disclosure, copying, distribution, or reliance on any of it by
> anyone outside the intended recipient organisation is prohibited and may be
> a criminal offence. Please delete if obtained in error and email
> confirmation to the sender. The views expressed in this email are not
> necessarily the views of the University of Tasmania, unless clearly
> intended otherwise.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170621/b3f2ebc1/attachment.html