gwfinterfacemodelmodule Module Reference

Data Types
type	gwfinterfacemodeltype
	The GWF Interface Model is a utility to calculate the solution's exchange coefficients from the interface between a GWF model and its GWF neighbors. The interface model itself will not be part of the solution, it is not being solved. Patching (a part of the) discretizations of two GWF models in a general way, e.g. DIS+DIS with refinement, requires the resulting. More...

Functions/Subroutines
subroutine	gwfifm_cr (this, name, iout, gridConn)
	set up the interface model, analogously to what More...
subroutine	gwfifm_df (this)
	Define, mostly DISU and the NPF package. More...
subroutine	gwfifm_ar (this)
	allocate and read the packages More...
subroutine	gwfifm_da (this)
	Clean up. More...
subroutine	setnpfoptions (this, npfOptions)
	Copy NPF options from the model owning the interface to the data structure. More...
subroutine	setbuydata (this, buyData)
	Sets the BUY input data from the models that make up this interface. We adopt everything from the owning model, but during validation it should be. More...

Function/Subroutine Documentation

gwfifm_ar()

subroutine gwfinterfacemodelmodule::gwfifm_ar ( class(gwfinterfacemodeltype)  this )

private

Parameters

this	the GWF interface model

Definition at line 119 of file GwfInterfaceModel.f90.

    class(GwfInterfaceModelType) :: this !< the GWF interface model
 
    call this%npf%npf_ar(this%ic, this%vsc, this%ibound, this%x)
    if (this%inbuy > 0) call this%buy%buy_ar(this%npf, this%ibound)
 

gwfifm_cr()

subroutine gwfinterfacemodelmodule::gwfifm_cr	(	class(gwfinterfacemodeltype)	this,
		character(len=*), intent(in)	name,
		integer(i4b), intent(in)	iout,
		class(gridconnectiontype), intent(in), pointer	gridConn
	)

Parameters

	this	the GWF interface model
[in]	name	the interface model's name
[in]	iout	the output unit
[in]	gridconn	the grid connection for creating a DISU

Definition at line 46 of file GwfInterfaceModel.f90.

    class(GwfInterfaceModelType) :: this !< the GWF interface model
    character(len=*), intent(in) :: name !< the interface model's name
    integer(I4B), intent(in) :: iout !< the output unit
    class(GridConnectionType), pointer, intent(in) :: gridConn !< the grid connection for creating a DISU
    ! local
    class(*), pointer :: modPtr
 
    this%memoryPath = create_mem_path(name)
    call this%allocate_scalars(name)
 
    this%iout = iout
 
    this%gridConnection => gridconn
    modptr => this%gridConnection%model
    this%owner => castasgwfmodel(modptr)
 
    this%innpf = huge(1_i4b)
    this%inewton = this%owner%inewton
    this%inewtonur = this%owner%inewtonur
 
    if (this%owner%inbuy > 0) then
      this%inbuy = huge(1_i4b)
    end if
 
    ! create discretization and packages
    call disu_cr(this%dis, this%name, '', -1, this%iout)
    call npf_cr(this%npf, this%name, '', -this%innpf, this%iout)
    call xt3d_cr(this%xt3d, this%name, -this%innpf, this%iout)
    call buy_cr(this%buy, this%name, this%inbuy, this%iout)
 

Here is the call graph for this function:

gwfifm_da()

subroutine gwfinterfacemodelmodule::gwfifm_da ( class(gwfinterfacemodeltype)  this )

private

Parameters

this	the GWF interface model

Definition at line 129 of file GwfInterfaceModel.f90.

    use memorymanagermodule, only: mem_deallocate
    class(GwfInterfaceModelType) :: this !< the GWF interface model
 
    ! -- Internal flow packages deallocate
    call this%dis%dis_da()
    call this%npf%npf_da()
    call this%xt3d%xt3d_da()
    call this%buy%buy_da()
    !
    ! -- Internal package objects
    deallocate (this%dis)
    deallocate (this%npf)
    deallocate (this%xt3d)
    !
    ! -- Scalars
    call mem_deallocate(this%inic)
    call mem_deallocate(this%inoc)
    call mem_deallocate(this%inobs)
    call mem_deallocate(this%innpf)
    call mem_deallocate(this%inbuy)
    call mem_deallocate(this%invsc)
    call mem_deallocate(this%insto)
    call mem_deallocate(this%incsub)
    call mem_deallocate(this%inmvr)
    call mem_deallocate(this%inhfb)
    call mem_deallocate(this%ingnc)
    call mem_deallocate(this%iss)
    call mem_deallocate(this%inewtonur)
    !
    ! -- NumericalModelType
    call this%NumericalModelType%model_da()
 

gwfifm_df()

subroutine gwfinterfacemodelmodule::gwfifm_df ( class(gwfinterfacemodeltype)  this )

private

Parameters

this	the GWF interface model

Definition at line 81 of file GwfInterfaceModel.f90.

    class(GwfInterfaceModelType) :: this !< the GWF interface model
    ! local
    type(GwfNpfOptionsType) :: npfOptions
    type(GwfBuyInputDataType) :: buyData
    class(*), pointer :: disPtr
 
    this%moffset = 0
 
    ! define DISU
    disptr => this%dis
    call this%gridConnection%getDiscretization(castasdisutype(disptr))
 
    ! define NPF package
    call npfoptions%construct()
    call this%setNpfOptions(npfoptions)
    call this%npf%npf_df(this%dis, this%xt3d, 0, 0, npfoptions)
    call npfoptions%destroy()
 
    ! define BUY package
    if (this%inbuy > 0) then
      call buydata%construct(this%owner%buy%nrhospecies)
      call this%setBuyData(buydata)
      call this%buy%buy_df(this%dis, buydata)
      call buydata%destruct()
    end if
 
    this%neq = this%dis%nodes
    this%nja = this%dis%nja
    this%ia => this%dis%con%ia
    this%ja => this%dis%con%ja
 
    call this%allocate_arrays()
 

Here is the call graph for this function:

setbuydata()

subroutine gwfinterfacemodelmodule::setbuydata ( class(gwfinterfacemodeltype)  this, type(gwfbuyinputdatatype)  buyData  )

private

Parameters

this	the interface model
buydata	the data for the buoyancy package

Definition at line 189 of file GwfInterfaceModel.f90.

    class(GwfInterfaceModelType) :: this !< the interface model
    type(GwfBuyInputDataType) :: buyData !< the data for the buoyancy package
    ! local
    integer(I4B) :: i
 
    buydata%denseref = this%owner%buy%denseref
    buydata%iform = this%owner%buy%iform
    buydata%nrhospecies = this%owner%buy%nrhospecies
 
    do i = 1, buydata%nrhospecies
      buydata%drhodc(i) = this%owner%buy%drhodc(i)
      buydata%crhoref(i) = this%owner%buy%crhoref(i)
      buydata%cmodelname(i) = this%owner%buy%cmodelname(i)
      buydata%cauxspeciesname(i) = this%owner%buy%cauxspeciesname(i)
    end do
 

setnpfoptions()

subroutine gwfinterfacemodelmodule::setnpfoptions	(	class(gwfinterfacemodeltype)	this,
		type(gwfnpfoptionstype)	npfOptions
	)

Parameters

this	the GWF interface model
npfoptions	the options data to be filled

Definition at line 167 of file GwfInterfaceModel.f90.

    class(GwfInterfaceModelType) :: this !< the GWF interface model
    type(GwfNpfOptionsType) :: npfOptions !< the options data to be filled
 
    ! for now, assuming full homogeneity, so just take
    ! the options from the owning model's npf package
    npfoptions%icellavg = this%owner%npf%icellavg
    npfoptions%ithickstrt = this%owner%npf%ithickstrt
    npfoptions%iperched = this%owner%npf%iperched
    npfoptions%ivarcv = this%owner%npf%ivarcv
    npfoptions%idewatcv = this%owner%npf%idewatcv
    npfoptions%irewet = this%owner%npf%irewet
    npfoptions%wetfct = this%owner%npf%wetfct
    npfoptions%iwetit = this%owner%npf%iwetit
    npfoptions%ihdwet = this%owner%npf%ihdwet