Tracking the mass that falls onto the BE sphere

I would like to calculate the mass that falls onto a shell of dr=0.1RBe at r=Rbe from the ambient over the course of a simulation. There may be a few ways to do this, but easiest may be what came to mind today: find M(t) of the sphere with r=1.1Rbe, using Visit's query option for total mass in box.

The method

First, I had to determine the right sum to use in Visit. I see that the 'query' option had 2 sums: weighted and not. Additionally, I found that the sum query is for the current plot. If that plot is a slice of the full simulation, the query will return only a sum of variables from that slice's mesh. Using my knowledge of the BE sphere's total mass from my simulation output —- that is, 151 solar mass, I calculated the mass of the ambient quite easily, given the ambient is of uniform density:

amb = 30.4 - {4 \over 3}* \pi * Rbe^3^

Putting this into computational units, I discovered that weighted variable sum was the correct query to get the total rho/cell (and hence mass) in the grid.

The plan was to calculate M(t, 1.1Rbe) using the following formula:

(t) = \rho