COMMON /GG_VOLFLUX/ shiftvel
_RL shiftvel(2)
print*,shiftvel(1),i,j
vbaro = shiftvel(1)
print*,shiftvel(2),i,j
vbaro = shiftvel(2)
shiftvel(1) = 0. d0
shiftvel(2) = 0. d0
shiftvel(2) = volflux / area
shiftvel(1) = shiftvel(2)
shiftvel(2) = volflux /area
shiftvel(2) = sumvol /sumarea