methodsubcellmodule Module Reference

Data Types
type	methodsubcelltype
	Abstract base type for subcell tracking methods. More...

Functions/Subroutines
subroutine	assess (this, particle, cell_defn, tmax)
	Assess conditions before tracking. More...
subroutine	subcellexit (this, particle)
	Particle exits a subcell. More...
integer(i4b) function	get_level (this)
	Get the subcell method's level. More...

Function/Subroutine Documentation

assess()

subroutine methodsubcellmodule::assess ( class(methodsubcelltype), intent(inout)  this, type(particletype), intent(inout), pointer  particle, type(celldefntype), intent(inout), pointer  cell_defn, real(dp), intent(in)  tmax  )

private

Definition at line 25 of file MethodSubcell.f90.

    ! dummy
    class(MethodSubcellType), intent(inout) :: this
    type(ParticleType), pointer, intent(inout) :: particle
    type(CellDefnType), pointer, intent(inout) :: cell_defn
    real(DP), intent(in) :: tmax
    ! noop

get_level()

integer(i4b) function methodsubcellmodule::get_level ( class(methodsubcelltype), intent(in)  this )

private

Definition at line 72 of file MethodSubcell.f90.

    class(MethodSubcellType), intent(in) :: this
    integer(I4B) :: level
    level = level_subfeature

subcellexit()

subroutine methodsubcellmodule::subcellexit ( class(methodsubcelltype), intent(inout)  this, type(particletype), intent(inout), pointer  particle  )

private

Definition at line 35 of file MethodSubcell.f90.

    class(MethodSubcellType), intent(inout) :: this
    type(ParticleType), pointer, intent(inout) :: particle
    class(ParticleEventType), pointer :: event
    ! local
    integer(I4B) :: i, nhist
    class(*), pointer :: prev
 
    allocate (subcellexiteventtype :: event)
    select type (event)
    type is (subcellexiteventtype)
      event%isc = particle%itrdomain(level_subfeature)
      event%exit_face = particle%iboundary(level_subfeature)
    end select
    call this%events%broadcast(particle, event)
    if (particle%icycwin == 0) then
      deallocate (event)
      return
    end if
    if (this%forms_cycle(particle, event)) then
      print *, "Cyclic subcell pathline detected"
      nhist = particle%history%Count()
      do i = 1, nhist
        prev => particle%history%GetItem(i)
        select type (prev)
        class is (particleeventtype)
          print *, "Back ", nhist - i + 1, ": ", prev%get_text()
        end select
      end do
      print *, "Current :", event%get_text()
      call pstop(1, 'Cyclic subcell pathline detected, aborting')
    else
      call this%store_event(particle, event)
    end if

Here is the call graph for this function: