0
|
1 |
#======================================================================
|
|
2 |
# M E R G E C T D + L A D C P . T I M E M A T C H I N G P R O B L E M S
|
|
3 |
# doc: Sun May 23 11:43:53 2010
|
|
4 |
# dlm: Sat Jun 5 12:09:43 2010
|
|
5 |
# (c) 2010 A.M. Thurnherr
|
|
6 |
# uE-Info: 23 22 NIL 0 0 72 0 2 4 NIL ofnI
|
|
7 |
#======================================================================
|
|
8 |
|
|
9 |
$sumCdt = $NCdt = $trgSumCdt = 0; # init stuff
|
|
10 |
|
|
11 |
sub detectTimeMatchingProblems($$$$)
|
|
12 |
{
|
|
13 |
|
|
14 |
local($ens,$depth,$otimoff,$timoff) = @_;
|
|
15 |
local($adtimoff) = $otimoff - $timoff;
|
|
16 |
|
|
17 |
warning("WARNING: Jump in CTD time offset at %dm (#$ens, %.1fs, %d CTD scans)",
|
|
18 |
$depth,$adtimoff,abs($adtimoff/$opt_z))
|
|
19 |
if (abs($adtimoff) >= 1);
|
|
20 |
|
|
21 |
if ($depth > 30) {
|
|
22 |
warning("WARNING: Jump in CTD time offset at %dm (#$ens, %.1fs, %d CTD scans)",
|
|
23 |
$depth,$adtimoff,abs($adtimoff/$opt_z))
|
|
24 |
if ($adtimoff >= 0.1);
|
|
25 |
$sumCdt += $adtimoff;
|
|
26 |
$trgSumCdt += 0.02;
|
|
27 |
$NCdt++;
|
|
28 |
if ($NCtd > 10) {
|
|
29 |
if ($sumCdt >= 5.0*$trgSumCdt) {
|
|
30 |
croak("$0: bad time match\n");
|
|
31 |
} elsif ($sumCdt >= 1.5*$trgSumCdt && !$warned) {
|
|
32 |
warning("WARNING: dodgy time match at ens #$ens --- mean absolute time-offset jigger is %g",$sumCdt/$NCdt);
|
|
33 |
}
|
|
34 |
}
|
|
35 |
}
|
|
36 |
}
|
|
37 |
|
|
38 |
1;
|