CellExitEvent.f90

Go to the documentation of this file.

module cellexiteventmodule
  use kindmodule, only: dp, i4b, lgp
  use constantsmodule, only: lenhugeline
  use errorutilmodule, only: pstop
  use particlemodule, only: particletype
  use particleeventmodule, only: featexit
  use featexiteventmodule, only: featexiteventtype
  implicit none
 
  private
  public :: cellexiteventtype
 
  type, extends(featexiteventtype) :: cellexiteventtype
    integer(I4B) :: exit_face ! face through which the particle exited
  contains
    procedure :: get_code
    procedure :: get_verb
    procedure :: get_text
  end type cellexiteventtype
 
contains
 
  function get_code(this) result(code)
    class(cellexiteventtype), intent(in) :: this
    integer(I4B) :: code
    code = featexit
  end function get_code
 
  function get_verb(this) result(verb)
    class(cellexiteventtype), intent(in) :: this
    character(len=:), allocatable :: verb
    verb = 'exited cell'
  end function get_verb
 
  function get_text(this) result(text)
    class(cellexiteventtype), intent(in) :: this
    character(len=:), allocatable :: text
    character(len=LENHUGELINE) :: temp
 
    write (temp, '(*(G0))') &
      'Particle from model ', this%imdl, &
      ', package ', this%iprp, &
      ', point ', this%irpt, &
      ', time ', this%trelease, &
      ' '//this%get_verb()// &
      ' in layer ', this%ilay, &
      ', cell ', this%icu, &
      ', zone ', this%izone, &
      ' through face ', this%exit_face, &
      ' at x ', this%x, &
      ', y ', this%y, &
      ', z ', this%z, &
      ', time ', this%ttrack, &
      ', period ', this%kper, &
      ', timestep ', this%kstp, &
      ' with status ', this%istatus
    text = trim(adjustl(temp))
  end function get_text
 
end module cellexiteventmodule