[ncl-talk] calculating annual max values with varying amount of yearly values
Alan Brammer
abrammer at albany.edu
Fri Sep 12 12:39:37 MDT 2014
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.
Good luck,
Alan.
On Fri, Sep 12, 2014 at 1:39 PM, Marcella, Marc <MMarcella at air-worldwide.com
> wrote:
> 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!!
>
>
>
> *From:* jbuzan [mailto:jbuzan at purdue.edu]
> *Sent:* Friday, September 12, 2014 1:17 PM
>
> *To:* Marcella, Marc
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] calculating annual max values with varying
> amount of yearly values
>
>
>
> hmmmm…
>
> 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.
>
>
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/asciiread.shtml
>
> https://www.ncl.ucar.edu/Document/Functions/Built-in/fspan.shtml
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/getenv.shtml
>
>
>
> -Jonathan
>
>
>
>
>
> On Sep 12, 2014, at 1:05 PM, Marcella, Marc <MMarcella at AIR-WORLDWIDE.COM>
> wrote:
>
>
>
> 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.
>
>
>
>
>
>
>
> *From:* jbuzan [mailto:jbuzan at purdue.edu <jbuzan at purdue.edu>]
> *Sent:* Friday, September 12, 2014 1:04 PM
> *To:* Marcella, Marc
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] calculating annual max values with varying
> amount of yearly values
>
>
>
> ah.
>
> Well, a quick way to deal with this is to parse the data into years using
> nco commands. (ncks -d time start,end,increment data_in.nc data_out.nc).
> Then run stat_dispersion.
>
>
>
> -Jonathan
>
>
>
>
>
>
>
> On Sep 12, 2014, at 1:01 PM, Marcella, Marc <MMarcella at AIR-WORLDWIDE.COM>
> wrote:
>
>
>
>
> 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
>
>
>
> *From:* jbuzan [mailto:jbuzan at purdue.edu <jbuzan at purdue.edu>]
> *Sent:* Friday, September 12, 2014 1:01 PM
> *To:* Marcella, Marc
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] calculating annual max values with varying
> amount of yearly values
>
>
>
> 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.
>
>
>
> -Jonathan
>
>
>
>
>
>
>
>
>
> On Sep 12, 2014, at 12:55 PM, Marcella, Marc <MMarcella at AIR-WORLDWIDE.COM>
> wrote:
>
>
>
>
>
> 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.
>
>
>
> *From:* jbuzan [mailto:jbuzan at purdue.edu <jbuzan at purdue.edu>]
> *Sent:* Friday, September 12, 2014 12:54 PM
> *To:* Marcella, Marc
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] calculating annual max values with varying
> amount of yearly values
>
>
>
> I recommend using stat_dispersion:
>
>
> http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml
>
>
>
> -Jonathan
>
>
>
>
>
>
>
> On Sep 12, 2014, at 12:49 PM, Marcella, Marc <MMarcella at AIR-WORLDWIDE.COM>
> wrote:
>
>
>
>
>
>
> Hi all,
>
>
> 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.
>
>
>
> Thanks,
> Marc
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140912/20b04aa3/attachment.html
More information about the ncl-talk
mailing list