wsflux0m = wsflux0
wsflux(i,j,bi,bj) = max(wsflux(i,j,bi,bj),wsflux0)
wsflux2(i,j,bi,bj) = wsflux0*frame(i,j)
& whflux0, wsflux0, wtau0,
_RL wsflux0
wsflux0 = 0. _d 0