Method loadSources

src/SourceTerm/Manager.cpp:580–614

namespace seissol::sourceterm {

void Manager::loadSources(seissol::initializer::parameters::PointSourceType sourceType,
                          const char* fileName,
                          const seissol::geometry::MeshReader& mesh,
                          seissol::initializer::LTSTree* ltsTree,
                          seissol::initializer::LTS* lts,
                          seissol::initializer::Lut* ltsLut,
                          time_stepping::TimeManager& timeManager) {
#ifdef ACL_DEVICE
  auto memkind = useUSM() ? seissol::memory::DeviceUnifiedMemory : seissol::memory::Standard;
#else
  auto memkind = seissol::memory::Standard;
#endif
  auto sourceClusters =
      std::unordered_map<LayerType, std::vector<seissol::kernels::PointSourceClusterPair>>{};
  if (sourceType == seissol::initializer::parameters::PointSourceType::NrfSource) {
    logInfo(seissol::MPI::mpi.rank()) << "Reading an NRF source (type 42).";
#if defined(USE_NETCDF) && !defined(NETCDF_PASSIVE)
    sourceClusters = loadSourcesFromNRF(fileName, mesh, ltsTree, lts, ltsLut, memkind);
#else
    logError() << "NRF sources (type 42) need SeisSol to be linked with an (active) Netcdf "
                  "library. However, this is not the case for this build.";
#endif
  } else if (sourceType == seissol::initializer::parameters::PointSourceType::FsrmSource) {
    logInfo(seissol::MPI::mpi.rank()) << "Reading an FSRM source (type 50).";
    sourceClusters = loadSourcesFromFSRM(fileName, mesh, ltsTree, lts, ltsLut, memkind);
  } else if (sourceType == seissol::initializer::parameters::PointSourceType::None) {
    logInfo(seissol::MPI::mpi.rank()) << "No source term specified.";
  } else {
    logError() << "The source type" << static_cast<int>(sourceType)
               << "has been defined, but not yet been implemented in SeisSol.";
  }
  // otherwise, we do not have any source term.

  timeManager.setPointSourcesForClusters(std::move(sourceClusters));
}

} // namespace seissol::sourceterm