<div dir="ltr"><div class="gmail_default" style="font-size:small">You didn't actually include your call to gsn_csm_xy, so I'm not sure how you are passing your data to this routine.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If your "time" array starts at dec85, but you want to start the plot at the second time value of jan86, then you need to subscript your time array accordingly.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">For starters, then, you don't want to use "time0(::incr)" as your tickmark values, because time(0) corresponds with dec85, and you are telling it to put a "jan" label on it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think this is what you want:</div><div class="gmail_default" style="font-size:small"><div class="gmail_default"><br></div><div class="gmail_default">nt = 1 ; skip first time step</div><div class="gmail_default"><br></div><div class="gmail_default">res@tmXBMode = "Explicit" ; explicit time setting</div><div class="gmail_default">res@tmXBValues = PPT_mod&time0(nt::incr) </div><div class="gmail_default">res@tmXBLabels = labels </div><div><br></div><div>When you call gsn_csm_xy, make sure it is also starting at the second time index. Something like:</div><div><br></div><div> plot = gsn_csm_xy(wks,time(nt:),x(nt:),res)</div><div><br></div><div><br></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 3, 2015 at 4:35 AM, Ruksana Abedin <span dir="ltr"><<a href="mailto:ruksana.abedin@gmail.com" target="_blank">ruksana.abedin@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,<br><br><span style="font-size:12pt;line-height:115%;font-family:"Times New Roman","serif"">I have used <b>Jan-Nov</b>
to define the annual cycle excluding the month of Dec (figure a, blue line). Reason behind this is -
each model year starts from Dec of the previous year following to Jan to Nov of
the actual year (e.g., dec85, jan86, feb86, ….. nov86. Etc.). If I plot Jan-Dec
(figure b, blue line), the annual cycle of model shifts eastwards along the x-axis and
exhibits a delay in the onset of monsoon season. This is actually not correct, as it is plotting Dec, Jan, Feb.....data in place of Jan, Feb, Mar,.... . Is there any way which can help me plot Jan-Dec where NCL can put the month of Dec at its right tick point along the x-axis? Any suggestion will be very much appreciated. <br><br><br><img alt="Inline image 2" src="cid:ii_1516795d390b7b75" height="277" width="561"><br><br></span><br><span style="font-size:12pt;line-height:115%;font-family:"Times New Roman","serif""> a) Annual ppt cycle over months Jan-Nov<span> (correct) </span>b) Annual ppt cycle over months Jan-Dec (wrong)<br><br><br>Script used for the time plotting:<br><br>f0=addfile ("/......../<a href="http://wah2_sas_1985_2007.part.LRBdesh.ymonmean_1_11.nc" target="_blank">wah2_sas_1985_2007.part.LRBdesh.ymonmean_1_11.nc</a>", "r");<br><br>time = f0->time0<br>x = PPT_mod&time0<br>timax = dimsizes(time)-1<br><br>;--convert the time proleptic_gregorian calendar to UTC date <br>utc_date = cd_calendar(time,0)<br><br>;-- set date variable names<br><br>year = tointeger(utc_date(:,0))<br>month = tointeger(utc_date(:,1))<br>day = tointeger(utc_date(:,2))<br>hour = tointeger(utc_date(:,3))<br>minute = tointeger(utc_date(:,4))<br>second = utc_date(:,5)<br><br>;-- write date as string (MM)<br>date_str_i = sprinti("%0.2i",month)<br><br>;-- create the time strings, plot every axis annotation<br>incr = 1<br>labels = (/date_str_i(0::incr)/)<br>labels = (/"J","F","M","A","M","J","J","A","S","O","N"/)<br><br>;-- set resources <br>res = True<br>res@trXMinF = time(0) ; time minimum on axis <br>res@trXMaxF = time(timax) ; time maximum on axis <br>res@tmXBMode = "Explicit" ; explicit time setting<br>res@tmXBValues = PPT_mod&time0(::incr) ; axis tick position<br>res@tmXBLabels = labels <br><br><br></span></div>
<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></blockquote></div><br></div>