search

rniswon / mfnwt

MODFLOW-NWT - Version: 1.1.0 Newton Formulation for MODFLOW-2005 For more information go to http://water.usgs.gov/ogw/modflow-nwt
4 stars 7 forks source link

smoothing of conductance in SFR uses stream head and should use depth #22

Closed rniswon closed 6 years ago

rniswon commented 6 years ago

From Peter Vermeulen "Though I know MF2005 is history, I discovered a bug in the FM-routine of the SFR of the latest released MF2005 version. Below is what was original stated at line 3400 in green. In black what it should be in my opinion. I came along this bug as my bottom heights and levels are negative (stupid Dutch negative-elevation J). Anyhow, it might be that you included this code in MF6 and therefore you might want to change this.

! IF ( icalc.EQ.1 ) cstrsmooth = cstr*

! + smooth(hstrave,dwdh)

        IF ( icalc.EQ.1 ) cstrsmooth = cstr*

 +                        smooth(hstrave-sbot,dwdh)"
rniswon commented 6 years ago

fixed by replacing line above with depthave = hstrave - strtop if ( depthave < 0.0 ) depthave = 0.0 cstrsmooth = cstr IF ( icalc.EQ.1 ) cstrsmooth = cstr*