In the example, in the context of “if plotIt:” in function “run”, the line (about 201):
dpredFD = prbFD.dpred(moptTD)
before ax1.loglog(…
I think it should be moptFD, but not moptTD, because they are inversion solutions from Frequency-Domain and Time-Domain respectively. Of cause, it can be right, but the purpose of the author(s) here is to show the different results of different methods.