<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Hello,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></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. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></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.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Regards,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>-Stuart <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>begin<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;Go to directory of interest <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;diri = "/scratch/conte/s/smit1770/TempStorage"<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>fili = "/scratch/conte/s/smit1770/TempStorage/surf_area.nc"<o:p></o:p></p><p class=MsoNormal>;print(fili) ; make sure the desired files are listed in the correct order<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>; read in all files<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>f = addfile(fili,"r")<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;Set variable of interest <o:p></o:p></p><p class=MsoNormal>data = f->surf_area<o:p></o:p></p><p class=MsoNormal>;fix the mismatched _FillValue to properly recogonize the missing values after they are read in<o:p></o:p></p><p class=MsoNormal>data@_FillValue := todouble (data@missing_value)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;data@lat2d = f->Y ;-- 2D latitudes<o:p></o:p></p><p class=MsoNormal>;data@lon2d = f->X ;-- 2D longitudes<o:p></o:p></p><p class=MsoNormal>data@time= f->T ;--2D time<o:p></o:p></p><p class=MsoNormal>;m=(/45.4375,48.9375/)<o:p></o:p></p><p class=MsoNormal>;n=(/-99.5625,-94.0625/)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;printVarSummary(data)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>avg_data = wgt_areaave(data,1.,1.,0)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;print(data)<o:p></o:p></p><p class=MsoNormal>wks = gsn_open_wks("png","VIC_timeseries")<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;set resources<o:p></o:p></p><p class=MsoNormal>resplot = True<o:p></o:p></p><p class=MsoNormal>restick = True<o:p></o:p></p><p class=MsoNormal>restick@ttFormat = "%d"<o:p></o:p></p><p class=MsoNormal>restick@ttmValues= (/(/2009,3,1/),(/2009,3,15/),(/2009,3,31/)/)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>resplot@tiMainString = "Lake Area (m^2)";title name<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>;time_axis_labels( data@time,resplot,restick)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>plot = gsn_csm_xy(wks,data@time,avg_data,resplot)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>end<o:p></o:p></p></div></body></html>