18 character(len=LENFTYPE) ::
ftype =
'ESL'
19 character(len=16) ::
text =
' ESL'
46 subroutine esl_create(packobj, id, ibcnum, inunit, iout, namemodel, pakname, &
49 class(
bndtype),
pointer :: packobj
50 integer(I4B),
intent(in) :: id
51 integer(I4B),
intent(in) :: ibcnum
52 integer(I4B),
intent(in) :: inunit
53 integer(I4B),
intent(in) :: iout
54 character(len=*),
intent(in) :: namemodel
55 character(len=*),
intent(in) :: pakname
65 call packobj%set_names(ibcnum, namemodel, pakname,
ftype)
69 call eslobj%allocate_scalars()
72 call packobj%pack_initialize()
74 packobj%inunit = inunit
77 packobj%ibcnum = ibcnum
84 eslobj%gwecommon => gwecommon
99 call this%BndType%bnd_da()
116 call this%BndType%allocate_scalars()
136 integer(I4B) :: i, node
140 if (this%nbound == 0)
return
143 do i = 1, this%nbound
144 node = this%nodelist(i)
146 if (this%ibound(node) <= 0)
then
162 subroutine esl_fc(this, rhs, ia, idxglo, matrix_sln)
165 real(DP),
dimension(:),
intent(inout) :: rhs
166 integer(I4B),
dimension(:),
intent(in) :: ia
167 integer(I4B),
dimension(:),
intent(in) :: idxglo
170 integer(I4B) :: i, n, ipos
173 if (this%imover == 1)
then
174 call this%pakmvrobj%fc()
178 do i = 1, this%nbound
180 rhs(n) = rhs(n) + this%rhs(i)
182 call matrix_sln%add_value_pos(idxglo(ipos), this%hcof(i))
186 if (this%imover == 1 .and. this%rhs(i) >
dzero)
then
187 call this%pakmvrobj%accumulate_qformvr(i, this%rhs(i))
205 this%listlabel = trim(this%filtyp)//
' NO.'
206 if (this%dis%ndim == 3)
then
207 write (this%listlabel,
'(a, a7)') trim(this%listlabel),
'LAYER'
208 write (this%listlabel,
'(a, a7)') trim(this%listlabel),
'ROW'
209 write (this%listlabel,
'(a, a7)') trim(this%listlabel),
'COL'
210 elseif (this%dis%ndim == 2)
then
211 write (this%listlabel,
'(a, a7)') trim(this%listlabel),
'LAYER'
212 write (this%listlabel,
'(a, a7)') trim(this%listlabel),
'CELL2D'
214 write (this%listlabel,
'(a, a7)') trim(this%listlabel),
'NODE'
216 write (this%listlabel,
'(a, a16)') trim(this%listlabel),
'STRESS RATE'
217 if (this%inamedbound == 1)
then
218 write (this%listlabel,
'(a, a16)') trim(this%listlabel),
'BOUNDARY NAME'
257 call this%obs%StoreObsType(
'esl', .true., indx)
262 call this%obs%StoreObsType(
'to-mvr', .true., indx)
279 integer(I4B) :: i, nlinks
282 nlinks = this%TsManager%boundtslinks%Count()
285 if (
associated(tslink))
then
286 if (tslink%JCol == 1)
then
287 tslink%Text =
'SMASSRATE'
This module contains block parser methods.
This module contains the base boundary package.
This module contains simulation constants.
real(dp), parameter dem1
real constant 1e-1
integer(i4b), parameter lenftype
maximum length of a package type (DIS, WEL, OC, etc.)
real(dp), parameter dzero
real constant zero
real(dp), parameter done
real constant 1
subroutine, public esl_create(packobj, id, ibcnum, inunit, iout, namemodel, pakname, gwecommon)
Create an energy source loading package.
subroutine esl_da(this)
Deallocate memory.
subroutine esl_cf(this)
Formulate the HCOF and RHS terms.
character(len=lenftype) ftype
subroutine define_listlabel(this)
Define list labels.
subroutine esl_df_obs(this)
Define observations.
subroutine esl_fc(this, rhs, ia, idxglo, matrix_sln)
Add matrix terms related to specified energy source loading.
logical function esl_obs_supported(this)
Support function for specified energy source loading observations.
subroutine esl_rp_ts(this)
Procedure related to time series.
subroutine esl_allocate_scalars(this)
Allocate scalars.
This module defines variable data types.
This module contains the derived type ObsType.
subroutine, public defaultobsidprocessor(obsrv, dis, inunitobs, iout)
@ brief Process IDstring provided for each observation
type(timeserieslinktype) function, pointer, public gettimeserieslinkfromlist(list, indx)
Get time series link from a list.