The 'and' operator in this line causes data points that have only top or botm to be discarded:
|
is_original_arrays[layer_name] = ts["is_original"].values & ds["is_original"].values |
Results in loss of many data points. An 'or' operator would fix this. Will test if this produces meaningful results and if it does create a PR.
The 'and' operator in this line causes data points that have only top or botm to be discarded:
data/src/nhflodata/data/mockup/bodemlagen_pwn_2024/v2.0.0/botm/botm.py
Line 296 in 882a24d