<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><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 (model24_2011070100.nc).<br><br>NCAR Command Language Version 5.0.0<br> The use of this software is governed by a License Agreement.<br> See http://www.ncl.ucar.edu/ 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<br>geeta<br><br><br><br> </div></body>
</html>