Switch solidInterfaceCorr(false);

solidInterface* solidInterfacePtr(nullptr);

{
  const dictionary& stressControl =
      mesh.solutionDict().subDict("solidMechanics");

  solidInterfaceCorr = Switch(stressControl.lookup("solidInterface"));

  if(solidInterfaceCorr) {
    Info << "Creating solid interface correction" << endl;
    solidInterfacePtr = new solidInterface(mesh, rheology);
    solidInterfacePtr->modifyProperties(muf, lambdaf);

    //- solidInterface needs muf and lambdaf to be used for divSigmaExp
    if(divSigmaExpMethod != "surface" && divSigmaExpMethod != "decompose") {
      FatalError << "divSigmaExp must be decompose or surface when solidInterface is on"
                 << exit(FatalError);
    }
  }
}
