vtk2grav: A computer program for gravity forward of 3D density distribution in VTK format

Gravity forward calculation of a 3D density distribution, which is stored in VTK format. This code is originally designed for computing the gravity effect of a 3D temperature field generated from a geodynamic model, e.g. ASPECT output file. It supports OpenMP for parallel computing. The gravity field of each cubic cell in the VTK file is calculated from an accurate gravity forward formula proposed by Nagy(1996). The supported VTK cell type is VTK_VOXEL.

The source code can also be found at Github: https://github.com/zguoch/vtk2grav 

Usage:

$ vtk2grav  -h
██╗   ██╗████████╗██╗  ██╗██████╗  ██████╗ ██████╗  █████╗ ██╗   ██╗
██║   ██║╚══██╔══╝██║ ██╔╝╚════██╗██╔════╝ ██╔══██╗██╔══██╗██║   ██║
██║   ██║   ██║   █████╔╝  █████╔╝██║  ███╗██████╔╝███████║██║   ██║
╚██╗ ██╔╝   ██║   ██╔═██╗ ██╔═══╝ ██║   ██║██╔══██╗██╔══██║╚██╗ ██╔╝
 ╚████╔╝    ██║   ██║  ██╗███████╗╚██████╔╝██║  ██║██║  ██║ ╚████╔╝
  ╚═══╝     ╚═╝   ╚═╝  ╚═╝╚══════╝ ╚═════╝ ╚═╝  ╚═╝╚═╝  ╚═╝  ╚═══╝

========================== vtk2grav ===========================
vtk2grav, calculate gravity of 3D thermal structure.
Author              Zhikui Guo
Locus               GEOMAR, Germany
Date                2021-04-13
Version             1.0.0-git-12bdcc1
Email               zguo@geomar.de
============================================================
Usage: vtk2grav [options]
options:
  -h                List descriptions of usage and available arguments
  -v                Print swEOS version number
  -i                Input vtk file of temperature field.
  -E                Only extract T and save to a new vtu file
  -p                Position file of calculation points (xyz).
  -D                Set reference density(kg/m3) which will be subtracted from the input density field, e.g. -D 3300. Default is 0
  -F                Specify field name of density in the vtu file, default is density.
  -o                Set output file, e.g. -o density.txt
  -t                Set number of thread for parallel computing.
  -V                Set output information level, e.g. -V 0 means silent mode.
Example: vtk2grav -i T_density.vtu -p sites.xyz -D 3300 -T 1300 -A 1E-5 -o grav.txt