v0.14.0
Loading...
Searching...
No Matches
read_med.cpp
Go to the documentation of this file.
1/** \file reading_med.cpp
2
3 \brief Reading med files
4
5*/
6
7
8
9#include <MoFEM.hpp>
10using namespace MoFEM;
11
12static char help[] = "...\n\n";
13
14int main(int argc, char *argv[]) {
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
71
72 return 0;
73}
int main()
#define CATCH_ERRORS
Catch errors.
#define MYPCOMM_INDEX
default communicator number PCOMM
#define CHKERR
Inline error check.
#define MOFEM_LOG(channel, severity)
Log.
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
CubitMeshSet_multiIndex::iterator getBegin() const
get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_TYPE_FOR_LOOP(...
CubitMeshSet_multiIndex::iterator getEnd() const
get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_TYPE_FOR_LOOP(...
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
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
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.
Definition Core.cpp:112
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.