Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

I've successfully tested your software without error, but there are two problems that can't be solved in the process of using it. #9

Open
Baylor96 opened this issue Mar 26, 2019 · 7 comments

Comments

@Baylor96
Copy link

Respected Lion Krischer:
When I execute py.test ,I got the same result as you(), indicating that my LASIF program was installed correctly.
image
The First Question:When i using lasif launch_model_gui, the images I get from your tutorial can only show part of the images.In the physical domain, data information is displayed only in the Mediterranean Sea, but not on European land. Why?
The Last Question:When I use the script that you update the model, although I adjust the script according to my own file, there is a problem that can not run. How can I solve it?
image
If you can get your help, we would be very grateful!!!

@Baylor96
Copy link
Author

image
this is the Frist Question.

@krischer
Copy link
Owner

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:

im = m.pcolormesh(x, y, depth_data, cmap=tomo_colormap, vmin=vmin,

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!

@Baylor96
Copy link
Author

Baylor96 commented Apr 12, 2019

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)
image
What can I do to get the gradient_csv file needed for the script?

@krischer
Copy link
Owner

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.

@Baylor96
Copy link
Author

Baylor96 commented Apr 12, 2019

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.
WX20190412-231736@2x

@krischer
Copy link
Owner

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.

@Baylor96
Copy link
Author

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants