Det J is negative

Yes it is the search length.
In my example (from a regional model) the mesh is large and element edge lengths range from 200m - 2000m . The method identifysubdomain didn’t work with s smaller than 1. It seems that s must fit somehow to your model.