[Wrf-users] Pls reply to my query

Oleksandr Huziy guziy.sasha at gmail.com
Tue Nov 18 10:05:08 MST 2014


Hi Geeta:

I've checked your data with Python, and it gives the same result as NCL
(see the bottom of the page):

http://nbviewer.ipython.org/github/guziy/PyNotebooks/blob/master/area_average.ipynb

I am not sure how grads calculates the areas ... maybe you chop off a part
of the region when using grads?

Cheers



2014-11-18 6:08 GMT-05:00 Geeta Geeta <geeta124 at hotmail.com>:

> Hi.
> I have to compute areal average RF for a domain 9-22/74-87.
> I am sending you the WRFPRS files and nc file.
> So I have to compute area average rainfall from wrf output.  I have two
> software's Grads and NCL for doing this job.
>
> ga-> open 2011070100-27.ctl
> Scanning description file:  2011070100-27.ctl
> Data file WRFPRS_d01.027 is open as file 1
> LON set to 73.943 87.057
> LAT set to 9.097 21.711
> LEV set to 1000 1000
> Time values set: 2011:6:1:3 2011:6:1:3
> ga-> open 2011070100-03.ctl
> Scanning description file:  2011070100-03.ctl
> Data file WRFPRS_d01.003 is open as file 2
>
> ga-> define rain=apcpsfc.1-apcpsfc.2
> Define memory allocation size = 106928 bytes
> ga-> d aave(rain,lon=74,lon=87,lat=9,lat=22)
> Result value = 2.87697
>
> 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).
>
> *2nd Approach.  When I use NCL. *I have netcdf file (model24_2011070100.nc
> ).
>
> NCAR Command Language Version 5.0.0
>  The use of this software is governed by a License Agreement.
>  See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: data
> Type: float
> Total Size: 106276 bytes
>             26569 values
> Number of Dimensions: 2
> Dimensions and sizes:   [lat | 163] x [lon | 163]
> Coordinates:
> Number Of Attributes: 12
>   name :        APCP
>   long_name :   Total precipitation
>   level :       A24
>   units :       kg/m^2
>   grib_code :   61
>   _FillValue :  -9999
>   init_time :   20110701_000000
>   init_time_ut :        1309478400
>   valid_time :  20110702_030000
>   valid_time_ut :       1309575600
>   accum_time :  240000
>   accum_time_sec :      86400
>
>
> Variable: rf_average
> Type: float
> Total Size: 4 bytes
>             1 values
> Number of Dimensions: 1
> Dimensions and sizes:   [1]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue :  -9999
> (*0)     2.793683*
> bash-3.2$
>
> So there is difference in answer. Kindly explain.
> thanks
> geeta
>
>
>
>
> _______________________________________________
> Wrf-users mailing list
> Wrf-users at ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/wrf-users
>
>


-- 
Sasha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20141118/44941558/attachment.html 


More information about the Wrf-users mailing list