Integrate 2D point loads into GRASS GIS: v.flexure - Githubissues

awickert / gFlex

Multiple methods to solve elastic plate flexure, designed for applications to Earth's lithosphere.

http://csdms.colorado.edu/wiki/Model:gFlex

GNU General Public License v3.0

21 stars 12 forks source link

Integrate 2D point loads into GRASS GIS: v.flexure #6

Closed awickert closed 9 years ago

awickert commented 9 years ago

Use a set of points with a loads attribute (or this set proportzional to z).
Pull x and y from geometry: no need to be in attribute table.

awickert commented 9 years ago

And add a recipe for generating voronoi cells around them to compute reasonable areas over which those loads are important.

The other option is to interpolate point loads into a gridded set of loads -- show how to do that too in help/tutorial file.

NOT doing this -- what to do with point loads on edges? Instead, just expecting users to give full mass.

awickert commented 9 years ago

It will become a vector command then! So must be v.flexure

awickert commented 9 years ago

Now at the point at which I pull x and y. Lat/lon check. If ll, then use a function to find great circle distances. This might have to be implemented directly inside f2d.py, with a flag or something, as well as a note on R_planet (wow, another possible variable!)

awickert commented 9 years ago

Complete! Works for both unprojected and lat/lon locations. The latter is really exciting: it uses great circle paths, and thus could be used for global calculations! Of course, spherical harmonics with FFT-style solutions would be better for a globe, but still nice to know that it is possible to operate around the poles.