MethodCellPassToBot.f90

Go to the documentation of this file.

module methodcellpasstobotmodule
 
  use kindmodule, only: dp, i4b
  use methodmodule, only: methodtype
  use celldefnmodule, only: celldefntype, create_defn
  use prtfmimodule, only: prtfmitype
  use basedismodule, only: disbasetype
  use particlemodule, only: particletype
  use cellmodule, only: celltype
  use subcellmodule, only: subcelltype
  use trackcontrolmodule, only: trackcontroltype
  implicit none
 
  private
  public :: methodcellpasstobottype
  public :: create_method_cell_ptb
 
  type, extends(methodtype) :: methodcellpasstobottype
    private
  contains
    procedure, public :: apply => apply_ptb
    procedure, public :: deallocate
  end type methodcellpasstobottype
 
contains
 
  !> @brief Create a new pass-to-bottom tracking method
  subroutine create_method_cell_ptb(method)
    type(methodcellpasstobottype), pointer :: method
    allocate (method)
    allocate (method%name)
    method%name = "passtobottom"
    method%delegates = .false.
  end subroutine create_method_cell_ptb
 
  !> @brief Deallocate the pass-to-bottom tracking method
  subroutine deallocate (this)
    class(methodcellpasstobottype), intent(inout) :: this
    deallocate (this%name)
  end subroutine deallocate
 
  !> @brief Pass particle vertically and instantaneously to the cell bottom
  subroutine apply_ptb(this, particle, tmax)
    ! dummy
    class(methodcellpasstobottype), intent(inout) :: this
    type(particletype), pointer, intent(inout) :: particle
    real(DP), intent(in) :: tmax
 
    ! Check termination/reporting conditions
    call this%check(particle, this%cell%defn)
    if (.not. particle%advancing) return
 
    ! Pass to bottom face
    particle%z = this%cell%defn%bot
    particle%iboundary(2) = this%cell%defn%npolyverts + 2
 
    ! Save datum
    call this%save(particle, reason=1)
  end subroutine apply_ptb
 
end module methodcellpasstobotmodule