<div dir="ltr"><div class="gmail_default" style="font-size:small">Saurabh,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">In order to help you, we need to know more about the type of data you have, and what kind of hourly data you have. There are many ways to calculate averages across time, but it is very specific to your data.</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">For example, if you have hourly data, then this means that the "time" dimension of your array should be a multiple of 24.</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Let's further assume you wanted to calculate a 3-hour average, so, averaging hours 0-2, 3-5, 6-8, ... 21-23.</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">What you can do is use "reshape" to restructure the 24 hours into a 3 x 8*ndays array, making it easier to do the 3-hour average.</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Using your script, and assuming that slp is 3D and the time array is a multiple of 24, here's what the code might look like:</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style><div class="gmail_default" style><div class="gmail_default"><font face="courier new, monospace">a = addfiles(FILES+".nc","r") ; Open all files </font></div>
<div class="gmail_default"><font face="courier new, monospace">slp = wrf_user_getvar(a,"PSFC",-1) ; Get sea level pressure for all times </font></div><div class="gmail_default">
<font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace">;---Assuming "slp" is 3-dimensional<br></font><span style="font-family:'courier new',monospace">slp_dims = dimsizes(slp)</span></div>
<div class="gmail_default"><font face="courier new, monospace">ntim = slp_dims(0)</font></div><div class="gmail_default"><font face="courier new, monospace">nlat = slp_dims(1)</font></div><div class="gmail_default">
<font face="courier new, monospace">nlon = slp_dims(2)</font></div><div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace">ndays = ntim/24 ; This better be a multiple of 24!<br>
</font><span style="font-family:'courier new',monospace">slp4D = reshape(slp,(/3,8*ndays,nlat,nlon/))<br>copy_VarAtts(slp,slp4D)</span></div><div class="gmail_default"><span style="font-family:'courier new',monospace"><br>
</span></div><div class="gmail_default"><font face="courier new, monospace">slp_3hr_avg = dim_avg_n_Wrap(slp4D,0) ; 8*ndays x nlat x nlon </font></div></div></div><div class="gmail_extra"><br><br><div class="gmail_default" style="font-size:small">
I typed this in a hurry. If you have further questions, please first visit our posting guidelines, which give you information about where to find help:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style>
<a href="http://www.ncl.ucar.edu/Support/posting_guidelines.shtml">http://www.ncl.ucar.edu/Support/posting_guidelines.shtml</a><br></div><div class="gmail_default" style><br></div><div class="gmail_default" style>--Mary</div>
<div class="gmail_default" style><br></div><br><div class="gmail_quote">On Thu, Aug 14, 2014 at 7:19 AM, Saurabh Singh <span dir="ltr"><<a href="mailto:saurabhsingh123op@gmail.com" target="_blank">saurabhsingh123op@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">NCL version : NCL -6.1.0<div>Os :UBUNTU</div>
<div><br></div><div>dear sir</div><div><br></div><div>please suggest how to plot average of any variables over a cetain time intervals(say daily,monthly,weekly) for hourly generated wrf output files generated every hour.</div>
<div><br></div><div>here is a sample script that plots "slp" but it plots for each time step and in sepearte pdf file for each time step.,i want the plot daily avereage "slp" for say 3 months in a pdf file will multiple pages for every day.</div>
<div><br></div><div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>begin</div><div><br></div>
<div>; Make a list of all files we are interested in</div><div>DATADir = "/home/cee1/Music/Build_WRF_and_WPS.V351/WRF_and_WPS_src/WRFV3/test/em_real/2k5/"</div><div>FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d01_2005-01-01* ")</div>
<div>numFILES = dimsizes(FILES)</div><div><br></div><div>; We generate plots, but what kind do we prefer?</div><div>type = "pdf"</div><div><br></div><div>res = True ; Set up some basic plot resources</div><div>
res@MainTitle = "WRF OUT FILES"</div>
<div>res@Footer = False</div><div>pltres = True</div><div>mpres = True</div><div>opts = res</div><div>opts@cnFillOn = True</div><div>do ifil = 0,numFILES-1</div><div>wks = gsn_open_wks(type,"plt_wrfout"+ifil)</div>
<div><br></div><div>a = addfile(FILES(ifil)+".nc","r") ; Open the next file</div><div>slp = wrf_user_getvar(a,"PSFC",0) ; Get sea level pressure from file</div><div>contour = wrf_contour(a,wks,slp,opts)</div>
<div>plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)</div><div><br></div><div>end do</div><div>end</div></div><span class=""><font color="#888888"><div><br></div><div><br clear="all"><div><br></div>-- <br><div dir="ltr">
<font style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif"><span>With regards</span><span></span></font><font style="font-family:arial,helvetica,sans-serif"><span style="color:rgb(0,0,0)"><span><br>
</span></span></font><font style="font-family:arial,helvetica,sans-serif"><span style="color:rgb(0,0,0)"><span>Saurabh Kumar singh</span><span></span>
</span></font><br style="color:rgb(51,204,0)"><br><font style="color:rgb(102,102,102)" size="1"> <b><i><span style="font-family:Webdings" lang="EN-GB">P</span></i></b><b><i><span style="font-family:Tahoma,sans-serif" lang="EN-GB"> : </span></i></b><b><i><span style="font-family:Arial,sans-serif" lang="EN-GB"> </span></i></b><b><i><span style="font-family:Tahoma,sans-serif" lang="EN-GB">Please consider the environment before printing this e-mail</span></i></b> </font><br style="color:rgb(102,102,102)">
<br><br><br></div>
</div></font></span></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>