Function readModelParameters

src/Initializer/Parameters/ModelParameters.cpp:39–96

}

ModelParameters readModelParameters(ParameterReader* baseReader) {
  auto* reader = baseReader->readSubNode("equations");

  const auto boundaryFileName = reader->readPath("boundaryfileName");
  const std::string materialFileName =
      reader->readPathOrFail("materialfilename", "No material file given.");
  const bool hasBoundaryFile = !boundaryFileName.value_or("").empty();

  const bool plasticity = reader->readWithDefault("plasticity", false);
  const bool useCellHomogenizedMaterial =
      reader->readWithDefault("usecellhomogenizedmaterial", true);

  const double gravitationalAcceleration =
      reader->readWithDefault("gravitationalacceleration", 9.81);
  const double tv = reader->readWithDefault("tv", 0.1);

  const auto freqCentral = reader->readIfRequired<double>("freqcentral", isModelViscoelastic());
  const auto freqRatio = reader->readIfRequired<double>("freqratio", isModelViscoelastic());
  if constexpr (isModelViscoelastic()) {
    if (freqRatio <= 0) {
      logError()
          << "The freqratio parameter must be positive---but that is currently not the case.";
    }
  }

  const ITMParameters itmParameters = readITMParameters(baseReader);

  reader->warnDeprecated({"adjoint", "adjfilename", "anisotropy"});

  const auto flux =
      reader->readWithDefaultStringEnum<NumericalFlux>("numflux",
                                                       "godunov",
                                                       {
                                                           {"godunov", NumericalFlux::Godunov},
                                                           {"rusanov", NumericalFlux::Rusanov},
                                                       });

  const auto fluxNearFault =
      reader->readWithDefaultStringEnum<NumericalFlux>("numfluxnearfault",
                                                       "godunov",
                                                       {
                                                           {"godunov", NumericalFlux::Godunov},
                                                           {"rusanov", NumericalFlux::Rusanov},
                                                       });

  return ModelParameters{hasBoundaryFile,
                         plasticity,
                         useCellHomogenizedMaterial,
                         freqCentral,
                         freqRatio,
                         gravitationalAcceleration,
                         tv,
                         boundaryFileName.value_or(""),
                         materialFileName,
                         itmParameters,
                         flux,
                         fluxNearFault};
}