ExitSolution.f90

Go to the documentation of this file.

module exitsolutionmodule
 
  use kindmodule, only: i4b, dp, lgp
  use constantsmodule, only: dzero
 
  implicit none
 
  !> @brief Exit status codes
  enum, bind(C)
    enumerator :: ok_exit = 0 !< exit found using velocity interpolation
    ! below only used for linear solution
    enumerator :: ok_exit_constant = 1 !< exit found, constant velocity
    enumerator :: no_exit_stationary = 2 !< no exit, zero velocity
    enumerator :: no_exit_no_outflow = 3 !< no exit, no outflow
  end enum
 
  !> @brief Base type for exit solutions
  type :: exitsolutiontype
    integer(I4B) :: status = -1 !< domain exit status code
    integer(I4B) :: iboundary = 0 !< boundary number
    real(dp) :: dt = 1.0d+20 !< time to exit
  end type exitsolutiontype
 
  !> @brief Linear velocity interpolation exit solution
  type, extends(exitsolutiontype) :: linearexitsolutiontype
    real(dp) :: v = dzero !< particle velocity
    real(dp) :: dvdx = dzero !< velocity gradient
  end type linearexitsolutiontype
 
end module exitsolutionmodule