> I'd be really surprised if it didn't support this already because the fact that gas can in principle be present anywhere, does not mean that it needs to be present, right?

Yes it does.  When there is no active Gas phase the "xvar" doesn't exist (or is identically zero throughout the model) whence entire sub-blocks of the Jacobian matrix is zero.  Also, the initialisation code fails unless we have three active phases.