Started Mon Apr 12 15:37:48 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/grav24jul03.txt with 501 data records
>>> Reading station parameter file
->> Removing stations in skip list
--> Removed 0 of 1 stations in skip list (0 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-24/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 25.
->>............................
->> Reading with s.d. > threshhold; check reading 30 of station 25.
->>........
->> Reading with s.d. > threshhold; check reading 39 of station 25.
->>.................................
->>............
->> Reading with s.d. > threshhold; check reading 1 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 28.
->>.
->> Reading with s.d. > threshhold; check reading 8 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 9 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 12 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 13 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 14 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 15 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 16 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 17 of station 28.
->>
->> Reading with s.d. > threshhold; check reading 1 of station 29.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 29.
->>.
->> Reading with s.d. > threshhold; check reading 4 of station 29.
->>........................
->> Reading with s.d. > threshhold; check reading 2 of station 30.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 30.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 30.
->>............
->>.......
->> Reading with s.d. > threshhold; check reading 1 of station 31.
->>.
->> Reading with s.d. > threshhold; check reading 3 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 7 of station 31.
->>.
->> Reading with s.d. > threshhold; check reading 9 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 31.
->>..
->> Reading with s.d. > threshhold; check reading 13 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 14 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 15 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 16 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 17 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 18 of station 31.
->>
->> Reading with s.d. > threshhold; check reading 1 of station 32.
->>.
->> Reading with s.d. > threshhold; check reading 3 of station 32.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 32.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 32.
->>...........
->> Reading with s.d. > threshhold; check reading 17 of station 32.
->>.....
->> Reading with s.d. > threshhold; check reading 1 of station 33.
->>.............
->> Reading with s.d. > threshhold; check reading 15 of station 33.
->>
->> Reading with s.d. > threshhold; check reading 16 of station 33.
->>.
->> Reading with s.d. > threshhold; check reading 18 of station 33.
->>..
->> Reading with s.d. > threshhold; check reading 21 of station 33.
->>....
->>......................
->> Reading with s.d. > threshhold; check reading 4 of station 35.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 35.
->>.............................................
->> Reading with s.d. > threshhold; check reading 1 of station 37.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 37.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 37.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 37.
->>.
->> Reading with s.d. > threshhold; check reading 6 of station 37.
->>.........
->> Reading with s.d. > threshhold; check reading 1 of station 38.
->>................
->> Reading with s.d. > threshhold; check reading 1 of station 39.
->>.
->> Reading with s.d. > threshhold; check reading 3 of station 39.
->>.....
->> Reading with s.d. > threshhold; check reading 9 of station 39.
->>.....
->> Reading with s.d. > threshhold; check reading 15 of station 39.
->>..
->> Reading with s.d. > threshhold; check reading 18 of station 39.
->>...
->> Reading with s.d. > threshhold; check reading 1 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 4 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 7 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 8 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 9 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 12 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 13 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 14 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 15 of station 40.
->>
->> Reading with s.d. > threshhold; check reading 1 of station 41.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 41.
->>....
->> Reading with s.d. > threshhold; check reading 7 of station 41.
->>.................
->> Reading with s.d. > threshhold; check reading 1 of station 42.
->>..
->> Reading with s.d. > threshhold; check reading 4 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 5 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 6 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 7 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 8 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 9 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 10 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 11 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 12 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 13 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 14 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 15 of station 42.
->>....
->> Reading with s.d. > threshhold; check reading 20 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 21 of station 42.
->>.
->> Reading with s.d. > threshhold; check reading 23 of station 42.
->>
->> Reading with s.d. > threshhold; check reading 1 of station 43.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 43.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 43.
->>..
->> Reading with s.d. > threshhold; check reading 6 of station 43.
->>.
->> Reading with s.d. > threshhold; check reading 8 of station 43.
->>.........
->> Reading with s.d. > threshhold; check reading 18 of station 43.
->>..............
->>....................................................
->> Reading with s.d. > threshhold; check reading 1 of station 47.
->>
->> Reading with s.d. > threshhold; check reading 2 of station 47.
->>
->> Reading with s.d. > threshhold; check reading 3 of station 47.
->>...
->> Reading with s.d. > threshhold; check reading 7 of station 47.
->>..........
--> checked 501 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)
>>> ---------- ------------ -------------- ------------ -------------- ---------------
->>        25.    -0.000001    4808.170196    -0.000000    4808.169963            -1.3
->>        26.     0.000001    4797.636758     0.000002    4797.636787             6.1
->>        27.    -0.000008    4816.522431    -0.000009    4816.521203           -31.0
->>        28.    -0.000015    4808.679559    -0.000010    4808.675787           -35.8
->>        29.    -0.000005    4809.620627    -0.000003    4809.618184           -10.8
->>        30.    -0.000006    4809.969069    -0.000005    4809.967448           -18.2
->>        31.    -0.000013    4809.234244     0.000003    4809.227362             9.5
->>        32.    -0.000028    4810.660523    -0.000027    4810.655065           -98.2 m>threshhold; detrended
->>        33.    -0.000008    4810.232137    -0.000008    4810.230521           -28.9
->>        34.    -0.000014    4809.957699    -0.000015    4809.955691           -54.9
->>        35.    -0.000003    4810.225089    -0.000003    4810.224300            -9.5
->>        36.    -0.000004    4816.476936    -0.000003    4816.475681           -10.1
->>        37.    -0.000011    4808.624414    -0.000011    4808.622419           -39.5
->>        38.    -0.000008    4809.580989    -0.000015    4809.581288           -54.8
->>        39.     0.000011    4809.927767     0.000012    4809.929558            42.3
->>        40.    -0.000018    4809.219458     0.000029    4809.202630           103.4 m>threshhold; detrended
->>        41.    -0.000019    4810.638926    -0.000019    4810.635196           -69.7
->>        42.     0.000017    4810.216746    -0.000005    4810.229594           -17.0
->>        43.     0.000004    4809.953561     0.000001    4809.955481             2.4
->>        44.    -0.000009    4810.237716    -0.000012    4810.237440           -41.9
->>        45.     0.000004    4816.497647     0.000001    4816.499778             3.3
->>        46.    -0.000006    4797.624887    -0.000008    4797.624318           -28.2
->>        47.     0.000008    4808.139351     0.000005    4808.141516            19.5
--> 2 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)
>>> ---------- ------------ -------------- ------------ -------------- ---------------
->>        25.    -0.000001    4808.170196    -0.000000    4808.169963            -1.3
->>        26.     0.000001    4797.636758     0.000002    4797.636787             6.1
->>        27.    -0.000008    4816.522431    -0.000009    4816.521203           -31.0
->>        28.    -0.000015    4808.679559    -0.000010    4808.675787           -35.8
->>        29.    -0.000005    4809.620627    -0.000003    4809.618184           -10.8
->>        30.    -0.000006    4809.969069    -0.000005    4809.967448           -18.2
->>        31.    -0.000013    4809.234244     0.000003    4809.227362             9.5
->>        32.    -0.000000    4810.655203    -0.000000    4810.655065            -0.0
->>        33.    -0.000008    4810.232137    -0.000008    4810.230521           -28.9
->>        34.    -0.000014    4809.957699    -0.000015    4809.955691           -54.9
->>        35.    -0.000003    4810.225089    -0.000003    4810.224300            -9.5
->>        36.    -0.000004    4816.476936    -0.000003    4816.475681           -10.1
->>        37.    -0.000011    4808.624414    -0.000011    4808.622419           -39.5
->>        38.    -0.000008    4809.580989    -0.000015    4809.581288           -54.8
->>        39.     0.000011    4809.927767     0.000012    4809.929558            42.3
->>        40.    -0.000046    4809.225058     0.000000    4809.202630             0.0
->>        41.    -0.000019    4810.638926    -0.000019    4810.635196           -69.7
->>        42.     0.000017    4810.216746    -0.000005    4810.229594           -17.0
->>        43.     0.000004    4809.953561     0.000001    4809.955481             2.4
->>        44.    -0.000009    4810.237716    -0.000012    4810.237440           -41.9
->>        45.     0.000004    4816.497647     0.000001    4816.499778             3.3
->>        46.    -0.000006    4797.624887    -0.000008    4797.624318           -28.2
->>        47.     0.000008    4808.139351     0.000005    4808.141516            19.5
--> 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.
--> 23 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.020018  L2: 0.008220   POST-DRIFT- L1:-0.020018  L2: 0.008220 mGal
	C[00]=-0.000000 
	.
	G=N/A (m == 0.0)
    ORDER:  1 PRE-DRIFT- L1:-0.020018  L2: 0.008220   POST-DRIFT- L1: 0.011133  L2: 0.008268 mGal
	C[00]= 0.036659 C[01]=-0.120256 
	.
	G= 1.000000
    ORDER:  2 PRE-DRIFT- L1:-0.020018  L2: 0.008220   POST-DRIFT- L1: 0.004426  L2: 0.004925 mGal
	C[00]= 0.379443 C[01]=-1.690532 C[02]= 1.462407 
	.
	G=23.396814
    ORDER:  3 PRE-DRIFT- L1:-0.020018  L2: 0.008220   POST-DRIFT- L1:-0.000915  L2: 0.002971 mGal
	C[00]= 0.165024 C[01]=-0.367525 C[02]=-1.042724 C[03]= 1.550023 
	.
	G=227.129755
--> [STAIRCASE DRIFT FUNCTION]
    STEPS: 22 PRE-DRIFT- L1:-0.020018  L2: 0.008220   POST-DRIFT- L1: 0.000000  L2: 0.000000 mGal
	.
	G=39.179645
--> [LINEAR DRIFT RATES (uGal/day)]
                25. ->             47. ==>      -58.787
                26. ->             46. ==>      -32.941
                27. ->             36. ==>     -265.509
                27. ->             45. ==>      -56.999
                28. ->             37. ==>     -276.504
                29. ->             38. ==>     -198.866
                30. ->             39. ==>     -164.438
                31. ->             40. ==>     -102.150
                32. ->             41. ==>     -137.463
                33. ->             42. ==>      -25.107
                34. ->             43. ==>       16.544
                35. ->             44. ==>       52.095
--> Instrument drift correction done.
--> Constants of drift function:
        C[0] ==          -0.006815
        C[1] ==           0.003272
        C[2] ==           0.004269
        C[3] ==          -0.009107
        C[4] ==          -0.003750
        C[5] ==          -0.007531
        C[6] ==           0.002001
        C[7] ==          -0.012152
        C[8] ==          -0.005176
        C[9] ==          -0.004497
        C[10] ==          -0.008026
        C[11] ==          -0.006513
        C[12] ==           0.006863
        C[13] ==           0.001506
        C[14] ==           0.005287
        C[15] ==          -0.004245
        C[16] ==           0.009908
        C[17] ==           0.002933
        C[18] ==           0.002253
        C[19] ==           0.005782
        C[20] ==           0.003272
        C[21] ==          -0.006815

*** COMPUTE FINAL GRAVITY ***

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

*** 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-24/grav.out
->> Control file written to reduce.cmd.
Completed Mon Apr 12 15:37:49 2004
