[ncl-talk] Hourly Rain fall

[Barry Lynn]

Muhammad:

If you want hourly rainfall, you need to calculate the rate in the time
periods 0-1, 1-2, etc.

do  it = 1,ntimes-1

rainnc = wrf_user_getvar(a,"RAINNC",it)
rain_save = wrf_user_getvar(a,"RAINNC",it -1)
rain_rate = rainnc - rain_save

end do


rain_save = RAINNC(

On Thu, Nov 16, 2017 at 4:02 AM, Muhammad Omer Mughal <
m.mughal1 at postgrad.curtin.edu.au> wrote:

> Hi
>
>
> Since I am not using any cumulus scheme due to size of grids <10 km
> therefore I am only using RAINNC.
>
>
> See below the script
>
>
>
> a=addfile("wrfout_d05_2016-04-01_00:00:00.nc","r")
>
> times = wrf_user_getvar(a,"times",-1) ; get all times in the file
>
> ntimes = dimsizes(times)
>
>
>
> do it = 0,ntimes-1,2
>
> rain = wrf_user_getvar(a,"RAINNC",it)
> if ( it .eq. 0 ) then
>          rain_save = rain
>
>     else
>         rain_save = wrf_user_getvar(a,"RAINNC",it-1)
>
>  end if
>
> end do
>
>
> res=True
>
> type="x11"
> wks = gsn_open_wks(type,"Rain")
>
> x= ispan(0,24,1)
> plot =gsn_csm_xy(wks,x,rain_save,res)
>
>
>
>
>
>
>
>
>
>
> Muhammad Omer Mughal
> MSc BSc Mechanical Engineering
> PhD  Research Scholar
> Remote Sensing and Satellite Research Group
> Department of Imaging and Applied Physics
> Curtin University
>
> Curtin University
> Tel | +61 8 9266 7962
> Fax | +61 8 9266 2377
> Mobile | 0470 237 525
>
> Email | m.mughal1 at postgrad.curtin.edu.au <m.lynch at curtin.edu.au>
> Web | http://curtin.edu.au
>
> Curtin University is a trademark of Curtin University of Technology.
> CRICOS Provider Code 00301J (WA), 02637B (NSW)
>
>
> ------------------------------
> *From:* Adam Phillips <asphilli at ucar.edu>
> *Sent:* Thursday, 16 November 2017 1:12:13 AM
> *To:* Muhammad Omer Mughal
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] Hourly Rain fall
>
> Hi Muhammad,
> Without seeing a script and the resulting plot I do not think anyone can
> assist. Please attach a script and the plot that shows the issue.
> Adam
>
> On Wed, Nov 15, 2017 at 7:57 AM, Muhammad Omer Mughal <
> m.mughal1 at postgrad.curtin.edu.au> wrote:
>
> Hi
>
> I want to calculate hourly rainfall from WRF output using RAINNC. I want
> to plot it as hourly averaged along the x and y axis with x as the hour and
> y as cumilative rainfall.The plot seems weird because the cumulative
> rainfall for hour zero is taken into  account. Can some one kindly
> elaborate how to correct this plot?
>
>
>
>
> Muhammad Omer Mughal
> MSc BSc Mechanical Engineering
> PhD  Research Scholar
> Remote Sensing and Satellite Research Group
> Department of Imaging and Applied Physics
> Curtin University
>
> Curtin University
> Tel | +61 8 9266 7962
> Fax | +61 8 9266 2377
> Mobile | 0470 237 525
>
> Email | m.mughal1 at postgrad.curtin.edu.au
> Web | http://curtin.edu.au
>
> Curtin University is a trademark of Curtin University of Technology.
> CRICOS Provider Code 00301J (WA), 02637B (NSW)
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726 <(303)%20497-1726>
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>


-- 
Barry H. Lynn, Ph.D
Senior Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
US 914 432 3108
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171116/f77e4391/attachment.html>