fixing the example

    # porosity = a * exp(-b * depth), with depth in km

    poro_a = 28.97

    poro_b = 0.21

    poro_std = 5  # 1% uncertainty

    poro_std = 1  # 1% uncertainty, if this goes larger it blows up massively...

    # ln(permeability) = a * poro + b

    perm_a = 0.334

    ax2.fill_betweenx(p50.depth, p10.power, p90.power, color="black", alpha=0.5)

    ax2.plot(p50.power, p50.power.depth, label="Power [MW]", color="black")

    # transmissivity Dm

    # transmissivity

    ax.fill_betweenx(

        p50.depth,

        p10.transmissivity_with_ntg,