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

issue with surface processes #696

Open
hanybit opened this issue May 20, 2024 · 0 comments
Open

issue with surface processes #696

hanybit opened this issue May 20, 2024 · 0 comments

Comments

@hanybit
Copy link

hanybit commented May 20, 2024

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

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

1 participant