Example of Parameter Input File

The following is an example of a parameter input file. It demonstrates how one can incorporate comments into the file, as well as how to specify variable values.

#  This is an example parameter file that can be used with the routine
#  RASTERTC.  This file also documents the values that may be used with
#  the routine.
#
#     ===  The following rules for constructing the file are used.  ===
#
#  1.  Any line beginning with a pound sign (#) is completely ignored.
#      This permits users to enter copious comments into the file. Also,
#      lines that are completely blank are ignored, as well.
#  2.  Any information on a line after a pound sign (#), including the
#      pound sign itself, is ignored.
#  3.  Lines, except for comment lines, should never exceed 80 characters
#      in length.
#  4.  An equals sign (=) MUST appear between the name of a variable and
#      the value of the variable.
#  5.  File names provided should be <= 64 characters in length.
#  6.  The case (upper, lower, or mixed) of a variable name provided is
#      irrelevant. The case of the variable values is also irrelevant,
#      except for the case of any file names provided.
#
#          Example: DENSITY=2.67 and density=2.67 are equivalent.
#
#  7.  If one wishes to use the debugging feature of RASTERTC, it is best
#      to always provide the value of the variable DEBUG as the first
#      variable provided in the input file.
#  8.  The BOOLEAN variables DEBUG and VERBOSE can be set TRUE by specifying
#      them as either "true", "on", or "yes" (without the quotation marks).
#      Any other values imply FALSE.
#  9.  The only variables for which one MUST provide a value are "gravfile"
#      and "demfile". The program will not execute if these values are
#      not provided. Defaults are taken for ALL other variable values.
# 10.  Blanks may appear on either side of the equals sign (=). For example,
#      "rmin = 50.0", "rmin=  50.0", and "rmin =50.0" are equivalent.
#
debug=off
verbose=on

# Specify whether terrain surfaces near are to be saved as Surfer
# gridfiles (1) or Geosoft gridfiles (2). If no such surfaces are
# to be saved (the default), specify 0.

savesurf=1   # Should be 0, 1, or 2.

gravfile=TestData.dat
auxelevfile=
listfile=RasterTC4.lst
outfile=RasterTC4.out
csvfile=RasterTC4.csv
rmin=10.0
rmed=300.0
angle=6   # specified in degrees (must divide evenly into 360)
rcurve=15000.0
gridtype=2  #  Note SURFER type = 1, GEOSOFT type = 2!
gridfile=GravElev.grd
rmaxgrid=2000.0  #  Must be specified in meters (not km)
density=2.67     #  In grams/cc

#  fileformat: Flag indicating the order of variables to be
#              read from the gravity data input file, as well
#              as the aux elevation file. The values are
#              interpreted as follows:
#
#                0 = StationID,x,y,elevation,obs gravity
#                1 = StationID,y,x,elevation,obs gravity
#                2 = x,y,StationID,elevation,obs gravity
#                3 = y,x,StationID,elevation,obs gravity
#                4 = x,y,elevation,obs gravity
#                5 = y,x,elevation,obs gravity
#                6 = x,y,elevation
#                7 = y,x,elevation

fileformat=6
elevunits=0  # 0 ==> feet; otherwise, meters


Table of Contents


Example of Listing File

RasterTC version 1.30[beta], Copyright 1998-2000, Geophysical Software.

Today is Tue Aug  8 12:43:09 2000

List of input and output file names.
====================================
Parameter File Name......... ParameterFile.dat
Gravity File Name........... TestData.dat
Aux Elev File Name.......... [no aux elev data provided]
Listing File Name........... RasterTC4.lst
Output File Name............ RasterTC4.out

Parameters controlling terrain correction calculation.
======================================================
Minimum radius for TC calculation .............     10.0.
Intermediate radius for TC calculation.........    300.0.
Angular increment for integration.............. 6 degrees.
Elevation units flag........................... 0 [0 ==> feet, otherwise, meters].
Elevation bias (added to grid elevevations)....      0.0 m.
Terrain density (assumed constant) ............    2.670 grams/cc.
Surface tension factor (between 0 and 100).....      0.0
Elevation grid format flag .................... 2.
Terrain grids will be saved for each station location.


List of Grid File(s) provided.
==============================
      Rmax (m)  GridFile Name
        2000.0  GravElev.grd


The above grid(s) are assumed to be in Geosoft Version 2 [16-bit] format.

**  Starting to read the gravity data file **

**  Finished reading gravity data file; 5 stations read  **

5 gravity station locations read from TestData.dat.
 Count StationID   Easting(m)   Northing(m)   Elev(m)   Elev(ft)  ObsGravity
     1        1     195000.0      99000.0     192.38     631.17       0.000
     2        2     195000.0      99100.0     184.61     605.66       0.000
     3        3     195000.0      99200.0     180.24     591.33       0.000
     4        4     195000.0      99300.0     172.69     566.56       0.000
     5        5     195000.0      99400.0     175.20     574.79       0.000

**  Geosoft Grid Parameters  **
===============================
   Geosoft grid name: GravElev.grd
   No. of elements per vector: 851.
   No. of vectors per grid:    891.
   Grid data element size:     2 [must be 2].
   Grid orientation flag:      -1 [must be 1 or -1].
   Grid increment between vector elements: 10.00.
   Grid increment between vectors        : 10.00.
   Grid origin (x0,y0):  (193000.00,97000.00).

Elevations for station # 1 (station ID 1).
   Interpolated Elevation at Station location:...     191.89
   Actual Elevation at Station location: ........     192.38
   Difference: ..................................      -0.49
   Note: Elevation grid saved for station 1.

Elevations for station # 2 (station ID 2).
   Interpolated Elevation at Station location:...     185.01
   Actual Elevation at Station location: ........     184.61
   Difference: ..................................       0.41
   Note: Elevation grid saved for station 2.

Elevations for station # 3 (station ID 3).
   Interpolated Elevation at Station location:...     179.32
   Actual Elevation at Station location: ........     180.24
   Difference: ..................................      -0.92
   Note: Elevation grid saved for station 3.

Elevations for station # 4 (station ID 4).
   Interpolated Elevation at Station location:...     173.02
   Actual Elevation at Station location: ........     172.69
   Difference: ..................................       0.33
   Note: Elevation grid saved for station 4.

Elevations for station # 5 (station ID 5).
   Interpolated Elevation at Station location:...     175.20
   Actual Elevation at Station location: ........     175.20
   Difference: ..................................       0.01
   Note: Elevation grid saved for station 5.


Results of RasterTC (version 1.30[beta]) Terrain Correction Calculation.  Today is Tue Aug  8 12:43:09 2000
Corrections calculated from a min radius of 10.0 m to a max radius of 2000.0 m.
An intermediate radius of 300.0 m was used (numerical integration between min and intermediate radius).
All terrain corrections below are in milliGals, calculated for a terrain density of 2.670 grams/cc.

 Count StationID   Easting   Northing   Elev, m   Elev, ft  Int Elev, m  Elev Diff  ObsGravity  TC-Near QF-Inner N-Evals Est.Error  TC-Far  QF-Outer  TC-Total
     1 1          195000.0    99000.0    192.38     631.17     191.89       0.49         0.000    0.097      1      1980 2.807e-06   0.069        0     0.166
     2 2          195000.0    99100.0    184.61     605.66     185.01      -0.41         0.000    0.097      0      1980 2.250e-06   0.076        0     0.173
     3 3          195000.0    99200.0    180.24     591.33     179.32       0.92         0.000    0.092      0      1980 1.375e-06   0.080        0     0.172
     4 4          195000.0    99300.0    172.69     566.56     173.02      -0.33         0.000    0.125      0      1980 4.272e-07   0.095        0     0.219
     5 5          195000.0    99400.0    175.20     574.79     175.20      -0.01         0.000    0.132      0      1980 2.216e-07   0.076        0     0.208



Table of Contents


Example of Output File

#Results of RasterTC (version 1.30[beta]) Terrain Correction Calculation.  Today is Tue Aug  8 12:43:09 2000
#Corrections calculated from a min radius of 10.0 m to a max radius of 2000.0 m.
#An intermediate radius of 300.0 m was used (numerical integration between min and intermediate radius).
#All terrain corrections below are in milliGals, calculated for a terrain density of 2.670 grams/cc.
#Count StationID   Easting   Northing   Elev, m   Elev, ft  Int Elev, m  Elev Diff  ObsGravity  TC-Near QF-Inner N-Evals Est.Error  TC-Far  QF-Outer  TC-Total
     1 1          195000.0    99000.0    192.38     631.17     191.89       0.49         0.000    0.097      1      1980 2.807e-06   0.069        0     0.166
     2 2          195000.0    99100.0    184.61     605.66     185.01      -0.41         0.000    0.097      0      1980 2.250e-06   0.076        0     0.173
     3 3          195000.0    99200.0    180.24     591.33     179.32       0.92         0.000    0.092      0      1980 1.375e-06   0.080        0     0.172
     4 4          195000.0    99300.0    172.69     566.56     173.02      -0.33         0.000    0.125      0      1980 4.272e-07   0.095        0     0.219
     5 5          195000.0    99400.0    175.20     574.79     175.20      -0.01         0.000    0.132      0      1980 2.216e-07   0.076        0     0.208

Table of Contents