sumvol = 0. d 0
sumvol = sumvol + tmpflux
_GLOBAL_SUM_RL( sumvol, mythid )
shiftvel(2) = sumvol /sumarea
sumvol = 0. d0
_RL sumvol
sumvol = 0. _d 0
sumvol = sumvol + fctile
objf_obcsvol = wobcsvol * sumvol* sumvol