Started Mon Apr 12 15:37:47 2004
########################################################
##      High-Precision Gravity Reduction Program      ##
##            Version 1.4.0 (Jan 2004)                ##
##                   Paul Gettings                    ##
##         Department of Geology & Geophysics         ##
##                University of Utah                  ##
##        http://thermal.gg.utah.edu/~gettings        ##
########################################################

*** READ RAW DATA ***

>>> Reading station data file
--> Processed CG-3m data file /home/gettings/projects/USGS_Nebraska/gravity/2003-07/grav23jul03.txt with 636 data records
>>> Reading station parameter file
->> Removing stations in skip list
--> Removed 1 of 1 stations in skip list (115 data records)
->> Checking for missing stations in parameter file
->> Removing parameter entries with no matching station data
->> Checking station repeat information for missing stations
--> station repeat information has 16 changes
--> Processed station parameter file /home/gettings/projects/USGS_Nebraska/gravity/2003-07/jul.par
>>> Reading tare data from file
!!! No tare data found in /home/gettings/projects/USGS_Nebraska/gravity/2003-07/07-23/grav.out.  Assuming no previous tare data.
--> done

*** VIEW UNCORRECTED DATA ***

>>> skipping view data

*** CORRECTIONS TO RAW DATA ***

>>> Pre-average corrections
->> Remove Longman Earth Tide Correction?
--> removing ETC applied by meter...done.
->> Remove Meter Drift Correction?
--> leaving drift correction applied by meter.
->> Atmospheric Pressure Correction?
--> no atmospheric corrections
->> Tamura Earth Tide Correction?
--> performing earth tide correction...done.
->> Tilt Correction?
--> not correcting meter tip/tilt corrections.

*** QUALITY CONTROL CHECK OF TIME SERIES ***
>>> Checking station time series for s.d. > threshhold=0.050000
->>
->> Reading with s.d. > threshhold; check reading 1 of station 1.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 1.
->>..
->> Reading with s.d. > threshhold; check reading 5 of station 1.
->>..............
->> Reading with s.d. > threshhold; check reading 1 of station 2.
->>....................................................
->>....................
->> Reading with s.d. > threshhold; check reading 1 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 5.
->>.
->> Reading with s.d. > threshhold; check reading 8 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 9 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 5.
->>.
->> Reading with s.d. > threshhold; check reading 13 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 14 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 15 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 16 of station 5.
->>.
->> Reading with s.d. > threshhold; check reading 18 of station 5.
->>
->> Reading with s.d. > threshhold; check reading 19 of station 5.
->>.
->> Reading with s.d. > threshhold; check reading 1 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 7 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 8 of station 6.
->>.
->> Reading with s.d. > threshhold; check reading 10 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 12 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 13 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 14 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 15 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 16 of station 6.
->>..
->> Reading with s.d. > threshhold; check reading 19 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 20 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 21 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 22 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 23 of station 6.
->>
->> Reading with s.d. > threshhold; check reading 24 of station 6.
->>........
->>..................
->> Reading with s.d. > threshhold; check reading 7 of station 8.
->>....
->> Reading with s.d. > threshhold; check reading 12 of station 8.
->>............................
->> Reading with s.d. > threshhold; check reading 1 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 7 of station 10.
->>.
->> Reading with s.d. > threshhold; check reading 9 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 12 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 13 of station 10.
->>...
->> Reading with s.d. > threshhold; check reading 17 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 18 of station 10.
->>
->> Reading with s.d. > threshhold; check reading 1 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 11.
->>.
->> Reading with s.d. > threshhold; check reading 4 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 7 of station 11.
->>.
->> Reading with s.d. > threshhold; check reading 9 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 11.
->>
->> Reading with s.d. > threshhold; check reading 12 of station 11.
->>......
->> Reading with s.d. > threshhold; check reading 19 of station 11.
->>.......................................................
->>.......
->> Reading with s.d. > threshhold; check reading 16 of station 14.
->>..
->> Reading with s.d. > threshhold; check reading 19 of station 14.
->>...............................
->> Reading with s.d. > threshhold; check reading 7 of station 16.
->>.........
->> Reading with s.d. > threshhold; check reading 17 of station 16.
->>...................
->>...............................
->> Reading with s.d. > threshhold; check reading 2 of station 19.
->>.................
->> Reading with s.d. > threshhold; check reading 1 of station 20.
->>......................
->>..................................................................
->> Reading with s.d. > threshhold; check reading 1 of station 24.
->>.
->> Reading with s.d. > threshhold; check reading 3 of station 24.
->>...
->>.................
--> checked 521 data points
>>> Time Series Linear Fits (y=m*x+b), threshhold=0.000027 mGal/sec (97.200000 uGal/hr)
>>>            <------- all data --------> <---- w/o 1st 180 sec ----->
>>> Station ID Linear slope Zero intercept Linear slope Zero intercept Slope (uGal/hr)
>>> ---------- ------------ -------------- ------------ -------------- ---------------
->>         1.    -0.000002    4808.144935    -0.000002    4808.144544            -8.9
->>         2.     0.000005    4808.085623     0.000005    4808.086695            18.1
->>         3.    -0.000004    4808.149525    -0.000001    4808.147195            -4.8
->>         4.    -0.000011    4816.497605    -0.000020    4816.498069           -72.0
->>         5.    -0.000002    4815.315988     0.000007    4815.312865            23.9
->>         6.     0.000027    4810.976613     0.000028    4810.981478           101.0 m>threshhold; detrended
->>         7.    -0.000009    4810.669408    -0.000006    4810.666280           -20.2
->>         8.    -0.000007    4811.117374     0.000006    4811.111861            19.9
->>         9.    -0.000005    4811.955389    -0.000004    4811.954134           -12.9
->>        10.     0.000012    4810.408587     0.000013    4810.410306            47.2
->>        11.    -0.000013    4809.807692    -0.000024    4809.810538           -87.4
->>        12.     0.000003    4807.931590     0.000004    4807.932006            13.2
->>        13.    -0.000004    4816.472408    -0.000002    4816.470582            -6.2
->>        14.    -0.000009    4815.288559    -0.000002    4815.285234            -5.6
->>        15.    -0.000010    4810.985323    -0.000013    4810.984065           -45.9
->>        16.    -0.000015    4810.663630    -0.000017    4810.661393           -60.5
->>        17.    -0.000005    4811.113456    -0.000008    4811.113474           -30.6
->>        18.     0.000004    4811.956394     0.000004    4811.956654            16.1
->>        19.    -0.000007    4810.422588    -0.000010    4810.422106           -36.0
->>        20.    -0.000012    4809.823725    -0.000010    4809.820692           -37.5
->>        21.    -0.000008    4807.950113    -0.000011    4807.949994           -39.0
->>        22.    -0.000012    4816.504379    -0.000013    4816.502818           -47.3
->>        23.    -0.000013    4797.632791    -0.000012    4797.630016           -43.9
->>        24.    -0.000007    4808.144972    -0.000005    4808.143168           -18.6
--> 1 stations detrended.
--> Rerunning with updated data.
>>> Time Series Linear Fits (y=m*x+b), threshhold=0.000027 mGal/sec (97.200000 uGal/hr)
>>>            <------- all data --------> <---- w/o 1st 180 sec ----->
>>> Station ID Linear slope Zero intercept Linear slope Zero intercept Slope (uGal/hr)
>>> ---------- ------------ -------------- ------------ -------------- ---------------
->>         1.    -0.000002    4808.144935    -0.000002    4808.144544            -8.9
->>         2.     0.000005    4808.085623     0.000005    4808.086695            18.1
->>         3.    -0.000004    4808.149525    -0.000001    4808.147195            -4.8
->>         4.    -0.000011    4816.497605    -0.000020    4816.498069           -72.0
->>         5.    -0.000002    4815.315988     0.000007    4815.312865            23.9
->>         6.    -0.000001    4810.982138     0.000000    4810.981478             0.0
->>         7.    -0.000009    4810.669408    -0.000006    4810.666280           -20.2
->>         8.    -0.000007    4811.117374     0.000006    4811.111861            19.9
->>         9.    -0.000005    4811.955389    -0.000004    4811.954134           -12.9
->>        10.     0.000012    4810.408587     0.000013    4810.410306            47.2
->>        11.    -0.000013    4809.807692    -0.000024    4809.810538           -87.4
->>        12.     0.000003    4807.931590     0.000004    4807.932006            13.2
->>        13.    -0.000004    4816.472408    -0.000002    4816.470582            -6.2
->>        14.    -0.000009    4815.288559    -0.000002    4815.285234            -5.6
->>        15.    -0.000010    4810.985323    -0.000013    4810.984065           -45.9
->>        16.    -0.000015    4810.663630    -0.000017    4810.661393           -60.5
->>        17.    -0.000005    4811.113456    -0.000008    4811.113474           -30.6
->>        18.     0.000004    4811.956394     0.000004    4811.956654            16.1
->>        19.    -0.000007    4810.422588    -0.000010    4810.422106           -36.0
->>        20.    -0.000012    4809.823725    -0.000010    4809.820692           -37.5
->>        21.    -0.000008    4807.950113    -0.000011    4807.949994           -39.0
->>        22.    -0.000012    4816.504379    -0.000013    4816.502818           -47.3
->>        23.    -0.000013    4797.632791    -0.000012    4797.630016           -43.9
->>        24.    -0.000007    4808.144972    -0.000005    4808.143168           -18.6
--> 0 stations detrended.
--> Done with quality control checking

*** VIEW PRE_AVERAGE DATA ***

>>> skipping view data

*** REDUCING CORRECTED GRAVITY TIME SERIES ***

>>> Using weighted averages
--> skipping first 3.000 minutes
->> averaging gravity stations with identical station IDs...done.
--> 24 unique stations in data file.
--> 0 records skipped due to bad std. deviations.
--> 0 stations dropped due to bad date (JD = 0)

*** CORRECTIONS TO AVERAGED DATA ***

>>> Post-average Corrections
->> Elevation Change Correction?
--> performing elevation change correction...done.
->> Instrument Drift Correction?
-----------------------
Performing drift correction.
-----------------------
--> Please wait; this may take a few minutes.
--> [POLYNOMIAL DRIFT FUNCTION]
    ORDER:  0 PRE-DRIFT- L1:-0.005593  L2: 0.005628   POST-DRIFT- L1:-0.005593  L2: 0.005628 mGal
	C[00]=-0.000000 
	.
	G=N/A (m == 0.0)
    ORDER:  1 PRE-DRIFT- L1:-0.005593  L2: 0.005628   POST-DRIFT- L1: 0.002238  L2: 0.005928 mGal
	C[00]= 0.006278 C[01]=-0.018913 
	.
	G= 1.000000
    ORDER:  2 PRE-DRIFT- L1:-0.005593  L2: 0.005628   POST-DRIFT- L1: 0.007819  L2: 0.008518 mGal
	C[00]= 0.039655 C[01]=-0.134457 C[02]= 0.045160 
	.
	G=17209.255561
    ORDER:  3 PRE-DRIFT- L1:-0.005593  L2: 0.005628   POST-DRIFT- L1: 0.005308  L2: 0.006089 mGal
	C[00]= 0.124718 C[01]=-0.522635 C[02]= 0.488413 C[03]=-0.138137 
	.
	G=11.285897
--> [STAIRCASE DRIFT FUNCTION]
    STEPS: 23 PRE-DRIFT- L1:-0.005593  L2: 0.005628   POST-DRIFT- L1: 0.004710  L2: 0.003345 mGal
	.
	G=37.297626
--> [LINEAR DRIFT RATES (uGal/day)]
                 1. ->              2. ==>     -111.595
                 1. ->              3. ==>        2.838
                 1. ->             24. ==>       -1.330
                 4. ->             13. ==>     -133.068
                 4. ->             22. ==>       14.264
                 5. ->             14. ==>     -153.282
                 6. ->             15. ==>       -3.169
                 7. ->             16. ==>      -44.683
                 8. ->             17. ==>      -20.145
                 9. ->             18. ==>       24.008
                10. ->             19. ==>       43.089
                11. ->             20. ==>       96.998
                12. ->             21. ==>       77.388
--> Instrument drift correction done.
--> Constants of drift function:
        C[0] ==          -0.056833
        C[1] ==           0.034283
        C[2] ==          -0.029118
        C[3] ==          -0.000470
        C[4] ==          -0.000048
        C[5] ==           0.000386
        C[6] ==          -0.000029
        C[7] ==           0.000006
        C[8] ==          -0.000340
        C[9] ==           0.000125
        C[10] ==           0.000288
        C[11] ==          -0.022218
        C[12] ==          -0.008106
        C[13] ==           0.029270
        C[14] ==          -0.007673
        C[15] ==           0.004752
        C[16] ==           0.008443
        C[17] ==           0.003258
        C[18] ==           0.009809
        C[19] ==          -0.003454
        C[20] ==          -0.009062
        C[21] ==           0.009768
        C[22] ==           0.009768

*** COMPUTE FINAL GRAVITY ***

->> Referencing gravity to station 4.
Final gravity computed relative to 4.

*** OUTPUT REDUCED DATA ***

-----------------------
Select or enter the filename of the output file.  Cancelling the dialog cancels all output.
-----------------------
--> Final gravity written to file /home/gettings/projects/USGS_Nebraska/gravity/2003-07/07-23/grav.out
->> Control file written to reduce.cmd.
Completed Mon Apr 12 15:37:48 2004
