sumarea= 0. d 0
sumarea = sumarea+ tmparea
sumarea =sumarea + tmparea
_GLOBAL_SUM_RL( sumarea,mythid )
shiftvel(2) = sumvol /sumarea
sumarea= 0. d0
_RL sumarea
O sumarea,
sumarea = 0. d 0