test_MOAB_GtoM Subroutine

subroutine test_MOAB_GtoM(correct, rc)
  logical :: correct
  integer :: rc

  type(ESMF_Grid) :: grid
  type(ESMF_Mesh) :: mesh
  type(ESMF_VM) :: vm
  integer :: localPet, petCount

  ! get global VM
  call ESMF_VMGetGlobal(vm, rc=rc)
  if (rc /= ESMF_SUCCESS) return
  call ESMF_VMGet(vm, localPet=localPet, petCount=petCount, rc=rc)
  if (rc /= ESMF_SUCCESS) return

  ! Init correct
  correct=.true.


  ! Create Grid
  grid=ESMF_GridCreateNoPeriDimUfrm( &
       maxIndex=(/4,4/), &
       minCornerCoord=(/0.0_ESMF_KIND_R8,0.0_ESMF_KIND_R8/), &
       maxCornerCoord=(/10.0_ESMF_KIND_R8,10.0_ESMF_KIND_R8/), &
       regDecomp=(/2,2/), & 
       staggerLocList=(/ESMF_STAGGERLOC_CENTER, ESMF_STAGGERLOC_CORNER/), &
       coordSys=ESMF_COORDSYS_CART, &
       rc=localrc)
  if (localrc /=ESMF_SUCCESS) then
    rc=ESMF_FAILURE
    return
  endif

  ! Turn on MOAB mesh creation
  call ESMF_MeshSetMOAB(.true., rc=localrc)
  if (localrc /=ESMF_SUCCESS) then
    rc=ESMF_FAILURE
    return
  endif

   ! Create Mesh structure in 1 step
   mesh=ESMF_MeshCreate(grid, rc=localrc)
   if (localrc /=ESMF_SUCCESS) then
      call ESMF_MeshSetMOAB(.false.)
      rc=ESMF_FAILURE
      return
   endif

  ! Turn off MOAB mesh creation
  call ESMF_MeshSetMOAB(.false., rc=localrc)
  if (localrc /=ESMF_SUCCESS) then
    rc=ESMF_FAILURE
    return
  endif

#if 0
   ! Output Mesh for debugging
   call ESMF_MeshWrite(mesh,"moab_gtom_mesh",rc=localrc)
   if (localrc /=ESMF_SUCCESS) then
      rc=ESMF_FAILURE
      return
   endif
#endif

   ! Get rid of Grid
   call ESMF_GridDestroy(grid, rc=localrc)
   if (localrc /=ESMF_SUCCESS) then
      rc=ESMF_FAILURE
      return
   endif

   ! Get rid of Mesh
   call ESMF_MeshDestroy(mesh, rc=localrc)
   if (localrc /=ESMF_SUCCESS) then
      rc=ESMF_FAILURE
      return
   endif


   ! Return success
   rc=ESMF_SUCCESS

end subroutine test_MOAB_GtoM