[teqc] another request for teqc +qc

Lou Estey lou at unavco.org
Mon Aug 4 14:31:02 MDT 2008


All,

Brief update, and your next-to-last chance for comments before this goes into
testing before release:

> Bug #1:
> To address Herb Dragert's first question (below), I still have to fiddle
> a bit so that the hrs on the new SUM line will give the expected result
> if there are one or more gaps of epochs in between the first and last
> data epochs in the window.

Fixed.

> Bug #2:
> If I do a qc on the full 24-hours of data at this site:
> 
> [1142] teqc +qcq full.obs 2> /dev/null | grep ^SUM
> SUM 08  1 28 00:00 08  1 28 23:59 24.00  15  48599  43273  89  0.90  0.79    231
> 
> you'll notice that the expected number of observations is 48599, rather
> than the 48583 observations determined above for the 3-hour file.

Fixed.

> Bug #3:
> The SV count on the "+<cutoff>" line is wrong in the windowed case, e.g.
> 
> [1144] teqc -st 00:00:00 -e 23:59:45 +qcq part.obs 2> /dev/null | grep ^+10
> +10|                              999999888888888888888888888888888888888888|+10
> 
> compared to:
> [1145] teqc +qcq full.obs 2> /dev/null | grep ^+10
> +10|aa99bba9999999988778999988789999999988898888867788888aa99bbbcca988aabbba|+10 

Not fixed (yet), but altered so that the line will populated in the time bins
bounded by where there is data and blank elsewhere, i.e.

[1424] teqc -st 00:00:00 -e 23:59:45 +qcq part.obs 2> /dev/null | grep ^+10
+10|                              999999888                                 |+10

The test case I've been using is file part.obs, which has data from 10:00:00 to
10:59:45 every 15 seconds and from 12:00:00 to 12:59:45 every 15 seconds.  This
is an edited subset of full.obs, which has complete data from 00:00:00 to 23:59:45
every 15 seconds.  Results of various runs, initially all using a "part.nav" which
is a copy of full.nav -- all the possible navigation messages from this site:

normal "qc full" run of the complete 24-hrs of data:
[1415] teqc +qcq full.obs 2>/dev/null | grep ^S
SUM 08  1 28 00:00 08  1 28 23:59 24.00  15  48599  43273  89  0.90  0.79    231

same as above, except windowed to the full extent of the data:
[1416] teqc -st 00:00:00 -e 23:59:45 +qcq full.obs 2>/dev/null | grep ^S
SUM 08  1 28 00:00 08  1 28 23:59 24.00  15  48599  43273  89  0.90  0.79    231

Results are identical.

normal "qc full" run of the 2-hrs of data spread over 3 hours:
[1417] teqc +qcq part.obs 2>/dev/null | grep ^S
SUM 08  1 28 10:00 08  1 28 12:59 2.000  15   5986   3754  63  0.74  0.69    375

Note the "hrs" value now shows "2.000" (rather than "2.996" as with the previous
qc results of teqc).

24-hour window of the 2-hrs of data:
[1418] teqc -st 00:00:00 -e 23:59:45 +qcq part.obs 2>/dev/null | grep ^S
SUM 08  1 28 10:00 08  1 28 12:59 2.000  15   5986   3754  63  0.74  0.69    375
SWN 08  1 28 00:00 08  1 28 23:59 24.00  15  48599   3754   8  0.74  0.69    375

The SUM line is identical to the previous non-windowed run.  The SWN line
gives the window bounds, the expected span of data over that window, and the
expected number of observations over that window, and the computed % observations
over that window that are present in the data.

+ssv results:
No change was made to this.

qc lite:

After removing the nav file part.gps -- thus forcing a "qc lite" run (i.e.
no SV or antenna positions can be calculated):

[1420] rm part.gps

no window specified:
[1421] teqc +qcq part.obs 2>/dev/null | grep ^S
SUM 08  1 28 10:00 08  1 28 12:59 2.000  15    n/a   3754 n/a  0.74  0.69    375

The start/end times on the SUM line show the bounds of the data.

24-hour window specified:
[1422] teqc -st 00:00:00 -e 23:59:45 +qcq part.obs 2>/dev/null | grep ^S
SUM 08  1 28 00:00 08  1 28 23:59 2.000  15    n/a   3754 n/a  0.74  0.69    375

The start/ends times on the SUM line show the bound of the window.

Question: should the last "qc lite" case _instead_ yield something like:

[1422] teqc -st 00:00:00 -e 23:59:45 +qcq part.obs 2>/dev/null | grep ^S
SUM 08  1 28 10:00 08  1 28 12:59 2.000  15    n/a   3754 n/a  0.74  0.69    375
SWN 08  1 28 00:00 08  1 28 23:59 24.00  15    n/a   3754 n/a  0.74  0.69    375

If so, then the meaning of the start/end values would be consistent across qc full
and qc lite runs for all SUM lines.

--lou

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Louis H. Estey, Ph.D.              office:  [+001] 303-381-7456
UNAVCO, 6350 Nautilus Drive           FAX:  [+001] 303-381-7451
Boulder, CO  80301-5554            e-mail:  lou  unavco.org
    WWW:  http://www.unavco.org   http://jules.unavco.org

"If the universe is the answer, what is the question?"
                                                -- Leon Lederman
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


More information about the teqc mailing list