src/service/itm_grid.f90

Go to the documentation of this file.

module itm_grid

  !use itm_types , ITM_R8 => R8, ITM_R4 => R4
  use itm_assert
  use euitm_schemas ! IGNORE
  use itm_string

  use itm_grid_common 
  use itm_grid_access
  use itm_grid_structured
  use itm_grid_object
  use itm_grid_objectlist
  use itm_grid_subgrid

  implicit none      

contains


  !> Check a grid description for possible errors

  subroutine gridSanityCheck( grid, doStop )
    type(type_complexgrid), intent(in) :: grid
    logical, intent(in), optional :: doStop
    
    ! integer
    integer :: is

    ! disable immediate stopping in assertion module
    call assertSetStopMode( doStop = .false. )   

    do is = 1, size( grid % spaces )
       call assertSetMsgPrefix( "Space "// itmInt2str( is ) )
       call gridSanityCheckSpace( grid % spaces( is ) ) 
    end do

    ! Output error message if any errors found, and stop (if requested)
    call assertStopOnFailed( "Grid contained errors", doStop )


    ! reset assertion module
    call assertReset()

    ! re-enable immediate stop for assertions
    call assertSetStopMode()   
    
    ! clear message prefix
    call assertSetMsgPrefix()

  end subroutine gridSanityCheck
  

  !> Check a grid space description for possible errors

  subroutine gridSanityCheckSpace( space )
    type(type_complexgrid_space), intent(in) :: space
    
    call assert( ( size( space % coordtype ) <= size( space % objects ) ), &
         & "coordtype, objects,have inconsistent size" )
        

  end subroutine gridSanityCheckSpace


end module itm_grid