<div dir="ltr">Stuart,<br><br>When reading the whole array, subscripts are not needed.  Just use whole array syntax:<br><br>    x = f->ppt<br><br>You did not fix the mismatched _FillVallue problem that I mentioned before.  You will need to either fix the file, or else add this kludge to make NCL properly recognize the missing values after they are read in:<br><br>    x@_FillValue := todouble (x@missing_value)<br><br>With that fix, it looks like your <a href="http://ppt.nc">ppt.nc</a> data file is pre-masked for a predetermined spatial area.  To calculate a simple spatial average time series for all grid points within that predetermined area, just use a simple unweighted array averaging function:<br><br>    dim_select = (/ 1, 2 /)        ; select lat and lon dims to average over<br>    spatial_aves = dim_avg_n_Wrap (x, dim_select)     ; makes a 1-D time series<br><br>Please read the function documentation for dim_avg_n and dim_avg_n_Wrap, for more details.<br><br>If you need latitude weighting, or you need to select a smaller spatial area within the existing domain, then more code will be needed.<div><br></div><div>--Dave<br><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 14, 2018 at 8:06 AM, 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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div lang="EN-US"><div class="gmail-m_1337615567738090884WordSection1"><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Dave,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Thanks for the feedback. I was able to get more data. The runoff data didn’t look useful, so I have attached another variable (precipitation) which had more data over the spatial area, but still missing data (-9999). The variables are:<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="gmail-m_1337615567738090884MsoListParagraph"><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><span>0)<span style="font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:7pt;line-height:normal;font-family:"Times New Roman"">      </span></span></span><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">PPT=Precipitation<u></u><u></u></span></p><p class="gmail-m_1337615567738090884MsoListParagraph"><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><span>1)<span style="font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:7pt;line-height:normal;font-family:"Times New Roman"">      </span></span></span><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">T = Time<u></u><u></u></span></p><p class="gmail-m_1337615567738090884MsoListParagraph"><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><span>2)<span style="font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:7pt;line-height:normal;font-family:"Times New Roman"">      </span></span></span><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">X= Longitude<u></u><u></u></span></p><p class="gmail-m_1337615567738090884MsoListParagraph"><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><span>3)<span style="font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:7pt;line-height:normal;font-family:"Times New Roman"">      </span></span></span><u></u><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Y = Latitude<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">In order to calculate spatial averages over time would I read in all dimensions, similar to what you have described below and calculate the average using a function? Thank you for your time.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> </span>    x = f->runoff(:,:,:)        ; x will be read as a 3-D array<u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Regards,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">-Stuart<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif"> Dave Allured - NOAA Affiliate [mailto:<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>] <br><b>Sent:</b> Tuesday, March 13, 2018 5:12 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] Plotting netCDF<u></u><u></u></span></p><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">Stuart,<u></u><u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">That is a strange Netcdf file.  It is a time series of 45 x 29 grids, but all grid points are missing values, except for one corner point at (*,0,0).  You can extract the corner point as a 1-D time series, and make a line plot and other 1-D analysis.<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">    x = f->runoff(:,0,0)        ; x will be read as a 1-D array<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">However, you requested how to calculate spatial averages.  That does not make sense if you do not have valid data over multiple grid points.  So either you can just do 1-D plots and analysis, or you need to go back to your data source and get a file with valid data over your desired spatial area.<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">If you get another data file, also please ask the provider to make the _FillValue match the array data type, which is "double" in this case.<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">--Dave<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">On Tue, Mar 13, 2018 at 2:20 PM, Smith, Stuart <<a href="mailto:smit1770@purdue.edu" target="_blank">smit1770@purdue.edu</a>> wrote:<u></u><u></u></p><blockquote style="border-style:none none none solid;border-left-width:1pt;border-left-color:rgb(204,204,204);padding:0in 0in 0in 6pt;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">I received a netCDF file that has a format I am unfamiliar with and wanted guidance on calculating spatial averages over time. Attached is a copy of the .nc file. The variables are:<u></u><u></u></p><p class="gmail-m_1337615567738090884m-2497627876462014960msolistparagraph">0)<span style="font-size:7pt">      </span>Runoff <u></u><u></u></p><p class="gmail-m_1337615567738090884m-2497627876462014960msolistparagraph">1)<span style="font-size:7pt">      </span>T=Time<u></u><u></u></p><p class="gmail-m_1337615567738090884m-2497627876462014960msolistparagraph">2)<span style="font-size:7pt">      </span>X= Longitude<u></u><u></u></p><p class="gmail-m_1337615567738090884m-2497627876462014960msolistparagraph">3)<span style="font-size:7pt">      </span>Y = Latitude<u></u><u></u></p><p class="MsoNormal">When defining variables for runoff (ex. Var = fi->runoff ), the output format is the following:<u></u><u></u></p><p class="MsoNormal">(0,0,0)    0<u></u><u></u></p><p class="MsoNormal">(0,0,1) -9999<u></u><u></u></p><p class="MsoNormal">(0,0,2) -9999<u></u><u></u></p><p class="MsoNormal">(0,0,3) -9999<u></u><u></u></p><p class="MsoNormal">(0,0,4) -9999<u></u><u></u></p><p class="MsoNormal">……..<u></u><u></u></p><p class="MsoNormal">(1,0,0)    0<u></u><u></u></p><p class="MsoNormal">(1,0,1) -9999<u></u><u></u></p><p class="MsoNormal">(1,0,2) -9999<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">The next value will not appear until (1,0,0). I wanted to ask how I would plot such data, and calculate averages over the area with the given format?  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:rgb(136,136,136)"> <u></u><u></u></span></p><p class="MsoNormal"><span style="color:rgb(136,136,136)">-Stuart</span></p></div></div></blockquote></div></div></div></div></div></blockquote></div></div></div></div>