bug: flopy.utils.get_transmissivities bug report

[Paszka1]

Please refer to cell 4 of the flopy3_get_transmissivities_example notebook. Is there a bug in the case of well 0? I.e., transmissivities should all be 0 for this well, but it is 2 in layer 3.

[Paszka1]

As a follow up on the above message, I checked the source code of the module and I found some unclear steps at the end of the function code. I.e., thick[closest, i] = 1.0 in the last line of the for loop doesn't make sense to me.

    not_in_layer = np.sum(thick < 0, axis=0)
    not_in_any_layer = not_in_layer == thick.shape[0]
    for i, n in enumerate(not_in_any_layer):
        if n:
            closest = np.argmax(thick[:, i])
            thick[closest, i] = 1.0
    thick[thick < 0] = 0
    thick[heads == nodata] = 0  # exclude nodata cells