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