gimli-org / gimli

Geophysical Inversion and Modeling Library :earth_africa:
https://www.pygimli.org
Other
365 stars 133 forks source link

IP data inversion with ERTIPManager #765

Closed HenryWHR closed 2 weeks ago

HenryWHR commented 2 weeks ago

Hi,

Hope you all had a nice summer!

I am using ERTIPManager for IP data inversion.

mgr = ert.ERTIPManager(data,verbose=True)
mgr.setMesh(mesh)
mgr.invert(zWeight=0.3, lam=20)

I have two question regarding the inversion and the output. 1), How can I define the lambda for IP inversion, similar to 'lam' for DC resistivity. 2), How can I get the response of the final IP model, like 'mgr.inv.response' for resistivity model.

Many thanks for your kind help!

Best regards, Henry

halbmy commented 2 weeks ago

Dear Henry, instead of mgr.invert which currently uses the same lam argument, you should better do it in two calls

mgr.invertDC(lam=...)
mgr.invertTDIP(lam=...)
halbmy commented 2 weeks ago

To the second question: I made the inversion object invIP available (2eec1c77)

HenryWHR commented 2 weeks ago

Thanks a lot, Thomas! Solved!