gwfgwfconnectionmodule Module Reference

Data Types
type	gwfgwfconnectiontype
	Connecting a GWF model to other models in space, implements NumericalExchangeType so the solution can used this object to determine the coefficients for the coupling between two adjacent models. More...

Functions/Subroutines
subroutine	gwfgwfconnection_ctor (this, model, gwfEx)
	Basic construction of the connection. More...
subroutine	gwfgwfcon_df (this)
	Define the connection. More...
subroutine	cfg_dist_vars (this)
	Configure distributed variables for this interface model. More...
subroutine	setgridextent (this)
	Set the required size of the interface grid from. More...
subroutine	allocatescalars (this)
	allocation of scalars in the connection More...
subroutine	gwfgwfcon_ar (this)
	Allocate and read the connection. More...
subroutine	gwfgwfcon_rp (this)
	Read time varying data when required. More...
subroutine	gwfgwfcon_ad (this)
	Advance this connection. More...
subroutine	gwfgwfcon_cf (this, kiter)
subroutine	gwfgwfcon_fc (this, kiter, matrix_sln, rhs_sln, inwtflag)
	Write the calculated coefficients into the global. More...
subroutine	validateconnection (this)
	Validate this connection This is called before proceeding to construct the interface model. More...
subroutine	validategwfexchange (this)
	Validate the exchange, intercepting those cases where two models have to be connected with an interface model, where the individual configurations don't allow this. More...
subroutine	gwfgwfcon_da (this)
	Deallocate all resources. More...
subroutine	gwfgwfcon_cq (this, icnvg, isuppress_output, isolnid)
	Calculate intra-cell flows The calculation will be dispatched to the interface model, and then mapped back to real-world cell ids. More...
subroutine	setflowtoexchange (this)
	Set the flows (flowja from interface model) to the. More...
subroutine	setflowtomodel (this)
	Set the flows (flowja from the interface model) to. More...
subroutine	setnpfedgeprops (this)
	Set flowja as edge properties in the model,. More...
subroutine	gwfgwfcon_bd (this, icnvg, isuppress_output, isolnid)
	Calculate the budget terms for this connection, this is dispatched to the GWF-GWF exchange. More...
subroutine	gwfgwfcon_ot (this)
	Write output for exchange (and calls. More...
class(gwfgwfconnectiontype) function, pointer, public	castasgwfgwfconnection (obj)
	Cast to GwfGwfConnectionType. More...

Function/Subroutine Documentation

allocatescalars()

subroutine gwfgwfconnectionmodule::allocatescalars ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 254 of file GwfGwfConnection.f90.

    use memorymanagermodule, only: mem_allocate
    class(GwfGwfConnectionType) :: this !< the connection
    ! local
 
    call mem_allocate(this%iXt3dOnExchange, 'IXT3DEXG', this%memoryPath)
 

castasgwfgwfconnection()

class(gwfgwfconnectiontype) function, pointer, public gwfgwfconnectionmodule::castasgwfgwfconnection

(

class(*), intent(inout), pointer

obj

)

Parameters

[in,out]

obj

object to be cast

Returns

the GwfGwfConnection

Definition at line 712 of file GwfGwfConnection.f90.

    implicit none
    class(*), pointer, intent(inout) :: obj !< object to be cast
    class(GwfGwfConnectionType), pointer :: res !< the GwfGwfConnection
 
    res => null()
    if (.not. associated(obj)) return
 
    select type (obj)
    class is (gwfgwfconnectiontype)
      res => obj
    end select

Here is the caller graph for this function:

cfg_dist_vars()

subroutine gwfgwfconnectionmodule::cfg_dist_vars ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 203 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< the connection
 
    call this%cfg_dv('X', '', sync_nds, &
                     (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
    call this%cfg_dv('IBOUND', '', sync_nds, &
                     (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
    call this%cfg_dv('XOLD', '', sync_nds, (/stg_bfr_exg_ad, stg_bfr_exg_cf/))
    call this%cfg_dv('ICELLTYPE', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    call this%cfg_dv('K11', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    call this%cfg_dv('K22', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    call this%cfg_dv('K33', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    if (this%gwfInterfaceModel%npf%iangle1 == 1) then
      call this%cfg_dv('ANGLE1', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    end if
    if (this%gwfInterfaceModel%npf%iangle2 == 1) then
      call this%cfg_dv('ANGLE2', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    end if
    if (this%gwfInterfaceModel%npf%iangle3 == 1) then
      call this%cfg_dv('ANGLE3', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    end if
    if (this%gwfInterfaceModel%npf%iwetdry == 1) then
      call this%cfg_dv('WETDRY', 'NPF', sync_nds, (/stg_bfr_con_ar/))
    end if
    call this%cfg_dv('TOP', 'DIS', sync_nds, (/stg_bfr_con_ar/))
    call this%cfg_dv('BOT', 'DIS', sync_nds, (/stg_bfr_con_ar/))
    call this%cfg_dv('AREA', 'DIS', sync_nds, (/stg_bfr_con_ar/))
    if (this%gwfInterfaceModel%inbuy > 0) then
      call this%cfg_dv('DENSE', 'BUY', sync_nds, (/stg_bfr_exg_cf/))
    end if
 

gwfgwfcon_ad()

subroutine gwfgwfconnectionmodule::gwfgwfcon_ad ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 306 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
 
    ! this triggers the BUY density calculation
    !if (this%gwfInterfaceModel%inbuy > 0) call this%gwfInterfaceModel%buy%buy_ad()
 
    if (this%owns_exchange) then
      call this%gwfExchange%exg_ad()
    end if
 

gwfgwfcon_ar()

subroutine gwfgwfconnectionmodule::gwfgwfcon_ar

(

class(gwfgwfconnectiontype)

this

)

Parameters

this	this connection

Definition at line 265 of file GwfGwfConnection.f90.

    use gridconnectionmodule, only: gridconnectiontype
    class(GwfGwfConnectionType) :: this !< this connection
    ! local
 
    ! check if we can construct an interface model
    ! NB: only makes sense after the models' allocate&read have been
    ! called, which is why we do it here
    call this%validateConnection()
 
    ! allocate and read base
    call this%spatialcon_ar()
 
    ! ... and now the interface model
    call this%gwfInterfaceModel%model_ar()
 
    ! AR the movers and obs through the exchange
    if (this%owns_exchange) then
      if (this%gwfExchange%inmvr > 0) then
        call this%gwfExchange%mvr%mvr_ar()
      end if
      if (this%gwfExchange%inobs > 0) then
        call this%gwfExchange%obs%obs_ar()
      end if
    end if
 

gwfgwfcon_bd()

subroutine gwfgwfconnectionmodule::gwfgwfcon_bd ( class(gwfgwfconnectiontype)  this, integer(i4b), intent(inout)  icnvg, integer(i4b), intent(in)  isuppress_output, integer(i4b), intent(in)  isolnid  )

private

Parameters

	this	this connection
[in,out]	icnvg	convergence flag
[in]	isuppress_output	suppress output when =1
[in]	isolnid	solution id

Definition at line 680 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    integer(I4B), intent(inout) :: icnvg !< convergence flag
    integer(I4B), intent(in) :: isuppress_output !< suppress output when =1
    integer(I4B), intent(in) :: isolnid !< solution id
    ! local
 
    ! call exchange budget routine, also calls bd
    ! for movers.
    if (this%owns_exchange) then
      call this%gwfExchange%exg_bd(icnvg, isuppress_output, isolnid)
    end if
 

gwfgwfcon_cf()

subroutine gwfgwfconnectionmodule::gwfgwfcon_cf ( class(gwfgwfconnectiontype)  this, integer(i4b), intent(in)  kiter  )

private

Parameters

	this	this connection
[in]	kiter	the iteration counter

Definition at line 318 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    integer(I4B), intent(in) :: kiter !< the iteration counter
 
    call this%SpatialModelConnectionType%spatialcon_cf(kiter)
 
    ! CF the movers through the exchange
    if (this%owns_exchange) then
      if (this%gwfExchange%inmvr > 0) then
        call this%gwfExchange%mvr%xmvr_cf()
      end if
    end if
 

gwfgwfcon_cq()

subroutine gwfgwfconnectionmodule::gwfgwfcon_cq	(	class(gwfgwfconnectiontype)	this,
		integer(i4b), intent(inout)	icnvg,
		integer(i4b), intent(in)	isuppress_output,
		integer(i4b), intent(in)	isolnid
	)

Parameters

	this	this connection
[in,out]	icnvg	convergence flag
[in]	isuppress_output	suppress output when =1
[in]	isolnid	solution id

Definition at line 498 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    integer(I4B), intent(inout) :: icnvg !< convergence flag
    integer(I4B), intent(in) :: isuppress_output !< suppress output when =1
    integer(I4B), intent(in) :: isolnid !< solution id
 
    call this%gwfInterfaceModel%model_cq(icnvg, isuppress_output)
 
    call this%setFlowToExchange()
 
    call this%setFlowToModel()
 
    !cdl Could we allow GwfExchange to do this instead, using
    !    simvals?
    ! if needed, we add the edge properties to the model's NPF
    ! package for its spdis calculation:
    if (this%gwfModel%npf%icalcspdis == 1) then
      call this%setNpfEdgeProps()
    end if
 
    ! Add exchange flows to each model flowja diagonal.  This used
    ! to be done in setNpfEdgeProps, but there was a sign issue
    ! and flowja was only updated if icalcspdis was 1 (it should
    ! always be updated.
    if (this%owns_exchange) then
      call this%gwfExchange%gwf_gwf_add_to_flowja()
    end if
 

gwfgwfcon_da()

subroutine gwfgwfconnectionmodule::gwfgwfcon_da

(

class(gwfgwfconnectiontype)

this

)

Parameters

this	this connection

Definition at line 468 of file GwfGwfConnection.f90.

    use kindmodule, only: lgp
    class(GwfGwfConnectionType) :: this !< this connection
    ! local
    logical(LGP) :: isOpen
 
    ! scalars
    call mem_deallocate(this%iXt3dOnExchange)
 
    call this%gwfInterfaceModel%model_da()
    deallocate (this%gwfInterfaceModel)
 
    call this%spatialcon_da()
 
    inquire (this%iout, opened=isopen)
    if (isopen) then
      close (this%iout)
    end if
 
    ! we need to deallocate the baseexchange we own:
    if (this%owns_exchange) then
      call this%gwfExchange%exg_da()
    end if
 

gwfgwfcon_df()

subroutine gwfgwfconnectionmodule::gwfgwfcon_df

(

class(gwfgwfconnectiontype)

this

)

This sets up the GridConnection (for creating the interface grid), creates and defines the interface

Parameters

this	this connection

Definition at line 142 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    ! local
    character(len=LENCOMPONENTNAME) :: imName !< the interface model's name
    integer(I4B) :: i
 
    ! this sets up the GridConnection
    call this%spatialcon_df()
 
    ! Now grid conn is defined, we create the interface model
    ! here, and the remainder of this routine is define.
    ! we basically follow the logic that is present in sln_df()
    if (this%prim_exchange%v_model1 == this%owner) then
      write (imname, '(a,i0)') 'GWFIM1_', this%gwfExchange%id
    else
      write (imname, '(a,i0)') 'GWFIM2_', this%gwfExchange%id
    end if
    call this%gwfInterfaceModel%gwfifm_cr(imname, this%iout, this%ig_builder)
    call this%gwfInterfaceModel%set_idsoln(this%gwfModel%idsoln)
    this%gwfInterfaceModel%npf%satomega = this%gwfModel%npf%satomega
    this%gwfInterfaceModel%npf%ixt3d = this%iXt3dOnExchange
    this%gwfInterfaceModel%npf%icellavg = this%gwfExchange%icellavg
    call this%gwfInterfaceModel%model_df()
 
    ! Take these settings from the owning model
    this%gwfInterfaceModel%npf%ik22 = this%gwfModel%npf%ik22
    this%gwfInterfaceModel%npf%ik33 = this%gwfModel%npf%ik33
    this%gwfInterfaceModel%npf%iwetdry = this%gwfModel%npf%iwetdry
    this%gwfInterfaceModel%npf%iangle1 = this%gwfModel%npf%iangle1
    this%gwfInterfaceModel%npf%iangle2 = this%gwfModel%npf%iangle2
    this%gwfInterfaceModel%npf%iangle3 = this%gwfModel%npf%iangle3
 
    call this%cfg_dist_vars()
 
    if (this%gwfInterfaceModel%npf%ixt3d > 0) then
      this%gwfInterfaceModel%npf%iangle1 = 1
      this%gwfInterfaceModel%npf%iangle2 = 1
      this%gwfInterfaceModel%npf%iangle3 = 1
    end if
 
    ! set defaults
    do i = 1, size(this%gwfInterfaceModel%npf%angle1)
      this%gwfInterfaceModel%npf%angle1(i) = 0.0_dp
    end do
    do i = 1, size(this%gwfInterfaceModel%npf%angle2)
      this%gwfInterfaceModel%npf%angle2(i) = 0.0_dp
    end do
    do i = 1, size(this%gwfInterfaceModel%npf%angle3)
      this%gwfInterfaceModel%npf%angle3(i) = 0.0_dp
    end do
 
    ! point X, RHS, IBOUND to connection
    call this%spatialcon_setmodelptrs()
 
    ! connect interface model to spatial connection
    call this%spatialcon_connect()
 

gwfgwfcon_fc()

subroutine gwfgwfconnectionmodule::gwfgwfcon_fc ( class(gwfgwfconnectiontype)  this, integer(i4b), intent(in)  kiter, class(matrixbasetype), pointer  matrix_sln, real(dp), dimension(:), intent(inout)  rhs_sln, integer(i4b), intent(in), optional  inwtflag  )

private

Parameters

	this	this connection
[in]	kiter	the iteration counter
	matrix_sln	global system matrix coefficients
[in,out]	rhs_sln	global right-hand-side
[in]	inwtflag	newton-raphson flag

Definition at line 335 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    integer(I4B), intent(in) :: kiter !< the iteration counter
    class(MatrixBaseType), pointer :: matrix_sln !< global system matrix coefficients
    real(DP), dimension(:), intent(inout) :: rhs_sln !< global right-hand-side
    integer(I4B), optional, intent(in) :: inwtflag !< newton-raphson flag
    ! local
 
    call this%SpatialModelConnectionType%spatialcon_fc( &
      kiter, matrix_sln, rhs_sln, inwtflag)
 
    ! FC the movers through the exchange; we cannot call
    ! exg_fc() directly because it calculates matrix terms
    if (this%owns_exchange) then
      if (this%gwfExchange%inmvr > 0) then
        call this%gwfExchange%mvr%mvr_fc()
      end if
    end if
 

gwfgwfcon_ot()

subroutine gwfgwfconnectionmodule::gwfgwfcon_ot ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 697 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    ! local
 
    ! Call exg_ot() here as it handles all output processing
    ! based on gwfExchange%simvals(:), which was correctly
    ! filled from gwfgwfcon
    if (this%owns_exchange) then
      call this%gwfExchange%exg_ot()
    end if
 

gwfgwfcon_rp()

subroutine gwfgwfconnectionmodule::gwfgwfcon_rp

(

class(gwfgwfconnectiontype)

this

)

Parameters

this	this connection

Definition at line 295 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
 
    ! Call exchange rp routines
    if (this%owns_exchange) then
      call this%gwfExchange%exg_rp()
    end if

gwfgwfconnection_ctor()

subroutine gwfgwfconnectionmodule::gwfgwfconnection_ctor	(	class(gwfgwfconnectiontype)	this,
		class(numericalmodeltype), pointer	model,
		class(disconnexchangetype), pointer	gwfEx
	)

Parameters

this	the connection
model	the model owning this connection, this must of course be a GwfModelType
gwfex	the exchange the interface model is created for

Definition at line 81 of file GwfGwfConnection.f90.

    use numericalmodelmodule, only: numericalmodeltype
    use inputoutputmodule, only: openfile
    class(GwfGwfConnectionType) :: this !< the connection
    class(NumericalModelType), pointer :: model !< the model owning this connection,
                                                !! this must of course be a GwfModelType
    class(DisConnExchangeType), pointer :: gwfEx !< the exchange the interface model is created for
    ! local
    character(len=LINELENGTH) :: fname
    character(len=LENCOMPONENTNAME) :: name
    class(*), pointer :: objPtr
    logical(LGP) :: write_ifmodel_listfile = .false.
 
    objptr => model
    this%gwfModel => castasgwfmodel(objptr)
    objptr => gwfex
    this%gwfExchange => castasgwfexchange(objptr)
 
    if (gwfex%v_model1%is_local .and. gwfex%v_model2%is_local) then
      this%owns_exchange = (gwfex%v_model1 == model)
    else
      this%owns_exchange = .true.
    end if
 
    if (gwfex%v_model1 == model) then
      write (name, '(a,i0)') 'GWFCON1_', gwfex%id
    else
      write (name, '(a,i0)') 'GWFCON2_', gwfex%id
    end if
 
    ! .lst file for interface model
    if (write_ifmodel_listfile) then
      fname = trim(name)//'.im.lst'
      call openfile(this%iout, 0, fname, 'LIST', filstat_opt='REPLACE')
      write (this%iout, '(4a)') 'Creating GWF-GWF connection for model ', &
        trim(this%gwfModel%name), ' from exchange ', &
        trim(gwfex%name)
    end if
 
    ! first call base constructor
    call this%SpatialModelConnectionType%spatialConnection_ctor(model, &
                                                                gwfex, &
                                                                name)
 
    call this%allocateScalars()
 
    this%typename = 'GWF-GWF'
 
    ! determine the required size of the interface grid
    call this%setGridExtent()
 
    allocate (this%gwfInterfaceModel)
    this%interface_model => this%gwfInterfaceModel
 

Here is the call graph for this function:

Here is the caller graph for this function:

setflowtoexchange()

subroutine gwfgwfconnectionmodule::setflowtoexchange ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 530 of file GwfGwfConnection.f90.

    use indexmapmodule
    class(GwfGwfConnectionType) :: this !< this connection
    ! local
    integer(I4B) :: i
    class(GwfExchangeType), pointer :: gwfEx
    type(IndexMapSgnType), pointer :: map
 
    if (this%owns_exchange) then
      gwfex => this%gwfExchange
      map => this%interface_map%exchange_maps(this%interface_map%prim_exg_idx)
 
      ! use (half of) the exchange map in reverse:
      do i = 1, size(map%src_idx)
        if (map%sign(i) < 0) cycle ! simvals is defined from exg%m1 => exg%m2
        gwfex%simvals(map%src_idx(i)) = &
          this%gwfInterfaceModel%flowja(map%tgt_idx(i))
      end do
    end if
 

setflowtomodel()

subroutine gwfgwfconnectionmodule::setflowtomodel

(

class(gwfgwfconnectiontype)

this

)

Parameters

this	this connection

Definition at line 554 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    ! local
    integer(I4B) :: n, m, ipos, iposLoc
    integer(I4B) :: nLoc, mLoc
    type(ConnectionsType), pointer :: imCon !< interface model connections
    type(GlobalCellType), dimension(:), pointer :: toGlobal !< map interface index to global cell
 
    ! for readability
    imcon => this%gwfInterfaceModel%dis%con
    toglobal => this%ig_builder%idxToGlobal
 
    do n = 1, this%neq
      if (.not. toglobal(n)%v_model == this%owner) then
        ! only add flows to own model
        cycle
      end if
 
      nloc = toglobal(n)%index
 
      do ipos = imcon%ia(n) + 1, imcon%ia(n + 1) - 1
        if (imcon%mask(ipos) < 1) cycle ! skip this connection, it's masked so not determined by us
 
        m = imcon%ja(ipos)
        mloc = toglobal(m)%index
        if (toglobal(m)%v_model == this%owner) then
 
          ! internal, need to set flowja for n-m
          iposloc = getcsrindex(nloc, mloc, this%gwfModel%ia, this%gwfModel%ja)
 
          ! update flowja with correct value
          this%gwfModel%flowja(iposloc) = this%gwfInterfaceModel%flowja(ipos)
 
        end if
      end do
    end do
 

Here is the call graph for this function:

setgridextent()

subroutine gwfgwfconnectionmodule::setgridextent ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 238 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< the connection
    ! local
 
    this%iXt3dOnExchange = this%gwfExchange%ixt3d
    if (this%iXt3dOnExchange > 0) then
      this%exg_stencil_depth = 2
      if (this%gwfModel%npf%ixt3d > 0) then
        this%int_stencil_depth = 2
      end if
    end if
 

setnpfedgeprops()

subroutine gwfgwfconnectionmodule::setnpfedgeprops ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 595 of file GwfGwfConnection.f90.

    class(GwfGwfConnectionType) :: this !< this connection
    ! local
    integer(I4B) :: n, m, ipos, isym
    integer(I4B) :: nLoc, mLoc
    integer(I4B) :: ihc
    real(DP) :: rrate
    real(DP) :: area
    real(DP) :: satThick
    real(DP) :: nx, ny, nz
    real(DP) :: cx, cy, cz
    real(DP) :: conLen
    real(DP) :: dist
    real(DP) :: cl
    logical :: nozee
    type(ConnectionsType), pointer :: imCon !< interface model connections
    class(GwfNpfType), pointer :: imNpf !< interface model npf package
    class(DisBaseType), pointer :: imDis !< interface model discretization
    type(GlobalCellType), dimension(:), pointer :: toGlobal !< map interface index to global cell
 
    ! for readability
    imdis => this%gwfInterfaceModel%dis
    imcon => this%gwfInterfaceModel%dis%con
    imnpf => this%gwfInterfaceModel%npf
    toglobal => this%ig_builder%idxToGlobal
 
    nozee = .false.
    if (imnpf%ixt3d > 0) then
      nozee = imnpf%xt3d%nozee
    end if
 
    ! loop over flowja in the interface model and set edge properties
    ! for flows crossing the boundary, and set flowja for internal
    ! flows affected by the connection.
    do n = 1, this%neq
      if (.not. toglobal(n)%v_model == this%owner) then
        ! only add flows to own model
        cycle
      end if
 
      nloc = toglobal(n)%index
 
      do ipos = imcon%ia(n) + 1, imcon%ia(n + 1) - 1
        if (imcon%mask(ipos) < 1) then
          ! skip this connection, it's masked so not determined by us
          cycle
        end if
 
        m = imcon%ja(ipos)
        mloc = toglobal(m)%index
 
        if (.not. toglobal(m)%v_model == this%owner) then
          ! boundary connection, set edge properties
          isym = imcon%jas(ipos)
          ihc = imcon%ihc(isym)
          area = imcon%hwva(isym)
          satthick = imnpf%calcSatThickness(n, m, ihc)
          rrate = this%gwfInterfaceModel%flowja(ipos)
 
          call imdis%connection_normal(n, m, ihc, nx, ny, nz, ipos)
          call imdis%connection_vector(n, m, nozee, imnpf%sat(n), imnpf%sat(m), &
                                       ihc, cx, cy, cz, conlen)
 
          if (ihc == 0) then
            ! check if n is below m
            if (nz > 0) rrate = -rrate
          else
            area = area * satthick
          end if
 
          cl = imcon%cl1(isym)
          if (m < n) then
            cl = imcon%cl2(isym)
          end if
          dist = conlen * cl / (imcon%cl1(isym) + imcon%cl2(isym))
          call this%gwfModel%npf%set_edge_properties(nloc, ihc, rrate, area, &
                                                     nx, ny, dist)
        end if
      end do
    end do
 

validateconnection()

subroutine gwfgwfconnectionmodule::validateconnection ( class(gwfgwfconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 360 of file GwfGwfConnection.f90.

    use simvariablesmodule, only: errmsg
    use simmodule, only: count_errors
    class(GwfGwfConnectionType) :: this !< this connection
    ! local
 
    ! base validation (geometry/spatial)
    call this%SpatialModelConnectionType%validateConnection()
    call this%validateGwfExchange()
 
    ! abort on errors
    if (count_errors() > 0) then
      write (errmsg, '(a)') 'Errors occurred while processing exchange(s)'
      call ustop(errmsg)
    end if
 

Here is the call graph for this function:

validategwfexchange()

subroutine gwfgwfconnectionmodule::validategwfexchange

(

class(gwfgwfconnectiontype)

this

)

Stops with error message on config mismatch

Parameters

this	this connection

Definition at line 384 of file GwfGwfConnection.f90.

    use simvariablesmodule, only: errmsg, simulation_mode
    use simmodule, only: store_error
    use gwfnpfmodule, only: gwfnpftype
    class(GwfGwfConnectionType) :: this !< this connection
    ! local
    class(GwfExchangeType), pointer :: gwfEx
    class(*), pointer :: modelPtr
    class(GwfModelType), pointer :: gwfModel1
    class(GwfModelType), pointer :: gwfModel2
    type(GwfBuyType), pointer :: buy1, buy2
    logical(LGP) :: compatible
 
    gwfex => this%gwfExchange
 
    ! GNC not allowed
    if (gwfex%ingnc /= 0 .and. gwfex%ixt3d /= 0) then
      write (errmsg, '(2a)') 'Ghost node correction not supported '// &
        'combined with XT3D for exchange ', trim(gwfex%name)
      call store_error(errmsg)
    end if
    if (gwfex%ingnc /= 0 .and. simulation_mode == 'PARALLEL') then
      write (errmsg, '(2a)') 'Ghost node correction not supported '// &
        'in parallel run for exchange ', trim(gwfex%name)
      call store_error(errmsg)
    end if
 
    ! we cannot validate the remainder (yet) in parallel mode
    if (.not. gwfex%v_model1%is_local) return
    if (.not. gwfex%v_model2%is_local) return
 
    modelptr => this%gwfExchange%model1
    gwfmodel1 => castasgwfmodel(modelptr)
    modelptr => this%gwfExchange%model2
    gwfmodel2 => castasgwfmodel(modelptr)
 
    if ((gwfmodel1%inbuy > 0 .and. gwfmodel2%inbuy == 0) .or. &
        (gwfmodel1%inbuy == 0 .and. gwfmodel2%inbuy > 0)) then
      write (errmsg, '(2a)') 'Buoyancy package should be enabled/disabled '// &
        'simultaneously in models connected with the '// &
        'interface model for exchange ', &
        trim(gwfex%name)
      call store_error(errmsg)
 
    end if
 
    if (gwfmodel1%inbuy > 0 .and. gwfmodel2%inbuy > 0) then
      ! does not work with XT3D
      if (this%iXt3dOnExchange > 0) then
        write (errmsg, '(2a)') 'Connecting models with BUY package not '// &
          'allowed with XT3D enabled on exchange ', &
          trim(gwfex%name)
        call store_error(errmsg)
      end if
 
      ! check compatibility of buoyancy
      compatible = .true.
      buy1 => gwfmodel1%buy
      buy2 => gwfmodel2%buy
      if (buy1%iform /= buy2%iform) compatible = .false.
      if (buy1%denseref /= buy2%denseref) compatible = .false.
      if (buy1%nrhospecies /= buy2%nrhospecies) compatible = .false.
      if (.not. all(buy1%drhodc == buy2%drhodc)) compatible = .false.
      if (.not. all(buy1%crhoref == buy2%crhoref)) compatible = .false.
      if (.not. all(buy1%cauxspeciesname == buy2%cauxspeciesname)) then
        compatible = .false.
      end if
 
      if (.not. compatible) then
        write (errmsg, '(6a)') 'Buoyancy packages in model ', &
          trim(gwfex%model1%name), ' and ', &
          trim(gwfex%model2%name), &
          ' should be equivalent to construct an '// &
          ' interface model for exchange ', &
          trim(gwfex%name)
        call store_error(errmsg)
      end if
 
    end if
 

Here is the call graph for this function: