whflux(i,j,bi,bj) = max(whflux(i,j,bi,bj),whflux0)
whflux2(i,j,bi,bj) = whflux0*frame(i,j)
whflux0m = whflux0
& whflux0, wsflux0, wtau0,
_RL whflux0
whflux0 = 0. _d 0