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