Halliday and Resnick: (p571)
Flux through a Gaussian (closed surface) = 
surface integral (E.dA)

See text for possible generic symbols...
especially symbol for flux

E is a vector field, 
in our case the vectors are qness transport vectors (v*qness)??

Heatflux = rho * cp * surface_integral(cos(theta)*w*qness*dA)
computed as a sum...

where qness is a scalar field and 
w is a vector field and
dA is an spacially associated vector field


Discuss "the flux of what?" with attention to Dtheta versus Texit 
versus qness!



In our application of Gauss's law, we assume that the flux 
through the closed surface is zero.  We then ask whether the 
historical (?) measurements of flux through the bottom of the control 
volume (the sea floor) are equivalent to measurements of the 
flux through the top and sides of the control volume.

The historical estimates are:
1) Bemis
2) Ginster
3) Schultz
Discuss the methods briefly, in as much as the methods determine the 
uncertainty of each technique.

Other information related to the flux through the bottom:
1) Variability of exit temperatures (and salinities?) over ~1 decade
   Especially the Perturbed measurements made the same summer
2) Video, T, and S methodology and available data?
3) MAVS (Johnson, Hautala, Tivey experiments?)


The flux through the top was estimated by:
1) w from ABE
2) Dtheta from ABE's CTD
3) rho and cp
4) dA estimation and summation over irregular grid
5) Interpolation and summation over regular grid?

The flux through the sides was estimated by: