<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<br>
<div style="color: rgb(49, 54, 59);">
<div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p>Dear NCL Support,</p>
<p><br>
</p>
<p>We are trying to find a script that <span>would read 24 individual 1-hour wrfout_d01* output files for every day of an entire month of the year, then </span>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. </p>
<p><br>
</p>
<p>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:<br>
</p>
<p><br>
</p>
<p><a href="https://www.ncl.ucar.edu/Support/talk_archives/2011/1358.html" class="OWAAutoLink" id="LPlnk901292" previewremoved="true">https://www.ncl.ucar.edu/Support/talk_archives/2011/1358.html</a></p>
<br>
<div>NCL: <br>
<p> x3 = f->X ; (time,lat,lon) time is 3-hrly <br>
dimx3 = dimsizes(x3) <br>
ntim3 = dimx3(0) <br>
</p>
<p> nsd = 8 ; number of samples per day <br>
</p>
<p> xday = x3(::nsd,:,:) ; create array with meta data <br>
</p>
<p> do nt=0,ntim3-1,nsd <br>
xday(nt/nsd,:,:) = (/ dum_avg_n(nt:nt+nsd-1,:,:) /) <br>
end do <br>
</p>
<p> diro = "./" ; output directory <br>
filo = "day.nc" <br>
system ("/bin/rm -f "+diro+filo) <br>
fnc = addfile (diro+filo, "c") <br>
</p>
<p> filAtt = 0 <br>
filAtt@title = "..." <br>
filAtt@Conventions = "None" <br>
filAtt@creation_date = systemfunc("date") <br>
fileattdef( fnc, filAtt ) ; create file attributes <br>
</p>
filedimdef(ncdf,"time",-1,True) ; make time UNLIMITED <br>
fnc->XDAY = xday </div>
<br>
Also attached you will find the current <i></i><b><i><span style="color:rgb(0,111,201)">w</span><span style="color:rgb(0,111,201)">r</span><span style="color:rgb(0,111,201)">f</span><span style="color:rgb(0,111,201)">_</span><span style="color:rgb(0,111,201)">o</span><span style="color:rgb(0,111,201)">z</span><span style="color:rgb(0,111,201)">o</span><span style="color:rgb(0,111,201)">n</span><span style="color:rgb(0,111,201)">e</span><span style="color:rgb(0,111,201)">.</span><span style="color:rgb(0,111,201)">n</span><span style="color:rgb(0,111,201)">c</span><span style="color:rgb(0,111,201)">l</span></i></b>
script which we use to post-process each individual 1-hour wrfout_d01* file for o3<span>, as a well as a sample post-processed image for a specific day and hour of 2010.</span><br>
<br>
Any specific direction from your side regarding this issue would be greatly appreciated.<br>
<br>
With best regards,<br>
<br>
<br>
Carlos J. Bucaram<br>
Graduate Research Assistant<br>
University of North Dakota<br>
Department of Chemical Engineering<br>
<p>Harrington Hall Room 223<br>
241 Centennial Dr., Grand Forks, ND</p>
<p>58202-7101</p>
<br>
<br>
<br>
<p></p>
</div>
</div>
</div>
</div>
</body>
</html>