added a check for volaqu in ssm qss calc for itype=15

MODFLOW-USGS / mt3d-usgs

MT3D-USGS Repository

23 stars 12 forks source link

added a check for volaqu in ssm qss calc for itype=15 #22

Closed jtwhite79 closed 6 years ago

jtwhite79 commented 6 years ago

I was getting divide by zero for itype 15 for a cell that had gone dry (volaqu = 0.0). So I added a check reset volaqu to 1.0e-5 for this case. Seems to work for me, but not sure if this has other implications...

vivekbedekar commented 6 years ago

Your changes make sense. I added some more code in this commit: c5bbfc2.