HenrikSpa

This implementation would be great! I came to the same conclusion and the same solution, so that b_new = np.log(b_old), and your suggestion of b = np.exp(b) for the step...

Maybe you also have to make a change in get_tmax where b is used?: https://github.com/pastas/pastas/blob/1cd9fd94502920de8ae1bb03a902fd0d5ab25354/pastas/rfunc.py#L346

Great! But at the same time not great that you also had the error:) In my specific case it was a very short well series that caused it. I think...

But! Changing this behaviour will affect calibrated models that depend on the "flash drain"-effect of large amounts of precipitation and snowmelt.

_My suggested solution doesn't completely solve even with above 1000 steps it seems._ _EDIT: It does seems that 1000 solve the issue. (I made an error in the other code....

An early "test"-version of this is now working in qgis3-branch.  I haven't made any gui-changes for it yet (not sure where to put those settings). The grading is done...

The plan is to make it generic so that any dem could be coloured by any polygon layer. (it uses geoshort and zz_stratigraphy_plots currently). The bedrock could be styled in...

It's in qgis3_secplot_dem_style, not in general qgis3-branch.

I've made a change to the way this works to better separate the bore hole stratigraphy from the surface layer stratigraphy. It now picks the color from the surface geology...

Hi! You need to select the line layer containing you profile line (layer Polyline?).