Add a new message to MNW2 for case when well is deactivated

[langevin-usgs]

Richard,

Good suggestion. To implement this, I inserted the following 6 lines of new code into Subroutine SMNW2SEEP. Attached is an updated version of the file GWF2MNW27.f with this change in it. 

c write message that inode was deactivated this time step [LFK; 3/9/18]
if(mnwprnt.gt.0) then
write(iout,212) ND,WELLID(iw)
212 FORMAT ('MNW2: Node no. ',I3,' of Multi-Node Well ', A20)
write(iout,*) ' deactivated this time step because cell is dry'
end if

Chris,
Please incorporate this change in the next update/release of Modflow2005. It should not affect any calculations--only write an additional information line if conditions warrant. So no rush in releasing it. Please double check that I made the change in the latest version of the gwf2mnw27.f file in mf2005.

Let me know if any of you have any questions about this.

Lenny

On Fri, Mar 9, 2018 at 4:17 PM, Winston, Richard [email protected] wrote:
At line 4060 in gwf2mnw27.f, MODFLOW can deactivate a node in a multinode well without printing a warning message. I suggest that a warning message be printed.

For context, here are lines 4040 to 4075

        if(IBOUND(ic,ir,il).ne.0) then
          Bottom = BOTM(ic,ir,LBOTM(il))

c
c only allow seepage or flow into a cell if the head in the cell is greater than the bottom of the cell
c
if(hnew(ic,ir,il).gt.bottom) then
c Second time through loop, now that we have a guess for hwell, check to see if there is a seepage face
c If so, use the bottom (instead of the low hwell) as the gradient driver in the qact calc
c Set the hwell_by_node to the bottom to flag this as a seepage face cell
c Sum the seepage (there may be more than one node's worth
if( kSeep.gt.1 .and. hwell.lt.Bottom )then
MNWNOD(4,INODE) = (bottom - hnew(ic,ir,il))
& * MNWNOD(14,INODE)
MNWNOD(15,INODE) = bottom
Qseep = Qseep + MNWNOD(4,INODE)
else
csum = csum + MNWNOD(14,INODE)
chsum = chsum + MNWNOD(14,INODE)hnew(ic,ir,il)
endif
else
MNWNOD(4,INODE) = 0.0D0 ! <- MNW node is deactivated here on line 4060
MNWNOD(14,INODE) = 0.0D0
MNWNOD(15,INODE) = hdry
end if
else
MNWNOD(4,INODE) = 0.0D0 ! <- A warning is printed here if the node is deactivated.
MNWNOD(15,INODE) = hdry
c write message that inode was deactivated this time step
if(kseep.gt.1.and.mnwprnt.gt.0) then
write(iout,210) ND,WELLID(iw)
C-LFK rev. format 11/2012
c 210 FORMAT ('Node no. ',I3,' of Multi-Node Well ', A20)
210 FORMAT ('MNW2: Node no. ',I3,' of Multi-Node Well ', A20)
write(iout,) ' deactivated this time step because IBOUND=0'
end if
end if

--