MODFLOW 6  version 6.7.0.dev2
USGS Modular Hydrologic Model
Data Types | Enumerations | Functions/Subroutines | Variables
particlemodule Module Reference

Data Types

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

Enumerations

enum  {
  active = 1 , term_boundary = 2 , term_weaksink = 3 , term_no_exits = 5 ,
  term_stopzone = 6 , term_inactive = 7 , term_unreleased = 8 , term_no_exits_sub = 9 ,
  term_timeout = 10
}
 Particle status enumeration. More...
 

Functions/Subroutines

subroutine create_particle (particle)
 Create a new particle. More...
 
subroutine create_particle_store (store, np, mempath)
 Allocate particle store. More...
 
subroutine destroy (this, mempath)
 Destroy particle store after use. More...
 
subroutine destroy_particle (particle)
 Destroy a particle after use. More...
 
subroutine resize (this, np, mempath)
 Reallocate particle storage to the given size. More...
 
subroutine get (this, particle, imdl, iprp, ip)
 Load a particle from the particle store. More...
 
subroutine put (this, particle, ip)
 Save a particle's state to the particle store. More...
 
subroutine transform_coords (this, xorigin, yorigin, zorigin, sinrot, cosrot, invert)
 Transform particle coordinates. More...
 
subroutine reset_transform (this)
 Reset particle coordinate transformation properties. More...
 
subroutine get_model_coords (this, x, y, z)
 Return the particle's model coordinates, inverting any applied transformation if needed. The particle's state is not altered. More...
 
integer function num_stored (this)
 Return the number of particles. More...
 

Variables

integer, parameter max_level = 4
 Tracking "levels" defined in method modules. Currently only 3 used. More...
 

Enumeration Type Documentation

◆ anonymous enum

anonymous enum

Particles begin in status 1 (active) at release time. Status may only increase over time. Status values greater than one imply termination. A particle may terminate for several reasons, all mutually exclusive. A particle's final tracking status will always be greater than one.

Status codes 0-3 and 5-8 correspond directly to MODPATH 7 status codes. Code 4 does not apply to PRT because PRT does not distinguish forwards from backwards tracking. Status code 9 provides more specific, subcell- level information about a particle which terminates due to no outflow. Code 10 distinguishes particles which have "timed out" upon reaching a user-specified stop time or the end of the simulation.

Enumerator
active 
term_boundary 

terminated at a boundary face

term_weaksink 

terminated in a weak sink cell

term_no_exits 

terminated in a cell with no exit face

term_stopzone 

terminated in a cell with a stop zone number

term_inactive 

terminated in an inactive cell

term_unreleased 

terminated permanently unreleased

term_no_exits_sub 

terminated in a subcell with no exit face

term_timeout 

terminated at stop time or end of simulation

Definition at line 28 of file Particle.f90.

Function/Subroutine Documentation

◆ create_particle()

subroutine particlemodule::create_particle ( type(particletype), pointer  particle)

Definition at line 145 of file Particle.f90.

146  type(ParticleType), pointer :: particle !< particle
147  allocate (particle)
148  allocate (particle%history)

◆ create_particle_store()

subroutine particlemodule::create_particle_store ( type(particlestoretype), pointer  store,
integer(i4b), intent(in)  np,
character(*), intent(in)  mempath 
)
Parameters
[in]npnumber of particles
[in]mempathpath to memory

Definition at line 152 of file Particle.f90.

153  type(ParticleStoreType), pointer :: store !< store
154  integer(I4B), intent(in) :: np !< number of particles
155  character(*), intent(in) :: mempath !< path to memory
156 
157  allocate (store)
158  call mem_allocate(store%imdl, np, 'PLIMDL', mempath)
159  call mem_allocate(store%irpt, np, 'PLIRPT', mempath)
160  call mem_allocate(store%iprp, np, 'PLIPRP', mempath)
161  call mem_allocate(store%name, lenboundname, np, 'PLNAME', mempath)
162  call mem_allocate(store%icu, np, 'PLICU', mempath)
163  call mem_allocate(store%ilay, np, 'PLILAY', mempath)
164  call mem_allocate(store%izone, np, 'PLIZONE', mempath)
165  call mem_allocate(store%izp, np, 'PLIZP', mempath)
166  call mem_allocate(store%istatus, np, 'PLISTATUS', mempath)
167  call mem_allocate(store%x, np, 'PLX', mempath)
168  call mem_allocate(store%y, np, 'PLY', mempath)
169  call mem_allocate(store%z, np, 'PLZ', mempath)
170  call mem_allocate(store%trelease, np, 'PLTRELEASE', mempath)
171  call mem_allocate(store%tstop, np, 'PLTSTOP', mempath)
172  call mem_allocate(store%ttrack, np, 'PLTTRACK', mempath)
173  call mem_allocate(store%istopweaksink, np, 'PLISTOPWEAKSINK', mempath)
174  call mem_allocate(store%istopzone, np, 'PLISTOPZONE', mempath)
175  call mem_allocate(store%idrymeth, np, 'PLIDRYMETH', mempath)
176  call mem_allocate(store%ifrctrn, np, 'PLIFRCTRN', mempath)
177  call mem_allocate(store%iexmeth, np, 'PLIEXMETH', mempath)
178  call mem_allocate(store%extol, np, 'PLEXTOL', mempath)
179  call mem_allocate(store%extend, np, 'PLIEXTEND', mempath)
180  call mem_allocate(store%icycwin, np, 'PLICYCWIN', mempath)
181  call mem_allocate(store%itrdomain, np, max_level, 'PLIDOMAIN', mempath)
182  call mem_allocate(store%iboundary, np, max_level, 'PLIBOUNDARY', mempath)
Here is the caller graph for this function:

◆ destroy()

subroutine particlemodule::destroy ( class(particlestoretype), intent(inout)  this,
character(*), intent(in)  mempath 
)
Parameters
[in,out]thisstore
[in]mempathpath to memory

Definition at line 186 of file Particle.f90.

187  class(ParticleStoreType), intent(inout) :: this !< store
188  character(*), intent(in) :: mempath !< path to memory
189 
190  call mem_deallocate(this%imdl, 'PLIMDL', mempath)
191  call mem_deallocate(this%iprp, 'PLIPRP', mempath)
192  call mem_deallocate(this%irpt, 'PLIRPT', mempath)
193  call mem_deallocate(this%name, 'PLNAME', mempath)
194  call mem_deallocate(this%icu, 'PLICU', mempath)
195  call mem_deallocate(this%ilay, 'PLILAY', mempath)
196  call mem_deallocate(this%izone, 'PLIZONE', mempath)
197  call mem_deallocate(this%izp, 'PLIZP', mempath)
198  call mem_deallocate(this%istatus, 'PLISTATUS', mempath)
199  call mem_deallocate(this%x, 'PLX', mempath)
200  call mem_deallocate(this%y, 'PLY', mempath)
201  call mem_deallocate(this%z, 'PLZ', mempath)
202  call mem_deallocate(this%trelease, 'PLTRELEASE', mempath)
203  call mem_deallocate(this%tstop, 'PLTSTOP', mempath)
204  call mem_deallocate(this%ttrack, 'PLTTRACK', mempath)
205  call mem_deallocate(this%istopweaksink, 'PLISTOPWEAKSINK', mempath)
206  call mem_deallocate(this%istopzone, 'PLISTOPZONE', mempath)
207  call mem_deallocate(this%idrymeth, 'PLIDRYMETH', mempath)
208  call mem_deallocate(this%ifrctrn, 'PLIFRCTRN', mempath)
209  call mem_deallocate(this%iexmeth, 'PLIEXMETH', mempath)
210  call mem_deallocate(this%extol, 'PLEXTOL', mempath)
211  call mem_deallocate(this%extend, 'PLIEXTEND', mempath)
212  call mem_deallocate(this%icycwin, 'PLICYCWIN', mempath)
213  call mem_deallocate(this%itrdomain, 'PLIDOMAIN', mempath)
214  call mem_deallocate(this%iboundary, 'PLIBOUNDARY', mempath)

◆ destroy_particle()

subroutine particlemodule::destroy_particle ( class(particletype), intent(inout)  particle)

Definition at line 218 of file Particle.f90.

219  class(ParticleType), intent(inout) :: particle !< particle
220  deallocate (particle%history)

◆ get()

subroutine particlemodule::get ( class(particlestoretype), intent(inout)  this,
class(particletype), intent(inout)  particle,
integer(i4b), intent(in)  imdl,
integer(i4b), intent(in)  iprp,
integer(i4b), intent(in)  ip 
)

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

Parameters
[in,out]thisparticle store
[in]imdlindex of model particle originated in
[in]iprpindex of particle release package particle originated in
[in]ipindex into the particle list

Definition at line 263 of file Particle.f90.

264  class(ParticleStoreType), intent(inout) :: this !< particle store
265  class(ParticleType), intent(inout) :: particle !< particle
266  integer(I4B), intent(in) :: imdl !< index of model particle originated in
267  integer(I4B), intent(in) :: iprp !< index of particle release package particle originated in
268  integer(I4B), intent(in) :: ip !< index into the particle list
269 
270  call particle%reset_transform()
271  call particle%history%Clear()
272  particle%imdl = imdl
273  particle%iprp = iprp
274  particle%irpt = this%irpt(ip)
275  particle%ip = ip
276  particle%name = this%name(ip)
277  particle%istopweaksink = this%istopweaksink(ip)
278  particle%istopzone = this%istopzone(ip)
279  particle%idrymeth = this%idrymeth(ip)
280  particle%icp = 0
281  particle%icu = this%icu(ip)
282  particle%ilay = this%ilay(ip)
283  particle%izone = this%izone(ip)
284  particle%izp = this%izp(ip)
285  particle%istatus = this%istatus(ip)
286  particle%x = this%x(ip)
287  particle%y = this%y(ip)
288  particle%z = this%z(ip)
289  particle%trelease = this%trelease(ip)
290  particle%tstop = this%tstop(ip)
291  particle%ttrack = this%ttrack(ip)
292  particle%advancing = .true.
293  particle%itrdomain(1:max_level) = &
294  this%itrdomain(ip, 1:max_level)
295  particle%itrdomain(1) = imdl
296  particle%iboundary(1:max_level) = &
297  this%iboundary(ip, 1:max_level)
298  particle%ifrctrn = this%ifrctrn(ip)
299  particle%iexmeth = this%iexmeth(ip)
300  particle%extol = this%extol(ip)
301  particle%iextend = this%extend(ip)
302  particle%icycwin = this%icycwin(ip)

◆ get_model_coords()

subroutine particlemodule::get_model_coords ( class(particletype), intent(in)  this,
real(dp), intent(out)  x,
real(dp), intent(out)  y,
real(dp), intent(out)  z 
)
Parameters
[in]thisparticle
[out]xx coordinate
[out]yy coordinate
[out]zz coordinate

Definition at line 390 of file Particle.f90.

391  use geomutilmodule, only: transform, compose
392  class(ParticleType), intent(in) :: this !< particle
393  real(DP), intent(out) :: x !< x coordinate
394  real(DP), intent(out) :: y !< y coordinate
395  real(DP), intent(out) :: z !< z coordinate
396 
397  if (this%transformed) then
398  call transform(this%x, this%y, this%z, x, y, z, &
399  this%xorigin, this%yorigin, this%zorigin, &
400  this%sinrot, this%cosrot, invert=.true.)
401  else
402  x = this%x
403  y = this%y
404  z = this%z
405  end if
geomutilmodule::transform
subroutine, public transform(xin, yin, zin, xout, yout, zout, xorigin, yorigin, zorigin, sinrot, cosrot, invert)
Apply a 3D translation and optional 2D rotation to coordinates.
Definition: GeomUtil.f90:183
geomutilmodule::compose
subroutine, public compose(xorigin, yorigin, zorigin, sinrot, cosrot, xorigin_new, yorigin_new, zorigin_new, sinrot_new, cosrot_new, invert)
Apply a 3D translation and 2D rotation to an existing transformation.
Definition: GeomUtil.f90:243
Here is the call graph for this function:

◆ num_stored()

integer function particlemodule::num_stored ( class(particlestoretype)  this)

Definition at line 409 of file Particle.f90.

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

◆ put()

subroutine particlemodule::put ( class(particlestoretype), intent(inout)  this,
class(particletype), intent(in)  particle,
integer(i4b), intent(in)  ip 
)
Parameters
[in,out]thisparticle storage
[in]ipparticle index

Definition at line 306 of file Particle.f90.

307  class(ParticleStoreType), intent(inout) :: this !< particle storage
308  class(ParticleType), intent(in) :: particle !< particle
309  integer(I4B), intent(in) :: ip !< particle index
310 
311  this%imdl(ip) = particle%imdl
312  this%iprp(ip) = particle%iprp
313  this%irpt(ip) = particle%irpt
314  this%name(ip) = particle%name
315  this%istopweaksink(ip) = particle%istopweaksink
316  this%istopzone(ip) = particle%istopzone
317  this%idrymeth(ip) = particle%idrymeth
318  this%icu(ip) = particle%icu
319  this%ilay(ip) = particle%ilay
320  this%izone(ip) = particle%izone
321  this%izp(ip) = particle%izp
322  this%istatus(ip) = particle%istatus
323  this%x(ip) = particle%x
324  this%y(ip) = particle%y
325  this%z(ip) = particle%z
326  this%trelease(ip) = particle%trelease
327  this%tstop(ip) = particle%tstop
328  this%ttrack(ip) = particle%ttrack
329  this%itrdomain( &
330  ip, &
331  1:max_level) = &
332  particle%itrdomain(1:max_level)
333  this%iboundary( &
334  ip, &
335  1:max_level) = &
336  particle%iboundary(1:max_level)
337  this%ifrctrn(ip) = particle%ifrctrn
338  this%iexmeth(ip) = particle%iexmeth
339  this%extol(ip) = particle%extol
340  this%extend(ip) = particle%iextend
341  this%icycwin(ip) = particle%icycwin

◆ reset_transform()

subroutine particlemodule::reset_transform ( class(particletype), intent(inout)  this)
Parameters
[in,out]thisparticle

Definition at line 375 of file Particle.f90.

376  class(ParticleType), intent(inout) :: this !< particle
377 
378  this%xorigin = dzero
379  this%yorigin = dzero
380  this%zorigin = dzero
381  this%sinrot = dzero
382  this%cosrot = done
383  this%cosrot = done
384  this%transformed = .false.

◆ resize()

subroutine particlemodule::resize ( class(particlestoretype), intent(inout)  this,
integer(i4b), intent(in)  np,
character(*), intent(in)  mempath 
)
Parameters
[in,out]thisparticle store
[in]npnumber of particles
[in]mempathpath to memory

Definition at line 224 of file Particle.f90.

225  ! dummy
226  class(ParticleStoreType), intent(inout) :: this !< particle store
227  integer(I4B), intent(in) :: np !< number of particles
228  character(*), intent(in) :: mempath !< path to memory
229 
230  ! resize arrays
231  call mem_reallocate(this%imdl, np, 'PLIMDL', mempath)
232  call mem_reallocate(this%iprp, np, 'PLIPRP', mempath)
233  call mem_reallocate(this%irpt, np, 'PLIRPT', mempath)
234  call mem_reallocate(this%name, lenboundname, np, 'PLNAME', mempath)
235  call mem_reallocate(this%icu, np, 'PLICU', mempath)
236  call mem_reallocate(this%ilay, np, 'PLILAY', mempath)
237  call mem_reallocate(this%izone, np, 'PLIZONE', mempath)
238  call mem_reallocate(this%izp, np, 'PLIZP', mempath)
239  call mem_reallocate(this%istatus, np, 'PLISTATUS', mempath)
240  call mem_reallocate(this%x, np, 'PLX', mempath)
241  call mem_reallocate(this%y, np, 'PLY', mempath)
242  call mem_reallocate(this%z, np, 'PLZ', mempath)
243  call mem_reallocate(this%trelease, np, 'PLTRELEASE', mempath)
244  call mem_reallocate(this%tstop, np, 'PLTSTOP', mempath)
245  call mem_reallocate(this%ttrack, np, 'PLTTRACK', mempath)
246  call mem_reallocate(this%istopweaksink, np, 'PLISTOPWEAKSINK', mempath)
247  call mem_reallocate(this%istopzone, np, 'PLISTOPZONE', mempath)
248  call mem_reallocate(this%idrymeth, np, 'PLIDRYMETH', mempath)
249  call mem_reallocate(this%ifrctrn, np, 'PLIFRCTRN', mempath)
250  call mem_reallocate(this%iexmeth, np, 'PLIEXMETH', mempath)
251  call mem_reallocate(this%extol, np, 'PLEXTOL', mempath)
252  call mem_reallocate(this%extend, np, 'PLIEXTEND', mempath)
253  call mem_reallocate(this%icycwin, np, 'PLICYCWIN', mempath)
254  call mem_reallocate(this%itrdomain, np, max_level, 'PLIDOMAIN', mempath)
255  call mem_reallocate(this%iboundary, np, max_level, 'PLIBOUNDARY', mempath)

◆ 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 
)

Apply a translation and/or rotation to particle coordinates. No rescaling. It's also possible to invert a transformation. Be sure to reset the transformation after using it.

Parameters
[in,out]thisparticle
[in]xoriginx coordinate of origin
[in]yoriginy coordinate of origin
[in]zoriginz coordinate of origin
[in]sinrotsine of rotation angle
[in]cosrotcosine of rotation angle
[in]invertwhether to invert

Definition at line 350 of file Particle.f90.

352  use geomutilmodule, only: transform, compose
353  class(ParticleType), intent(inout) :: this !< particle
354  real(DP), intent(in), optional :: xorigin !< x coordinate of origin
355  real(DP), intent(in), optional :: yorigin !< y coordinate of origin
356  real(DP), intent(in), optional :: zorigin !< z coordinate of origin
357  real(DP), intent(in), optional :: sinrot !< sine of rotation angle
358  real(DP), intent(in), optional :: cosrot !< cosine of rotation angle
359  logical(LGP), intent(in), optional :: invert !< whether to invert
360 
361  call transform(this%x, this%y, this%z, &
362  this%x, this%y, this%z, &
363  xorigin, yorigin, zorigin, &
364  sinrot, cosrot, invert)
365 
366  call compose(this%xorigin, this%yorigin, this%zorigin, &
367  this%sinrot, this%cosrot, &
368  xorigin, yorigin, zorigin, &
369  sinrot, cosrot, invert)
370 
371  this%transformed = .true.
Here is the call graph for this function:

Variable Documentation

◆ max_level

integer, parameter particlemodule::max_level = 4

Definition at line 12 of file Particle.f90.

12  integer, parameter :: MAX_LEVEL = 4