[CIG-SHORT] PyLith Question

[Brad Aagaard]

Cecep,

You should not use a mesh in lon/lat/elev coordinates (CSGeo), because 
it does not accurately describe distances between points. You need to 
use a geographic projection for the mesh or a conventional Cartesian 
coordinate system. See examples/twocells/twotet4-geoproj for an example 
that uses a mesh in a geographic projection.

For OutputSolnPoints, after reading in the list of points from the file, 
the coordinates are converted to the mesh coordinate system. The 
interpolation of the solution fields for output is then done in the mesh 
coordinate system.

Regards,
Brad


On 12/02/2015 09:59 PM, Cecep Pratama wrote:
> Dear Mr. Brad,
>
> Let me introduce myself simply. My name, Cecep from Indonesia. Currently studying Crustal Deformation.
>
> I use Pylith, the software that you developed, for simulation. I am stuck with OutputSolnPoints.
>
> I already put the coordinate in output_points.txt. The coordinate related with GPS site using Lat Long and Height dimension. And I set the coordsys for OutputSolnPoints in CSGeo also.  But, I am afraid Pylith convert point list in ouput_points.txt file and assuming that list points are 3D cartesian.
>
> How do you solve this issue?
>
> Note:
> My mesh also in CSGeo coordinate
>
> Thank you
>
> Regards,
> Cecep
>