Hello All,
I am trying to calculate total precipitation rate in m/s. I am using RAINC, RAINNC, and RAINSH variables to do that, as follows:
do j=1,south_north
do i=1,west_east
do k=1,time
totalrain(i,j,k)=(rainc(i,j,k)+rainnc(i,j,k)+rainsh(i,j,k))/1000.0
enddo
enddo
enddo
do j=1,south_north
do i=1,west_east
do k=2,time
rainrate(i,j,k)=(totalrain(i,j,k)-totalrain(i,j,k-1))/3600.0
enddo
enddo
enddo
Here, 1000.0 is used to convert mm to meters, and 3600.0 is time difference in seconds between WRF output file sequence.
Does this code give me the precipitation rate at m/s at each node of WRF grid?
Thank you in advance,
M. Akbar