[CIG-SHORT] 答复: question about coordinate transform

Sat Aug 19 00:28:40 PDT 2017

Thank you for your reply.

Now, I have check the demo and found some errors. I corrected the errors, then, run the demo with Pylith 2.2.0 version. But it still report some errors:

[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[0]PETSC ERROR: to get more information on the crash.
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Signal received
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.7.5-3660-g9beae12  GIT Date: 2017-03-23 13:26:44 -0500
[0]PETSC ERROR: /home/tx/pylith/pylith-2.2.0-linux-x86_64/bin/mpinemesis on a arch-pylith named insar.geodyn by tx Sat Aug 19 07:25:43 2017
[0]PETSC ERROR: Configure options --prefix=/home/brad/pylith-binary/dist --with-c2html=0 --with-x=0 --with-clanguage=C --with-mpicompilers=1 --with-debugging=0 --with-shared-libraries=1 --with-64-bit-points=1 --with-large-file-io=1 --download-chaco=1 --download-ml=1 --download-f2cblaslapack=1 --with-hwloc=0 --with-ssl=0 --with-x=0 --with-c2html=0 --with-lgrind=0 --with-hdf5=1 --with-hdf5-dir=/home/brad/pylith-binary/dist --with-zlib=1 --LIBS=-lz --with-fc=0 CPPFLAGS="-I/home/brad/pylith-binary/dist/include " LDFLAGS="-L/home/brad/pylith-binary/dist/lib " CFLAGS="-g -O2" CXXFLAGS="-g -O2 -DMPICH_IGNORE_CXX_SEEK" FCFLAGS= PETSC_DIR=/home/brad/pylith-binary/build/petsc-pylith PETSC_ARCH=arch-pylith
[0]PETSC ERROR: #1 User provided function() line 0 in  unknown file
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0
/home/tx/pylith/pylith-2.2.0-linux-x86_64/bin/nemesis: mpirun: exit 59
/home/tx/pylith/pylith-2.2.0-linux-x86_64/bin/pylith: /home/tx/pylith/pylith-2.2.0-linux-x86_64/bin/nemesis: exit 1

________________________________
发件人: tu xiang
发送时间: 2017年8月18日 7:25:44
收件人: cig-short at geodynamics.org
主题: 答复: question about coordinate transform

Thank you for your help about coordinate transform.

But I have another two questions.

question 1

Yes, I am importing the coordsys.py from the example, the question is z-coordinates are greater than zero. I want to bury the fault into the mesh,  but z-coordinates of the mesh are smaller than zero. What do I need to do?

change the cs.inventory.datumVert = "mean sea level" or the cs.inventory.projector.inventory.projOptions = "+lat_0=45.5231 +lon_0=-122.6765 +k=0.9996"

and how to change?

question2
I have construct a very simple 3D demo mesh, and I want to do dynamic simulation with pylith 1.9.0 version, then I set all the configure file. but it output errors. I attach the demo. Hope someone would like to help me check the demo.

output errors
Traceback (most recent call last):
  File "/home/xzhou/pylith-1.9.0-linux-x86_64/bin/pylith", line 39, in <module>
    start(applicationClass=PyLithApp)
  File "/home/xzhou/pylith-1.9.0-linux-x86_64/lib/python2.7/site-packages/pythia-0.8.1.14-py2.7.egg/pyre/applications/__init__.py", line 41, in start
    shell.run(**kwds)
  File "/home/xzhou/pylith-1.9.0-linux-x86_64/lib/python2.7/site-packages/pythia-0.8.1.14-py2.7.egg/pyre/applications/Shell.py", line 125, in run
    app.applyConfiguration(context)
  File "/home/xzhou/pylith-1.9.0-linux-x86_64/lib/python2.7/site-packages/pythia-0.8.1.14-py2.7.egg/pyre/inventory/Configurable.py", line 84, in applyConfiguration
    raise ValueError("%s\nBacktrace - Component %s" % (err.message, aliases))
ValueError: Error while configuring spatial database reader (simpleioascii, iohandler):
Filename for spatial database not specified.
Backtrace - Component simpleioascii, iohandler
Backtrace - Component simpledb, slip
Backtrace - Component stepslipfn, slip_function
Backtrace - Component eqkinsrc, rupture
Backtrace - Component singlerupture, eq_srcs
Backtrace - Component faultcohesivekin, demo_fault
Backtrace - Component interfaces
Backtrace - Component timedependent, problem
Backtrace - Component pylithapp

________________________________
发件人: tu xiang
发送时间: 2017年8月15日 10:59:45
收件人: cig-short at geodynamics.org
主题: question about coordinate transform

Hello, everyone

Does the geographical coordinates (+lat_0=45.5231 +lon_0=-122.6765 +k=0.9996) correspond to the origin of Cartesian coordinates for mesh in the 3D subduction example?

I want to try to transform some geographical coordinates, then, I write next code:
#!/usr/bin/env nemesis

import coordsys
import numpy as np
point=np.array([[45.5231,-122.6765,0.9996],[46.0,-122.0,1.0],[45.0,-122.0,0.5]])
print point
print type(point)
coordsys.geoToMesh(point)

But it output errors!

Traceback (most recent call last):
  File "demo.py", line 9, in <module>
    coordsys.geoToMesh(point)
  File "/home/bxu/covert/coordsys.py", line 54, in geoToMesh
    convert(xyz, cs_mesh(), cs_geo3D())
  File "/home/bxu/pylith/pylith-2.2.1rc1-linux-x86_64/lib/python2.7/site-packages/spatialdata/geocoords/Converter.py", line 33, in convert
    geocoords.Converter.convert(coords, csDest, csSrc)
RuntimeError: Error while projecting location.
  latitude or longitude exceeded limits
  projection: tmerc
  units: m
  proj options: +lat_0=45.5231 +lon_0=-122.6765 +k=0.9996
  lon: 0.794528
  lat: -2.14111

It said the latitude or longitude exceeded limits, I do not know what does it mean?

Looking forward to your reply.
Best regards,
Tu Xiang

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170819/5390da1c/attachment-0001.html>