search

underworldcode / underworld2

underworld2: A parallel, particle-in-cell, finite element code for Geodynamics.
http://www.underworldcode.org/
Other
168 stars 58 forks source link

issue with surface processes #696

Open hanybit opened 4 months ago

hanybit commented 4 months ago

hello all,

I have a small issue with my surface processes when I run on cluster on high resolution i got this error, but when I lowered the resolution of my model I got another error I added at bottom, any ideas how to solve this??

File "/monfs01/scratch/hkha0013/./2D_base.py", line 246, in Model.run_for(30.0 u.megayear, checkpoint_interval=0.5 u.megayear) File "/opt/venv/lib/python3.9/site-packages/UWGeodynamics/_model.py", line 1713, in run_for self._update() File "/opt/venv/lib/python3.9/site-packages/UWGeodynamics/_model.py", line 1829, in _update self.surfaceProcesses.solve(dt) File "/opt/venv/lib/python3.9/site-packages/UWGeodynamics/surfaceProcesses.py", line 830, in solve tracer_velocity = self.Model.velocityField.evaluate(self.surface_data_local) File "/opt/venv/lib/python3.9/site-packages/underworld/function/_function.py", line 938, in evaluate return _cfn.Query(self._fncself).query(_cfn.NumpyInput(inputData,inputType)) ValueError: Issue utilising function of class 'MeshVariable'

For full func debug info on all processes, set the UW_WORLD_FUNC_MESSAGES environment variable.

Error message: FeVariable interpolation at location (0.766373, -0.0365942) does not appear to be valid. Location is probably outside local domain. application called MPI_Abort(comm=0x84000002, 1) - process 11

but when I lowered the resolution of my model I got this error

An uncaught exception was encountered on processor 0. Traceback (most recent call last): File "/monfs01/scratch/hkha0013/./2D_base.py", line 246, in Model.run_for(30.0 u.megayear, checkpoint_interval=0.5 u.megayear) File "/opt/venv/lib/python3.9/site-packages/UWGeodynamics/_model.py", line 1713, in run_for self._update() File "/opt/venv/lib/python3.9/site-packages/UWGeodynamics/_model.py", line 1829, in _update self.surfaceProcesses.solve(dt) File "/opt/venv/lib/python3.9/site-packages/UWGeodynamics/surfaceProcesses.py", line 900, in solve f = interp1d(x2, z2, kind='cubic', fill_value='extrapolate') File "/opt/venv/lib/python3.9/site-packages/scipy/interpolate/_interpolate.py", line 564, in init self._spline = make_interp_spline(xx, yy, k=order, File "/opt/venv/lib/python3.9/site-packages/scipy/interpolate/_bsplines.py", line 1295, in make_interp_spline raise ValueError("Expect x to not have duplicates") ValueError: Expect x to not have duplicates application called MPI_Abort(comm=0x84000004, 1) - process 0