README.TimeLagging
author A.M. Thurnherr <athurnherr@yahoo.com>
Fri, 06 Mar 2015 15:51:10 -0500
changeset 32 9b972ce37c3b
parent 17 343ba18738cf
permissions -rw-r--r--
merged 1
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
17
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     1
======================================================================
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     2
                    R E A D M E . T I M E L A G G I N G 
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     3
                    doc: Fri Oct 19 10:08:19 2012
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     4
                    dlm: Fri Oct 19 12:13:59 2012
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     5
                    (c) 2012 A.M. Thurnherr
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     6
                    uE-Info: 106 0 NIL 0 0 72 3 2 8 NIL ofnI
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     7
======================================================================
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     8
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     9
=Introduction=
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    10
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    11
In order to derive velocity profiles the data from the CTD and LADCP
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    12
instruments need to be merged. This is accomplished by calculating lag
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    13
correlations between the two corresponding time series of vertical
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    14
velocities calculated from the two instruments. In this software, the
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    15
time lagging is accomplished WITHOUT regard of the clock time reported
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    16
by the instruments, i.e. the instrument clocks do not have to be
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    17
synchronized. Instead of clock time, elapsed time in seconds is used. In
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    18
case of the CTD data, an elapsed time field can be supplied by the user
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    19
(see [README.ProcessData]); in case of the LADCP data, the
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    20
"elapsed-time" field is calculated by the software. The "elapsed-time"
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    21
fields in the processing output are always consistent with the CTD
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    22
elapsed times. While the time-lagging algorithm implemented in the
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    23
software is fairly robust, it has been known to fail. Possible reasons
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    24
include:
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    25
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    26
1) CTD PRESSURE SPIKES. Significant pressure spikes must be removed
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    27
prior to processing, *without* adding or removing CTD time-series
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    28
records.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    29
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    30
2) LACK OF SURFACE VESSEL MOTION. If there is no surface-wave motion
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    31
affecting the vessel, time lagging is much more difficult. In rare
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    32
cases, time lagging must be carried out manually (see below).
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    33
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    34
3) MISSING CTD SCANS. For SeaBird 911 systems, if the connection
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    35
between the CTD and the deck box is not clean CTD scans will be
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    36
dropped. For the software, this looks like the CTD clock running faster
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    37
than the ADCP clock. There are cases where the CTD clock appears to
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    38
have gained more than 5 seconds during a 2000m-deep cast.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    39
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    40
4) MULTIPLE CTD FILES. When CTD acquisition is restarted during a cast,
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    41
multiple files are created. In order to process the LADCP data from such
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    42
a cast, a CTD time-series file without any missing records must be
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    43
constructed manually.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    44
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    45
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    46
=Solving Time-Lagging Problems=
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    47
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    48
While there are several run-time options that can be used to help the
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    49
time-lagging algorithm, detailed knowledge of the algorithm is required
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    50
to understand when and how to use these options, i.e. the user is
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    51
referred to the code and comments in [LADCPproc.bestLag]. However, the
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    52
following method can always be used to solve time-lagging problems, as
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    53
long as the CTD time series does not have any gaps.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    54
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    55
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    56
-Step 1: Produce and Plot a Combined CTD/LADCP Time-Series File-
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    57
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    58
This is accomplished by processing the data with the "-l 0" option and
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    59
using "-t <time-series file>" to produce the file. Plot the resulting
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    60
time series of CTD_w and LADCP_w in the same panel. The plot should
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    61
show immediately whether there are problems with the CTD pressure data
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    62
(spikes). Often, standard processing works after setting any bad
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    63
pressure values to nan.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    64
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    65
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    66
-Step 2: Manually Determine an Approximate Time Lag-
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    67
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    68
Use the output file generated in step 1 to determine how many seconds
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    69
have to be added to the elapsed field when plotting LADCP_w to bring
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    70
the two time series into approximate (a few seconds accuracy)
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    71
agreement. Often, the data can now be processed normally by using "-i
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    72
<estimated lag>".
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    73
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    74
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    75
-Step 3: Manually Determine an Accurate Time Lag-
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    76
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    77
If preprocessing with the -i option still does not succeed, time lagging
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    78
must be carried out manually. If this happens, there is most likely a
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    79
serious problem with either the CTD or LADCP data that should be solved
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    80
before proceeding. This is done exactly as in step 2 but to higher
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    81
accuracy (as high as you can). Once the best lag has been determined
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    82
manually, the data can be reprocessed with the "-l <manually determined
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    83
lag>" option.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    84
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    85
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    86
After solving any time-lagging problems the results should be checked by
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    87
creating a time-series file (with -t) during final processing and
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    88
overplotting the LADCP_w and CTD_w time series. If there is still a
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    89
visible lag between the time series time lagging was not carried out
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    90
correctly. 
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    91
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    92
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    93
=Patching Together CTD Time-Series Files=
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    94
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    95
The LADCP processing software requires the CTD data to be supplied as a
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    96
single time series file with a constant sampling interval. When CTD
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    97
data acquisition is restarted during a cast, multiple files are
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    98
produced. The resulting files cannot simply be pasted together because
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    99
the resulting time series would have gaps. The only way to solve this
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   100
problem is to determine separate time lags for each of the CTD files
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   101
manually (using the method described above). The difference between the
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   102
resulting time lags is equal to the length of the gap between the two
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   103
files. The user can now create a dummy (all fields set to nan) CTD file
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   104
with required number of records that must be added between the CTD
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   105
files to create a single continuous regularly-space time series.
343ba18738cf V1.0beta8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   106
Fractional seconds can be ignored.