I've successfully tested your software without error, but there are two problems that can't be solved in the process of using it. about lasif HOT 7 OPEN

this is the Frist Question.

from lasif.

Hi!

regarding the first question: There has been some change in matplotlib or basemap with the default z-order (i.e. what get's drawn on top of what). Please play around a bit with the zorder arguments in the domain.py code as well as here:

The argument here has no zorder argument but it might need one. Basically you want to sort the zorders in a sense that only the continent outlines are above the model. Please consider sending a pull request if you end up solving this!

Your second question seems to a be a general python/numpy question. The first argument in np.array() is clearly a list of strings and it seems to contains some non-numeric entries which you want to set to dtype=float. This does not work but hard to help you via remote debugging. Just carefully look at all the pieces of your code and execute it step by step to understand where it goes wrong. This is also a very helpful community if you are just starting out to learn Python: https://www.reddit.com/r/learnpython/

Cheers!

from lasif.

Thank you very much for your help.
Let me first change the first question according to your suggestion.
The second problem is that I didn't make it clear.
Is the gradient_csv file in the script directly obtained by SES3D?
All the gradient files I get from SES3D are grad_csv_*.And I don't find any part of the output gradient_csv in the naming format of the output file in theses3d_output.f90file.(please below)

What can I do to get the gradient_csv file needed for the script?

from lasif.

FYI there is newer C based SES3D version but the inputs and outputs are almost identical: http://www.cos.ethz.ch/software/production/ses3d.html

Where does the model update script come from? If it is part of SES3D then it operates on the blockfiles so you'd first have to convert the grad_* files to the blockfiles with one of the fortran programs (forgot which one). But I'm not entirely sure I understand your question to be honest.

from lasif.

The ses3d-master comes from https://github.com/echolite/ses3d.
The model update script comes from http://lasif.net/tutorial/14_model_update.html
Is the grad file structure obtained by using ses3d of f90 and F70 different?
Because I do not know the specific technical details of the model update section in step 15 of the LASIF tutorial(please below), I can only rely on my own guess.When I adjust the script path and provide updates to the required files, I don't know if I'm doing the right thing. I hope I can get your help.
In addition, I'm reading articles about lasif and obspy from you and your team, convinced by the power of these two software. I admire you and your team for this great achievement.
Thank you very much for your patient guidance and help.

from lasif.

Yea okay - this script takes files in the blockfile format. These files here can convert form the binary files to the blockfile format: https://github.com/echolite/ses3d/tree/master/TOOLS/SOURCE

Yes - this is the link to the github repository, but the link I sent still contains a more recent version that has been rewritten in C. To my knowledge there is no repository for this code.

from lasif.

Okay, I'll try it first according to the method you provided.
Thank you for your guidance.

from lasif.