Closed cneyens closed 3 years ago
vivekbedekar
commented
3 years ago @cneyens - did you develop an example problem that will help us test the code changes you made? Could you share the example problem including the MODFLOW files? I will be QAing the code changes and merging the code this week. Thanks!
cneyens
commented
3 years ago @cneyens - did you develop an example problem that will help us test the code changes you made? Could you share the example problem including the MODFLOW files? I will be QAing the code changes and merging the code this week. Thanks!
Please find an example model in attachment. It is based on the 2ED5EA example problem shipped with the MT3D-USGS code. I've adjusted the example to include a 6th EA (MnO2) simulated as a solid phase EA in addition to a Fe-oxide solid phase. Only the BTN and reaction module files were adjusted (+ renaming filename prefixes to 2ed6ea). The folder also includes a README file where this is explained in more detail.
I ran the flow model with MODFLOW-NWT v1.2.0 and the solute transport with the MT3D-USGS code as presented in this pull request which I compiled for Windows through pymake using the default compiler settings with gfortran.
We've also applied this approach to a field-case where we simulated DOC degradation in the presence of multiple EA's including pyrolusite and goethite. If you want, I can email you some of those results as well.
vivekbedekar
commented
3 years ago Thank you, Cas!
Vivek
Sent from Outlookhttp://aka.ms/weboutlook
From: Cas Neyens @.> Sent: Wednesday, September 8, 2021 7:57 AM To: MODFLOW-USGS/mt3d-usgs @.> Cc: vivekbedekar @.>; Review requested @.> Subject: Re: [MODFLOW-USGS/mt3d-usgs] Feature/edea multiple solids (#87)
@cneyenshttps://github.com/cneyens - did you develop an example problem that will help us test the code changes you made? Could you share the example problem including the MODFLOW files? I will be QAing the code changes and merging the code this week. Thanks!
Please find an example model in attachment. It is based on the 2ED5EA example problem shipped with the MT3D-USGS code. I've adjusted the example to include a 6th EA (MnO2) simulated as a solid phase EA in addition to a Fe-oxide solid phase. Only the BTN and reaction module files were adjusted (+ renaming filename prefixes to 2ed6ea). The folder also includes a README file where this is explained in more detail.
I ran the flow model with MODFLOW-NWT v1.2.0 and the solute transport with the MT3D-USGS code as presented in this pull request which I compiled for Windows through pymakehttps://github.com/modflowpy/pymake using the default compiler settings with gfortran.
We've also applied this approach to a field-case where we simulated DOC degradation in the presence of multiple EA's including pyrolusite and goethite. If you want, I can email you some of those results as well.
2ED6EAs.ziphttps://github.com/MODFLOW-USGS/mt3d-usgs/files/7128610/2ED6EAs.zip
— You are receiving this because your review was requested. Reply to this email directly, view it on GitHubhttps://github.com/MODFLOW-USGS/mt3d-usgs/pull/87#issuecomment-915170148, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AA5YRSFNT3AGNOTG6EL64G3UA5FRTANCNFSM4UWPHVQQ. Triage notifications on the go with GitHub Mobile for iOShttps://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Androidhttps://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.
Per #86.
Summary:
'SOLID'can be specified for EA's corresponding to these immobile phases (recordSPECIAL(ISPEC)). The first SOLID keyword case is linked to the first immobile species, the second SOLID keyword to the second immobile species and so forth. Previously, only 1 immobile phase and 1 mobile EA could be linked. All input pertaining to 'solid phase Fe' can also be used for other solid phases (e.g. the IFESLD flag in the RCT input file).Changes to input:
None. A specific ordering of the simulated species is assumed however, as was already the case previously:
Changes to code:
ICOMP==NCOMPcheck where the immobile phase was assumed to be the NCOMPth species. Using ISLDPH also allows more immobile species to be specified than are used by the ED/EA reactions, e.g. for use by other packages.