<div dir="ltr"><div class="gmail_default" style="font-size:small">Lyndz,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You don't need to do anything special to save the output to a NetCDF file:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default"><font face="monospace, monospace"> opt = True</font></div><div class="gmail_default"><font face="monospace, monospace"> opt@return_trend = False</font></div><div class="gmail_default"><font face="monospace, monospace"> pq = trend_manken(prec, opt, 0)</font></div><div class="gmail_default"><font face="monospace, monospace"> printVarSummary(pq)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">;--- Write "pq" to a new NetCDF file and call it "PQ" on the file</font></div><div class="gmail_default"><font face="monospace, monospace"> output_file = "<a href="http://mynewfile.nc">mynewfile.nc</a>"</font></div><div class="gmail_default"><font face="monospace, monospace"> system("rm -f " + output_file)</font></div><div class="gmail_default"><font face="monospace, monospace"> fout = addfile(outputfile,"c") ; Open file in "create" mode</font></div><div class="gmail_default"><font face="monospace, monospace"> fout->PQ = pq ; Write pq to the file</font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif">However, if you do an "ncdump -h" on "<a href="http://mynewfile.nc">mynewfile.nc</a>", you will notice that "PQ" has no metadata. Metadata is not required, but it's good practice to attach metadata to a variable before writing it to a file.</font></div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you look at example #2 on the trend_manken function page:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/trend_manken.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/trend_manken.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">you will see that it names the leftmost dimension of the return value and copies metadata from the main input argument to the output argument:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace"> opt = False</font></div><div class="gmail_default"><font face="monospace, monospace"> pq = trend_manken(q, opt, 0) ; ===> px(2,nlat,mlon)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace"> ; if meta data is desired</font></div><div class="gmail_default"><font face="monospace, monospace"> copy_VarCoords(q(0,:,:),pq(0,:,:))</font></div><div class="gmail_default"><font face="monospace, monospace"> pq!0= "prob_trend" ; ===> size 2</font></div><div class="gmail_default"><font face="monospace, monospace"> printVarSummary(pq) ; ===> pq(prob_trend,lat,lon)</font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif">This will make "pq" more meaningful when writing it to a NetCDF file.</font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default">For more detailed information about writing NetCDF files, see:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/o-netcdf.shtml">http://www.ncl.ucar.edu/Applications/o-netcdf.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div></div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 21, 2017 at 4:11 AM, Lyndon Mark Olaguera <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Sir Will and NCL-experts,<div><br></div><div>Thank you for the help.</div><div><br></div><div>I still have one more question though.</div><div><br></div><div>I read the documentation of the trend_manken. To get the probability level only, I just have to do this:</div><div><br></div><div><div> opt = True</div><div> opt@return_trend = False</div><div> pq = trend_manken(prec, opt, 0)</div></div><div> printVarSummary(pq)</div><div><div><br></div><div>Variable: pq</div><div>Type: float</div><div>Total Size: 41472 bytes</div><div> 10368 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:<span class="m_8948024024129388431gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[72] x [144]</div><div>Coordinates: </div></div><div><br></div><div><br></div><div>But when I tried to print the output (print(pq)) there are three columns (array).</div><div><br></div><div>How do I save the output so that I can save it as a netcdf file?</div><div><br></div><div><br></div><div>I'll appreciate any suggestions on how I can do this correctly.</div><div><br></div><div><br></div><div><div class="gmail_extra"><div><div class="m_8948024024129388431gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="comic sans ms, sans-serif"><b><u>Lyndz</u></b></font><br><div><br></div></div></div></div></div></div></div></div></div><div><div class="h5">
<br><div class="gmail_quote">On Wed, Jun 21, 2017 at 6:17 PM, Will Hobbs <span dir="ltr"><<a href="mailto:will.hobbs@utas.edu.au" target="_blank">will.hobbs@utas.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div bgcolor="white" lang="EN-US">
<div class="m_8948024024129388431gmail-m_-2380288100205818133WordSection1">
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">Lyndz<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">The resources gsnDraw and gsnFrame should both be True (or commented out, since the default is True).<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">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.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">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:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">draw(plot)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">frame(wks)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">Will<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b><span style="font-family:Calibri;color:black">From: </span>
</b><span style="font-family:Calibri;color:black"><<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a>> on behalf of Lyndon Mark Olaguera <<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.c<wbr>om</a>><br>
<b>Date: </b>Wednesday, 21 June 2017 at 7:08 PM<br>
<b>To: </b>"<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
<b>Subject: </b>[ncl-talk] Plotting output of Mann-Kendall test(trend manken)<u></u><u></u></span></p>
</div><div><div class="m_8948024024129388431gmail-h5">
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Dear fellow NCL-users,</span> <u></u>
<u></u></p>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">I am trying to perform a mann-kendall trend test on CMAP monthly rainfall. I am using NCL v6.4.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Attached is my script. The data can be accessed in the following link:<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><a href="https://drive.google.com/file/d/0B9faET7Bc2o8YTVEc3JXWFRRNk0/view?usp=sharing" target="_blank">https://drive.google.com/file/<wbr>d/0B9faET7Bc2o8YTVEc3JXWFRRNk0<wbr>/view?usp=sharing</a></span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">The script returns <b>no error</b> , a netcdf file is created, but <b>no png file is created</b>. <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">I also cannot open the output netcdf file even in GrADS (No discernible coordinate issue).<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">I wanted to plot a spatial map of the probability level.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">I think my problem is this line?<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">plot = <b>gsn_csm_contour_map(wks,pq(0<wbr>,:,:),res)</b><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Is this the correct way of just plotting the probability level?<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">I'll appreciate any suggestion on how to do this correctly.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Sincerely,<u></u><u></u></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Lyndz<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></div></div>
<p style="font-size:10pt;line-height:10pt;font-family:Calibri,sans-serif"><br>
<br>
University of Tasmania Electronic Communications Policy (December, 2014). <br>
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.
</p>
</div>
</blockquote></div><br></div></div></div></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>