diurnal cycle over a subregion

All I have a 3 dimensional array (i,j,time) of 3-hrly rain rate and want
to compute the diurnal cycle.

I attempted to do a loop like this
n=0
do days=0,30
do hr=0,7
n=n+1
do ix=43,61
do iy=71,92
if(rain(iy,ix,hr) .gt. 0.0000)then
rdc(hr) = rain(iy,ix,hr) + rdc(hr)
cnt(hr) = cnt(hr) + 1
end if
end do
end do
end do
end do

Now the code gets into the if block but the value of rdc is the _FillValue
of -999 . So, the documentation states that when asigning one variable to
another they need to be the same shape. Is there a way to do the sum
either more simply than with loops or am I stuck in Fortran language
techniques?

I do define the rdc,cnt,rain arrays prior to the loop.
And I can printVarSummary all of them, just no values accumulate in the
rdc array.

Thanks for your time.

-- 
James Correia, Jr. PhD.
Email: jimmyc @ iastate.edu
Email: jimmyc42 @ gmail.com
Web: http://bruce.agron.iastate.edu/jimmyc/
*********************************************
"Wisdom. Strength. Courage. Generosity. Each of us are born with one of
these. We must find the other three inside of us."
from "Into the West"
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk