gwegweconnectionmodule Module Reference

Data Types
type	gwegweconnectiontype
	Connects a GWE model to other GWE models in space. Derives from NumericalExchangeType so the solution can use it to fetch the coefficients for this connection. More...

Functions/Subroutines
subroutine	gwegweconnection_ctor (this, model, gweEx)
	Basic construction of the connection. More...
subroutine	allocate_scalars (this)
	Allocate scalar variables for this connection. More...
subroutine	gwegwecon_df (this)
	define the GWE-GWE connection More...
subroutine	cfg_dist_vars (this)
	Configure distributed variables for this interface model. More...
subroutine	allocate_arrays (this)
	Allocate array variables for this connection. More...
subroutine	setgridextent (this)
	Set required extent of the interface grid from. More...
subroutine	gwegwecon_ar (this)
	allocate and read/set the connection's data structures More...
subroutine	validateconnection (this)
	validate this connection prior to constructing More...
subroutine	gwegwecon_rp (this)
subroutine	gwegwecon_ad (this)
	Advance this connection. More...
subroutine	gwegwecon_cf (this, kiter)
subroutine	gwegwecon_fc (this, kiter, matrix_sln, rhs_sln, inwtflag)
subroutine	gwegwecon_cq (this, icnvg, isuppress_output, isolnid)
subroutine	setflowtoexchange (this)
	Set the flows (flowja from interface model) to the. More...
subroutine	gwegwecon_bd (this, icnvg, isuppress_output, isolnid)
subroutine	gwegwecon_ot (this)
subroutine	validategweexchange (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	gwegwecon_da (this)
	Deallocate all resources. More...
class(gwegweconnectiontype) function, pointer, public	castasgwegweconnection (obj)
	Cast to GweGweConnectionType. More...

Function/Subroutine Documentation

allocate_arrays()

subroutine gwegweconnectionmodule::allocate_arrays ( class(gwegweconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 248 of file GweGweConnection.f90.

    class(GweGweConnectionType) :: this !< the connection
 
    call mem_allocate(this%exgflowjaGwe, this%ig_builder%nrOfBoundaryCells, &
                      'EXGFLOWJAGWT', this%memoryPath)
 

allocate_scalars()

subroutine gwegweconnectionmodule::allocate_scalars

(

class(gwegweconnectiontype)

this

)

Parameters

this	the connection

Definition at line 143 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
 
    call mem_allocate(this%iIfaceAdvScheme, 'IADVSCHEME', this%memoryPath)
    call mem_allocate(this%iIfaceXt3d, 'IXT3D', this%memoryPath)
    call mem_allocate(this%exgflowSign, 'EXGFLOWSIGN', this%memoryPath)
 

castasgwegweconnection()

class(gwegweconnectiontype) function, pointer, public gwegweconnectionmodule::castasgwegweconnection

(

class(*), intent(inout), pointer

obj

)

Parameters

[in,out]

obj

object to be cast

Returns

the GweGweConnection

Definition at line 639 of file GweGweConnection.f90.

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

Here is the caller graph for this function:

cfg_dist_vars()

subroutine gwegweconnectionmodule::cfg_dist_vars ( class(gwegweconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 213 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: 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/))
    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%gweInterfaceModel%cnd%idisp > 0) then
      call this%cfg_dv('ALH', 'CND', sync_nds, (/stg_bfr_con_ar/))
      call this%cfg_dv('ALV', 'CND', sync_nds, (/stg_bfr_con_ar/))
      call this%cfg_dv('ATH1', 'CND', sync_nds, (/stg_bfr_con_ar/))
      call this%cfg_dv('ATH2', 'CND', sync_nds, (/stg_bfr_con_ar/))
      call this%cfg_dv('ATV', 'CND', sync_nds, (/stg_bfr_con_ar/))
      call this%cfg_dv('KTW', 'CND', sync_nds, (/stg_bfr_con_ar/))
      call this%cfg_dv('KTS', 'CND', sync_nds, (/stg_bfr_con_ar/))
    end if
    call this%cfg_dv('GWFHEAD', 'FMI', sync_nds, (/stg_bfr_exg_ad/))
    call this%cfg_dv('GWFSAT', 'FMI', sync_nds, (/stg_bfr_exg_ad/))
    call this%cfg_dv('GWFSPDIS', 'FMI', sync_nds, (/stg_bfr_exg_ad/))
    call this%cfg_dv('GWFFLOWJA', 'FMI', sync_con, (/stg_bfr_exg_ad/))
    call this%cfg_dv('GWFFLOWJA', 'FMI', sync_exg, (/stg_bfr_exg_ad/), &
                     exg_var_name='GWFSIMVALS')
    ! fill porosity from est packages, needed for cnd
    if (this%gweModel%incnd > 0 .and. this%gweModel%inest > 0) then
      call this%cfg_dv('POROSITY', 'EST', sync_nds, (/stg_aft_con_ar/))
    end if
 

gwegwecon_ad()

subroutine gwegweconnectionmodule::gwegwecon_ad ( class(gwegweconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 389 of file GweGweConnection.f90.

 
    class(GweGweConnectionType) :: this !< this connection
 
    ! recalculate conduction ellipse
    if (this%gweInterfaceModel%incnd > 0) call this%gweInterfaceModel%cnd%cnd_ad()
 
    if (this%owns_exchange) then
      call this%gweExchange%exg_ad()
    end if
 

gwegwecon_ar()

subroutine gwegweconnectionmodule::gwegwecon_ar ( class(gwegweconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 291 of file GweGweConnection.f90.

    class(GweGweConnectionType) :: this !< the connection
 
    ! 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%gweInterfaceModel%model_ar()
 
    ! set a pointer in the interface model to the gwecommon data
    if (this%gweModel%inest > 0) then
      this%gweInterfaceModel%gwecommon%gwecpw => this%gweModel%gwecommon%gwecpw
      this%gweInterfaceModel%gwecommon%gwerhow => this%gweModel%gwecommon%gwerhow
    end if
 
    ! set the equation scaling factor in the interface model to that of
    !   underlying GWE model
    if (this%gweModel%incnd > 0) then
      this%gweInterfaceModel%ieqnsclfac = this%gweModel%cnd%eqnsclfac
    end if
 
    ! AR the movers and obs through the exchange
    if (this%owns_exchange) then
      !cdl implement this when MVT is ready
      !cdl if (this%gweExchange%inmvt > 0) then
      !cdl   call this%gweExchange%mvt%mvt_ar()
      !cdl end if
      if (this%gweExchange%inobs > 0) then
        call this%gweExchange%obs%obs_ar()
      end if
    end if
 

gwegwecon_bd()

subroutine gwegweconnectionmodule::gwegwecon_bd	(	class(gwegweconnectiontype)	this,
		integer(i4b), intent(inout)	icnvg,
		integer(i4b), intent(in)	isuppress_output,
		integer(i4b), intent(in)	isolnid
	)

Parameters

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

Definition at line 495 of file GweGweConnection.f90.

    ! modules
    use budgetmodule, only: rate_accumulator
    ! dummy
    class(GweGweConnectionType) :: this !< the 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 exchange budget routine, also calls _bd
    !    for movers.
    if (this%owns_exchange) then
      call this%gweExchange%exg_bd(icnvg, isuppress_output, isolnid)
    end if
 

Here is the call graph for this function:

gwegwecon_cf()

subroutine gwegweconnectionmodule::gwegwecon_cf ( class(gwegweconnectiontype)  this, integer(i4b), intent(in)  kiter  )

private

Parameters

	this	this connection
[in]	kiter	the iteration counter

Definition at line 402 of file GweGweConnection.f90.

    class(GweGweConnectionType) :: this !< this connection
    integer(I4B), intent(in) :: kiter !< the iteration counter
    ! local
    real(DP), dimension(:), pointer, contiguous :: x_m1, x_m2
 
    call this%SpatialModelConnectionType%spatialcon_cf(kiter)
 
    ! CF the movers in the exchange
    if (this%owns_exchange) then
      if (this%gweExchange%inmvt > 0) then
        x_m1 => null()
        x_m2 => null()
        if (associated(this%gweExchange%gwemodel1)) then
          x_m1 => this%gweExchange%gwemodel1%x
        end if
        if (associated(this%gweExchange%gwemodel2)) then
          x_m2 => this%gweExchange%gwemodel2%x
        end if
        call this%gweExchange%mvt%xmvt_cf(x_m1, x_m2)
      end if
    end if
 

gwegwecon_cq()

subroutine gwegweconnectionmodule::gwegwecon_cq ( class(gwegweconnectiontype)  this, integer(i4b), intent(inout)  icnvg, integer(i4b), intent(in)  isuppress_output, integer(i4b), intent(in)  isolnid  )

private

Parameters

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

Definition at line 457 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the 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%gweInterfaceModel%model_cq(icnvg, isuppress_output)
    call this%setFlowToExchange()
 

gwegwecon_da()

subroutine gwegweconnectionmodule::gwegwecon_da

(

class(gwegweconnectiontype)

this

)

Parameters

this	the connection

Definition at line 604 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
    ! local
    logical(LGP) :: isOpen
 
    ! scalars
    call mem_deallocate(this%iIfaceAdvScheme)
    call mem_deallocate(this%iIfaceXt3d)
    call mem_deallocate(this%exgflowSign)
 
    ! arrays
    call mem_deallocate(this%exgflowjaGwe)
 
    ! interface model
    call this%gweInterfaceModel%model_da()
    deallocate (this%gweInterfaceModel)
 
    ! dealloc base
    call this%spatialcon_da()
 
    inquire (this%iout, opened=isopen)
    if (isopen) then
      close (this%iout)
    end if
 
    ! we need to deallocate the exchange we own:
    if (this%owns_exchange) then
      call this%gweExchange%exg_da()
    end if
 

gwegwecon_df()

subroutine gwegweconnectionmodule::gwegwecon_df ( class(gwegweconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 155 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
    ! local
    character(len=LENCOMPONENTNAME) :: imName
 
    ! determine advection scheme (the GWE-GWE exchange
    !    has been read at this point)
    this%iIfaceAdvScheme = this%gweExchange%iAdvScheme
 
    ! determine xt3d setting on interface
    this%iIfaceXt3d = this%gweExchange%ixt3d
 
    ! turn off when off in the owning model
    if (this%gweModel%incnd > 0) then
      this%iIfaceXt3d = this%gweModel%cnd%ixt3d
    end if
 
    ! determine the required size of the interface model grid
    call this%setGridExtent()
 
    ! now set up the GridConnection
    call this%spatialcon_df()
 
    ! we have to 'catch up' and create the interface model
    !    here, then the remainder of this routine will be define
    if (this%prim_exchange%v_model1 == this%owner) then
      write (imname, '(a,i0)') 'GWEIM1_', this%gweExchange%id
    else
      write (imname, '(a,i0)') 'GWEIM2_', this%gweExchange%id
    end if
    call this%gweInterfaceModel%gweifmod_cr(imname, &
                                            this%iout, &
                                            this%ig_builder)
    call this%gweInterfaceModel%set_idsoln(this%gweModel%idsoln)
    this%gweInterfaceModel%iAdvScheme = this%iIfaceAdvScheme
    this%gweInterfaceModel%ixt3d = this%iIfaceXt3d
    call this%gweInterfaceModel%model_df()
 
    call this%cfg_dist_vars()
 
    call this%allocate_arrays()
    call this%gweInterfaceModel%allocate_fmi()
 
    ! connect X, RHS, IBOUND, and flowja
    call this%spatialcon_setmodelptrs()
 
    ! connect pointers (used by BUY)
    this%conc => this%gweInterfaceModel%x
    this%icbound => this%gweInterfaceModel%ibound
 
    ! add connections from the interface model to solution matrix
    call this%spatialcon_connect()
 

gwegwecon_fc()

subroutine gwegweconnectionmodule::gwegwecon_fc ( class(gwegweconnectiontype)  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	the connection
[in]	kiter	the iteration counter
	matrix_sln	the system matrix
[in,out]	rhs_sln	global right-hand-side
[in]	inwtflag	newton-raphson flag

Definition at line 427 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
    integer(I4B), intent(in) :: kiter !< the iteration counter
    class(MatrixBaseType), pointer :: matrix_sln !< the system matrix
    real(DP), dimension(:), intent(inout) :: rhs_sln !< global right-hand-side
    integer(I4B), optional, intent(in) :: inwtflag !< newton-raphson flag
    ! local
    real(DP), dimension(:), pointer, contiguous :: x_m1, x_m2
 
    call this%SpatialModelConnectionType%spatialcon_fc( &
      kiter, matrix_sln, rhs_sln, inwtflag)
 
    ! _fc the movers through the exchange
    if (this%owns_exchange) then
      if (this%gweExchange%inmvt > 0) then
        x_m1 => null()
        x_m2 => null()
        if (associated(this%gweExchange%gwemodel1)) then
          x_m1 => this%gweExchange%gwemodel1%x
        end if
        if (associated(this%gweExchange%gwemodel2)) then
          x_m2 => this%gweExchange%gwemodel2%x
        end if
        call this%gweExchange%mvt%mvt_fc(x_m1, x_m2)
      end if
    end if
 

gwegwecon_ot()

subroutine gwegweconnectionmodule::gwegwecon_ot

(

class(gwegweconnectiontype)

this

)

Parameters

this	the connection

Definition at line 512 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
 
    ! call exg_ot() here as it handles all output processing
    !    based on gweExchange%simvals(:), which was correctly
    !    filled from gwegwecon
    if (this%owns_exchange) then
      call this%gweExchange%exg_ot()
    end if
 

gwegwecon_rp()

subroutine gwegweconnectionmodule::gwegwecon_rp

(

class(gwegweconnectiontype)

this

)

Parameters

this	the connection

Definition at line 376 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
 
    ! call exchange rp routines
    if (this%owns_exchange) then
      call this%gweExchange%exg_rp()
    end if
 

gwegweconnection_ctor()

subroutine gwegweconnectionmodule::gwegweconnection_ctor	(	class(gwegweconnectiontype)	this,
		class(numericalmodeltype), pointer	model,
		class(disconnexchangetype), pointer	gweEx
	)

Parameters

this	the connection
model	the model owning this connection, this must be a GweModelType
gweex	the GWE-GWE exchange the interface model is created for

Definition at line 85 of file GweGweConnection.f90.

    ! modules
    use inputoutputmodule, only: openfile
    ! dummy
    class(GweGweConnectionType) :: this !< the connection
    class(NumericalModelType), pointer :: model !< the model owning this connection,
                                                !! this must be a GweModelType
    class(DisConnExchangeType), pointer :: gweEx !< the GWE-GWE 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%gweModel => castasgwemodel(objptr)
    objptr => gweex
    this%gweExchange => castasgweexchange(objptr)
 
    if (gweex%v_model1%is_local .and. gweex%v_model2%is_local) then
      this%owns_exchange = associated(model, gweex%model1)
    else
      this%owns_exchange = .true.
    end if
 
    if (gweex%v_model1 == model) then
      write (name, '(a,i0)') 'GWECON1_', gweex%id
    else
      write (name, '(a,i0)') 'GWECON2_', gweex%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 GWE-GWE connection for model ', &
        trim(this%gweModel%name), 'from exchange ', &
        trim(gweex%name)
    end if
 
    ! first call base constructor
    call this%SpatialModelConnectionType%spatialConnection_ctor(model, &
                                                                gweex, &
                                                                name)
 
    call this%allocate_scalars()
    this%typename = 'GWE-GWE'
    this%iIfaceAdvScheme = 0
    this%iIfaceXt3d = 0
    this%exgflowSign = 1
 
    allocate (this%gweInterfaceModel)
    this%interface_model => this%gweInterfaceModel
 

Here is the call graph for this function:

Here is the caller graph for this function:

setflowtoexchange()

subroutine gwegweconnectionmodule::setflowtoexchange ( class(gwegweconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 471 of file GweGweConnection.f90.

    ! modules
    use indexmapmodule
    ! dummy
    class(GweGweConnectionType) :: this !< this connection
    ! local
    integer(I4B) :: i
    class(GweExchangeType), pointer :: gweEx
    type(IndexMapSgnType), pointer :: map
 
    if (this%owns_exchange) then
      gweex => this%gweExchange
      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
        gweex%simvals(map%src_idx(i)) = &
          this%gweInterfaceModel%flowja(map%tgt_idx(i))
      end do
    end if
 

setgridextent()

subroutine gwegweconnectionmodule::setgridextent ( class(gwegweconnectiontype)  this )

private

Parameters

this	the connection

Definition at line 258 of file GweGweConnection.f90.

    ! dummy
    class(GweGweConnectionType) :: this !< the connection
    ! local
    logical(LGP) :: hasAdv, hasCnd
 
    hasadv = this%gweModel%inadv > 0
    hascnd = this%gweModel%incnd > 0
 
    if (hasadv) then
      if (this%iIfaceAdvScheme == adv_scheme_tvd .or. &
          this%iIfaceAdvScheme == adv_scheme_utvd) then
        this%exg_stencil_depth = 2
        if (this%gweModel%adv%iadvwt == adv_scheme_tvd .or. &
            this%gweModel%adv%iadvwt == adv_scheme_utvd) then
          this%int_stencil_depth = 2
        end if
      end if
    end if
 
    if (hascnd) then
      if (this%iIfaceXt3d > 0) then
        this%exg_stencil_depth = 2
        if (this%gweModel%cnd%ixt3d > 0) then
          this%int_stencil_depth = 2
        end if
      end if
    end if
 

validateconnection()

subroutine gwegweconnectionmodule::validateconnection ( class(gwegweconnectiontype)  this )

private

Parameters

this	this connection

Definition at line 332 of file GweGweConnection.f90.

    use simvariablesmodule, only: errmsg
    use simmodule, only: count_errors, store_error
    class(GweGweConnectionType) :: this !< this connection
 
    ! base validation, the spatial/geometry part
    call this%SpatialModelConnectionType%validateConnection()
 
    ! we cannot validate this (yet) in parallel mode
    if (.not. this%gweExchange%v_model1%is_local) return
    if (.not. this%gweExchange%v_model2%is_local) return
 
    ! check specific cross-interface options/values that should be the same
    call this%validateGweExchange()
 
    ! GWE related matters
    if ((this%gweExchange%gwemodel1%inadv > 0 .and. &
         this%gweExchange%gwemodel2%inadv == 0) .or. &
        (this%gweExchange%gwemodel2%inadv > 0 .and. &
         this%gweExchange%gwemodel1%inadv == 0)) then
      write (errmsg, '(1x,a,a,a)') 'Cannot connect GWE models in exchange ', &
        trim(this%gweExchange%name), ' because one model is configured with ADV &
        &and the other one is not'
      call store_error(errmsg)
    end if
 
    if ((this%gweExchange%gwemodel1%incnd > 0 .and. &
         this%gweExchange%gwemodel2%incnd == 0) .or. &
        (this%gweExchange%gwemodel2%incnd > 0 .and. &
         this%gweExchange%gwemodel1%incnd == 0)) then
      write (errmsg, '(1x,a,a,a)') 'Cannot connect GWE models in exchange ', &
        trim(this%gweExchange%name), ' because one model is configured with CND &
        &and the other one is not'
      call store_error(errmsg)
    end if
 
    ! 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:

validategweexchange()

subroutine gwegweconnectionmodule::validategweexchange ( class(gwegweconnectiontype)  this )

private

Stops with error message on config mismatch

Parameters

this	this connection

Definition at line 531 of file GweGweConnection.f90.

    ! modules
    use simvariablesmodule, only: errmsg
    use simmodule, only: store_error
    use gweestmodule, only: gweesttype
 
    ! dummy
    class(GweGweConnectionType) :: this !< this connection
 
    ! local
    class(GweExchangeType), pointer :: gweEx
    class(*), pointer :: modelPtr
    class(GweModelType), pointer :: gweModel1
    class(GweModelType), pointer :: gweModel2
    type(GweEstType), pointer :: est1, est2
    logical(LGP) :: compatible
 
    gweex => this%gweExchange
 
    ! we cannot validate the remainder (yet) in parallel mode
    if (.not. gweex%v_model1%is_local) return
    if (.not. gweex%v_model2%is_local) return
 
    modelptr => this%gweExchange%model1
    gwemodel1 => castasgwemodel(modelptr)
    modelptr => this%gweExchange%model2
    gwemodel2 => castasgwemodel(modelptr)
 
    ! check that EST package usage is the same on both side of the interface
    if ((gwemodel1%inest > 0 .and. gwemodel2%inest == 0) .or. &
        (gwemodel1%inest == 0 .and. gwemodel2%inest > 0)) then
      write (errmsg, '(2a)') 'Energy Storage and Transfer package should '// &
        'be enabled/disabled simultaneously in models connected with the '// &
        'interface model for exchange ', &
        trim(gweex%name)
      call store_error(errmsg)
    end if
 
    ! conduction options need to be the same in both model
    if ((gwemodel1%cnd%ixt3d > 0 .and. gwemodel2%cnd%ixt3d == 0) .or. &
        (gwemodel1%cnd%ixt3d == 0 .and. gwemodel2%cnd%ixt3d > 0)) then
      write (errmsg, '(2a)') 'Use of XT3D to calculate conduction should '// &
        'be the same in both models, either both use XT3D or neither for '// &
        ' exchange '//trim(gweex%name)
      call store_error(errmsg)
    end if
 
    ! check compatibility of Energy Storage and Transfer (EST) package
    compatible = .true.
    est1 => gwemodel1%est
    est2 => gwemodel2%est
    if (est1%rhow /= est2%rhow) compatible = .false.
    if (est1%cpw /= est2%cpw) compatible = .false.
    ! if (est1%nrhospecies /= est2%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)') 'Energy storage and transfer (EST) packages ', &
        'in model '//trim(gweex%model1%name), ' and ', &
        trim(gweex%model2%name), &
        ' should be equivalent to construct an '// &
        ' interface model for exchange ', &
        trim(gweex%name)
      call store_error(errmsg)
    end if
 

Here is the call graph for this function: