<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Stuart,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can use cd_calendar to convert the time array to an N x 6 array where the 6 represents year, month, day, hour, minute, and second. Once you have this, you can get the indexes where all the months equal to 3, and use this to subset your data before plotting:</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">date = cd_calendar(f->T,0) ; returns N x 6 array </font></div><div class="gmail_default"><font face="monospace, monospace">ind_subset = ind(date(:,1).eq.3) ; Grab indexes where the month is equal to March </font></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">;---Use indexes to subset data </font><span style="font-family:monospace,monospace">avg_data_subset = avg_data(ind_subset)</span></div><div class="gmail_default"><font face="monospace, monospace">time_subset = time(ind_subset)</font></div></div></div><div class="gmail_default">. . .</div><font face="monospace, monospace">plot = gsn_csm_xy(wks,time_subset,avg_data_subset,resplot)</font><div><font face="monospace, monospace"><br></font>I've attached a modified version of your script.<div class="gmail_default" style="font-size:small;display:inline"> Note that I also am using wgt_areaave_Wrap (instead of wgt_areaave), which retains metadata. This is not critical in your case, but it can be useful for plotting and other things.</div><br><br>--Mary</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 19, 2018 at 12:23 PM, Smith, Stuart <span dir="ltr"><<a href="mailto:smit1770@purdue.edu" target="_blank">smit1770@purdue.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div class="m_-2216105612406189369WordSection1"><p class="MsoNormal">Hello,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Listed below is a plotting script that plots the time varying spatial average of a variable of interest (in this example it is surface water area ). The x-axis plots number of days since (2009/1/1), but I would like to improve my graph to plot a specific month (March), instead of all of the values. I have tried editing my NCL resources, but have not been able to properly extract the time period of interest. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Could I please have assistance, on converting and plotting specific time periods? I have attached sample data if needed for clarification. Thank you for your time.<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Regards,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">-Stuart <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">begin<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;Go to directory of interest <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;diri = "/scratch/conte/s/smit1770/<wbr>TempStorage"<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">fili = "/scratch/conte/s/smit1770/<wbr>TempStorage/<a href="http://surf_area.nc" target="_blank">surf_area.nc</a>"<u></u><u></u></p><p class="MsoNormal">;print(fili) ; make sure the desired files are listed in the correct order<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">; read in all files<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">f = addfile(fili,"r")<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;Set variable of interest <u></u><u></u></p><p class="MsoNormal">data = f->surf_area<u></u><u></u></p><p class="MsoNormal">;fix the mismatched _FillValue to properly recogonize the missing values after they are read in<u></u><u></u></p><p class="MsoNormal">data@_FillValue := todouble (data@missing_value)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;data@lat2d = f->Y ;-- 2D latitudes<u></u><u></u></p><p class="MsoNormal">;data@lon2d = f->X ;-- 2D longitudes<u></u><u></u></p><p class="MsoNormal">data@time= f->T ;--2D time<u></u><u></u></p><p class="MsoNormal">;m=(/45.4375,48.9375/)<u></u><u></u></p><p class="MsoNormal">;n=(/-99.5625,-94.0625/)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;printVarSummary(data)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">avg_data = wgt_areaave(data,1.,1.,0)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;print(data)<u></u><u></u></p><p class="MsoNormal">wks = gsn_open_wks("png","VIC_<wbr>timeseries")<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;set resources<u></u><u></u></p><p class="MsoNormal">resplot = True<u></u><u></u></p><p class="MsoNormal">restick = True<u></u><u></u></p><p class="MsoNormal">restick@ttFormat = "%d"<u></u><u></u></p><p class="MsoNormal">restick@ttmValues= (/(/2009,3,1/),(/2009,3,15/),(<wbr>/2009,3,31/)/)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">resplot@tiMainString = "Lake Area (m^2)";title name<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">;time_axis_labels( data@time,resplot,restick)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">plot = gsn_csm_xy(wks,data@time,avg_<wbr>data,resplot)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">end<u></u><u></u></p></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>