jkennedy-usgs / heavy-mirror

This github.com site is a mirror of the official code release at code.usgs.gov. Issues can be opened on either platform (opening issues on code.usgs.gov requires an account on that platform).

Heavy

A forward-gravity program for MODFLOW-2005 and MODFLOW-NWT

Download Heavy Windows executable and source code (updated 2023-08-09)

Heavy is a Fortran program for calculating the change in gravity caused by the change in groundwater storage as predicted by a MODFLOW groundwater-flow model.

Gravity change is calculated relative to the initial time step at user-specified locations (e.g., measurement locations). At each location, the gravitational attraction is the sum of the attraction caused by storage change (confined or unconfined) within each model cell and layer.

Heavy is intended to be incorporated into workflows for parameter estimation such as PEST++.

Running Heavy on Windows:

$> Heavy -q
 Heavy debug routine
 Forsberg gravity, single 100000 x 100000-m layer,   1.0-m thick prism, g =  41.927 uGal
 Horizontal infinite slab (2*pi*G), g = 41.928 uGal
 Forsberg gravity, single layer 10000 x 10000-m layer composed of 10x10x1-m prisms, g = 41.85 uGal

 Gravitational attraction of a 20x20x20-m prism, 100 m deep, vs. an equivalent mass sphere:
 Heavy point-mass gravity, g = 6.590618 uGal; Newton's Law (G*m/r^2) = 6.553962 uGal

Process returned 0 (0x0)   execution time : 0.585 s
Press any key to continue.

(Copy /bin/Heavy.exe to /test/aac)

$> Heavy model1_hvy.nam

                                  HEAVY-0.0.2 2021-09-10
                         U.S. GEOLOGICAL SURVEY MODELING TOOL FOR
                    THE FORWARD MODELING OF MICROGRAVITY WITH MODFLOW

 Using NAME file: model1_hvy.nam
 Run start date and time (yyyy/mm/dd hh:mm:ss): 2021/12/17 13:33:49

 CALCULATING GRAVITY FOR ALL LAYERS
 Single Precision Binary Head file
 READ HEAD FOR PERIOD     1 STEP     1  PERTIM =   1.000000E+00  TOTIM =   1.000000E+00
 READ HEAD FOR PERIOD     1 STEP     1  PERTIM =   1.000000E+00  TOTIM =   1.000000E+00
...
 READ HEAD FOR PERIOD    20 STEP     1  PERTIM =   1.000000E+00  TOTIM =   7.600000E+01
 READ HEAD FOR PERIOD    21 STEP     1  PERTIM =   1.000000E+00  TOTIM =   7.700000E+01
Elapsed Run Time:      6.859 seconds

See /doc/User Guide.docx for information on running Heavy

Compiling Heavy:

Windows There is a compiled executable in the /bin/Release directory.

Heavy has been compiled on Windows using the Code::Blocks IDE with the gfortran compiler, and using Visual Studio 2019 with the Intel oneAPI Fortran compiler. The latter is somewhat faster but more difficult to compile a standalone executable.

Linux

1. Install pymake: pip install --user https://github.com/modflowpy/pymake/zipball/master

2. Compile: python -m pymake ./src ./bin/heavy -mc

Citation

Kennedy, J.R., and Larson, J.D., 2023, Heavy: software for forward-modeling gravity change from MODFLOW output. U.S. Geological Survey Software Release, https://doi.org/10.5066/P9IIHXN3

This software has been approved for release by the U.S. Geological Survey (USGS). Although the software has been subjected to rigorous review, the USGS reserves the right to update the software as needed pursuant to further analysis and review. No warranty, expressed or implied, is made by the USGS or the U.S. Government as to the functionality of the software and related material nor shall the fact of release constitute any such warranty. Furthermore, the software is released on condition that neither the USGS nor the U.S. Government shall be held liable for any damages resulting from its authorized or unauthorized use.

IP-137201

heavy-mirror's Issues

Obtaining a gravity time series formed by zeros

I am using HEAVY to calculate the gravitational signal observed from a point within a very simplified hydrogeological model.
I have created and filled the requiered files (nam and hvy) properly, and executed the program successfully with the default options.
The lst and out files are created well within Resultados folder, nevertheless, the created gravity signal is completeley formed by zeroes.

I do not know what could be happening. I attached all my files as a zip file (Data.zip), in case you could help.

For the hydrogeological model I am running MODFLOW 2005 withModelMuse, from a Windows 10 x64 bits machine.
Thank you very much in advance.

Is possible to obtain a non-absolute values gravity changes time-series?

Hello again,

I am using Heavy to estimate the gravimetry time-series generated by changes in aquifer storage. In the hydrogeological model I am working on (ExampleData.zip), the piezometric time series of the observation point fluctuates generating a concave curve upwards: it decreases and then increases exceeding the initial piezometric value.

I have realized that the gravity time-series that Heavy calculates for me is made up of absolute values, not including sign (negative or positive) of the gravity change. Since Heavy calculates gravity changes with respect to the value of the first time-step, this generates in the gravimetric time series a fictitious fluctuation when (during water table recovery) the piezometric value exceeds the initial value (reference point, ie. zero gravimetric difference).

I open the issue then to confirm that the output is intended (not a execution problem of mine) and to know if there is the option of obtaining a non-absolute gravimetric time-series, in which the changes include the direction of the difference (positive or negative sign). I think this could be very useful when comparing the estimated gravimetric signal with any other variable.

You can find attached the data I am using.

Thank you very much!