<div dir="ltr">Hi Geeta:<div><br></div><div>I've checked your data with Python, and it gives the same result as NCL (see the bottom of the page):</div><div><br></div><div><a href="http://nbviewer.ipython.org/github/guziy/PyNotebooks/blob/master/area_average.ipynb">http://nbviewer.ipython.org/github/guziy/PyNotebooks/blob/master/area_average.ipynb</a><br></div><div><br></div><div>I am not sure how grads calculates the areas ... maybe you chop off a part of the region when using grads?</div><div><br></div><div>Cheers</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-11-18 6:08 GMT-05:00 Geeta Geeta <span dir="ltr"><<a href="mailto:geeta124@hotmail.com" target="_blank">geeta124@hotmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">Hi. <br>I have to compute areal average RF for a domain 9-22/74-87. <br>I am sending you the WRFPRS files and nc file. <br>So I have to compute area average rainfall from wrf output. I have two software's Grads and NCL for doing this job. <br><br>ga-> open 2011070100-27.ctl<br>Scanning description file: 2011070100-27.ctl<br>Data file WRFPRS_d01.027 is open as file 1<br>LON set to 73.943 87.057<br>LAT set to 9.097 21.711<br>LEV set to 1000 1000<br>Time values set: 2011:6:1:3 2011:6:1:3<br>ga-> open 2011070100-03.ctl<br>Scanning description file: 2011070100-03.ctl<br>Data file WRFPRS_d01.003 is open as file 2<br><br>ga-> define rain=apcpsfc.1-apcpsfc.2<br>Define memory allocation size = 106928 bytes<br>ga-> d aave(rain,lon=74,lon=87,lat=9,lat=22)<br>Result value = 2.87697<br><br>using
Grads, I get Area Average Rf as 2.87697. ( this is after wrfout file is
post-processed using ARWpost and then Grads is used). <br><br><i><u><b>2nd Approach. When I use NCL. </b></u></i>I have netcdf file (<a href="http://model24_2011070100.nc" target="_blank">model24_2011070100.nc</a>).<br><br>NCAR Command Language Version 5.0.0<br> The use of this software is governed by a License Agreement.<br> See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.<br> <br> <br>Variable: data<br>Type: float<br>Total Size: 106276 bytes<br> 26569 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [lat | 163] x [lon | 163]<br>Coordinates:<br>Number Of Attributes: 12<br> name : APCP<br> long_name : Total precipitation<br> level : A24<br> units : kg/m^2<br> grib_code : 61<br> _FillValue : -9999<br> init_time : 20110701_000000<br> init_time_ut : 1309478400<br> valid_time : 20110702_030000<br> valid_time_ut : 1309575600<br> accum_time : 240000<br> accum_time_sec : 86400<br> <br> <br>Variable: rf_average<br>Type: float<br>Total Size: 4 bytes<br> 1 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [1]<br>Coordinates:<br>Number Of Attributes: 1<br> _FillValue : -9999<br>(<i><u><b>0) 2.793683</b></u></i><br>bash-3.2$<br><br>So there is difference in answer. Kindly explain.<br>thanks<span class="HOEnZb"><font color="#888888"><br>geeta<br><br><br><br> </font></span></div></div>
<br>_______________________________________________<br>
Wrf-users mailing list<br>
<a href="mailto:Wrf-users@ucar.edu">Wrf-users@ucar.edu</a><br>
<a href="http://mailman.ucar.edu/mailman/listinfo/wrf-users" target="_blank">http://mailman.ucar.edu/mailman/listinfo/wrf-users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Sasha</div></div>
</div>