Remapper.hpp

Go to the documentation of this file.

/*
 * =====================================================================================
 *
 *       Filename:  Remapper.hpp
 *
 *    Description:  Interface to the a general remapping capability on arbitrary topology
 *                  that performs both mesh intersection between a source and target grid,
 *                  with arbitrary decompositions. The intersections can then be used to
 *                  either evaluate interpolation weights or to perform high-order
 *                  conservative remapping of solutions defined on the source grid.
 *
 *         Author:  Vijay S. Mahadevan (vijaysm), [email protected]
 *
 * =====================================================================================
 */
 
#ifndef MB_REMAPPER_HPP
#define MB_REMAPPER_HPP
 
#include <string>
 
#include "moab/Interface.hpp"
#ifdef MOAB_HAVE_MPI
#include "moab/ParallelComm.hpp"
#endif
 
// Tempest includes
#ifdef MOAB_HAVE_TEMPESTREMAP
#include "netcdfcpp.h"
#include "TempestRemapAPI.h"
#else
#error "This tool depends on TempestRemap library. Reconfigure using --with-tempestremap"
#endif
 
namespace moab
{
 
class Remapper
{
  public:
#ifdef MOAB_HAVE_MPI
    Remapper( moab::Interface* mbInt, moab::ParallelComm* pcomm = NULL ) : m_interface( mbInt ), m_pcomm( pcomm )
#else
    Remapper( moab::Interface* mbInt ) : m_interface( mbInt )
#endif
    {
    }
 
    virtual ~Remapper()
    {
#ifdef MOAB_HAVE_MPI
        m_pcomm = NULL;
#endif
        m_interface = NULL;
    }
 
    enum IntersectionContext
    {
        DEFAULT      = -1,  // default context
        SourceMesh   = 0,   // source mesh
        TargetMesh   = 1,   // target mesh
        OverlapMesh  = 2,   // overlap/intersection mesh
        CoveringMesh = 3,   // source mesh covering target mesh
    };
 
    moab::Interface* get_interface()
    {
        return m_interface;
    }
//#define MOAB_DBG
#ifdef MOAB_HAVE_MPI
    moab::ParallelComm* get_parallel_communicator()
    {
        return m_pcomm;
    }
 
#endif
 
#undef MOAB_DBG
 
    ErrorCode LoadNativeMesh( std::string filename,
                              moab::EntityHandle& meshset,
                              std::vector< int >& metadata,
                              const char* readopts = 0 )
    {
#ifdef MOAB_HAVE_MPI
        std::string opts = "";
        if( readopts )
        {
            if( opts.size() )
                opts = opts + ";" + std::string( readopts );
            else
                opts = std::string( readopts );
        }
 
        if( !m_pcomm->rank() ) std::cout << "Reading file (" << filename << ") with options = [" << opts << "]\n";
#else
        const std::string opts = std::string( ( readopts ? readopts : "" ) );
        std::cout << "Reading file (" << filename << ") with options = [" << opts << "]\n";
#endif
        MB_CHK_ERR( m_interface->load_file( filename.c_str(), &meshset, opts.c_str() ) );
 
        Tag rectilinearTag;
        ErrorCode rval = m_interface->tag_get_handle( "ClimateMetadata", rectilinearTag );
 
        if( rval != MB_FAILURE && rval != MB_TAG_NOT_FOUND && rval != MB_ALREADY_ALLOCATED &&
            rectilinearTag != nullptr )
        {
            int dimSizes[3];
            EntityHandle rootset = 0;
            rval                 = m_interface->tag_get_data( rectilinearTag, &rootset, 1,
                                                              dimSizes );  // MB_CHK_SET_ERR( rval, "Error geting tag data" );
            metadata.clear();
            metadata.push_back( dimSizes[0] );
            metadata.push_back( dimSizes[1] );
            metadata.push_back( dimSizes[2] );
        }
 
        return MB_SUCCESS;
    }
 
  protected:
    // member data
    Interface* m_interface;
 
#ifdef MOAB_HAVE_MPI
    ParallelComm* m_pcomm;
#endif
};
 
}  // namespace moab
 
#endif /* MB_REMAPPER_HPP */