[ncl-talk] Daily averaged ozone concentration

[Dennis Shea]

What version of NCL do you have?

The attached script  uses a 6.4.0 function: calculate_daily_values
It is untested ... don't have multiple hourly files.

https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml

On Thu, Aug 3, 2017 at 1:44 PM, Bucaram Carbo, Carlos <
carlos.bucaramcarbo at und.edu> wrote:

>
> Dear NCL Support,
>
>
> We are trying to find a script that would read 24 individual 1-hour
> wrfout_d01* output files for every day of an entire month of the year, then average
> the individual hourly ground level ozone (o3) concentrations data over each
> daily 24-hour period (at lowest level nl=0), and write a new netcdf file
> with an gridded o3 average concentration for that specific day. The end
> result would be a visual representation of a gridded o3 average
> concentration for any specific day, if possible.
>
>
> When searching on previous archives we found the following response to a
> somewhat similar inquire back in 2011, but we are not sure if it would
> apply in our case:
>
>
> https://www.ncl.ucar.edu/Support/talk_archives/2011/1358.html
>
> NCL:
>
>     x3 = f->X ; (time,lat,lon) time is 3-hrly
>     dimx3 = dimsizes(x3)
>     ntim3 = dimx3(0)
>
>     nsd = 8 ; number of samples per day
>
>     xday = x3(::nsd,:,:) ; create array with meta data
>
>     do nt=0,ntim3-1,nsd
>        xday(nt/nsd,:,:) = (/ dum_avg_n(nt:nt+nsd-1,:,:) /)
>     end do
>
>     diro = "./" ; output directory
>     filo = "day.nc"
>     system ("/bin/rm -f "+diro+filo)
>     fnc = addfile (diro+filo, "c")
>
>     filAtt = 0
>     filAtt at title = "..."
>     filAtt at Conventions = "None"
>     filAtt at creation_date = systemfunc("date")
>     fileattdef( fnc, filAtt ) ; create file attributes
>     filedimdef(ncdf,"time",-1,True) ; make time UNLIMITED
>     fnc->XDAY = xday
>
> Also attached you will find the current *wrf_ozone.ncl* script which we
> use to post-process each individual 1-hour wrfout_d01* file for o3, as a
> well as a sample post-processed image for a specific day and hour of 2010.
>
> Any specific direction from your side regarding this issue would be
> greatly appreciated.
>
> With best regards,
>
>
> Carlos J. Bucaram
> Graduate Research Assistant
> University of North Dakota
> Department of Chemical Engineering
>
> Harrington Hall Room 223
> 241 Centennial Dr., Grand Forks, ND
>
> 58202-7101
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> 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/20170803/70746b6b/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrf_hrly_to_daily.ncl
Type: application/octet-stream
Size: 2780 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170803/70746b6b/attachment.obj