tspexchangemovermodule Module Reference

Data Types
type	tspexchangemovertype

Functions/Subroutines
subroutine, public	xmvt_cr (mvt, name_exg, tsp_model1, tsp_model2, gwfmodelname1, gwfmodelname2, inunit, iout)
subroutine	xmvt_rp (this)
	Read and prepare mover transport object. More...
subroutine	xmvt_cf (this, cnew1, cnew2)
subroutine	xmvt_fc (this, cnew1, cnew2)
subroutine	xmvt_da (this)

Function/Subroutine Documentation

xmvt_cf()

subroutine tspexchangemovermodule::xmvt_cf	(	class(tspexchangemovertype)	this,
		real(dp), dimension(:), intent(in), target, contiguous	cnew1,
		real(dp), dimension(:), intent(in), target, contiguous	cnew2
	)

Definition at line 106 of file TspExchangeMover.f90.

    class(TspExchangeMoverType) :: this
    real(DP), intent(in), dimension(:), contiguous, target :: cnew1
    real(DP), intent(in), dimension(:), contiguous, target :: cnew2
    ! local
    integer(I4B) :: i, n, idx
    real(DP), dimension(:), pointer, contiguous :: prov_array
    logical(LGP) :: prov_is_local
 
    call this%mvt_fill_mvrterm(cnew1, cnew2)
 
    ! add mvrterm data into synchronization array
    idx = 1
    do i = 1, size(this%mvrterm)
 
      ! point to provider array
      prov_is_local = .false.
      if (this%mvrterm(i)%prov_is_m1) then
        prov_is_local = this%model1%is_local
        prov_array => this%quantity_m1
      else
        prov_is_local = this%model2%is_local
        prov_array => this%quantity_m2
      end if
 
      if (prov_is_local) then
        do n = 1, size(this%mvrterm(i)%qty)
          prov_array(idx) = this%mvrterm(i)%qty(n)
          idx = idx + 1
        end do
      end if
    end do
 

xmvt_cr()

subroutine, public tspexchangemovermodule::xmvt_cr	(	type(tspexchangemovertype), pointer	mvt,
		character(len=*), intent(in)	name_exg,
		class(transportmodeltype), pointer	tsp_model1,
		class(transportmodeltype), pointer	tsp_model2,
		character(len=*), intent(in)	gwfmodelname1,
		character(len=*), intent(in)	gwfmodelname2,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout
	)

Parameters

[in]	name_exg	name of the exchange
	tsp_model1	gwt model 1, can be null()
	tsp_model2	gwt model 2, can be null()

Definition at line 31 of file TspExchangeMover.f90.

    type(TspExchangeMoverType), pointer :: mvt
    character(len=*), intent(in) :: name_exg !< name of the exchange
    class(TransportModelType), pointer :: tsp_model1 !< gwt model 1, can be null()
    class(TransportModelType), pointer :: tsp_model2 !< gwt model 2, can be null()
    character(len=*), intent(in) :: gwfmodelname1
    character(len=*), intent(in) :: gwfmodelname2
    integer(I4B), intent(in) :: inunit
    integer(I4B), intent(in) :: iout
    ! local
    type(TspFmiType), pointer :: fmi1, fmi2
    real(DP), pointer :: eqnsclfac !< governing equation scale factor
    character(len=LENVARNAME) :: depvartype !< dependent variable type ('concentration' or 'temperature')
 
    allocate (mvt)
 
    fmi1 => null()
    if (associated(tsp_model1)) then
      fmi1 => tsp_model1%fmi
      eqnsclfac => tsp_model1%eqnsclfac
      depvartype = tsp_model1%depvartype
    end if
    fmi2 => null()
    if (associated(tsp_model2)) then
      fmi2 => tsp_model2%fmi
      eqnsclfac => tsp_model2%eqnsclfac
      depvartype = tsp_model2%depvartype
    end if
 
    mvt%model1 => get_virtual_model(gwfmodelname1)
    mvt%model2 => get_virtual_model(gwfmodelname2)
 
    call mvt%mvt_init(name_exg, inunit, iout, fmi1, &
                      eqnsclfac, depvartype, gwfmodelname1, &
                      gwfmodelname2, fmi2)
 
    call mem_allocate(mvt%quantity_m1, 0, "QUANTITY_M1", mvt%memoryPath)
    call mem_allocate(mvt%quantity_m2, 0, "QUANTITY_M2", mvt%memoryPath)
    call mem_allocate(mvt%maxmvt, "MAXMVT", mvt%memoryPath)
    mvt%maxmvt = 0
 

Here is the caller graph for this function:

xmvt_da()

subroutine tspexchangemovermodule::xmvt_da ( class(tspexchangemovertype)  this )

private

Definition at line 176 of file TspExchangeMover.f90.

    class(TspExchangeMoverType) :: this
 
    call this%TspMvtType%mvt_da()
 
    call mem_deallocate(this%quantity_m1)
    call mem_deallocate(this%quantity_m2)
    call mem_deallocate(this%maxmvt)
 

xmvt_fc()

subroutine tspexchangemovermodule::xmvt_fc ( class(tspexchangemovertype)  this, real(dp), dimension(:), intent(in), target, contiguous  cnew1, real(dp), dimension(:), intent(in), target, contiguous  cnew2  )

private

Definition at line 141 of file TspExchangeMover.f90.

    class(TspExchangeMoverType) :: this
    real(DP), intent(in), dimension(:), contiguous, target :: cnew1
    real(DP), intent(in), dimension(:), contiguous, target :: cnew2
    ! local
    integer(I4B) :: i, n, idx
    real(DP), dimension(:), pointer, contiguous :: prov_array ! after synchronization when parallel
    logical(LGP) :: prov_is_remote
 
    ! get mvrterm data from synchronization array
    idx = 1
    do i = 1, size(this%mvrterm)
 
      ! point to provider array
      if (this%mvrterm(i)%prov_is_m1) then
        prov_is_remote = .not. this%model1%is_local
        prov_array => this%quantity_m1
      else
        prov_is_remote = .not. this%model2%is_local
        prov_array => this%quantity_m2
      end if
 
      if (prov_is_remote) then
        ! remote provider, so we copy it in here after sync'ing
        do n = 1, size(this%mvrterm(i)%qty)
          this%mvrterm(i)%qty(n) = prov_array(idx)
          idx = idx + 1
        end do
      end if
    end do
 
    call this%mvt_update_qmfrommvr()
 

xmvt_rp()

subroutine tspexchangemovermodule::xmvt_rp ( class(tspexchangemovertype)  this )

private

Definition at line 77 of file TspExchangeMover.f90.

    use tdismodule, only: kper, kstp
    class(TspExchangeMoverType) :: this
    ! local
    integer(I4B) :: i, max_array_size
 
    call this%TspMvtType%mvt_rp()
 
    ! only on first timestep
    if (kstp * kper == 1) then
      max_array_size = 0
      do i = 1, this%mvrbudobj%nbudterm
        max_array_size = max_array_size + this%mvrbudobj%budterm(i)%maxlist
      end do
 
      if (max_array_size > 0) then
        call mem_reallocate(this%quantity_m1, max_array_size, &
                            "QUANTITY_M1", this%memoryPath)
        call mem_reallocate(this%quantity_m2, max_array_size, &
                            "QUANTITY_M2", this%memoryPath)
        this%maxmvt = max_array_size
      end if
    end if
 
    this%quantity_m1 = dnodata
    this%quantity_m2 = dnodata