Simple prism model in Cartesian coordinatesΒΆ
The simple_prism.sh
script calculates
the gravitational potential, gravitational attraction,
and gravity gradient tensor due to a simple prism model
in Cartesian coordinates:
#!/bin/bash
# Generate a regular grid, pipe it to all the computation programs,
# and write the result to output.txt
tessgrd -r0/20000/0/20000 -b50/50 -z1000 | \
prismpot model.txt | \
prismgx model.txt | prismgy model.txt | prismgz model.txt | \
prismgxx model.txt | prismgxy model.txt | \
prismgxz model.txt | prismgyy model.txt | \
prismgyz model.txt | prismgzz model.txt > output.txt
The model file looks like this:
# Test prism model file
2000 5000 2000 15000 0 5000 1000
10000 18000 10000 18000 0 5000 -1000
The result should look like the following (“column” means the column of the output file).