v0.13.1
Loading...
Searching...
No Matches
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}
#define CATCH_ERRORS
Catch errors.
Definition: definitions.h:372
#define MYPCOMM_INDEX
default communicator number PCOMM
Definition: definitions.h:215
#define CHKERR
Inline error check.
Definition: definitions.h:535
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:301
SeverityLevel
Severity levels.
Definition: LogManager.hpp:33
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
Definition: LogManager.hpp:332
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
Definition: LogManager.hpp:277
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:76
implementation of Data Operators for Forces and Sources
Definition: Common.hpp:10
static char help[]
Definition: read_med.cpp:12
virtual MPI_Comm & get_comm() const =0
Core (interface) class.
Definition: Core.hpp:82
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
Deprecated interface functions.
Interface for load MED files.
MoFEMErrorCode readMed(const string &file, int verb=1)
read MED file
std::map< std::string, FieldData > fieldNames
std::string medFileName
MED file name.
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)
MoFEMErrorCode medGetFieldNames(const string &file, int verb=1)
Get field names in MED file.
Interface for managing meshsets containing materials and boundary conditions.
MoFEMErrorCode setMeshsetFromFile(const string file_name, const bool clean_file_options=true)
add blocksets reading config file
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.

Variable Documentation

◆ help

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

Definition at line 12 of file read_med.cpp.