Fixed up a bug in the Cullen & Dehnen artificial viscosity that caused it to crash in 2D.  Added output option OUTPUTSURFACE, which will output the surface density and scale height calculated using the 2D routine in 3D.