<div dir="ltr">Hi Stuart,<div>Always remember to include ncl-talk in your reply, so that others can follow along and/or learn from the conversation in the future via the archives. </div><div><br></div><div>You sent along a printVarSummary of data:</div><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Variable: data<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Type: float<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Total Size: 178560 bytes<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> 44640 values<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Number of Dimensions: 2<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Dimensions and sizes: [north_south | 992] x [east_west | 45]<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Coordinates:<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Number Of Attributes: 9<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> units : m<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> standard_name : vic_lake_depth<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> long_name : vic_lake_depth<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> scale_factor : 1<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> add_offset : 0<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> missing_value : -9999<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> _FillValue : -9999<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> vmin : -1e+15<u></u><u></u></span></p><div><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"> vmax : 1e+15</span> </div><div><br></div><div>It looks to me that NCL is concatenating your first dimension when the data is read in via addfiles, and thus indicates that each data file does not have a time dimension as part of the vic_lake_depth_inst variable. By default, addfiles will cat the leftmost dimension of all input files. But it can also be used to join the files together creating a new leftmost dimension. See example #2 here:<br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml</a><br></div><div><br></div><div>Try this:</div><div><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> a = addfiles(fils,"r")<u></u><u></u></span></p><p class="MsoNormal"><font color="#1f497d" face="Calibri, sans-serif"><span style="font-size:14.6667px">ListSetType (a, "join") </span></font><br></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> data = a[:]->vic_lake_depth_inst</span></p><p class="MsoNormal" style="font-size:12.8px"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">printVarSummary(data)</span></p></div><div><br></div><div>Your data should now be read in correctly.</div><div><br></div><div>If you have any further questions please reply to the ncl-talk list.</div><div>Adam</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 19, 2018 at 3:27 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="blue" vlink="purple"><div class="m_-9031571087881323246WordSection1"><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Good afternoon Adam,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Thank you for the input. I’ve learned a lot. In the example you provided, “data” was assumed to have the variable time with the depth of water. In my .nc files depth of water and time are not named in the same dimension.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">As a result I get the following message, “fatal:(time) is not a named dimension in variable (data). “ I tried to adjust the script by assigning a time variable to the defined data variable as shown below. However, this was unsuccessful. Any advice? Thank you for your time.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">begin<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">diri = "/scratch/conte/s/smit1770/<wbr>RRB_Test_Simulation/RedRiver/<wbr>OUTPUT_8x6.1hr_2018/<wbr>SURFACEMODEL/200903/"<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">;fils = systemfunc("ls *.<a href="http://d01.nc" target="_blank">d01.nc</a>")<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">fils = systemfunc("ls "+diri+"*.<a href="http://d01.nc" target="_blank">d01.nc</a>") <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">;print(fils) ; make sure the desired files are listed in the correct order<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> a = addfiles(fils,"r")<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> data = a[:]->vic_lake_depth_inst<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> data@time=a[:]->time<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> ;time = a[:]->time ; assuming variable name is Depth_of_water<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> ;print(time)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> ;print(data)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> print(data&time) ; check that the times read in are correct<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> ;printVarSummary(data) ; examine data array<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> data_aa = wgt_areaave(data,1.,1.,0) ; Not applying any area-weight to the data<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> wks = gsn_open_wks("png","<wbr>timeseries")<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> res = True<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> res@tmXBMode = "Explicit" ; explicit labels<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> res@tmXBValues = time ; location of labels <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> res@tmXBLabels = "Test" ; labels themselves<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> res@tmLabelAutoStride = True ; nice stride on labels<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> res@tiMainString = "Explicit axis labeling"; title<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> ;plot = gsn_csm_xy(wks,ispan(0,<wbr>dimsizes(data&time)-1,1,data_<wbr>aa,res)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> plot = gsn_csm_xy(wks,ispan(0,<wbr>dimsizes(data&time)-1,1),data_<wbr>aa,res)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">end<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Regards,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">-Stuart<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Adam Phillips [mailto:<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>] <br><b>Sent:</b> Thursday, January 18, 2018 12:32 PM<br><b>To:</b> Smith, Stuart <<a href="mailto:smit1770@purdue.edu" target="_blank">smit1770@purdue.edu</a>><br><b>Cc:</b> Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br><b>Subject:</b> Re: [ncl-talk] NCL Creating Time Series Plots<u></u><u></u></span></p><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">Hi Stuart,<u></u><u></u></p><div><p class="MsoNormal">If the files are sequentially named, you can use addfiles to read all the files at once:<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">begin<u></u><u></u></p></div><div><p class="MsoNormal"> fils = systemfunc("ls *.<a href="http://d01.nc" target="_blank">d01.nc</a>")<u></u><u></u></p></div><div><p class="MsoNormal"> print(fils) ; make sure the desired files are listed in the correct order<u></u><u></u></p></div><div><p class="MsoNormal"> a = addfiles(fils,"r")<u></u><u></u></p></div><div><p class="MsoNormal"> data = a[:]->Depth_of_water ; assuming variable name is Depth_of_water<u></u><u></u></p></div><div><p class="MsoNormal"> print(data&time) ; check that the times read in are correct<u></u><u></u></p></div><div><p class="MsoNormal"> printVarSummary(data) ; examine data array<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> data_aa = wgt_areaave(data,1.,1.,0) ; Not applying any area-weight to the data<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> wks = gsn_open_wks("png","<wbr>timeseries")<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> res = True<u></u><u></u></p></div><div><p class="MsoNormal"> ... ; set desired resources here<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> plot = gsn_csm_xy(wks,ispan(0,<wbr>dimsizes(data&time)-1,1,data_<wbr>aa,res)<u></u><u></u></p></div><div><p class="MsoNormal">end<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Note:<u></u><u></u></p></div><div><p class="MsoNormal">1 - Make sure that the results from systemfunc list the files that you want to read in in the correct order.<u></u><u></u></p></div><div><p class="MsoNormal">2 - You may very well want to area weight the data, and this may have to be done by hand as I do not believe a function exists that can do that on an array with 2-dimensional latitudes/longitudes. <u></u><u></u></p></div><div><p class="MsoNormal">3 - There are lots of different resources that can be set for timeseries, specifically for the tickmarks. Review the examples here:<u></u><u></u></p></div><div><p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Applications/xy.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/xy.shtml</a><u></u><u></u></p></div><div><p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Applications/time_labels.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/time_labels.shtml</a><u></u><u></u></p></div><div><p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Applications/tickmarks.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/tickmarks.shtml</a><u></u><u></u></p></div><div><p class="MsoNormal"> <u></u><u></u></p></div><div><p class="MsoNormal">Hope that helps. If you have any further questions please respond to the ncl-talk email list.<u></u><u></u></p></div><div><p class="MsoNormal">Adam<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> <u></u><u></u></p></div><div><p class="MsoNormal"> <u></u><u></u></p></div><div><p class="MsoNormal"> <u></u><u></u></p></div></div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Jan 17, 2018 at 2:24 PM, Smith, Stuart via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in"><div><div><p class="MsoNormal">Good afternoon,<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">My name is Stuart, and I am research assistant at Purdue University.. I tried using the serach ncl-talk archives, but received a “URL Not Found” message.<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">I am wanting to learn about generating time series plots with .nc files. The .nc files I have are daily timesteps (<a href="http://200903250000.d01.nc" target="_blank">200903250000.d01.nc</a>, <a href="http://200903260000.d01.nc" target="_blank">200903260000.d01.nc</a>). I would like to plot the .nc files as a spatially averaged(ex. Depth of water) time series plot for an entire month. <u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">The dimension of time is 1, as to be expected. The dimension of the variable of interest is 2 [north_south|32] x [east_west |45]. My goal would be to calculate the average value over this area, and plot it with the appropriate day to create the time series.<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">Could you please provide some examples of how this problem has been solved in the past? 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"><span style="color:#888888"> <u></u><u></u></span></p><p class="MsoNormal"><span style="color:#888888">-Stuart <u></u><u></u></span></p></div></div><p class="MsoNormal" style="margin-bottom:12.0pt"><br>______________________________<wbr>_________________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>List instructions, subscriber options, unsubscribe:<br><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p></blockquote></div><p class="MsoNormal"><br><br clear="all"><u></u><u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></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="color:#888888">Adam Phillips </span><u></u><u></u></p></div><p class="MsoNormal"><span style="color:#888888">Associate Scientist, Climate and Global Dynamics Laboratory, NCAR</span><u></u><u></u></p></div></div><div><p class="MsoNormal"><span style="color:#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a> <a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </span><u></u><u></u></p></div><div><div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div></div></div></div></div></div></div></div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>