particlemodule Module Reference

Data Types
type	particletype
	Particle tracked by the PRT model. More...
type	particlestoretype
	Structure of arrays to store particles. More...

Functions/Subroutines
subroutine, public	create_particle (particle)
	Create a new particle. More...
subroutine, public	allocate_particle_store (this, np, mempath)
	Create a new particle store. More...
subroutine	deallocate (this, mempath)
	Deallocate particle arrays. More...
subroutine	resize (this, np, mempath)
	Reallocate particle arrays. More...
subroutine	load_particle (this, store, imdl, iprp, ip)
	Load a particle from the particle store. More...
subroutine	save_particle (this, particle, ip)
	Save a particle's state to the particle store. More...
subroutine	transform_coords (this, xorigin, yorigin, zorigin, sinrot, cosrot, invert, reset)
	Apply the given global-to-local transformation to the particle. More...
subroutine	get_model_coords (this, x, y, z)
	Return the particle's model (global) coordinates. More...
integer function	num_stored (this)

Variables
integer, parameter, public	levelmax = 4

Function/Subroutine Documentation

allocate_particle_store()

subroutine, public particlemodule::allocate_particle_store	(	type(particlestoretype), pointer	this,
		integer(i4b), intent(in)	np,
		character(*), intent(in)	mempath
	)

Parameters

	this	store
[in]	np	number of particles
[in]	mempath	path to memory

Definition at line 119 of file Particle.f90.

    type(ParticleStoreType), pointer :: this !< store
    integer(I4B), intent(in) :: np !< number of particles
    character(*), intent(in) :: mempath !< path to memory
 
    allocate (this)
    call mem_allocate(this%imdl, np, 'PLIMDL', mempath)
    call mem_allocate(this%irpt, np, 'PLIRPT', mempath)
    call mem_allocate(this%iprp, np, 'PLIPRP', mempath)
    call mem_allocate(this%name, lenboundname, np, 'PLNAME', mempath)
    call mem_allocate(this%icu, np, 'PLICU', mempath)
    call mem_allocate(this%ilay, np, 'PLILAY', mempath)
    call mem_allocate(this%izone, np, 'PLIZONE', mempath)
    call mem_allocate(this%istatus, np, 'PLISTATUS', mempath)
    call mem_allocate(this%x, np, 'PLX', mempath)
    call mem_allocate(this%y, np, 'PLY', mempath)
    call mem_allocate(this%z, np, 'PLZ', mempath)
    call mem_allocate(this%trelease, np, 'PLTRELEASE', mempath)
    call mem_allocate(this%tstop, np, 'PLTSTOP', mempath)
    call mem_allocate(this%ttrack, np, 'PLTTRACK', mempath)
    call mem_allocate(this%istopweaksink, np, 'PLISTOPWEAKSINK', mempath)
    call mem_allocate(this%istopzone, np, 'PLISTOPZONE', mempath)
    call mem_allocate(this%ifrctrn, np, 'PLIFRCTRN', mempath)
    call mem_allocate(this%iexmeth, np, 'PLIEXMETH', mempath)
    call mem_allocate(this%extol, np, 'PLEXTOL', mempath)
    call mem_allocate(this%extend, np, 'PLIEXTEND', mempath)
    call mem_allocate(this%idomain, np, levelmax, 'PLIDOMAIN', mempath)
    call mem_allocate(this%iboundary, np, levelmax, 'PLIBOUNDARY', mempath)

Here is the caller graph for this function:

create_particle()

subroutine, public particlemodule::create_particle

(

type(particletype), pointer

particle

)

Definition at line 111 of file Particle.f90.

    type(ParticleType), pointer :: particle !< particle
    allocate (particle)
    allocate (particle%idomain(levelmax))
    allocate (particle%iboundary(levelmax))

Here is the caller graph for this function:

deallocate()

subroutine particlemodule::deallocate ( class(particlestoretype), intent(inout)  this, character(*), intent(in)  mempath  )

private

Parameters

[in,out]	this	store
[in]	mempath	path to memory

Definition at line 150 of file Particle.f90.

    class(ParticleStoreType), intent(inout) :: this !< store
    character(*), intent(in) :: mempath !< path to memory
 
    call mem_deallocate(this%imdl, 'PLIMDL', mempath)
    call mem_deallocate(this%iprp, 'PLIPRP', mempath)
    call mem_deallocate(this%irpt, 'PLIRPT', mempath)
    call mem_deallocate(this%name, 'PLNAME', mempath)
    call mem_deallocate(this%icu, 'PLICU', mempath)
    call mem_deallocate(this%ilay, 'PLILAY', mempath)
    call mem_deallocate(this%izone, 'PLIZONE', mempath)
    call mem_deallocate(this%istatus, 'PLISTATUS', mempath)
    call mem_deallocate(this%x, 'PLX', mempath)
    call mem_deallocate(this%y, 'PLY', mempath)
    call mem_deallocate(this%z, 'PLZ', mempath)
    call mem_deallocate(this%trelease, 'PLTRELEASE', mempath)
    call mem_deallocate(this%tstop, 'PLTSTOP', mempath)
    call mem_deallocate(this%ttrack, 'PLTTRACK', mempath)
    call mem_deallocate(this%istopweaksink, 'PLISTOPWEAKSINK', mempath)
    call mem_deallocate(this%istopzone, 'PLISTOPZONE', mempath)
    call mem_deallocate(this%ifrctrn, 'PLIFRCTRN', mempath)
    call mem_deallocate(this%iexmeth, 'PLIEXMETH', mempath)
    call mem_deallocate(this%extol, 'PLEXTOL', mempath)
    call mem_deallocate(this%extend, 'PLIEXTEND', mempath)
    call mem_deallocate(this%idomain, 'PLIDOMAIN', mempath)
    call mem_deallocate(this%iboundary, 'PLIBOUNDARY', mempath)

get_model_coords()

subroutine particlemodule::get_model_coords	(	class(particletype), intent(inout)	this,
		real(dp), intent(out)	x,
		real(dp), intent(out)	y,
		real(dp), intent(out)	z
	)

Parameters

[in,out]	this	particle
[out]	x	x coordinate
[out]	y	y coordinate
[out]	z	z coordinate

Definition at line 337 of file Particle.f90.

    use geomutilmodule, only: transform, compose
    class(ParticleType), intent(inout) :: this !< particle
    real(DP), intent(out) :: x !< x coordinate
    real(DP), intent(out) :: y !< y coordinate
    real(DP), intent(out) :: z !< z coordinate
 
    if (this%transformed) then
      ! Transform back from local to model coordinates
      call transform(this%x, this%y, this%z, x, y, z, &
                     this%xorigin, this%yorigin, this%zorigin, &
                     this%sinrot, this%cosrot, .true.)
    else
      ! Already in model coordinates
      x = this%x
      y = this%y
      z = this%z
    end if

Here is the call graph for this function:

load_particle()

subroutine particlemodule::load_particle ( class(particletype), intent(inout)  this, type(particlestoretype), intent(in)  store, integer(i4b), intent(in)  imdl, integer(i4b), intent(in)  iprp, integer(i4b), intent(in)  ip  )

private

This routine is used to initialize a particle for tracking. The advancing flag and coordinate transformation are reset.

Parameters

[in,out]	this	particle
[in]	store	particle storage
[in]	imdl	index of model particle originated in
[in]	iprp	index of particle release package particle originated in
[in]	ip	index into the particle list

Definition at line 215 of file Particle.f90.

    class(ParticleType), intent(inout) :: this !< particle
    type(ParticleStoreType), intent(in) :: store !< particle storage
    integer(I4B), intent(in) :: imdl !< index of model particle originated in
    integer(I4B), intent(in) :: iprp !< index of particle release package particle originated in
    integer(I4B), intent(in) :: ip !< index into the particle list
 
    call this%transform(reset=.true.)
    this%imdl = imdl
    this%iprp = iprp
    this%irpt = store%irpt(ip)
    this%ip = ip
    this%name = store%name(ip)
    this%istopweaksink = store%istopweaksink(ip)
    this%istopzone = store%istopzone(ip)
    this%icu = store%icu(ip)
    this%ilay = store%ilay(ip)
    this%izone = store%izone(ip)
    this%istatus = store%istatus(ip)
    this%x = store%x(ip)
    this%y = store%y(ip)
    this%z = store%z(ip)
    this%trelease = store%trelease(ip)
    this%tstop = store%tstop(ip)
    this%ttrack = store%ttrack(ip)
    this%advancing = .true.
    this%idomain(1:levelmax) = &
      store%idomain(ip, 1:levelmax)
    this%idomain(1) = imdl
    this%iboundary(1:levelmax) = &
      store%iboundary(ip, 1:levelmax)
    this%ifrctrn = store%ifrctrn(ip)
    this%iexmeth = store%iexmeth(ip)
    this%extol = store%extol(ip)
    this%iextend = store%extend(ip)

num_stored()

integer function particlemodule::num_stored

(

class(particlestoretype)

this

)

Definition at line 357 of file Particle.f90.

    class(ParticleStoreType) :: this
    n = size(this%imdl)

resize()

subroutine particlemodule::resize ( class(particlestoretype), intent(inout)  this, integer(i4b), intent(in)  np, character(*), intent(in)  mempath  )

private

Parameters

[in,out]	this	particle store
[in]	np	number of particles
[in]	mempath	path to memory

Definition at line 179 of file Particle.f90.

    ! -- dummy
    class(ParticleStoreType), intent(inout) :: this !< particle store
    integer(I4B), intent(in) :: np !< number of particles
    character(*), intent(in) :: mempath !< path to memory
 
    ! resize arrays
    call mem_reallocate(this%imdl, np, 'PLIMDL', mempath)
    call mem_reallocate(this%iprp, np, 'PLIPRP', mempath)
    call mem_reallocate(this%irpt, np, 'PLIRPT', mempath)
    call mem_reallocate(this%name, lenboundname, np, 'PLNAME', mempath)
    call mem_reallocate(this%icu, np, 'PLICU', mempath)
    call mem_reallocate(this%ilay, np, 'PLILAY', mempath)
    call mem_reallocate(this%izone, np, 'PLIZONE', mempath)
    call mem_reallocate(this%istatus, np, 'PLISTATUS', mempath)
    call mem_reallocate(this%x, np, 'PLX', mempath)
    call mem_reallocate(this%y, np, 'PLY', mempath)
    call mem_reallocate(this%z, np, 'PLZ', mempath)
    call mem_reallocate(this%trelease, np, 'PLTRELEASE', mempath)
    call mem_reallocate(this%tstop, np, 'PLTSTOP', mempath)
    call mem_reallocate(this%ttrack, np, 'PLTTRACK', mempath)
    call mem_reallocate(this%istopweaksink, np, 'PLISTOPWEAKSINK', mempath)
    call mem_reallocate(this%istopzone, np, 'PLISTOPZONE', mempath)
    call mem_reallocate(this%ifrctrn, np, 'PLIFRCTRN', mempath)
    call mem_reallocate(this%iexmeth, np, 'PLIEXMETH', mempath)
    call mem_reallocate(this%extol, np, 'PLEXTOL', mempath)
    call mem_reallocate(this%extend, np, 'PLIEXTEND', mempath)
    call mem_reallocate(this%idomain, np, levelmax, 'PLIDOMAIN', mempath)
    call mem_reallocate(this%iboundary, np, levelmax, 'PLIBOUNDARY', mempath)

save_particle()

subroutine particlemodule::save_particle ( class(particlestoretype), intent(inout)  this, type(particletype), intent(in)  particle, integer(i4b), intent(in)  ip  )

private

Parameters

[in,out]	this	particle storage
[in]	ip	particle index

Definition at line 253 of file Particle.f90.

    class(ParticleStoreType), intent(inout) :: this !< particle storage
    type(ParticleType), intent(in) :: particle !< particle
    integer(I4B), intent(in) :: ip !< particle index
 
    this%imdl(ip) = particle%imdl
    this%iprp(ip) = particle%iprp
    this%irpt(ip) = particle%irpt
    this%name(ip) = particle%name
    this%istopweaksink(ip) = particle%istopweaksink
    this%istopzone(ip) = particle%istopzone
    this%icu(ip) = particle%icu
    this%ilay(ip) = particle%ilay
    this%izone(ip) = particle%izone
    this%istatus(ip) = particle%istatus
    this%x(ip) = particle%x
    this%y(ip) = particle%y
    this%z(ip) = particle%z
    this%trelease(ip) = particle%trelease
    this%tstop(ip) = particle%tstop
    this%ttrack(ip) = particle%ttrack
    this%idomain( &
      ip, &
      1:levelmax) = &
      particle%idomain(1:levelmax)
    this%iboundary( &
      ip, &
      1:levelmax) = &
      particle%iboundary(1:levelmax)
    this%ifrctrn(ip) = particle%ifrctrn
    this%iexmeth(ip) = particle%iexmeth
    this%extol(ip) = particle%extol
    this%extend(ip) = particle%iextend

transform_coords()

subroutine particlemodule::transform_coords ( class(particletype), intent(inout)  this, real(dp), intent(in), optional  xorigin, real(dp), intent(in), optional  yorigin, real(dp), intent(in), optional  zorigin, real(dp), intent(in), optional  sinrot, real(dp), intent(in), optional  cosrot, logical(lgp), intent(in), optional  invert, logical(lgp), intent(in), optional  reset  )

private

Parameters

[in,out]	this	particle
[in]	xorigin	x coordinate of origin
[in]	yorigin	y coordinate of origin
[in]	zorigin	z coordinate of origin
[in]	sinrot	sine of rotation angle
[in]	cosrot	cosine of rotation angle
[in]	invert	whether to invert
[in]	reset	whether to reset

Definition at line 289 of file Particle.f90.

    use geomutilmodule, only: transform, compose
    class(ParticleType), intent(inout) :: this !< particle
    real(DP), intent(in), optional :: xorigin !< x coordinate of origin
    real(DP), intent(in), optional :: yorigin !< y coordinate of origin
    real(DP), intent(in), optional :: zorigin !< z coordinate of origin
    real(DP), intent(in), optional :: sinrot !< sine of rotation angle
    real(DP), intent(in), optional :: cosrot !< cosine of rotation angle
    logical(LGP), intent(in), optional :: invert !< whether to invert
    logical(LGP), intent(in), optional :: reset !< whether to reset
 
    ! Reset if requested
    if (present(reset)) then
      if (reset) then
        this%xorigin = dzero
        this%yorigin = dzero
        this%zorigin = dzero
        this%sinrot = dzero
        this%cosrot = done
        this%cosrot = done
        this%transformed = .false.
        return
      end if
    end if
 
    ! Otherwise, transform coordinates
    call transform(this%x, this%y, this%z, &
                   this%x, this%y, this%z, &
                   xorigin, yorigin, zorigin, &
                   sinrot, cosrot, invert)
 
    ! Modify transformation from model coordinates to particle's new
    ! local coordinates by incorporating this latest transformation
    call compose(this%xorigin, this%yorigin, this%zorigin, &
                 this%sinrot, this%cosrot, &
                 xorigin, yorigin, zorigin, &
                 sinrot, cosrot, invert)
 
    ! Set isTransformed flag to true. Note that there is no check
    ! to see whether the modification brings the coordinates back
    ! to model coordinates (in which case the origin would be very
    ! close to zero and sinrot and cosrot would be very close to 0.
    ! and 1., respectively, allowing for roundoff error).
    this%transformed = .true.

Here is the call graph for this function:

Variable Documentation

levelmax

integer, parameter, public particlemodule::levelmax = 4

Definition at line 14 of file Particle.f90.

  integer, parameter, public :: levelmax = 4