subroutine ESMF_MeshFindPnt(mesh, unmappedaction, &
pntDim, pntCount, pntList, &
petList, rc)
!
! !ARGUMENTS:
type(ESMF_Mesh), intent(in) :: mesh
type(ESMF_UnmappedAction_Flag), intent(in), optional :: unmappedaction
integer, intent(in) :: pntDim
integer, intent(in) :: pntCount
real(ESMF_KIND_R8), pointer :: pntList(:)
integer, pointer :: petList(:)
integer, intent(out),optional :: rc
!
! !DESCRIPTION:
! Write a mesh to VTK file.
!
! \begin{description}
! \item [mesh]
! The mesh.
! \item [{[unmappedaction]}]
! Specifies what should happen if there are destination points that
! can't be mapped to a source cell. Options are
! {\tt ESMF\_UNMAPPEDACTION\_ERROR} or
! {\tt ESMF\_UNMAPPEDACTION\_IGNORE}. If not specified, defaults
! to {\tt ESMF\_UNMAPPEDACTION\_ERROR}.
! \item [pntDim]
! The dimension of the points in pntList.
! \item [pntNum]
! The number of points in pntList
! \item [petList]
! The generated list of pets for each point.
! \item [{[rc]}]
! Return code; equals {\tt ESMF\_SUCCESS} if there are no errors.
! \end{description}
!
!EOPI
!------------------------------------------------------------------------------
integer :: localrc
type(ESMF_Logical) :: isfree
type(ESMF_MeshStatus_Flag) :: status
type(ESMF_UnmappedAction_Flag) :: localunmappedaction
! initialize return code; assume routine not implemented
localrc = ESMF_RC_NOT_IMPL
if (present(rc)) rc = ESMF_RC_NOT_IMPL
! Check init status of arguments
ESMF_INIT_CHECK_DEEP(ESMF_MeshGetInit, mesh, rc)
call C_ESMC_MeshGetIsFree(mesh, isfree)
if (isfree == ESMF_TRUE) then
call ESMF_LogSetError(rcToCheck=ESMF_RC_OBJ_WRONG, &
msg="- the mesh internals have been freed", &
ESMF_CONTEXT, rcToReturn=rc)
return
endif
call C_ESMC_MeshGetStatus(mesh, status)
if (status .ne. ESMF_MESHSTATUS_COMPLETE) then
call ESMF_LogSetError(rcToCheck=ESMF_RC_OBJ_WRONG, &
msg="- the mesh has not been fully created", &
ESMF_CONTEXT, rcToReturn=rc)
return
endif
! Set default vale for unmappedaction
if (present(unmappedaction)) then
localunmappedaction=unmappedaction
else
localunmappedaction=ESMF_UNMAPPEDACTION_ERROR
endif
! Call into mesh find point subroutine
! TODO: ADD GIDS TO THIS INTERFACE AS THEY'RE AVAILABLE FROM C++ METHOD
call C_ESMC_MeshFindPnt(mesh%this, localunmappedaction, pntDim, pntCount, &
pntList, petList, localrc)
if (ESMF_LogFoundError(localrc, ESMF_ERR_PASSTHRU, &
ESMF_CONTEXT, rcToReturn=rc)) return
! return success
if (present(rc)) rc = ESMF_SUCCESS
end subroutine ESMF_MeshFindPnt