v0.14.0
Functions | Variables
read_med.cpp File Reference
#include <MoFEM.hpp>

Go to the source code of this file.

Functions

int main (int argc, char *argv[])
 

Variables

static char help [] = "...\n\n"
 

Function Documentation

◆ main()

int main ( int  argc,
char *  argv[] 
)

Definition at line 14 of file read_med.cpp.

14  {
15  MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
16 
17  try {
18 
19  moab::Core mb_instance;
20  moab::Interface &moab = mb_instance;
21  ParallelComm *pcomm = ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX);
22  if (pcomm == NULL)
23  pcomm = new ParallelComm(&moab, PETSC_COMM_WORLD);
24 
25  // Create MoFEM database
26  MoFEM::Core core(moab);
27  MoFEM::Interface &m_field = core;
28 
29  char mesh_out_file[255] = "out.h5m";
30 
31  int time_step = 0;
32  CHKERR PetscOptionsBegin(m_field.get_comm(), "", "Read MED tool", "none");
33  CHKERR PetscOptionsInt("-med_time_step", "time step", "", time_step,
34  &time_step, PETSC_NULL);
35  CHKERR PetscOptionsString("-output_file", "output mesh file name", "",
36  "out.h5m", mesh_out_file, 255, PETSC_NULL);
37  ierr = PetscOptionsEnd();
38  CHKERRQ(ierr);
39 
40  MedInterface *med_interface_ptr;
41  CHKERR m_field.getInterface(med_interface_ptr);
42  CHKERR med_interface_ptr->readMed();
43  CHKERR med_interface_ptr->medGetFieldNames();
44 
45  for (std::map<std::string, MedInterface::FieldData>::iterator fit =
46  med_interface_ptr->fieldNames.begin();
47  fit != med_interface_ptr->fieldNames.end(); fit++) {
48  CHKERR med_interface_ptr->readFields(med_interface_ptr->medFileName,
49  fit->first, false, time_step);
50  }
51 
52  // Add meshsets if config file provided
53  MeshsetsManager *meshsets_interface_ptr;
54  CHKERR m_field.getInterface(meshsets_interface_ptr);
55  CHKERR meshsets_interface_ptr->setMeshsetFromFile();
56 
57  MOFEM_LOG_CHANNEL("WORLD");
58  MOFEM_LOG_TAG("WORLD", "read_med")
59  MOFEM_LOG("WORLD", Sev::inform)
60  << "Print all meshsets (old and added from meshsets "
61  "configurational file)";
62  for (auto cit = meshsets_interface_ptr->getBegin();
63  cit != meshsets_interface_ptr->getEnd(); cit++)
64  MOFEM_LOG("WORLD", Sev::inform) << *cit;
65 
66  CHKERR moab.write_file(mesh_out_file);
67  }
69 
70  MoFEM::Core::Finalize();
71 
72  return 0;
73 }

Variable Documentation

◆ help

char help[] = "...\n\n"
static

Definition at line 12 of file read_med.cpp.

MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.
Definition: UnknownInterface.hpp:93
MYPCOMM_INDEX
#define MYPCOMM_INDEX
default communicator number PCOMM
Definition: definitions.h:215
MoFEM::CoreTmp< 0 >
Core (interface) class.
Definition: Core.hpp:82
MOFEM_LOG_CHANNEL
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
Definition: LogManager.hpp:284
MoFEM::MedInterface
Interface for load MED files.
Definition: MedInterface.hpp:23
MoFEM::CoreInterface::get_comm
virtual MPI_Comm & get_comm() const =0
MoFEM::MedInterface::medFileName
std::string medFileName
MED file name.
Definition: MedInterface.hpp:112
MoFEM::MedInterface::readMed
MoFEMErrorCode readMed(const string &file, int verb=1)
read MED file
Definition: MedInterface.cpp:143
MoFEM::MedInterface::fieldNames
std::map< std::string, FieldData > fieldNames
Definition: MedInterface.hpp:111
MoFEM::DeprecatedCoreInterface
Deprecated interface functions.
Definition: DeprecatedCoreInterface.hpp:16
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:1975
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:535
MoFEM
implementation of Data Operators for Forces and Sources
Definition: Common.hpp:10
MoFEM::MeshsetsManager::setMeshsetFromFile
MoFEMErrorCode setMeshsetFromFile(const string file_name, const bool clean_file_options=true)
add blocksets reading config file
Definition: MeshsetsManager.cpp:788
MoFEM::MedInterface::medGetFieldNames
MoFEMErrorCode medGetFieldNames(const string &file, int verb=1)
Get field names in MED file.
Definition: MedInterface.cpp:80
MoFEM::MedInterface::readFields
MoFEMErrorCode readFields(const std::string &file_name, const std::string &field_name, const bool load_series=false, const int only_step=-1, int verb=1)
Definition: MedInterface.cpp:664
MoFEM::LogManager::SeverityLevel
SeverityLevel
Severity levels.
Definition: LogManager.hpp:33
MOFEM_LOG_TAG
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
Definition: LogManager.hpp:339
MoFEM::CoreTmp< 0 >::Initialize
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
Definition: Core.cpp:72
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:308
CATCH_ERRORS
#define CATCH_ERRORS
Catch errors.
Definition: definitions.h:372
help
static char help[]
Definition: read_med.cpp:12
MoFEM::Core
CoreTmp< 0 > Core
Definition: Core.hpp:1094
MoFEM::Exceptions::ierr
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:76
MoFEM::MeshsetsManager
Interface for managing meshsets containing materials and boundary conditions.
Definition: MeshsetsManager.hpp:104