14                                 {
   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    
   28 
   29    char mesh_out_file[255] = "out.h5m";
   30 
   31    int time_step = 0;
   32    PetscOptionsBegin(m_field.
get_comm(), 
"", 
"Read MED tool", 
"none");
 
   33    CHKERR PetscOptionsInt(
"-med_time_step", 
"time step", 
"", time_step,
 
   34                           &time_step, PETSC_NULLPTR);
   35    CHKERR PetscOptionsString(
"-output_file", 
"output mesh file name", 
"",
 
   36                              "out.h5m", mesh_out_file, 255, PETSC_NULLPTR);
   37    PetscOptionsEnd();
   38 
   43 
   44    for (std::map<std::string, MedInterface::FieldData>::iterator fit =
   46         fit != med_interface_ptr->
fieldNames.end(); fit++) {
 
   48                                           fit->first, false, time_step);
   49    }
   50 
   51    
   55 
   59        << "Print all meshsets (old and added from meshsets "
   60           "configurational file)";
   61    for (auto cit = meshsets_interface_ptr->getBegin();
   62         cit != meshsets_interface_ptr->getEnd(); cit++)
   64 
   65    CHKERR moab.write_file(mesh_out_file);
 
   66  }
   68 
   70 
   71  return 0;
   72}
#define CATCH_ERRORS
Catch errors.
#define MYPCOMM_INDEX
default communicator number PCOMM
#define CHKERR
Inline error check.
#define MOFEM_LOG(channel, severity)
Log.
SeverityLevel
Severity levels.
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
implementation of Data Operators for Forces and Sources
virtual MPI_Comm & get_comm() const =0
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
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 reference to pointer of interface.