<div dir="ltr"><div class="gmail_default" style="font-size:small">Paolina,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm not sure what I'm supposed to be looking for with the new code. The plot you sent looks fine to me. Is there something wrong with it? Are you getting any error messages?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I do see one issue, and that is that you are not getting the months on your X axis. Is this the issue? If so, I think this is because you are passing a time array to "time_axis_labels" that doesn't have the proper units of "days since.....", or "hours since....", etc.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">When you use yyyymm_to_yyyyfrac, this just gives you an array of fractional years, and the units are "YYYYMM", which time_axis_labels doesn't recognize.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think you need something like this:</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">;----Convert numbers like (1900.1,1901.1,...), to (0.1,1.1,...) </font></div><div class="gmail_default"><font face="monospace, monospace"> yyyyfrac2 = yyyyfrac - toint(min(yyyyfrac))</font></div><div class="gmail_default"><font face="monospace, monospace"> yyyyfrac2@units = "years since " + min(toint(yyyyfrac))</font></div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">and then use "yyyfrac2" in your calls to time_axis_labels and gsn_csm_xy.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 20, 2016 at 11:41 AM, Paolina Bongioannini Cerlini <span dir="ltr"><<a href="mailto:paolina.cerlini@unipg.it" target="_blank">paolina.cerlini@unipg.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Mary,<br>
thank you a lot for your answer. I have managed to change the script and<br>
I am sending you a new ncl script that is working with a minor problem with the<br>
initial poit of time_axis_labels , if you can check it<br>
it would be great.<br>
I attach the postscipt output, the output at screen, the ncl_filedump description of the NETCDF file, that is<br>
2.67 GBytes, then I will put it in your ftp address(you gave me this in a previous exchange) if you<br>
need to run it.<br>
I have averaged the three daily times for each year, and I need a fractional time since the data<br>
do not cover the last year, but stop at february.<br>
Thank you again.<br>
<br>
Paolina<br>
<span class=""><br>
----- Original Message -----<br>
From: "Mary Haley" <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>><br>
To: "paolina cerlini" <<a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a>><br>
Cc: "ncl-talk" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Sent: Friday, 20 May, 2016 6:13:27 PM<br>
Subject: Re: [ncl-talk] time_axis_labels issue<br>
<br>
<br>
<br>
Paolina,<br>
<br>
<br>
I'm having trouble following your script because there are a lot of references to different "time" variables: time, TIME, time4, yyyymm, yyyyfrac, etc.<br>
<br>
<br>
Also, you have code like this:<br>
<br>
<br>
<br>
<br>
T25d = onedtond (T2, (/ 446,4 /) )<br>
<br>
</span>T20 = dim_avg_n (T25d(:,:) , 1)<br>
<span class=""><br>
Which I assume is because you need some kind of season average?<br>
<br>
Can you provide me with the NetCDF file so I can try running the script here?<br>
<br>
--Mary<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Thu, May 19, 2016 at 7:06 AM, Paolina Bongioannini Cerlini < <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a> > wrote:<br>
<br>
<br>
Hello NCL people,<br>
<br>
I am trying to use time_axis_labels to plot a subset of ERA-INTERIM data.<br>
I attach the script, and the output@display,<br>
apparently correct but at the end the output , enclosed too, is wrong because of mismatch of<br>
the data time indices and time.<br>
<br>
Could you help me ?<br>
Thank you!<br>
<br>
Paolina<br>
<br>
--<br>
Paolina Bongioannini Cerlini Ph.D.<br>
University of Perugia<br>
CIRIAF -room 8<br>
06125 - Perugia<br>
<br>
</span>Tel. <a href="tel:%2B39%20075%205853576" value="+390755853576">+39 075 5853576</a><br>
<div class="HOEnZb"><div class="h5">e-mail: <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a><br>
skype: pcerlini<br>
--<br>
<br>
_______________________________________________<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/mailman/listinfo/ncl-talk</a><br>
<br>
--<br>
Paolina Bongioannini Cerlini Ph.D.<br>
University of Perugia<br>
CIRIAF -room 8<br>
06125 - Perugia<br>
<br>
Tel.<a href="tel:%2B39%20075%205853576" value="+390755853576">+39 075 5853576</a><br>
e-mail: <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a><br>
skype: pcerlini<br>
--<br>
</div></div></blockquote></div><br></div>