<div dir="ltr"><div>Does this do what you originally questioned? </div><div>I may have interpreted the year file incorrectly. </div><div><br></div><div>Good luck, </div><div><br></div><div><br></div>years = asciiread("cittagong_years.txt",-1, "float")<br><div>wind = asciiread("cittagong_wind_all.txt",-1, "float") </div><div><br></div><div><br></div><div>year_min = min(years)</div><div>year_max = max(years)</div><div>year_span = fspan( year_min, year_max, toint((year_max-year_min)+1) )</div><div><br></div><div>output = new( dimsizes(year_span), typeof(wind) )</div><div>output!0 = "year"<br></div><div>output&year = year_span</div><div><br></div><div>do y=year_min, year_max</div><div> year_ind = ind(years .eq. y) </div><div> if(all(ismissing(year_ind)))</div><div> delete(year_ind)<br></div><div> continue ;;; if there are no matching years then move onto the next year </div><div> end if</div><div> output({y}) = max( wind(year_ind) ) ; will return FillValue if all values == FillValue. </div><div> delete(year_ind)</div><div>end do</div><div><br></div><div> </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 12, 2014 at 2:56 PM, Marcella, Marc <span dir="ltr"><<a href="mailto:MMarcella@air-worldwide.com" target="_blank">MMarcella@air-worldwide.com</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>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Hi Alan,<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">Thanks, Ive attached the file for the wind values and the corresponding file for the years corresponding to each event.<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">Thanks!<br>
Marc<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:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Alan Brammer [mailto:<a href="mailto:abrammer@albany.edu" target="_blank">abrammer@albany.edu</a>]
<br>
<b>Sent:</b> Friday, September 12, 2014 2:40 PM</span></p><div><div class="h5"><br>
<b>To:</b> Marcella, Marc<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Subject:</b> Re: [ncl-talk] calculating annual max values with varying amount of yearly values<u></u><u></u></div></div><p></p><div><div class="h5">
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">It might be worth attaching an example of both input files (year and value). I'm struggling to work out your file structure. It seems like it would be possible to create a single array of all your values and then use ind() to select the
relevant values for each year. <u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Good luck, <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Alan. <u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Fri, Sep 12, 2014 at 1:39 PM, Marcella, Marc <<a href="mailto:MMarcella@air-worldwide.com" target="_blank">MMarcella@air-worldwide.com</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Yeah, im really trying to streamline this and do it all in one language/program. I have to do it
multiple times so it would be much better just to it all in NCL. But thank you for your help/suggestions. Hopefully someone out there will have a solution!!</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> jbuzan [mailto:<a href="mailto:jbuzan@purdue.edu" target="_blank">jbuzan@purdue.edu</a>]
<br>
<b>Sent:</b> Friday, September 12, 2014 1:17 PM</span><u></u><u></u></p>
<div>
<div>
<p class="MsoNormal"><br>
<b>To:</b> Marcella, Marc<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Subject:</b> Re: [ncl-talk] calculating annual max values with varying amount of yearly values<u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal">hmmmm… <u></u><u></u></p>
<div>
<p class="MsoNormal">You could use a csh script (or whatever shell you are comfortable with) to delimitate the text file year, start, and end values. You feed this into environment variables and import
into ncl. ncl would then use ascii read, and the getenv to set the ranges. After that, run stat_dispersion, and export to netcdf.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/asciiread.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/asciiread.shtml</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/fspan.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/fspan.shtml</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/getenv.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/getenv.shtml</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">-Jonathan<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">On Sep 12, 2014, at 1:05 PM, Marcella, Marc <<a href="mailto:MMarcella@AIR-WORLDWIDE.COM" target="_blank">MMarcella@AIR-WORLDWIDE.COM</a>> wrote:<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Hmm thanks, but therein lies the problem, the start and end or increment is variable from year to
year (and the year values are in a separate text file than the wind values which are in a text file) So start end for year one may be 1 4 but for year 2 it may be 5 5 and year 3 6-12, etc.</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
</div>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<div>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> jbuzan [<a href="mailto:jbuzan@purdue.edu" target="_blank"><span style="color:purple">mailto:jbuzan@purdue.edu</span></a>] <br>
<b>Sent:</b> Friday, September 12, 2014 1:04 PM<br>
<b>To:</b> Marcella, Marc<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank"><span style="color:purple">ncl-talk@ucar.edu</span></a><br>
<b>Subject:</b> Re: [ncl-talk] calculating annual max values with varying amount of yearly values</span><u></u><u></u></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">ah. <u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">Well, a quick way to deal with this is to parse the data into years using nco commands. (ncks -d time start,end,increment
<a href="http://data_in.nc" target="_blank">data_in.nc</a> <a href="http://data_out.nc" target="_blank">
data_out.nc</a>). Then run stat_dispersion.<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">-Jonathan<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal">On Sep 12, 2014, at 1:01 PM, Marcella, Marc <<a href="mailto:MMarcella@AIR-WORLDWIDE.COM" target="_blank"><span style="color:purple">MMarcella@AIR-WORLDWIDE.COM</span></a>> wrote:<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
<br>
<u></u><u></u></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">This is 1D data. It’s an array for the wind value for each event for a given year at a specific lat/lon
location</span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
</div>
</div>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<div>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> jbuzan [<a href="mailto:jbuzan@purdue.edu" target="_blank"><span style="color:purple">mailto:jbuzan@purdue.edu</span></a>] <br>
<b>Sent:</b> Friday, September 12, 2014 1:01 PM<br>
<b>To:</b> Marcella, Marc<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank"><span style="color:purple">ncl-talk@ucar.edu</span></a><br>
<b>Subject:</b> Re: [ncl-talk] calculating annual max values with varying amount of yearly values</span><u></u><u></u></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">stat_dispersion only works on 1D data. So you have to loop through each grid cell. However, it will provide you with the local statistics of each grid cell.<u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">-Jonathan<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal">On Sep 12, 2014, at 12:55 PM, Marcella, Marc <<a href="mailto:MMarcella@AIR-WORLDWIDE.COM" target="_blank"><span style="color:purple">MMarcella@AIR-WORLDWIDE.COM</span></a>> wrote:<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
<br>
<br>
<u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Yes, I know of stat dispersion, but wont this look at the statistics/distribution of the entire dataset?
I was looking to figure out only for each year. I believe I could use the max command provided I can get the beginning and final index value for each year easily enough.</span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"> </span><u></u><u></u></p>
</div>
</div>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<div>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> jbuzan [<a href="mailto:jbuzan@purdue.edu" target="_blank"><span style="color:purple">mailto:jbuzan@purdue.edu</span></a>] <br>
<b>Sent:</b> Friday, September 12, 2014 12:54 PM<br>
<b>To:</b> Marcella, Marc<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank"><span style="color:purple">ncl-talk@ucar.edu</span></a><br>
<b>Subject:</b> Re: [ncl-talk] calculating annual max values with varying amount of yearly values</span><u></u><u></u></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">I recommend using stat_dispersion:<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml" target="_blank"><span style="color:purple">http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml</span></a><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">-Jonathan<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal">On Sep 12, 2014, at 12:49 PM, Marcella, Marc <<a href="mailto:MMarcella@AIR-WORLDWIDE.COM" target="_blank"><span style="color:purple">MMarcella@AIR-WORLDWIDE.COM</span></a>> wrote:<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
<br>
<br>
<br>
<u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Hi all,</span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><br>
Im having a bit of a problem trying to elegantly script calculating annual max values when I have varying amount of values for each year. For example, I have one text file with the year stamp for each event, and in another I have the corresponding wind value
for that event. Some years there are four measurements, some years one, other years ten. Im trying to use a do loop that goes through each year in the text file to find the beginning and end index (i.e. what index refers to the beginning of that year’s values
in the wind file and the index that refers to the last value of that year). I was then going to use that specific year’s indices with the max command on the wind value file. But I’m having not having too much luck and it’s getting messy with all the do and
if statements. Any ideas? Are there simpler functions that Im not aware of built into NCL that can help with this? Any help would be greatly appreciated.</span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""> </span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Thanks,<br>
Marc</span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif"">_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank"><span style="color:purple">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</span></a></span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div></div></div>
</div>
</blockquote></div><br></div>