This module contains the WEL package methods. More...

Data Types
type	weltype

Functions/Subroutines
subroutine, public	wel_create (packobj, id, ibcnum, inunit, iout, namemodel, pakname, mempath)
	@ brief Create a new package object More...

subroutine	wel_da (this)
	@ brief Deallocate package memory More...

subroutine	wel_allocate_scalars (this)
	@ brief Allocate scalars More...

subroutine	wel_allocate_arrays (this, nodelist, auxvar)
	@ brief Allocate arrays More...

subroutine	wel_options (this)
	@ brief Source additional options for package More...

subroutine	log_wel_options (this, found)
	@ brief Log WEL specific package options More...

subroutine	wel_rp (this)
	@ brief WEL read and prepare More...

subroutine	wel_cf (this)
	@ brief Formulate the package hcof and rhs terms. More...

subroutine	wel_fc (this, rhs, ia, idxglo, matrix_sln)
	@ brief Copy hcof and rhs terms into solution. More...

subroutine	wel_fn (this, rhs, ia, idxglo, matrix_sln)
	@ brief Add Newton-Raphson terms for package into solution. More...

subroutine	wel_afr_csv_init (this, fname)
	Initialize the auto flow reduce csv output file. More...

subroutine	wel_afr_csv_write (this)
	Write out auto flow reductions only when & where they occur. More...

subroutine	define_listlabel (this)
	@ brief Define the list label for the package More...

logical function	wel_obs_supported (this)
	Determine if observations are supported. More...

subroutine	wel_df_obs (this)
	Define the observation types available in the package. More...

subroutine	wel_bd_obs (this)
	Save observations for the package. More...

real(dp) function	q_mult (this, row)

real(dp) function	wel_bound_value (this, col, row)
	@ brief Return a bound value More...

Variables
character(len=lenftype)	ftype = 'WEL'
	package ftype More...

character(len=16)	text = ' WEL'
	package flow text string More...

Detailed Description

This module contains the overridden methods for the standard WEL package. Several methods need to be overridden because of the AUTO_FLOW_REDUCE option. Overridden methods include:

bnd_cf (AUTO_FLOW_REDUCE)
bnd_fc (AUTO_FLOW_REDUCE)
bnd_fn (AUTO_FLOW_REDUCE Newton-Raphson terms)
bnd_ot_package_flows (write AUTO_FLOW_REDUCE terms to csv file)
bnd_da (deallocate AUTO_FLOW_REDUCE variables)
bnd_bd_obs (wel-reduction observation added)

Function/Subroutine Documentation

◆ define_listlabel()

subroutine welmodule::define_listlabel ( class(weltype), intent(inout) this )

Method defined the list label for the WEL package. The list label is the heading that is written to iout when PRINT_INPUT option is used.

Parameters

[in,out] this WelType object

Definition at line 510 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType), intent(inout) :: this !< WelType object
     !
     ! -- create the header list label
     this%listlabel = trim(this%filtyp)//' NO.'
     if (this%dis%ndim == 3) then
       write (this%listlabel, '(a, a7)') trim(this%listlabel), 'LAYER'
       write (this%listlabel, '(a, a7)') trim(this%listlabel), 'ROW'
       write (this%listlabel, '(a, a7)') trim(this%listlabel), 'COL'
     elseif (this%dis%ndim == 2) then
       write (this%listlabel, '(a, a7)') trim(this%listlabel), 'LAYER'
       write (this%listlabel, '(a, a7)') trim(this%listlabel), 'CELL2D'
     else
       write (this%listlabel, '(a, a7)') trim(this%listlabel), 'NODE'
     end if
     write (this%listlabel, '(a, a16)') trim(this%listlabel), 'STRESS RATE'
     if (this%inamedbound == 1) then
       write (this%listlabel, '(a, a16)') trim(this%listlabel), 'BOUNDARY NAME'
     end if
     !
     ! -- return
     return

◆ log_wel_options()

subroutine welmodule::log_wel_options	(	class(weltype), intent(inout)	this,
		type(gwfwelparamfoundtype), intent(in)	found
	)

Parameters

[in,out] this BndExtType object

Definition at line 249 of file gwf-wel.f90.

     ! -- modules
     use gwfwelinputmodule, only: gwfwelparamfoundtype
     ! -- dummy variables
     class(WelType), intent(inout) :: this !< BndExtType object
     type(GwfWelParamFoundType), intent(in) :: found
     ! -- local variables
     ! -- format
     character(len=*), parameter :: fmtflowred = &
       &"(4x, 'AUTOMATIC FLOW REDUCTION OF WELLS IMPLEMENTED.')"
     character(len=*), parameter :: fmtflowredv = &
       &"(4x, 'AUTOMATIC FLOW REDUCTION FRACTION (',g15.7,').')"
     !
     ! -- log found options
     write (this%iout, '(/1x,a)') 'PROCESSING '//trim(adjustl(this%text)) &
       //' OPTIONS'
     !
     if (found%flowred) then
       if (this%iflowred > 0) &
         write (this%iout, fmtflowred)
       write (this%iout, fmtflowredv) this%flowred
     end if
     !
     if (found%afrcsvfile) then
       ! -- currently no-op
     end if
     !
     if (found%mover) then
       write (this%iout, '(4x,A)') 'MOVER OPTION ENABLED'
     end if
     !
     ! -- close logging block
     write (this%iout, '(1x,a)') &
       'END OF '//trim(adjustl(this%text))//' OPTIONS'
     !
     ! -- return
     return

◆ q_mult()

real(dp) function welmodule::q_mult	(	class(weltype), intent(inout)	this,
		integer(i4b), intent(in)	row
	)

private

Parameters

[in,out] this BndExtType object

Definition at line 644 of file gwf-wel.f90.

     ! -- modules
     use constantsmodule, only: dzero
     ! -- dummy variables
     class(WelType), intent(inout) :: this !< BndExtType object
     integer(I4B), intent(in) :: row
     ! -- result
     real(DP) :: q
     !
     if (this%iauxmultcol > 0) then
       q = this%q(row) * this%auxvar(this%iauxmultcol, row)
     else
       q = this%q(row)
     end if
     !
     ! -- return
     return

◆ wel_afr_csv_init()

subroutine welmodule::wel_afr_csv_init	(	class(weltype), intent(inout)	this,
		character(len=*), intent(in)	fname
	)

private

Parameters

[in,out] this WelType object

Definition at line 456 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType), intent(inout) :: this !< WelType object
     character(len=*), intent(in) :: fname
     ! -- format
     character(len=*), parameter :: fmtafrcsv = &
       "(4x, 'AUTO FLOW REDUCE INFORMATION WILL BE SAVED TO FILE: ', a, /4x, &
       &'OPENED ON UNIT: ', I0)"
  
     this%ioutafrcsv = getunit()
     call openfile(this%ioutafrcsv, this%iout, fname, 'CSV', &
                   filstat_opt='REPLACE')
     write (this%iout, fmtafrcsv) trim(adjustl(fname)), &
       this%ioutafrcsv
     write (this%ioutafrcsv, '(a)') &
       'time,period,step,boundnumber,cellnumber,rate-requested,&
       &rate-actual,wel-reduction'
     return

Here is the call graph for this function:

◆ wel_afr_csv_write()

subroutine welmodule::wel_afr_csv_write ( class(weltype), intent(inout) this )

private

Parameters

[in,out] this WelType object

Definition at line 477 of file gwf-wel.f90.

     ! -- modules
     use tdismodule, only: totim, kstp, kper
     ! -- dummy variables
     class(WelType), intent(inout) :: this !< WelType object
     ! -- local
     integer(I4B) :: i
     integer(I4B) :: nodereduced
     integer(I4B) :: nodeuser
     real(DP) :: v
     ! -- format
     do i = 1, this%nbound
       nodereduced = this%nodelist(i)
       !
       ! -- test if node is constant or inactive
       if (this%ibound(nodereduced) <= 0) then
         cycle
       end if
       v = this%q_mult(i) + this%rhs(i)
       if (v < dzero) then
         nodeuser = this%dis%get_nodeuser(nodereduced)
         write (this%ioutafrcsv, '(*(G0,:,","))') &
           totim, kper, kstp, i, nodeuser, this%q_mult(i), this%simvals(i), v
       end if
     end do

◆ wel_allocate_arrays()

subroutine welmodule::wel_allocate_arrays	(	class(weltype)	this,
		integer(i4b), dimension(:), optional, pointer, contiguous	nodelist,
		real(dp), dimension(:, :), optional, pointer, contiguous	auxvar
	)

Allocate and initialize arrays for the WEL package

Definition at line 169 of file gwf-wel.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate, mem_setptr, mem_checkin
     ! -- dummy
     class(WelType) :: this
     integer(I4B), dimension(:), pointer, contiguous, optional :: nodelist
     real(DP), dimension(:, :), pointer, contiguous, optional :: auxvar
     ! -- local
     !
     ! -- call BndExtType allocate scalars
     call this%BndExtType%allocate_arrays(nodelist, auxvar)
     !
     ! -- set constant head array input context pointer
     call mem_setptr(this%q, 'Q', this%input_mempath)
     !
     ! -- checkin constant head array input context pointer
     call mem_checkin(this%q, 'Q', this%memoryPath, &
                      'Q', this%input_mempath)
     !
     ! -- return
     return

◆ wel_allocate_scalars()

subroutine welmodule::wel_allocate_scalars ( class(weltype) this )

Allocate and initialize scalars for the WEL package. The base model allocate scalars method is also called.

Parameters

this	WelType object

Definition at line 141 of file gwf-wel.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate
     ! -- dummy variables
     class(WelType) :: this !< WelType object
     !
     ! -- call base type allocate scalars
     call this%BndExtType%allocate_scalars()
     !
     ! -- allocate the object and assign values to object variables
     call mem_allocate(this%iflowred, 'IFLOWRED', this%memoryPath)
     call mem_allocate(this%flowred, 'FLOWRED', this%memoryPath)
     call mem_allocate(this%ioutafrcsv, 'IOUTAFRCSV', this%memoryPath)
     !
     ! -- Set values
     this%iflowred = 0
     this%ioutafrcsv = 0
     this%flowred = dzero
     !
     ! -- return
     return

◆ wel_bd_obs()

subroutine welmodule::wel_bd_obs ( class(weltype) this )

private

Method to save simulated values for the WEL package.

Parameters

this	WelType object

Definition at line 590 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType) :: this !< WelType object
     ! -- local variables
     integer(I4B) :: i
     integer(I4B) :: n
     integer(I4B) :: jj
     real(DP) :: v
     type(ObserveType), pointer :: obsrv => null()
     !
     ! -- clear the observations
     call this%obs%obs_bd_clear()
     !
     ! -- Save simulated values for all of package's observations.
     do i = 1, this%obs%npakobs
       obsrv => this%obs%pakobs(i)%obsrv
       if (obsrv%BndFound) then
         do n = 1, obsrv%indxbnds_count
           v = dnodata
           jj = obsrv%indxbnds(n)
           select case (obsrv%ObsTypeId)
           case ('TO-MVR')
             if (this%imover == 1) then
               v = this%pakmvrobj%get_qtomvr(jj)
               if (v > dzero) then
                 v = -v
               end if
             end if
           case ('WEL')
             v = this%simvals(jj)
           case ('WEL-REDUCTION')
             if (this%iflowred > 0) then
               v = this%q_mult(jj) + this%rhs(jj)
             end if
           case default
             errmsg = 'Unrecognized observation type: '//trim(obsrv%ObsTypeId)
             call store_error(errmsg)
           end select
           call this%obs%SaveOneSimval(obsrv, v)
         end do
       else
         call this%obs%SaveOneSimval(obsrv, dnodata)
       end if
     end do
     !
     ! -- Write the auto flow reduce csv file entries for this step
     if (this%ioutafrcsv > 0) then
       call this%wel_afr_csv_write()
     end if
     !
     ! -- return
     return

Here is the call graph for this function:

◆ wel_bound_value()

real(dp) function welmodule::wel_bound_value	(	class(weltype), intent(inout)	this,
		integer(i4b), intent(in)	col,
		integer(i4b), intent(in)	row
	)

Return a bound value associated with an ncolbnd index and row.

Parameters

[in,out] this BndExtType object

Definition at line 669 of file gwf-wel.f90.

     ! -- modules
     use constantsmodule, only: dzero
     ! -- dummy variables
     class(WelType), intent(inout) :: this !< BndExtType object
     integer(I4B), intent(in) :: col
     integer(I4B), intent(in) :: row
     ! -- result
     real(DP) :: bndval
     !
     select case (col)
     case (1)
       bndval = this%q_mult(row)
     case default
       errmsg = 'Programming error. WEL bound value requested column '&
                &'outside range of ncolbnd (1).'
       call store_error(errmsg)
       call store_error_filename(this%input_fname)
     end select
     !
     ! -- return
     return

Here is the call graph for this function:

◆ wel_cf()

subroutine welmodule::wel_cf ( class(weltype) this )

Formulate the hcof and rhs terms for the WEL package that will be added to the coefficient matrix and right-hand side vector.

Parameters

this	WelType object

Definition at line 317 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType) :: this !< WelType object
     ! -- local variables
     integer(I4B) :: i, node, ict
     real(DP) :: qmult
     real(DP) :: q
     real(DP) :: tp
     real(DP) :: bt
     real(DP) :: thick
     !
     ! -- Return if no wells
     if (this%nbound == 0) return
     !
     ! -- Calculate hcof and rhs for each well entry
     do i = 1, this%nbound
       node = this%nodelist(i)
       this%hcof(i) = dzero
       if (this%ibound(node) <= 0) then
         this%rhs(i) = dzero
         cycle
       end if
       q = this%q_mult(i)
       if (this%iflowred /= 0 .and. q < dzero) then
         ict = this%icelltype(node)
         if (ict /= 0) then
           tp = this%dis%top(node)
           bt = this%dis%bot(node)
           thick = tp - bt
           tp = bt + this%flowred * thick
           qmult = sqsaturation(tp, bt, this%xnew(node))
           q = q * qmult
         end if
       end if
       this%rhs(i) = -q
     end do
     !
     return

Here is the call graph for this function:

◆ wel_create()

subroutine, public welmodule::wel_create	(	class(bndtype), pointer	packobj,
		integer(i4b), intent(in)	id,
		integer(i4b), intent(in)	ibcnum,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout,
		character(len=*), intent(in)	namemodel,
		character(len=*), intent(in)	pakname,
		character(len=*), intent(in)	mempath
	)

Create a new WEL Package object

Parameters

	packobj	pointer to default package type
[in]	id	package id
[in]	ibcnum	boundary condition number
[in]	inunit	unit number of WEL package input file
[in]	iout	unit number of model listing file
[in]	namemodel	model name
[in]	pakname	package name
[in]	mempath	input mempath

Definition at line 73 of file gwf-wel.f90.

     ! -- dummy variables
     class(BndType), pointer :: packobj !< pointer to default package type
     integer(I4B), intent(in) :: id !< package id
     integer(I4B), intent(in) :: ibcnum !< boundary condition number
     integer(I4B), intent(in) :: inunit !< unit number of WEL package input file
     integer(I4B), intent(in) :: iout !< unit number of model listing file
     character(len=*), intent(in) :: namemodel !< model name
     character(len=*), intent(in) :: pakname !< package name
     character(len=*), intent(in) :: mempath !< input mempath
     ! -- local variables
     type(WelType), pointer :: welobj
     !
     ! -- allocate the object and assign values to object variables
     allocate (welobj)
     packobj => welobj
     !
     ! -- create name and memory path
     call packobj%set_names(ibcnum, namemodel, pakname, ftype, mempath)
     packobj%text = text
     !
     ! -- allocate scalars
     call welobj%allocate_scalars()
     !
     ! -- initialize package
     call packobj%pack_initialize()
  
     packobj%inunit = inunit
     packobj%iout = iout
     packobj%id = id
     packobj%ibcnum = ibcnum
     packobj%ictMemPath = create_mem_path(namemodel, 'NPF')
     !
     ! -- return
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ wel_da()

subroutine welmodule::wel_da ( class(weltype) this )

private

Deallocate WEL package scalars and arrays.

Parameters

this	WelType object

Definition at line 116 of file gwf-wel.f90.

     ! -- modules
     use memorymanagermodule, only: mem_deallocate
     ! -- dummy variables
     class(WelType) :: this !< WelType object
     !
     ! -- Deallocate parent package
     call this%BndExtType%bnd_da()
     !
     ! -- scalars
     call mem_deallocate(this%iflowred)
     call mem_deallocate(this%flowred)
     call mem_deallocate(this%ioutafrcsv)
     call mem_deallocate(this%q, 'Q', this%memoryPath)
     !
     ! -- return
     return

◆ wel_df_obs()

subroutine welmodule::wel_df_obs ( class(weltype) this )

private

Method to define the observation types available in the WEL package.

Parameters

this	WelType object

Definition at line 561 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType) :: this !< WelType object
     ! -- local variables
     integer(I4B) :: indx
     !
     ! -- initialize observations
     call this%obs%StoreObsType('wel', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => defaultobsidprocessor
     !
     ! -- Store obs type and assign procedure pointer
     !    for to-mvr observation type.
     call this%obs%StoreObsType('to-mvr', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => defaultobsidprocessor
     !
     ! -- Store obs type and assign procedure pointer
     !    for wel-reduction observation type.
     call this%obs%StoreObsType('wel-reduction', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => defaultobsidprocessor
     !
     ! -- return
     return

Here is the call graph for this function:

◆ wel_fc()

subroutine welmodule::wel_fc	(	class(weltype)	this,
		real(dp), dimension(:), intent(inout)	rhs,
		integer(i4b), dimension(:), intent(in)	ia,
		integer(i4b), dimension(:), intent(in)	idxglo,
		class(matrixbasetype), pointer	matrix_sln
	)

private

Add the hcof and rhs terms for the WEL package to the coefficient matrix and right-hand side vector.

Parameters

	this	WelType object
[in,out]	rhs	right-hand side vector for model
[in]	ia	solution CRS row pointers
[in]	idxglo	mapping vector for model (local) to solution (global)
	matrix_sln	solution coefficient matrix

Definition at line 363 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType) :: this !< WelType object
     real(DP), dimension(:), intent(inout) :: rhs !< right-hand side vector for model
     integer(I4B), dimension(:), intent(in) :: ia !< solution CRS row pointers
     integer(I4B), dimension(:), intent(in) :: idxglo !< mapping vector for model (local) to solution (global)
     class(MatrixBaseType), pointer :: matrix_sln !< solution coefficient matrix
     ! -- local variables
     integer(I4B) :: i
     integer(I4B) :: n
     integer(I4B) :: ipos
     !
     ! -- pakmvrobj fc
     if (this%imover == 1) then
       call this%pakmvrobj%fc()
     end if
     !
     ! -- Copy package rhs and hcof into solution rhs and amat
     do i = 1, this%nbound
       n = this%nodelist(i)
       rhs(n) = rhs(n) + this%rhs(i)
       ipos = ia(n)
       call matrix_sln%add_value_pos(idxglo(ipos), this%hcof(i))
       !
       ! -- If mover is active and this well is discharging,
       !    store available water (as positive value).
       if (this%imover == 1 .and. this%rhs(i) > dzero) then
         call this%pakmvrobj%accumulate_qformvr(i, this%rhs(i))
       end if
     end do
     !
     ! -- return
     return

◆ wel_fn()

subroutine welmodule::wel_fn	(	class(weltype)	this,
		real(dp), dimension(:), intent(inout)	rhs,
		integer(i4b), dimension(:), intent(in)	ia,
		integer(i4b), dimension(:), intent(in)	idxglo,
		class(matrixbasetype), pointer	matrix_sln
	)

private

Calculate and add the Newton-Raphson terms for the WEL package to the coefficient matrix and right-hand side vector.

Parameters

	this	WelType object
[in,out]	rhs	right-hand side vector for model
[in]	ia	solution CRS row pointers
[in]	idxglo	mapping vector for model (local) to solution (global)
	matrix_sln	solution coefficient matrix

Definition at line 404 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType) :: this !< WelType object
     real(DP), dimension(:), intent(inout) :: rhs !< right-hand side vector for model
     integer(I4B), dimension(:), intent(in) :: ia !< solution CRS row pointers
     integer(I4B), dimension(:), intent(in) :: idxglo !< mapping vector for model (local) to solution (global)
     class(MatrixBaseType), pointer :: matrix_sln !< solution coefficient matrix
     ! -- local variables
     integer(I4B) :: i
     integer(I4B) :: node
     integer(I4B) :: ipos
     integer(I4B) :: ict
     real(DP) :: drterm
     real(DP) :: q
     real(DP) :: tp
     real(DP) :: bt
     real(DP) :: thick
     !
     ! -- Copy package rhs and hcof into solution rhs and amat
     do i = 1, this%nbound
       node = this%nodelist(i)
       !
       ! -- test if node is constant or inactive
       if (this%ibound(node) <= 0) then
         cycle
       end if
       !
       ! -- well rate is possibly head dependent
       ict = this%icelltype(node)
       if (this%iflowred /= 0 .and. ict /= 0) then
         ipos = ia(node)
         q = -this%rhs(i)
         if (q < dzero) then
           ! -- calculate derivative for well
           tp = this%dis%top(node)
           bt = this%dis%bot(node)
           thick = tp - bt
           tp = bt + this%flowred * thick
           drterm = sqsaturationderivative(tp, bt, this%xnew(node))
           drterm = drterm * this%q_mult(i)
           !--fill amat and rhs with newton-raphson terms
           call matrix_sln%add_value_pos(idxglo(ipos), drterm)
           rhs(node) = rhs(node) + drterm * this%xnew(node)
         end if
       end if
     end do
     !
     ! -- return
     return

Here is the call graph for this function:

◆ wel_obs_supported()

logical function welmodule::wel_obs_supported ( class(weltype) this )

private

Function to determine if observations are supported by the WEL package. Observations are supported by the WEL package.

Returns: wel_obs_supported boolean indicating if observations are supported

Parameters

this	WelType object

Definition at line 545 of file gwf-wel.f90.

     ! -- dummy variables
     class(WelType) :: this !< WelType object
     !
     ! -- set boolean
     wel_obs_supported = .true.
     !
     ! -- return
     return

◆ wel_options()

subroutine welmodule::wel_options ( class(weltype), intent(inout) this )

Source additional options for WEL package.

Parameters

[in,out] this WelType object

Definition at line 197 of file gwf-wel.f90.

     ! -- modules
     use inputoutputmodule, only: urword
     use memorymanagerextmodule, only: mem_set_value
     use gwfwelinputmodule, only: gwfwelparamfoundtype
     ! -- dummy variables
     class(WelType), intent(inout) :: this !< WelType object
     ! -- local variables
     character(len=LINELENGTH) :: fname
     type(GwfWelParamFoundType) :: found
     ! -- formats
     character(len=*), parameter :: fmtflowred = &
       &"(4x, 'AUTOMATIC FLOW REDUCTION OF WELLS IMPLEMENTED.')"
     character(len=*), parameter :: fmtflowredv = &
       &"(4x, 'AUTOMATIC FLOW REDUCTION FRACTION (',g15.7,').')"
     !
     ! -- source base BndExtType options
     call this%BndExtType%source_options()
     !
     ! -- source well options from input context
     call mem_set_value(this%flowred, 'FLOWRED', this%input_mempath, found%flowred)
     call mem_set_value(fname, 'AFRCSVFILE', this%input_mempath, found%afrcsvfile)
     call mem_set_value(this%imover, 'MOVER', this%input_mempath, found%mover)
     !
     if (found%flowred) then
       !
       this%iflowred = 1
       !
       if (this%flowred <= dzero) then
         this%flowred = dem1
       else if (this%flowred > done) then
         this%flowred = done
       end if
     end if
     !
     if (found%afrcsvfile) then
       call this%wel_afr_csv_init(fname)
     end if
     !
     if (found%mover) then
       this%imover = 1
     end if
     !
     ! -- log WEL specific options
     call this%log_wel_options(found)
     !
     ! -- return
     return

Here is the call graph for this function:

◆ wel_rp()

subroutine welmodule::wel_rp ( class(weltype), intent(inout) this )

Definition at line 291 of file gwf-wel.f90.

     use tdismodule, only: kper
     ! -- dummy
     class(WelType), intent(inout) :: this
     ! -- local
     !
     if (this%iper /= kper) return
     !
     ! -- Call the parent class read and prepare
     call this%BndExtType%bnd_rp()
     !
     ! -- Write the list to iout if requested
     if (this%iprpak /= 0) then
       call this%write_list()
     end if
     !
     ! -- return
     return

Variable Documentation

◆ ftype

character(len=lenftype) welmodule::ftype = 'WEL'

private

Definition at line 36 of file gwf-wel.f90.

36 character(len=LENFTYPE) :: ftype = 'WEL' !< package ftype

◆ text

character(len=16) welmodule::text = ' WEL'

private

Definition at line 37 of file gwf-wel.f90.

37 character(len=16) :: text = ' WEL' !< package flow text string

Data Types

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ define_listlabel()

◆ log_wel_options()

◆ q_mult()

◆ wel_afr_csv_init()

◆ wel_afr_csv_write()

◆ wel_allocate_arrays()

◆ wel_allocate_scalars()

◆ wel_bd_obs()

◆ wel_bound_value()

◆ wel_cf()

◆ wel_create()

◆ wel_da()

◆ wel_df_obs()

◆ wel_fc()

◆ wel_fn()

◆ wel_obs_supported()

◆ wel_options()

◆ wel_rp()

Variable Documentation

◆ ftype

◆ text