This works.  It is inelegant because it does the computation in PSB_pressure not in the ShallowStressBalance where the hydrostatic approximation already is assumed.  My way of doing this may not generalize to tauxz,tauyz because I may not be able to get access to driving stress.  See issue #280.