Currently, to find the highest active head (e.g. the water table) the call is like this wt0 = get_water_table(hds0, nodata=hdry). This is calling get_water_table from the flopy.utils.postprocessing. However, with MF6, should also test for hnoflo. In my application, there are no dry cells, but there are inactive cells. If the hnoflo value is used for nodata arg, the results make sense, but with hdry the result is a patchwork of values.
I'm not sure if this can be fixed in this package or if flopy.utils.postprocessing should be updated to accomodate multiple nodata values. That's probably the fix, but it would also require passing both nodata values along from here, so wanted to start the discussion here.
Currently, to find the highest active head (e.g. the water table) the call is like this
wt0 = get_water_table(hds0, nodata=hdry). This is callingget_water_tablefrom theflopy.utils.postprocessing. However, withMF6, should also test forhnoflo. In my application, there are no dry cells, but there are inactive cells. If thehnoflovalue is used fornodataarg, the results make sense, but withhdrythe result is a patchwork of values.I'm not sure if this can be fixed in this package or if
flopy.utils.postprocessingshould be updated to accomodate multiplenodatavalues. That's probably the fix, but it would also require passing bothnodatavalues along from here, so wanted to start the discussion here.