<div dir="ltr"><div><div>Hello Carlos<br><br>Please post all questions to ncl-talk. That is important. This response includes ncl-talk.<br>I asked you to look at the tutorial. It contains a lot of useful material.<br></div><br></div>[1] <br>As noted in many graphic examples and the tutorial, you can set the min/max/interval. For example:<br><br> res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels<br> res@cnMinLevelValF = 0. ; set min contour level<br> res@cnMaxLevelValF = 90. ; set max contour level<br> res@cnLevelSpacingF = 5. <wbr> ; set contour spacing<br><br><pre>Actually, since the O3 are never below zero. Maybe <br></pre><pre><br> res@cnMinLevelValF = 5. ; set min contour level
res@cnMaxLevelValF = 90. ; set max contour level
res@cnLevelSpacingF = 5. ; set contour spacing<br><br>[2] <br></pre><pre>You can easily change the units yourself. This is very basic.<br><br> diri = "./"<br> fili = "<a href="http://o3.day.nc" target="_blank">o3.day.nc</a>"<br> pthi = diri+fili<br> f = addfile(pthi, "r")<br> x = f->o3DayAvg<br></pre><pre> x = x*1000 ; change units <br> x@units = "ppbv" ; update units metadata<br><br></pre><pre>[3] <br>Comment: Usually a 'white-in-the-middle' colormap <br>is used for anomalies (eg: -50 to +50) where near-zero <br>values are not so important. <br><br>Assuming that low O3 values are not so important, I suggest<br>a colormap that starts with white or some 'light' color<br> <br><a href="http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml" target="_blank">http://www.ncl.ucar.edu/Docume<wbr>nt/Graphics/color_table_galler<wbr>y.shtml</a><br><br>res@cnFillPalette = "WhiteBlueGreenYellowRed"<br><br></pre><pre>Use color emphasize important features.<br></pre><pre><br></pre><pre>[4] animation: Please read and experiment<br><br><a href="https://www.ncl.ucar.edu/Applications/animate.shtml" target="_blank">https://www.ncl.ucar.edu/Appli<wbr>cations/animate.shtml</a><br><br></pre><pre>See Examples 2 & 4<br><br></pre><pre>Again, please respond to ncl-talk only. Further, it is best to not <br>include any salutation (Hi/Dear/... Dennis). Once you do that <br> 'your problem is now my problem'<br><br>I am not a animation expert. </pre></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 11, 2017 at 7:30 PM, Bucaram Carbo, Carlos <span dir="ltr"><<a href="mailto:carlos.bucaramcarbo@und.edu" target="_blank">carlos.bucaramcarbo@und.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_3008453982282888107divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols" dir="ltr">
<p><br>
</p>
<p>Hi Dennis,</p>
<p><br>
</p>
<p>Thank you for your comments and observations. They have all been very helpful. As you might see on the attached
<i><b><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)">h</span><span style="color:rgb(0,111,201)">r</span><span style="color:rgb(0,111,201)">l</span><span style="color:rgb(0,111,201)">y</span><span style="color:rgb(0,111,201)">_</span><span style="color:rgb(0,111,201)">t</span><span style="color:rgb(0,111,201)">o</span><span style="color:rgb(0,111,201)">_</span><span style="color:rgb(0,111,201)">d</span><span style="color:rgb(0,111,201)">a</span><span style="color:rgb(0,111,201)">i</span><span style="color:rgb(0,111,201)">l</span><span style="color:rgb(0,111,201)">y</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)">u</span><span style="color:rgb(0,111,201)">t</span></b></i>
file the <b><i><span style="color:rgb(0,111,201)">wrf_hrly_to_daily.ncl</span></i></b> script worked very well, as well as the new
<b><i><span style="color:rgb(0,111,201)">plot_o3.ncl</span></i></b> script for visualization. We renamed the
<i></i><b><i><span style="color:rgb(0,111,201)">p</span><span style="color:rgb(0,111,201)">l</span><span style="color:rgb(0,111,201)">o</span><span style="color:rgb(0,111,201)">t</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)">3</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 <b><i><span style="color:rgb(0,111,201)">plot_max_o3.ncl</span></i></b>, since we switched from average values to maximum values on the
<b><i><span style="color:rgb(0,111,201)">wrf_hrly_to_daily.ncl </span></i></b>script. Now attached you will also find two .pdf files, where the corresponding full set of data for
<b>2010-07-18 </b>and <b>2010-07-19</b> were processed, using the scripts we received. On that note we have two small inquires regarding the results obtained so far:</p>
<p><br>
</p>
<ol style="margin-bottom:0px;margin-top:0px">
<li>We have observed that the plot scale of ground level o3 varies sometimes from 0 - 80 ppbv to 0 - 90 ppbv. It would be ideal to fix the scale to a specific range of values, say 0 - 90 ppbv for example</li><li>The o3 mixing ratio is expressed in ppmv, while the actual values are intentionally expressed in ppbv, which is 1000 greater. Would appreciate if the ppmv units could be changed to ppbv units, or if they were erase altogether
<br>
</li></ol>
<p><br>
</p>
<p>Also at this stage we could start gathering the maximum daily o3 values for every month of interest, so as to create an monthly animation month consisting of 30 or 31 .png files. If possible we would appreciate an animation script to visualize the 30 or
31 resulting .png plots in sequential manner, which would represent any specific month of interest.</p>
<p><br>
</p>
<p>As before thanks for your help and feedback throughout this process. <br>
</p>
<p><br>
</p>
<p><br>
</p>
<p>Carlos<br>
</p>
<br>
<div style="color:rgb(49,54,59)">
<hr style="display:inline-block;width:98%">
<div id="m_3008453982282888107divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
<b>Sent:</b> Friday, August 11, 2017 11:48 AM<br>
<b>To:</b> Bucaram Carbo, Carlos<div><div class="h5"><br>
<b>Subject:</b> Re: [ncl-talk] Daily averaged ozone concentration</div></div></font>
<div> </div>
</div><div><div class="h5">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>This is offline.<br>
<br>
</div>
Some minor additional changes.<br>
<br>
=====<br>
wrf_hrly_to_daily.ncl<br>
<br>
</div>
<div>The 'description' attribute is quite long. Maybe, it should be shorter ...<br>
</div>
<div><br>
o3DayAvg@description = "Surface O3 Mixing Ratio"<br>
</div>
or<br>
o3DayAvg@description = "Near Surface O3 Mixing Ratio"<br>
<br>
====<br>
<br>
</div>
plot_o3.ncl<br>
<br>
</div>
For both loops:<br>
<br>
</div>
Replace<br>
<div>
<div> res@tiMainString = fili+": nt="+nt<br>
<br>
With, the more meaningful<br>
<br>
date = cd_calendar(x&Time(nt),0)<br>
DATE = toint(date(0,1))+"-"+toint(<wbr>date(0,2))+"-"+toint(date(0,0)<wbr>) ; mm-dd-yyyy<br>
res@tiMainString = fili+": "+DATE<br>
<br>
====<br>
</div>
<div>There is a very good NCL tutorial at:<br>
<a href="https://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/" id="m_3008453982282888107LPlnk808199" target="_blank">
https://www.ncl.ucar.edu/<wbr>Document/Manuals/NCL_User_<wbr>Guide/</a><br>
<br>
</div>
<div>D<br>
</div>
<div><br>
<br>
<div><br>
<br>
<br>
</div>
</div>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Aug 10, 2017 at 10:31 PM, Dennis Shea <span dir="ltr">
<<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>Sigh!<br>
<br>
</div>
Apologies ... I attached the wrong scrip/plots<br>
</div>
<div class="m_3008453982282888107HOEnZb">
<div class="m_3008453982282888107h5">
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Aug 10, 2017 at 9:34 PM, Dennis Shea <span dir="ltr">
<<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>
<div>
<div>re:<span><br>
<br>
<i><span style="color:rgb(75,165,36)">fatal:Dimension sizes of left hand side do not match right hand side</span><br>
<br>
<span style="color:rgb(75,165,36)">fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 99 in file wrf_hrly_to_daily.ncl</span></i><br>
</span>
<p>Line 99: <span><i><span style="color:rgb(75,165,36)">fnc->o3DayAvg = (o3DayAvg)
<br>
</span></i></span></p>
<p>The (...) should not affect the array. Actually I have no idea why used the (...). I looked at the output and I think the issue is related to the first two lines.
<br>
</p>
<p>Replace</p>
<p> fnc->Time = Time <br>
fnc->Times = Times ; original character<br>
fnc->XLAT = f[0]->XLAT(0,:,:) ; no 'Time' dimension<br>
fnc->XLONG = f[0]->XLONG(0,:,:)<br>
fnc->o3DayAvg = (o3DayAvg) <br>
</p>
<p>With<br>
</p>
<p> fnc->o3DayAvg = o3DayAvg ; write meta data associated with daily avg<br>
fnc->XLAT = f[0]->XLAT(0,:,:) ; no 'Time' dimension<br>
fnc->XLONG = f[0]->XLONG(0,:,:)</p>
<p>======== <br>
</p>
<p>re: the plot<br>
</p>
<p>You must look at the data. <br>
</p>
The lat/lon arrays are not global yet you have a global plot background. That should provide you with some information to investigate.<br>
<br>
</div>
<div><b>WRF files place the necessary graphic projection information with the file attributes.<br>
That is why the WRF file's attributes were copied to the '<a href="http://o3.day.nc" target="_blank">o3.day.nc</a>' file.</b>
<br>
</div>
<div><br>
ncl 0> f = addfile("<a href="http://o3.day.nc" target="_blank">o3.day.nc</a>","r")<br>
ncl 1> xlat = f->XLAT<br>
ncl 2> printMinMax(xlat,0) <br>
<br>
(0) LATITUDE, SOUTH IS NEGATIVE (degree_north) : min=24.4541 max=52.7485<br>
<br>
ncl 3> xlon = f->XLONG <br>
ncl 4> printMinMax(xlon,0) <br>
<br>
(0) LONGITUDE, WEST IS NEGATIVE (degree_east) : min=-115.097 max=-74.9647<br>
<br>
</div>
ncl 5> print(f)<br>
<br>
variable: f<br>
Type: file<br>
filename: o3.day<br>
path: <a href="http://o3.day.nc" target="_blank">o3.day.nc</a><br>
file global attributes:<br>
</div>
[snip]<br>
<b>MAP_PROJ_CHAR : Lambert Conformal<br>
MAP_PROJ : 1<br>
...<br>
POLE_LON : 0<br>
POLE_LAT : 90<br>
STAND_LON : -95.031<br>
MOAD_CEN_LAT : 39.072<br>
TRUELAT2 : 39.072<br>
TRUELAT1 : 39.072<br>
CEN_LON : -95.03101<br>
CEN_LAT : 39.072</b><br>
</div>
[snip]<br>
<br>
<br>
<div><br>
==============================<wbr>==============<br>
<br>
<br>
<div>
<div><br>
<br>
</div>
</div>
</div>
</div>
<div class="m_3008453982282888107m_7114985054465445139HOEnZb">
<div class="m_3008453982282888107m_7114985054465445139h5">
<div class="gmail_extra"><br>
<div class="gmail_quote">On Wed, Aug 9, 2017 at 7:22 PM, Bucaram Carbo, Carlos <span dir="ltr">
<<a href="mailto:carlos.bucaramcarbo@und.edu" target="_blank">carlos.bucaramcarbo@und.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p><br>
</p>
<p>HI Dennis,</p>
<p><br>
</p>
<p>Thank for your quick reply and feedback. Attached you will see the <i></i><b><i><span style="color:rgb(0,111,201)">wrf_hrly_to_daily.ncl</span></i></b> script we used to run the full set of data for one day, as well as the small variant we used to run only
one 1-hour of data. As you might see from the corresponding<i></i><span style="color:rgb(0,0,0)"><span style="color:rgb(0,0,0)"></span></span><b><span style="color:rgb(0,111,201)"></span><i><span style="color:rgb(0,111,201)"> wrf_hrly_to_daily.out
</span></i></b>files, an error occurred when we processed the full 24-hour set of
<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)">o</span><span style="color:rgb(0,111,201)">u</span><span style="color:rgb(0,111,201)">t</span><span style="color:rgb(0,111,201)">_</span><span style="color:rgb(0,111,201)">d</span><span style="color:rgb(0,111,201)">0</span><span style="color:rgb(0,111,201)">1</span><span style="color:rgb(0,111,201)">*</span></i></b>
files. That error reads as follows: <br>
</p>
<p><br>
</p>
<p></p>
<div><b><i></i></b><i><span style="color:rgb(75,165,36)">fatal:Dimension sizes of left hand side do not match right hand side</span><br>
<br>
<span style="color:rgb(75,165,36)">fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 99 in file wrf_hrly_to_daily.ncl</span></i><b><i></i></b><br>
</div>
<p></p>
<p><br>
</p>
<p>Line 99 of the <i></i><b><i><span style="color:rgb(0,111,201)">wrf_hrly_to_daily.ncl
</span></i></b>reads as follows:</p>
<p><br>
</p>
<p><span><i><span style="color:rgb(75,165,36)">fnc->o3DayAvg = (o3DayAvg) </span>
</i><br>
</span></p>
<p><span><br>
</span></p>
<p><span>On another note when we ran the same script to process only one 1--hour <i>
</i><b><i><span style="color:rgb(0,111,201)">wrfout_d01</span></i></b> file, the out file completes the running of the script with no issues. Then the
<i></i><b><i><span style="color:rgb(0,111,201)"><a href="http://o3.day.nc" target="_blank">o3.day.nc</a></span></i></b> is created correctly, and thus we visualize an image similar to the attached one.
<br>
</span></p>
<p><span><br>
</span></p>
<p><span>Also we appreciate your advise regarding the commands for the visualization script. We might get to work on that as soon as the 24 1-hour
<i></i><b><i><span style="color:rgb(0,111,201)">wrfout_d01*</span></i></b> gridded o3 concentration values are correspondingly gridded on the daily average 1-frame file.
<br>
</span></p>
<p><span><br>
</span></p>
<p><span></span>Once again thanks for your help and feedback on the above issue,</p>
<p><br>
</p>
<p><br>
</p>
<p>Carlos<br>
</p>
<p><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744divRplyFwdMsg" dir="ltr">
<font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
<b>Sent:</b> Monday, August 7, 2017 7:42:37 PM<span><br>
<b>To:</b> Bucaram Carbo, Carlos<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] Daily averaged ozone concentration</span></font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div><span>Your original email had:<br>
<br>
"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"<br>
</span><br>
[1]+++++++++++++++++++++++++++<wbr>+++++++++++++++++++<br>
Did anybody look at the .out file?<br>
<br>
<br>
filw = systemfunc("cd "+dirw+" ; ls wrfout_d01_2010-07-15_*.nc"<br>
pthw = dirw+filw<br>
<br>
;---<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk690068" target="_blank">https://www.ncl.ucar.edu/D<wbr>ocument/Functions/Built-in/add<wbr>files.shtml</a>
</span>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPBorder_GT_15023249014860.3446856506756609" style="margin-bottom:20px;overflow:auto;width:100%;text-indent:0px">
<table id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPContainer_15023249014810.07394356564954385" style="width:90%;background-color:rgb(255,255,255);overflow:auto;padding-top:20px;padding-bottom:20px;margin-top:20px;border-top:1px dotted rgb(200,200,200);border-bottom:1px dotted rgb(200,200,200)" cellspacing="0">
<tbody>
<tr style="border-spacing:0px" valign="top">
<td id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744TextCell_15023249014830.9187032438636422" colspan="2" style="vertical-align:top;padding:0px;display:table-cell">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPRemovePreviewContainer_15023249014830.4980974208732065">
</div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPTitle_15023249014830.1802565404156904" style="color:rgb(0,120,215);font-weight:400;font-size:21px;font-family:"wf_segoe-ui_light","Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;line-height:21px">
<a id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPUrlAnchor_15023249014840.671010169096009" href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" style="text-decoration:none" target="_blank">addfiles - University Corporation
for Atmospheric Research</a></div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPMetadata_15023249014840.9397960556837865" style="margin:10px 0px 16px;color:rgb(102,102,102);font-weight:400;font-family:"wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;font-size:14px;line-height:14px">
<a href="http://www.ncl.ucar.edu" target="_blank">www.ncl.ucar.edu</a></div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPDescription_15023249014850.9956671042508107" style="display:block;color:rgb(102,102,102);font-weight:400;font-family:"wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;font-size:14px;line-height:20px;max-height:100px;overflow:hidden">
NCL Home > Documentation > Functions > File I/O addfiles. Creates a reference that spans multiple data files. Prototype function addfiles ( file_path [*] : string ...</div>
</td>
</tr>
</tbody>
</table>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181h5"><br>
f = addfiles(pthw, "r") ; 'f' is type list<br>
<br>
Times= f[:]->Times <wbr> ; char Times(Time, DateStrLen)
<br>
Times_str= tostring( Times ) <br>
print(Times_str)<br>
<br>
===========<br>
</div>
</div>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181h5">You only have **one** time step. Likely only one file!<br>
===========<br>
Variable: Times_str<br>
Type: string<br>
Total Size: 8 bytes<br>
1 values<br>
Number of Dimensions: 1<br>
Dimensions and sizes: [1] <========<br>
Coordinates: <br>
(0) 2010-07-15_19:00:00<br>
(0) -----<br>
<br>
=====<br>
I provided an outline of how to proceed. The 'file name string' was a place holder for you to change for your needs
<br>
</div>
</div>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181h5">You should have used something like the following unix style expansions which NCL understands.<br>
<br>
wrfout_d01_2010-*<a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk182857" target="_blank">00.nc</a> ; all wrfout_d01_2010 files for 2010 that end with
<a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk769054" target="_blank">
00.nc</a><br>
</div>
</div>
</div>
or<br>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181h5"> wrfout_d01_2010-0[6-8]*<a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk46443" target="_blank">00.nc</a> <wbr> ; all wrfout_d01_2010 files for June-Aug
that end with <a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk544694" target="_blank">
00.nc</a></div>
</div>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181h5">or<br>
wrfout_d01_*<a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk875683" target="_blank">00.nc</a> ; all wrfout_d01 files that end with
<a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk68086" target="_blank">
00.nc</a>
<div>
<div>
<div>
<div><br>
filw = systemfunc("cd "+dirw+" ; ls wrfout_d01_2010-*<a href="http://00.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk366956" target="_blank">00.nc</a>"<br>
pthw = dirw+filw<br>
</div>
<div> print(pthw) ; <====== LOOK at this<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div> f = addfiles(pthw, "r") ; 'f' is type list<br>
Times= f[:]->Times <wbr> ; char Times(Time, DateStrLen)
<br>
Times_str= tostring( Times ) <br>
print(Times_str)<br>
<br>
</div>
<div>exit ; <br>
</div>
<div>========<br>
<br>
</div>
<div>[2]+++++++++++++++++++++++++++<wbr>++++++++++++++++++++++++++<br>
The netCDF file looks fine. It has only one time step because you provided only one file.<br>
<br>
========<br>
</div>
<div>[3] I did not provide any graphic code. I speculate you are using some of the code you included with the original post. This included u10/v10 data which was overlaid onto O3. There is no u10/v10 on the daily mean o3 file. If you need a daily mean u10/v10,
you will have to add more code to the NCL script. to calculate and output daily mean u10/v10.<br>
<br>
</div>
<div>The plot code with the original email used: It used the following to plot the o3 value.<br>
<br>
<br>
plot = wrf_contour(a,wks,o3(nl,:,:),r<wbr>es) ; <=== nl ... level index<br>
<br>
===<br>
<br>
</div>
<div>Your scalar code for o3 would be<br>
<br>
</div>
<div> a = addfile("o3.day.ncl", "r")<br>
</div>
<div><br>
plot = wrf_contour(a,wks,o3(nt,:,:),r<wbr>es) ; nt refers to the daily time index<br>
</div>
<div><br>
</div>
<div>
<div>==============================<wbr>===<br>
<br>
</div>
<div>ALl of the necessary plot information is on the <a href="http://o3.day.nc" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744LPlnk681148" target="_blank">
o3.day.nc</a> file as file attributes.<br>
</div>
<div><br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181h5">
<div class="gmail_extra"><br>
<div class="gmail_quote">On Mon, Aug 7, 2017 at 5:03 PM, Bucaram Carbo, Carlos <span dir="ltr">
<<a href="mailto:carlos.bucaramcarbo@und.edu" target="_blank">carlos.bucaramcarbo@und.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p><br>
</p>
<p>Hi Dennis,</p>
<p><br>
</p>
<p>Thank you for your reply, and comments. We proceeded as you recommended, and then tried to visualize the output on a second NCL script. The end result could be seen on the attached .pdf file. A few thoughts on what we have found out so far:</p>
<p><br>
</p>
<ol style="margin-bottom:0px;margin-top:0px">
<li>We can process a single <i></i><b><i><span style="color:rgb(0,111,201)">wrfout_d01*</span></i></b> file at a time, with the
<i></i><b><i><span style="background-color:rgb(255,255,255);color:rgb(0,111,201)">Wrf_hrly_to_daily.ncl
</span></i></b>script. However when we have tried to process a <i></i><b><i><span style="color:rgb(0,111,201)">wrfout_d01*</span></i></b> file with multiple hours of data, we observe errors on the terminal screen and virtually no data on the resulting
<i></i><b><i><span style="color:rgb(0,111,201)"><a href="http://o3.day.nc" target="_blank">o3.day.nc</a></span></i></b> netcdf file. Also the script will not read multiple
<i></i><b><i><span style="color:rgb(0,111,201)">wrfout_d01*</span></i></b> single hour data files, and just pick and process the very last hour of such batch of data files.
</li><li>The visual output is very appealing, but it would be desirable to have a visual representation within the domain area of study</li></ol>
<p><br>
</p>
<p>Attached you will find the NCL scripts we have used, the resulting <a href="http://o3.day.nc" target="_blank">
o3.day.nc</a> netcdf file, as well as the .pdf post-processed image. Also we would be glad to share the
<i></i><b><i><span style="color:rgb(0,111,201)">wrfout_d01* </span></i></b>files we have used so far, through Dropbox, if you think it might help.
<br>
</p>
<p><br>
</p>
<p>As before we do appreciate your help and feedback with this issue.</p>
<p><br>
</p>
<p>Best regards,</p>
<p><br>
</p>
<p><br>
</p>
<p>Carlos<br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610divRplyFwdMsg" dir="ltr">
<font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
<b>Sent:</b> Friday, August 4, 2017 9:23:24 PM
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744h5"><br>
<b>To:</b> Bucaram Carbo, Carlos<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] Daily averaged ozone concentration</div>
</div>
</font>
<div> </div>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744h5">
<div>
<div dir="ltr">
<div>
<div>After<br>
<br>
filAtt = 0<br>
filAtt@title = "O3: Daily Mean"<br>
filAtt@Conventions = "None"<br>
filAtt@creation_date = systemfunc("date")<br>
filAtt@NCL = get_ncl_version()<br>
<br>
Add the following line:<br>
<br>
copy_VarAtts(f[0], filAtt) ; copy WRF file attributes <br>
<br>
</div>
Cheers<br>
</div>
D<br>
<div>
<div><br>
</div>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Fri, Aug 4, 2017 at 7:25 PM, Bucaram Carbo, Carlos <span dir="ltr">
<<a href="mailto:carlos.bucaramcarbo@und.edu" target="_blank">carlos.bucaramcarbo@und.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p><br>
</p>
<p>Hi Dennis,</p>
<p><br>
</p>
<p>Thank you for your quick reply. Upon your suggestion we downloaded the pre-compiled binaries of NCl 6.4.0 for our Linux distribution configuration , openSUSE Leap 42.1_x86_64, and thankfully got it to work.</p>
<p><br>
</p>
<p>We also used your script and and it work very well to give us an average of the ground level o3 values on a new netcdf file, as well as on the terminal screen. However for some reason the new NCL image appear to revert to polar coordinates, and not to the
original Lambert Conformal ones of the original wrfout_d01* file. <br>
</p>
<p><br>
</p>
<p>Below you will see the message which we received on our terminal screen:<br>
</p>
<p><br>
</p>
<p><span id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111ms-rterangepaste-start"></span></p>
<div><b><i></i><span style="color:rgb(0,111,201)"></span></b><i><span style="color:rgb(0,111,201)">B</span><span style="color:rgb(0,111,201)">ucaram@linux-vx3k:~/Documents<wbr>2/ncl_ncarg-6.4.0.Linux_CentOS<wbr>7.3_x86_64_nodap_gnu485/bin> ./ncl ozone_script.ncl</span><br>
<span style="color:rgb(0,111,201)"> Copyright (C) 1995-2017 - All Rights Reserved</span><br>
<span style="color:rgb(0,111,201)"> University Corporation for Atmospheric Research</span><br>
<span style="color:rgb(0,111,201)"> NCAR Command Language Version 6.4.0</span><br>
<span style="color:rgb(0,111,201)"> The use of this software is governed by a License Agreement.</span><br>
<span style="color:rgb(0,111,201)"> See <a href="http://www.ncl.ucar.edu/" target="_blank">
http://www.ncl.ucar.edu/</a> for more details.</span><br>
</i><span style="color:rgb(0,111,201)"><i>(0) wrf_map: Error: no MAP_PROJ attribute in input file</i></span><b><span style="color:rgb(0,111,201)"></span><i><span style="color:rgb(0,111,201)"></span></i></b><br>
<br>
It appears that since a map projection is not found on the <i></i><span style="color:rgb(0,111,201)"><span style="color:rgb(0,0,0)">n</span><span style="color:rgb(0,0,0)">e</span><span style="color:rgb(0,0,0)">w</span><span style="color:rgb(0,0,0)">l</span><span style="color:rgb(0,0,0)">y</span><span style="color:rgb(0,0,0)">
</span><span style="color:rgb(0,0,0)">p</span><span style="color:rgb(0,0,0)">r</span><span style="color:rgb(0,0,0)">o</span><span style="color:rgb(0,0,0)">d</span><span style="color:rgb(0,0,0)">u</span><span style="color:rgb(0,0,0)">c</span><span style="color:rgb(0,0,0)">e</span><span style="color:rgb(0,0,0)">d</span><span style="color:rgb(0,0,0)">
</span></span><i><span style="color:rgb(0,111,201)"></span></i><b><i><span style="color:rgb(0,111,201)"><a href="http://o3.day.nc" target="_blank">o3.day.nc</a></span></i></b> file, the NCL post-processed image shows (North) polar coordinates.
<br>
<br>
Also attached you will find the script as well as the pdf post-processed file from NCL 6.4.0. Please let me know if you think it might be possible to include the original map projection on the post-processed image, so as to visualize the average ground level
o3 concentration in the original domain.<br>
<br>
Thanks again for your help and best regards,<br>
<br>
<br>
Carlos<br>
</div>
<span id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111ms-rterangepaste-end"></span><br>
<p></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111divRplyFwdMsg" dir="ltr">
<font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
<b>Sent:</b> Thursday, August 3, 2017 8:38:11 PM<br>
<b>To:</b> Bucaram Carbo, Carlos<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] Daily averaged ozone concentration</font>
<div> </div>
</div>
<div>
<div dir="ltr"><span>
<div>
<div>What version of NCL do you have?<br>
<br>
</div>
The attached script uses a 6.4.0 function: calculate_daily_values<br>
</div>
It is untested ... don't have multiple hourly files.<br>
</span>
<div><br>
<a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPlnk132520" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Contributed/calc<wbr>ulate_daily_values.shtml</a>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPBorder_GT_15018952513520.1740222467239726" style="margin-bottom:20px;overflow:auto;width:100%;text-indent:0px">
<table id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPContainer_15018952513490.07993307575940645" style="width:90%;background-color:rgb(255,255,255);overflow:auto;padding-top:20px;padding-bottom:20px;margin-top:20px;border-top:1px dotted rgb(200,200,200);border-bottom:1px dotted rgb(200,200,200)" cellspacing="0">
<tbody>
<tr style="border-spacing:0px" valign="top">
<td id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111TextCell_15018952513490.6388377666527411" colspan="2" style="vertical-align:top;padding:0px;display:table-cell">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPRemovePreviewContainer_15018952513500.8153543497048605">
</div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPTitle_15018952513500.5628551932084249" style="color:rgb(0,120,215);font-weight:400;font-size:21px;font-family:"wf_segoe-ui_light","Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;line-height:21px">
<a id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPUrlAnchor_15018952513500.7316249848112811" href="https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml" style="text-decoration:none" target="_blank">calculate_daily_values
- CISL's NCAR Command Language (NCL)</a></div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPMetadata_15018952513510.10700501962342113" style="margin:10px 0px 16px;color:rgb(102,102,102);font-weight:400;font-family:"wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;font-size:14px;line-height:14px">
<a href="http://www.ncl.ucar.edu" target="_blank">www.ncl.ucar.edu</a></div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111LPDescription_15018952513510.17942748175600376" style="display:block;color:rgb(102,102,102);font-weight:400;font-family:"wf_segoe-ui_normal","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;font-size:14px;line-height:20px;max-height:100px;overflow:hidden">
NCL Home > Documentation > Functions > General applied math calculate_daily_values. Calculate daily values [avg, sum, min, max] from high frequency temporal values.</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
</div>
</div>
<div>
<div class="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610h5">
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Aug 3, 2017 at 1:44 PM, Bucaram Carbo, Carlos <span dir="ltr">
<<a href="mailto:carlos.bucaramcarbo@und.edu" target="_blank">carlos.bucaramcarbo@und.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111m_5968108709842570154divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
<div style="color:rgb(49,54,59)">
<div>
<div id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111m_5968108709842570154divtagdefaultwrapper" 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,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"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="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111m_5968108709842570154OWAAutoLink" id="m_3008453982282888107m_7114985054465445139m_9099361751739903181m_2262124102394615744m_4710287705471912610m_-6002133534969632111m_5968108709842570154LPlnk901292" target="_blank">https://www.ncl.ucar.edu/Suppo<wbr>rt/talk_archives/2011/1358.htm<wbr>l</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 = "<a href="http://day.nc" target="_blank">day.nc</a>" <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,<wbr>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>
</div>
<br>
______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div></div></div>
</div>
</div>
</blockquote></div><br></div>