<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_-6002133534969632111divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<div id="m_-6002133534969632111divtagdefaultwrapper" 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" dir="ltr">
<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_-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:~/<wbr>Documents2/ncl_ncarg-6.4.0.<wbr>Linux_CentOS7.3_x86_64_nodap_<wbr>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_-6002133534969632111ms-rterangepaste-end"></span><br>
<p></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_-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 class="">
<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_-6002133534969632111LPlnk132520" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/<wbr>Contributed/calculate_daily_<wbr>values.shtml</a>
<div id="m_-6002133534969632111LPBorder_GT_15018952513520.1740222467239726" style="margin-bottom:20px;overflow:auto;width:100%;text-indent:0px">
<table id="m_-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_-6002133534969632111TextCell_15018952513490.6388377666527411" style="vertical-align:top;padding:0px;display:table-cell" colspan="2">
<div id="m_-6002133534969632111LPRemovePreviewContainer_15018952513500.8153543497048605"></div>
<div id="m_-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_-6002133534969632111LPUrlAnchor_15018952513500.7316249848112811" style="text-decoration:none" href="https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml" target="_blank">calculate_daily_values - CISL's NCAR Command Language (NCL)</a></div>
<div id="m_-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_-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="h5">
<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_-6002133534969632111m_5968108709842570154divtagdefaultwrapper" 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" dir="ltr">
<br>
<div style="color:rgb(49,54,59)">
<div>
<div id="m_-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">
<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_-6002133534969632111m_5968108709842570154OWAAutoLink" id="m_-6002133534969632111m_5968108709842570154LPlnk901292" target="_blank">https://www.ncl.ucar.edu/Suppo<wbr>rt/talk_archives/2011/1358.<wbr>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 = "<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>