what should be the defaults for porosity?  these are a revised guess; the sum phi+phi_0 controls the timescale contrast between CFL and the pressure diffusion update