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.

Sign up for GitHub

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

Jump to bottom

Error when closing first body and refreshing #14

Open
matty-434 opened this issue Apr 14, 2021 · 8 comments
Open

Error when closing first body and refreshing #14

matty-434 opened this issue Apr 14, 2021 · 8 comments

Comments

@matty-434
Copy link

matty-434 commented Apr 14, 2021 •
edited
Loading

Hi,
I am getting this error as soon as I close my first body and refresh:

You have selected following color
#ffaa7f
This is the coverted color
[u'#ffaa7f']
addig this body
Total number of bodies
1
body number count
1
Total number of envelpos
1
Working on Body number
1
start @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
Start Index:
0
Start node
(0.0, 0.024386118)
end #################################################################################
End Index:
346
End node
(692.0, 0.13966578000000002)
Fixing start
Fixing End
here is the problme _ master polygon
[(0.0, -28.614538012747033), (10.0, -30.12426637562254), (40.0, -33.79074954260588), (74.0, -39.182636552875522), (118.0, -42.160852133325363), (188.0, -38.480641493626933), (280.0, -39.094009933576672), (364.0, -38.480641493626933), (436.0, -37.253904613727514), (514.0, -34.187062413978822), (594.0, -34.187062413978822), (660.0, -36.027167733828037), (692.0, -41.547483693375625), (692.0, -41.547483693375625), (692.0, -400), (0.0, -400)]
here is the problem
after closing current body
Total number of bodies
1
Body count number
1
Total number of envelpos
2
Exception in Tkinter callback
Traceback (most recent call last):
File "/usr/lib/python2.7/lib-tk/Tkinter.py", line 1550, in call
return self.func(*args)
File "/usr/local/lib/python2.7/dist-packages/matplotlib/backends/backend_tkagg.py", line 420, in button_release_event
FigureCanvasBase.button_release_event(self, x, y, num, guiEvent=event)
File "/usr/local/lib/python2.7/dist-packages/matplotlib/backend_bases.py", line 1932, in button_release_event
self.callbacks.process(s, event)
File "/usr/local/lib/python2.7/dist-packages/matplotlib/cbook.py", line 549, in process
proxy(*args, **kwargs)
File "/usr/local/lib/python2.7/dist-packages/matplotlib/cbook.py", line 416, in call
return mtd(*args, **kwargs)
File "/usr/local/lib/python2.7/dist-packages/matplotlib/widgets.py", line 229, in _release
func(event)
File "/mnt/d/LitMod2D_2.0_package_dist_users-master/GUI/my_lib.py", line 954, in plot_model
if(temp[0]!='>'):
IndexError: list index out of range

I have no idea how to solve this. :(
@matty-434 matty-434 changed the title Unable to import the wx backend for pyface Error when closing first body and refreshing Apr 14, 2021
@matty-434
Copy link
Author

Also, I noticed the bounding box between (x0,x_end,z0,400) around the model is not there anymore (as it was before).

@ajay6763
Copy link
Owner

I thing it related to the backend. It is looking for tkagg and should look for qt4

You can try a small hack below.
Go to the following file and change self.toolkit to qt4
/usr/lib/python/dist-packages/traits/etsconfig/etsconfig.py

in this file change the

change self._toolkit = "qt4"

let me know if it works

@matty-434
Copy link
Author

Hi, yeah it worked. I don't see that error message anymore.
Rather, There's this one:

Total no of bodies saved
3
[u'Crust', u'SCLM_Mantle', 'asthenosphere']
running the model
INTERCEPTION NOT FOUND: BODY 3 NODE 33279
VERTICES: X1,Z1: 692000.00000000000 -163607.00000000000
X2,Z2: 692000.00000000000 -163607.00000000000

...READING THERMODYNAMIC INFORMATION...
Num of different mantle + crust bodies = 2
Type of mantle and crust
(if < 10) bodies =
23
99

(...)

...DONE!!...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG

:/

@ajay6763
Copy link
Owner

This is simply related to model you made where you have ovellaping in the bodies. Go to these nodes in the model and make sure interfaces of neighbouring bodies are not overlapping

@matty-434
Copy link
Author

Thank you very much for your suggestions!
Unfortunately, they don't seem to overlap. I am sending litmod.inp, bodies.out and bodies_GUI_envelops.out files from the test model that I was doing. Interestingly, that always happens to the model for material types equal to 0 (i.e. my crustal bodies). The mantle materials (e.g. 23, 99, etc.) seem to be considered in the calculation no matter what. The same does not happen to the crustal ones (i.e. 0).
bodies_out.txt
bodies_GUI_envelops_out.txt
litmod_inp.txt

@ajay6763
Copy link
Owner

ajay6763 commented Apr 21, 2021 via email

@matty-434
Copy link
Author

Thank you!

This is how the model that you sent looks when I open it in the GUI:
Model_GUI

This is what I get after running the above test model.

['topo.inp', 'ELEVin']
2
Types of bodies read:
['normal', 'normal', 'normal', 'normal', 'anomaly', 'anomaly', 'normal']
Noumber of bodies read = 5
Start and End of the profile
(0.0, 1000.0)
No of envelops read = 6
Anomalous bodies information
['Vp(%)', 'Thermal']
['2.0', '100.0']
[90, 99]
running the model

Anomalous composition of material 5 is set to PUM(99) just for this 1st run
In the 2nd run effects of composition are included

INTERCEPTION NOT FOUND: BODY 3 NODE 7261
VERTICES: X1,Z1: 1000000.0000000000 -150000.00000000000
X2,Z2: 743000.00000000000 -150000.00000000000
INTERCEPTION NOT FOUND: BODY 4 NODE 4185
VERTICES: X1,Z1: 525000.00000000000 -79641.000000000000
X2,Z2: 428000.00000000000 -137632.00000000000
INTERCEPTION NOT FOUND: BODY 5 NODE 5135
VERTICES: X1,Z1: 430000.00000000000 -139402.00000000000
X2,Z2: 527000.00000000000 -81250.000000000000
INTERCEPTION NOT FOUND: BODY 7 NODE 9661
VERTICES: X1,Z1: 1000000.0000000000 -150000.00000000000
X2,Z2: 1000000.0000000000 -150000.00000000000

YOU ENTERED SUB-LITHOSPHERIC ANOMALIES, SO LITMOD
WILL RUN TWICE TO CALCULATE DECOUPLED TOPOGRAPHY

LITMOD INITIATES THE FIRST RUN

...READING THERMODYNAMIC INFORMATION...
Num of different mantle + crust bodies = 3
Type of mantle and crust
(if < 10) bodies =
90
99
97
...THERMAL STRUCTURE AND HEAT FLOW...
Number of non zero elements in matrix: 31771
Iteration no. 0: Max DT=********
Iteration no. 1: Max DT= 177.2
Iteration no. 2: Max DT= 13.9
Iteration no. 3: Max DT= 2.0
Iteration no. 4: Max DT= 0.4

...DENSITY, ELEVATION & SEISMIC VELOCITIES...

LITMOD ENDS THE FIRST RUN

Then second run...

INTERCEPTION NOT FOUND: BODY 3 NODE 7261
VERTICES: X1,Z1: 1000000.0000000000 -150000.00000000000
X2,Z2: 743000.00000000000 -150000.00000000000
INTERCEPTION NOT FOUND: BODY 4 NODE 4185
VERTICES: X1,Z1: 525000.00000000000 -79641.000000000000
X2,Z2: 428000.00000000000 -137632.00000000000
INTERCEPTION NOT FOUND: BODY 5 NODE 5135
VERTICES: X1,Z1: 430000.00000000000 -139402.00000000000
X2,Z2: 527000.00000000000 -81250.000000000000
INTERCEPTION NOT FOUND: BODY 7 NODE 9661
VERTICES: X1,Z1: 1000000.0000000000 -150000.00000000000
X2,Z2: 1000000.0000000000 -150000.00000000000

THIS IS THE SECOND RUN OF LITMOD

...READING THERMODYNAMIC INFORMATION...
Num of different mantle + crust bodies = 3
Type of mantle and crust
(if < 10) bodies =
90
99
97
...THERMAL STRUCTURE AND HEAT FLOW...
Number of non zero elements in matrix: 31771
Iteration no. 0: Max DT=********
Iteration no. 1: Max DT= 177.2
Iteration no. 2: Max DT= 13.9
Iteration no. 3: Max DT= 2.0
Iteration no. 4: Max DT= 0.4

Material 5 with litho= 90 corresponds to:
DV referred to Vp(%)= 2.0000000000000000

Material 6 with litho= 99 corresponds to:
Material with DT= 100.00000000000000 K

Your model includes 2 anomalous materials
-320.39199389648434 8.1099999999999994 1325.9320000000000 4419699923.2849998 90
-331.39999389648438 8.1199999999999992 1331.7500000000000 4515386878.1230001 90
-318.24199389648447 8.0899999999999999 1334.1420000000001 4241272320.4619999 90
-333.23199389648448 8.1099999999999994 1349.1320000000001 4487444293.9910002 90
-311.59999389648419 8.1099999999999994 1358.6099999999999 4643317864.6630001 90
-303.51799389648431 8.0600000000000005 1329.7880000000000 3997053999.6900001 90
-298.24599389648438 8.0700000000000003 1340.0660000000000 4165424148.9759998 90
-301.91399389648427 8.0800000000000001 1348.9240000000000 4310798255.2259998 90
-321.84999389648442 8.0999999999999996 1358.4900000000000 4467866963.4499998 90
-305.49599389648438 8.0999999999999996 1368.0560000000000 4625109008.6610003 90
-299.78999389648448 8.109999999
etc...

...DENSITY, ELEVATION & SEISMIC VELOCITIES...

...GRAVITY AND GEOID ANOMALIES...

...OUTPUT OF RESULTS...

                    ...DONE!!...

END OF THE SECOND RUN OF LITMOD

Note: The following floating-point exceptions are signalling: IEEE_DENORMAL

And the results:
Figure 1
Figure 2

@ajay6763
Copy link
Owner

ajay6763 commented Apr 22, 2021 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@ajay6763 @matty-434