vuribe1221 / gravity_modeling

Gravity Modeling For Jeff Kennedy

Introduction

Code Written by Victor Uribe

*This code assumes that you have access to Futuregrid or some other HTC/HPC and have cctools installed properly.

Purpose - Create a script that automatically generates a Makeflow file and uses Futuregrid (or other HTC systems) to simulate the gravity data collected by a gravity detector.

Result

Final Result: In red you can see the water (largest gravity) recede into the water table. About 7 seconds of simulation amounts to almost 7 hours of calculations on futuregrid. This simulations took several days to put together on a laptop workstation where it took only 7 hours on Futuregrid.

Code Samples

$ sh makeMakeflow.sh density_file N gravity_pos_file
   Where: N = Number of peices to chop the density file into

Output files:

      OUT
      -Found in the present working directory

Installation

Files required to run code:

1. grav.py
         *used to calculate the gravity given values in density_grid.txt.
2. prism.py
         *dependency used by grav.py to calculate gravity in each prism.
3. density_grid.txt
         *Predictions of subsurface water content (ie density) change.
4. grav_pos.txt
         *Measurement/prediction locations.
5. makeMakeflow.sh
         *Main Script file used to create a custom Makeflow file then runs torque on the data outputing .out files with 3 values.
6. addFiles.sh
         *Dependency of makeMakeflow.sh in order to make more modular and simplify the whole process.