Comments (4)
emorway-usgs
commented on September 28, 2024
@DarkForceOne If you have a small example problem demonstrating the problem, that would be helpful as we recall why this code was added? I can see that the block of code you are referring to (shown below) undoes a constant concentration condition specified in BTN when NTMP is specified as something great than 1 in SSM (otherwise, this code won't be triggered).
C-------RESET ICBUND TO ABS(ICBUND) FOR CONSTANT CONCENTRATION CELLS
DO INDEX=1,NCOMP
DO KK=1,NLAY
DO II=1,NROW
DO JJ=1,NCOL
IF(ICBUND(JJ,II,KK,INDEX).LT.0) THEN
ICBUND(JJ,II,KK,INDEX)=ABS(ICBUND(JJ,II,KK,INDEX))
ENDIF
ENDDO
ENDDO
ENDDO
ENDDO
@vivekbedekar Do you recall why this code was added? It looks like if first showed up with a commit you made in 2015:
from mt3d-usgs.
DarkForceOne
commented on September 28, 2024
Hello,
In a project to test a GSFLOW-MT3D-USGS interface, I used a model with 192 stress periods and set a few cells to constant concentration condition with icbund=-1. This block of code was triggered when the simulation advanced to the second period, thus removing the constant concentration condition.
Best regards,
Eric Chiang
Simcore Software
… On Aug 24, 2020, at 5:16 PM, Eric Morway ***@***.***> wrote:
@DarkForceOne If you have a small example problem demonstrating the problem, that would be helpful as we recall why this code was added? I can see that the block of code you are referring to (shown below) undoes a constant concentration condition specified in BTN when NTMP is specified as something great than 1 in SSM (otherwise, this code won't be triggered).
C-------RESET ICBUND TO ABS(ICBUND) FOR CONSTANT CONCENTRATION CELLS
DO INDEX=1,NCOMP
DO KK=1,NLAY
DO II=1,NROW
DO JJ=1,NCOL
IF(ICBUND(JJ,II,KK,INDEX).LT.0) THEN
ICBUND(JJ,II,KK,INDEX)=ABS(ICBUND(JJ,II,KK,INDEX))
ENDIF
ENDDO
ENDDO
ENDDO
ENDDO
@vivekbedekar Do you recall why this code was added? It looks like if first showed up with a commit you made in 2015:
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub, or unsubscribe.
from mt3d-usgs.
vivekbedekar
commented on September 28, 2024
This code was specifically added to reset ICBUND at every stress period. In earlier versions, once a cell was defined ICBUND=-1 (prescribed/constant concentration boundary), it would remain ICBUND=-1 for the remainder of the simulation. There was no way to reset that. For instance, if a transient simulation needs the removal of a prescribed/constant concentration boundary, say due to removal of a source by dredging, it was not possible to do with the earlier versions.
This change to reset ICBUND at every stress period is also consistent with all other boundary conditions. If the user wants to continue using a boundary condition (including prescribed/constant concentration boundary), they need to actively do that.
Vivek
Sent from Outlook<http://aka.ms/weboutlook>
________________________________
From: DarkForceOne <[email protected]>
Sent: Tuesday, August 25, 2020 10:35 AM
To: MODFLOW-USGS/mt3d-usgs <[email protected]>
Cc: vivekbedekar <[email protected]>; Mention <[email protected]>
Subject: Re: [MODFLOW-USGS/mt3d-usgs] Constant concentration condition (ICBUND < 0) is not persistent (#82)
Hello,
In a project to test a GSFLOW-MT3D-USGS interface, I used a model with 192 stress periods and set a few cells to constant concentration condition with icbund=-1. This block of code was triggered when the simulation advanced to the second period, thus removing the constant concentration condition.
Best regards,
Eric Chiang
Simcore Software
On Aug 24, 2020, at 5:16 PM, Eric Morway ***@***.***> wrote:
@DarkForceOne If you have a small example problem demonstrating the problem, that would be helpful as we recall why this code was added? I can see that the block of code you are referring to (shown below) undoes a constant concentration condition specified in BTN when NTMP is specified as something great than 1 in SSM (otherwise, this code won't be triggered).
C-------RESET ICBUND TO ABS(ICBUND) FOR CONSTANT CONCENTRATION CELLS
DO INDEX=1,NCOMP
DO KK=1,NLAY
DO II=1,NROW
DO JJ=1,NCOL
IF(ICBUND(JJ,II,KK,INDEX).LT.0) THEN
ICBUND(JJ,II,KK,INDEX)=ABS(ICBUND(JJ,II,KK,INDEX))
ENDIF
ENDDO
ENDDO
ENDDO
ENDDO
@vivekbedekar Do you recall why this code was added? It looks like if first showed up with a commit you made in 2015:
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub, or unsubscribe.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub<#82 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AA5YRSH6EIHNDJCHRQJLS5TSCPD3NANCNFSM4QKATPEA>.
from mt3d-usgs.
emorway-usgs
commented on September 28, 2024
@DarkForceOne Hi Eric, I think a work around for your case could be to set ICBUND=1 for the constant concentration cells and then in SSM input, you could set the ITYPE for your constant concentration cells equal to -1. Using this approach, the constant concentration condition will persist in multiple stress period simulations. If this solution works for you, would you mind closing this issue?
from mt3d-usgs.
Related Issues (20)
- LKNODE in LKT is not set properly with MODFLOW 6 models HOT 3
- Setting of LKT terms when using MF6-generated flow fields HOT 2
- Code should stop when UZT, LKT, or SFT package is used with FT6 linker input
- MT3D and MF6 HOT 1
- Bad input on GWT test problem HOT 1
- LKT example output file shows time value higher than simulation time
- Incorrect CTS unit number in documentation HOT 1
- GWT passes locally, but fails on Travis HOT 2
- help with input error? HOT 7
- MT3D-USGS crashes with the attached model. HOT 9
- MT3D-USGS version is missing from LIST output file HOT 1
- sft_crnkNic fails on Windows and macOS HOT 2
- MT3D-USGS input instructions show NROW/NCOL for some arrays HOT 1
- Crash as soon as MODFLOWSTYLEARRAYS is added HOT 2
- End Loops HOT 1
- Multiple solid phases for kinetic ED/EA reactions HOT 4
- MT3D-USG - GWV8 HOT 3
- MT3D-USGS doesn't handle recharge and EVT properly with MODFLOW 6 models.
- Intermittent floating over/underflow
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from mt3d-usgs.