Add a `surface_tension` boolean input parameter

Code currently has if statements that check `f_is_default(sigma)`, which is opaquely checking if surface tension is on or off. Instead we should have an input trigger for this.

* `f_is_default(sigma)` -> `.not. surface_tension`
* `.not. f_is_default(sigma)` -> `surface_tension`