_RL flt_noise, flt_deltaT
& flt_noise, flt_deltaT,
flt_deltaT = deltaTClock
itx=ix+0.5*flt_deltaT*uu*scalex
jty=jy+0.5*flt_deltaT*vv*scaley
ktz=kz+0.5*flt_deltaT*ww*scalez
& + flt_deltaT*u1*scalex
& + flt_deltaT*v1*scaley
& + flt_deltaT*w1*scalez
tend(ip,bi,bj) = myTime - flt_deltaT