Contrained TDIP inversion

[aurorec]

Question

Hello,

I am trying to performed a 3D TDIP constrained inversion. it seems that the ERT part of the constrained inversion worked but not the TDIP part.
I first build a 3D mesh with markers and a resistivity initial model. Then I performed the ERT constrained inversion and after try to implement constrained in the TDIP inversion. It appears that cheking the verbose the ERT constrained have been successfully applied but not the TDIP constraints.
I used this type of implementation :
mgr.invIP.setRegularization(2, limits=[0, 1], cType=1, startModel=0.01)

Is TDIP constrained inversion possible in pygimli ? If yes what would be the correct way to implement it ?

Thank you very much for your help

Best regards

Aurore

I join you the code and data.

ForPG.zip