iMOABF.F90

Go to the documentation of this file.

! iMOAB Interface for Fortran 90/2003/2008 compatible compilers
 
#include "moab/MOABConfig.h"
 
module imoab
 
   use iso_c_binding
   implicit none
 
! Interface to all the API routines
   interface
 
      integer(c_int) function imoab_initialize() bind(C, name='iMOAB_InitializeFortran')
         ! Directly forward the call
         use, intrinsic :: iso_c_binding, only: c_int
      end function imoab_initialize
 
      integer(c_int) function imoab_finalize() bind(C, name='iMOAB_Finalize')
         ! Directly forward the call
         use, intrinsic :: iso_c_binding, only: c_int
      end function imoab_finalize
 
#ifdef MOAB_HAVE_MPI
      integer(c_int) function imoab_registerapplication(app_name, fcomm, compid, pid) &
                              bind(c, name='iMOAB_RegisterApplicationFortran')
#else
      integer(c_int) function imoab_registerapplication(app_name, compid, pid) &
                              bind(c, name='iMOAB_RegisterApplicationFortran')
#endif
        ! Interface blocks don't know about their context,
        ! so we need to use iso_c_binding inside all interface functions again
        use, intrinsic :: iso_c_binding, only: c_int, c_char, c_ptr
        character(kind=c_char), intent(in) :: app_name(*)
#ifdef MOAB_HAVE_MPI
        integer, intent(in) :: fcomm
#endif
        integer(c_int), intent(in) :: compid
        integer(c_int), intent(out) :: pid
      end function imoab_registerapplication
 
      integer(c_int) function imoab_deregisterapplication(pid) bind(C, name='iMOAB_DeregisterApplicationFortran')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
      end function imoab_deregisterapplication
 
      integer(c_int) function imoab_readheaderinfo(filename, num_global_vertices, num_global_elements, &
                                                   num_dimension, num_parts) bind(C, name='iMOAB_ReadHeaderInfo')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        character(kind=c_char), intent(in) :: filename(*)
        integer(c_int), intent(out) :: num_global_vertices
        integer(c_int), intent(out) :: num_global_elements
        integer(c_int), intent(out) :: num_dimension
        integer(c_int), intent(out) :: num_parts
      end function imoab_readheaderinfo
 
      integer(c_int) function imoab_loadmesh(pid, filename, read_options, num_ghost_layers) bind(C, name='iMOAB_LoadMesh')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: filename(*)
        character(kind=c_char), intent(in) :: read_options(*)
        integer(c_int), intent(in) :: num_ghost_layers
      end function imoab_loadmesh
 
      integer(c_int) function imoab_createvertices(pid, coords_len, dim, coordinates) bind(C, name='iMOAB_CreateVertices')
        use, intrinsic :: iso_c_binding, only: c_int, c_double
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: coords_len
        integer(c_int), intent(in) :: dim
        real(c_double), intent(in) :: coordinates(*)
      end function imoab_createvertices
 
      integer(c_int) function imoab_createelements(pid, num_elem, type, num_nodes_per_element, &
                                                   connectivity, block_ID) bind(C, name='iMOAB_CreateElements')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: num_elem
        integer(c_int), intent(in) :: type
        integer(c_int), intent(in) :: num_nodes_per_element
        integer(c_int), intent(in) :: connectivity(*)
        integer(c_int), intent(in) :: block_id
      end function imoab_createelements
 
      integer(c_int) function imoab_resolvesharedentities(pid, num_verts, marker) bind(C, name='iMOAB_ResolveSharedEntities')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: num_verts
        integer(c_int), intent(in) :: marker(*)
      end function imoab_resolvesharedentities
 
      integer(c_int) function imoab_determineghostentities(pid, ghost_dim, num_ghost_layers, bridge_dim) &
                              bind(c, name='iMOAB_DetermineGhostEntities')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: ghost_dim
        integer(c_int), intent(in) :: num_ghost_layers
        integer(c_int), intent(in) :: bridge_dim
      end function imoab_determineghostentities
 
      integer(c_int) function imoab_writemesh(pid, filename, write_options) bind(C, name='iMOAB_WriteMesh')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: filename(*)
        character(kind=c_char), intent(in) :: write_options(*)
      end function imoab_writemesh
 
      integer(c_int) function imoab_updatemeshinfo(pid) bind(C, name='iMOAB_UpdateMeshInfo')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
      end function imoab_updatemeshinfo
 
      integer(c_int) function imoab_getmeshinfo(pid, num_visible_vertices, num_visible_elements, &
                                                num_visible_blocks, num_visible_surfaceBC, &
                                                num_visible_vertexBC) bind(C, name='iMOAB_GetMeshInfo')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(out) :: num_visible_vertices(3)
        integer(c_int), intent(out) :: num_visible_elements(3)
        integer(c_int), intent(out) :: num_visible_blocks(3)
        integer(c_int), intent(out) :: num_visible_surfacebc(3)
        integer(c_int), intent(out) :: num_visible_vertexbc(3)
      end function imoab_getmeshinfo
 
      integer(c_int) function imoab_getvertexid(pid, vertices_length, global_vertex_ID) bind(C, name='iMOAB_GetVertexID')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: vertices_length
        integer(c_int), intent(out) :: global_vertex_id(*)
      end function imoab_getvertexid
 
      integer(c_int) function imoab_getvertexownership(pid, vertices_length, visible_global_rank_ID) &
                              bind(c, name='iMOAB_GetVertexOwnership')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: vertices_length
        integer(c_int), intent(out) :: visible_global_rank_id(*)
      end function imoab_getvertexownership
 
      integer(c_int) function imoab_getvisibleverticescoordinates(pid, coords_length, coordinates) &
                              bind(c, name='iMOAB_GetVisibleVerticesCoordinates')
        use, intrinsic :: iso_c_binding, only: c_int, c_double
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: coords_length
        real(c_double), intent(out) :: coordinates(*)
      end function imoab_getvisibleverticescoordinates
 
      integer(c_int) function imoab_getblockid(pid, block_length, global_block_IDs) bind(C, name='iMOAB_GetBlockID')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: block_length
        integer(c_int), intent(out) :: global_block_ids(*)
      end function imoab_getblockid
 
      integer(c_int) function imoab_getblockinfo(pid, global_block_ID, vertices_per_element, &
                                                num_elements_in_block) bind(C, name='iMOAB_GetBlockInfo')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: global_block_id
        integer(c_int), intent(out) :: vertices_per_element
        integer(c_int), intent(out) :: num_elements_in_block
      end function imoab_getblockinfo
 
      integer(c_int) function imoab_getvisibleelementsinfo(pid, num_visible_elements, element_global_IDs, &
                                                          ranks, block_IDs) bind(C, name='iMOAB_GetVisibleElementsInfo')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(out) :: num_visible_elements
        integer(c_int), intent(out) :: element_global_ids(*)
        integer(c_int), intent(out) :: ranks(*)
        integer(c_int), intent(out) :: block_ids(*)
      end function imoab_getvisibleelementsinfo
 
      integer(c_int) function imoab_getblockelementconnectivities(pid, global_block_ID, connectivity_length, &
                                                                  element_connectivity) &
                                                                  bind(c, name='iMOAB_GetBlockElementConnectivities')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: global_block_id
        integer(c_int), intent(in) :: connectivity_length
        integer(c_int), intent(out) :: element_connectivity(*)
      end function imoab_getblockelementconnectivities
 
      integer(c_int) function imoab_getelementconnectivity(pid, elem_index, connectivity_length, &
                                                          element_connectivity) bind(C, name='iMOAB_GetElementConnectivity')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: elem_index
        integer(c_int), intent(in) :: connectivity_length
        integer(c_int), intent(out) :: element_connectivity(*)
      end function imoab_getelementconnectivity
 
      integer(c_int) function imoab_getelementownership(pid, global_block_ID, num_elements_in_block, &
                                                        element_ownership) bind(C, name='iMOAB_GetElementOwnership')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: global_block_id
        integer(c_int), intent(out) :: num_elements_in_block
        integer(c_int), intent(out) :: element_ownership(*)
      end function imoab_getelementownership
 
      integer(c_int) function imoab_getelementid(pid, global_block_ID, num_elements_in_block, &
                                                global_element_ID, local_element_ID) bind(C, name='iMOAB_GetElementID')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: global_block_id
        integer(c_int), intent(in) :: num_elements_in_block
        integer(c_int), intent(in) :: global_element_id
        integer(c_int), intent(out) :: local_element_id
      end function imoab_getelementid
 
      integer(c_int) function imoab_getpointertosurfacebc(pid, surface_BC_length, local_element_ID, &
                                                          reference_surface_ID, boundary_condition_value) &
                                                          bind(c, name='iMOAB_GetPointerToSurfaceBC')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: surface_bc_length
        integer(c_int), intent(in) :: local_element_id
        integer(c_int), intent(in) :: reference_surface_id
        integer(c_int), intent(out) :: boundary_condition_value
      end function imoab_getpointertosurfacebc
 
      integer(c_int) function imoab_getpointertovertexbc(pid, vertex_BC_length, local_vertex_ID, &
                                                        boundary_condition_value) bind(C, name='iMOAB_GetPointerToVertexBC')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: vertex_bc_length
        integer(c_int), intent(in) :: local_vertex_id
        integer(c_int), intent(out) :: boundary_condition_value
      end function imoab_getpointertovertexbc
 
      integer(c_int) function imoab_definetagstorage(pid, tag_storage_name, tag_type, components_per_entity, &
                                                    tag_index) bind(C, name='iMOAB_DefineTagStorage')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer(c_int), intent(in) :: tag_type
        integer(c_int), intent(in) :: components_per_entity
        integer(c_int), intent(in) :: tag_index
      end function imoab_definetagstorage
 
      integer(c_int) function imoab_setinttagstorage(pid, tag_storage_name, num_tag_storage_length, entity_type, &
                                                    tag_storage_data) bind(C, name='iMOAB_SetIntTagStorage')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer(c_int), intent(in) :: num_tag_storage_length
        integer(c_int), intent(in) :: entity_type
        integer(c_int), intent(in) :: tag_storage_data(*)
      end function imoab_setinttagstorage
 
      integer(c_int) function imoab_getinttagstorage(pid, tag_storage_name, num_tag_storage_length, entity_type, &
                                                     tag_storage_data) bind(C, name='iMOAB_GetIntTagStorage')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer(c_int), intent(in) :: num_tag_storage_length
        integer(c_int), intent(in) :: entity_type
        integer(c_int), intent(out) :: tag_storage_data(*)
      end function imoab_getinttagstorage
 
      integer(c_int) function imoab_setdoubletagstorage(pid, tag_storage_name, num_tag_storage_length, entity_type, &
                                                        tag_storage_data) bind(C, name='iMOAB_SetDoubleTagStorage')
        use, intrinsic :: iso_c_binding, only: c_int, c_char, c_double
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer(c_int), intent(in) :: num_tag_storage_length
        integer(c_int), intent(in) :: entity_type
        real(c_double), intent(in) :: tag_storage_data(*)
      end function imoab_setdoubletagstorage
 
      integer(c_int) function imoab_setdoubletagstoragewithgid(pid, tag_storage_name, num_tag_storage_length, entity_type, &
                                         tag_storage_data, globalIds) bind(C, name='iMOAB_SetDoubleTagStorageWithGid')
        use, intrinsic :: iso_c_binding, only: c_int, c_char, c_double
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer(c_int), intent(in) :: num_tag_storage_length
        integer(c_int), intent(in) :: entity_type
        real(c_double), intent(in) :: tag_storage_data(*)
        integer(c_int), intent(in) :: globalids(*)
      end function imoab_setdoubletagstoragewithgid
 
 
      integer(c_int) function imoab_getdoubletagstorage(pid, tag_storage_name, num_tag_storage_length, entity_type, &
                                                        tag_storage_data) bind(C, name='iMOAB_GetDoubleTagStorage')
        use, intrinsic :: iso_c_binding, only: c_int, c_char, c_double
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer(c_int), intent(in) :: num_tag_storage_length
        integer(c_int), intent(in) :: entity_type
        real(c_double), intent(out) :: tag_storage_data(*)
      end function imoab_getdoubletagstorage
 
      integer(c_int) function imoab_synchronizetags(pid, num_tag, tag_indices, entity_type) bind(C, name='iMOAB_SynchronizeTags')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: num_tag
        integer(c_int), intent(in) :: tag_indices(*)
        integer(c_int), intent(in) :: entity_type(*)
      end function imoab_synchronizetags
 
      integer(c_int) function imoab_reducetagsmax(pid, tag_index, entity_type) bind(C, name='iMOAB_ReduceTagsMax')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: tag_index
        integer(c_int), intent(in) :: entity_type
      end function imoab_reducetagsmax
 
      integer(c_int) function imoab_getneighborelements(pid, local_index, num_adjacent_elements, adjacent_element_IDs) &
                              bind(c, name='iMOAB_GetNeighborElements')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: local_index
        integer(c_int), intent(out) :: num_adjacent_elements
        integer(c_int), intent(out) :: adjacent_element_ids
      end function imoab_getneighborelements
 
      integer(c_int) function imoab_getneighborvertices(pid, local_index, num_adjacent_vertices, adjacent_vertex_IDs) &
                              bind(c, name='iMOAB_GetNeighborVertices')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: local_index
        integer(c_int), intent(out) :: num_adjacent_vertices
        integer(c_int), intent(out) :: adjacent_vertex_ids
      end function imoab_getneighborvertices
 
      integer(c_int) function imoab_setglobalinfo(pid, num_global_verts, num_global_elems) bind(C, name='iMOAB_SetGlobalInfo')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: num_global_verts
        integer(c_int), intent(in) :: num_global_elems
      end function imoab_setglobalinfo
 
      integer(c_int) function imoab_getglobalinfo(pid, num_global_verts, num_global_elems) bind(C, name='iMOAB_GetGlobalInfo')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(out) :: num_global_verts
        integer(c_int), intent(out) :: num_global_elems
      end function imoab_getglobalinfo
 
  integer(c_int) function imoab_writelocalmesh(pid, prefix) bind(C, name='iMOAB_WriteLocalMesh')
    use, intrinsic :: iso_c_binding, only : c_int, c_char
    integer(c_int), intent(in) :: pid
    character(kind=c_char), intent(in) :: prefix(*)
  end function imoab_writelocalmesh
 
#ifdef MOAB_HAVE_MPI
 
      integer(c_int) function imoab_sendmesh(pid, joint_comm, receivingGroup, rcompid, method) bind(C, name='iMOAB_SendMesh')
        use, intrinsic :: iso_c_binding, only: c_int, c_ptr
        integer(c_int), intent(in) :: pid
        integer, intent(in) :: joint_comm ! MPI_Comm
        integer, intent(in) :: receivinggroup ! MPI_Group
        integer(c_int), intent(in) :: rcompid
        integer(c_int), intent(in) :: method
      end function imoab_sendmesh
 
      integer(c_int) function imoab_receivemesh(pid, joint_comm, sendingGroup, scompid) bind(C, name='iMOAB_ReceiveMesh')
        use, intrinsic :: iso_c_binding, only: c_int, c_ptr
        integer(c_int), intent(in) :: pid
        integer, intent(in) :: joint_comm ! MPI_Comm
        integer, intent(in) :: sendinggroup ! MPI_Group
        integer(c_int), intent(in) :: scompid
      end function imoab_receivemesh
 
      integer(c_int) function imoab_freesenderbuffers(pid, context_id) bind(C, name='iMOAB_FreeSenderBuffers')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: context_id
      end function imoab_freesenderbuffers
 
      integer(c_int) function imoab_sendelementtag(pid, tag_storage_name, joint_comm, context_id) &
                              bind(c, name='iMOAB_SendElementTag')
        use, intrinsic :: iso_c_binding, only: c_int, c_char, c_ptr
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer, intent(in) :: joint_comm ! MPI_Comm
        integer(c_int), intent(in) :: context_id
      end function imoab_sendelementtag
 
      integer(c_int) function imoab_receiveelementtag(pid, tag_storage_name, joint_comm, context_id) &
                              bind(c, name='iMOAB_ReceiveElementTag')
        use, intrinsic :: iso_c_binding, only: c_int, c_char, c_ptr
        integer(c_int), intent(in) :: pid
        character(kind=c_char), intent(in) :: tag_storage_name(*)
        integer, intent(in) :: joint_comm ! MPI_Comm
        integer(c_int), intent(in) :: context_id
      end function imoab_receiveelementtag
 
      integer(c_int) function imoab_computecommgraph(pid1, pid2, joint_comm, group1, group2, type1, type2, &
                                                     comp1, comp2) bind(C, name='iMOAB_ComputeCommGraph')
        use, intrinsic :: iso_c_binding, only: c_int, c_ptr
        integer(c_int), intent(in) :: pid1
        integer(c_int), intent(in) :: pid2
        integer, intent(in) :: joint_comm ! MPI_Comm
        integer, intent(in) :: group1     ! MPI_Group
        integer, intent(in) :: group2     ! MPI_Group
        integer(c_int), intent(in) :: type1
        integer(c_int), intent(in) :: type2
        integer(c_int), intent(in) :: comp1
        integer(c_int), intent(in) :: comp2
      end function imoab_computecommgraph
 
      integer(c_int) function imoab_coveragegraph(joint_comm, pid_source, pid_migration, pid_intx, source_id, &
                                                  migration_id, context_id) bind(C, name='iMOAB_CoverageGraph')
        use, intrinsic :: iso_c_binding, only: c_int, c_ptr
        integer, intent(in) :: joint_comm ! MPI_Comm
        integer(c_int), intent(in) :: pid_source
        integer(c_int), intent(in) :: pid_migration
        integer(c_int), intent(in) :: pid_intx
        integer(c_int), intent(in) :: source_id
        integer(c_int), intent(in) :: migration_id
        integer(c_int), intent(in) :: context_id
      end function imoab_coveragegraph
 
      integer(c_int) function imoab_dumpcommgraph(pid, context_id, is_sender, prefix) bind(C, name='iMOAB_DumpCommGraph')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: context_id
        integer(c_int), intent(in) :: is_sender
        character(kind=c_char), intent(in) :: prefix(*)
      end function imoab_dumpcommgraph
 
      integer(c_int) function imoab_mergevertices(pid) bind(C, name='iMOAB_MergeVertices')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
      end function imoab_mergevertices
 
      integer(c_int) function imoab_migratemapmesh( pid1, pid2, pid3, jointcomm, groupA, groupB, type, comp1, comp2, &
                                                    direction) bind(C, name='iMOAB_MigrateMapMesh')
        use, intrinsic :: iso_c_binding, only : c_int
        integer(c_int), intent(in) :: pid1
        integer(c_int), intent(in) :: pid2
        integer(c_int), intent(in) :: pid3
        integer, intent(in) :: jointcomm  ! MPI_Comm
        integer, intent(in) :: groupa     ! MPI_Group
        integer, intent(in) :: groupb     ! MPI_Group
        integer(c_int), intent(in) :: type
        integer(c_int), intent(in) :: comp1
        integer(c_int), intent(in) :: comp2
        integer(c_int), intent(in) :: direction
      end function imoab_migratemapmesh
 
      integer(c_int) function imoab_setghostlayers(pid, num_layers) bind(C, name='iMOAB_SetGhostLayers')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid
        integer(c_int), intent(in) :: num_layers
      end function imoab_setghostlayers
 
! closing endif: MOAB_HAVE_MPI
#endif
 
#ifdef MOAB_HAVE_TEMPESTREMAP
 
      integer(c_int) function imoab_computemeshintersectiononsphere(pid_source, pid_target, pid_intersection) &
                                                                  bind(c, name='iMOAB_ComputeMeshIntersectionOnSphere')
        use, intrinsic :: iso_c_binding, only: c_int, c_double
        integer(c_int), intent(in) :: pid_source
        integer(c_int), intent(in) :: pid_target
        integer(c_int), intent(in) :: pid_intersection
      end function imoab_computemeshintersectiononsphere
 
      integer(c_int) function imoab_computepointdofintersection(pid_source, pid_target, pid_intersection) &
                                                                bind(c, name='iMOAB_ComputePointDoFIntersection')
        use, intrinsic :: iso_c_binding, only: c_int
        integer(c_int), intent(in) :: pid_source
        integer(c_int), intent(in) :: pid_target
        integer(c_int), intent(in) :: pid_intersection
      end function imoab_computepointdofintersection
 
#ifdef MOAB_HAVE_NETCDF
 
      integer(c_int) function imoab_loadmappingweightsfromfile(pid_intersection, pid_cpl, column_or_row, mtype, &
                                                               solution_weights_identifier, remap_weights_filename) &
                                                                  bind(c, name='iMOAB_LoadMappingWeightsFromFile')
            use, intrinsic :: iso_c_binding, only : c_int, c_char
            integer(c_int), intent(in) :: pid_intersection
            integer(c_int), intent(in) :: pid_cpl
            integer(c_int), intent(in) :: column_or_row
            integer(c_int), intent(in) :: mtype
            character(kind=c_char), intent(in) :: solution_weights_identifier(*)
            character(kind=c_char), intent(in) :: remap_weights_filename(*)
      end function imoab_loadmappingweightsfromfile
 
      integer(c_int) function imoab_writemappingweightstofile(pid_intersection, solution_weights_identifier,  &
                                                  remap_weights_filename) bind(C, name='iMOAB_WriteMappingWeightsToFile')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid_intersection
        character(kind=c_char), intent(in) :: solution_weights_identifier(*)
        character(kind=c_char), intent(in) :: remap_weights_filename(*)
      end function imoab_writemappingweightstofile
 
! closing endif: MOAB_HAVE_NETCDF
#endif
 
      integer(c_int) function imoab_computescalarprojectionweights(pid_intersection, solution_weights_identifier, &
                                                                  disc_method_source, disc_order_source, &
                                                                  disc_method_target, disc_order_target, &
                                                                  fv_methods, fNoBubble, fMonotoneTypeID, fVolumetric, &
                                                                  fInverseDistanceMap, fNoConservation, fValidate, &
                                                                  source_solution_tag_dof_name, target_solution_tag_dof_name) &
                                                                  bind(c, name='iMOAB_ComputeScalarProjectionWeights')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid_intersection
        character(kind=c_char), intent(in) :: solution_weights_identifier(*)
        character(kind=c_char), intent(in) :: disc_method_source(*)
        integer(c_int), intent(in) :: disc_order_source
        character(kind=c_char), intent(in) :: disc_method_target(*)
        character(kind=c_char), intent(in) :: fv_methods(*)
        integer(c_int), intent(in) :: disc_order_target
        integer(c_int), intent(in) :: fnobubble
        integer(c_int), intent(in) :: fmonotonetypeid
        integer(c_int), intent(in) :: fvolumetric
        integer(c_int), intent(in) :: finversedistancemap
        integer(c_int), intent(in) :: fnoconservation
        integer(c_int), intent(in) :: fvalidate
        character(kind=c_char), intent(in) :: source_solution_tag_dof_name(*)
        character(kind=c_char), intent(in) :: target_solution_tag_dof_name(*)
      end function imoab_computescalarprojectionweights
 
      integer(c_int) function imoab_applyscalarprojectionweights(pid_intersection, filter_type, solution_weights_identifier, &
                                                                source_solution_tag_name, target_solution_tag_name) &
                                                                bind(c, name='iMOAB_ApplyScalarProjectionWeights')
        use, intrinsic :: iso_c_binding, only: c_int, c_char
        integer(c_int), intent(in) :: pid_intersection
        integer(c_int), intent(in) :: filter_type
        character(kind=c_char), intent(in) :: solution_weights_identifier(*)
        character(kind=c_char), intent(in) :: source_solution_tag_name(*)
        character(kind=c_char), intent(in) :: target_solution_tag_name(*)
      end function imoab_applyscalarprojectionweights
 
! closing endif: MOAB_HAVE_TEMPESTREMAP
#endif
 
#ifdef MOAB_HAVE_MPI
 
! Add some helper MPI converter functions
 
      ! Convert a C-based MPI-Comm handle to fortran
      integer function imoab_mpi_comm_c2f(c_handle) bind(C, name="MOAB_MPI_Comm_c2f")
        use, intrinsic :: iso_c_binding, only: c_ptr
        type(c_ptr), value :: c_handle
      end function
 
      ! Convert a Fortran-based MPI-integer handle to MPI-Comm pointer
      type(c_ptr) function imoab_mpi_comm_f2c(f_handle) bind(C, name="MOAB_MPI_Comm_f2c")
        use, intrinsic :: iso_c_binding, only: c_ptr
        integer, intent(in) :: f_handle
      end function
 
      ! Convert a C-based MPI-Group handle to fortran
      integer function imoab_mpi_group_c2f(c_handle) bind(C, name="MOAB_MPI_Group_c2f")
        use, intrinsic :: iso_c_binding, only: c_ptr
        type(c_ptr), value :: c_handle
      end function
 
      ! Convert a Fortran-based MPI-integer handle to MPI-Group pointer
      type(c_ptr) function imoab_mpi_group_f2c(f_handle) bind(C, name="MOAB_MPI_Group_f2c")
        use, intrinsic :: iso_c_binding, only: c_ptr
        integer, intent(in) :: f_handle
      end function
 
#endif
 
   end interface
 
end module imoab