Stormwater-Management-Model icon indicating copy to clipboard operation
Stormwater-Management-Model copied to clipboard

Modified Horton Infiltration Bug

Open LRossman opened this issue 1 year ago • 3 comments

Line 545 in the modHorton_getInfil() function of infil.c should be changed from

        if ( infil->Fmax > 0.0 ) infil->Fe = MAX(infil->Fe, infil->Fmax);

to

        if ( infil->Fmax > 0.0 ) infil->Fe = MIN(infil->Fe, infil->Fmax);

LRossman avatar Jan 18 '24 00:01 LRossman

Recall that the modified Horton infiltration method was updated to limit the maximum infiltration in issue #59 to reflect changes that failed to make it into the main branch as part of PR #56 for the v5.1.15 update .

While the rationale for the fix suggested in this issue is clear enough to infer, it will be useful to provide some clarification as to why these changes are being suggested for our understanding as folks tasked with maintaining and advancing SWMM and for community understanding.

If the previous fix addresses the issue raised here, I will go ahead and close this issue.

cbuahin avatar Jan 18 '24 00:01 cbuahin

No, this is a different bug not addressed in PR #56. As seen on page 103 of the Hydrology Reference Manual,

image

the updated cumulative infiltration should be the smaller of the newly computed value and the max. limit, not the larger of the two.

LRossman avatar Jan 18 '24 02:01 LRossman

Understood. Thanks Lew.

cbuahin avatar Jan 18 '24 03:01 cbuahin