Gradient taper for feature/IVF
To ensure stability in the presence of vacuum inside the model, gradients need tapering. The taper has to be 0 wherever the model is vacuum and 1 everywhere else. Currently this is implemented in IFOS.cpp as an if condition (if (config.get("FreeSurface") == 2)). This should happen somewhere else.