README.ProcessData
changeset 16 ff157f4445ce
child 17 343ba18738cf
equal deleted inserted replaced
15:7896932f0206 16:ff157f4445ce
       
     1 ======================================================================
       
     2                     R E A D M E . P R O C E S S D A T A 
       
     3                     doc: Tue May 15 18:49:00 2012
       
     4                     dlm: Fri Jun 15 07:24:53 2012
       
     5                     (c) 2012 A.M. Thurnherr
       
     6                     uE-Info: 144 38 NIL 0 0 72 75 2 8 NIL ofnI
       
     7 ======================================================================
       
     8 
       
     9 =Overview=
       
    10 
       
    11 This README describes how to obtain profiles of vertical shear and
       
    12 velocity from CTD/LADCP data. It assumes that all of the required
       
    13 software has been installed (see [README.Install]).
       
    14 
       
    15 The re-implemented shear method software provides two commands:
       
    16 
       
    17 [LADCPproc]	This utility produces LADCP shear data from a raw ADCP
       
    18 		data file and the corresponding CTD time series.
       
    19 		Additionally, it can create profiles of acoustic
       
    20 		backscatter, as well as BT-referenced velocity
       
    21 		profiles near the seabed from downlooking ADCPs.
       
    22 
       
    23 [LADCPintsh]	This utility produces profiles of horizontal velocity
       
    24 		from the [LADCPproc] shear output. BT profiles (from
       
    25 		[LADCPproc] or from the LDEO_IX inversion software) or
       
    26 		SADCP profiles (manually constructed) can be used to
       
    27 		reference the velocity profiles.
       
    28 
       
    29 For non-standard processing, the shear output from [LADCPproc] can be
       
    30 post-edited before gridding, e.g. in order to filter data collected at
       
    31 very shallow depths when the ADCP may be affected by the magnetic field
       
    32 of the surface vessel (see [README.PostEdit] for details).
       
    33 
       
    34 
       
    35 =DATA REQUIREMENTS=
       
    36 
       
    37 ADCP DATA: The software reads binary RDI BB ADCP files from both down-
       
    38 and upward-looking ADCPs. Clock setting of the ADCP is not important.
       
    39 
       
    40 CTD DATA: LADCP processing requires a CTD-derived time series of
       
    41 elapsed time, pressure, temperature and salinity. A time resolution of
       
    42 1Hz is recommended. The software is capable of reading both binary and
       
    43 ASCII SeaBird .cnv files with lat/lon information in the header and
       
    44 with the following fields: timeS, prDM, t090C and/or t190C, sal00
       
    45 and/or sal11. Alternatively, the CTD time series can be supplied as an
       
    46 arbitrary headerless ASCII CTD file with the same information, as
       
    47 described in [LADCPproc.defaults].
       
    48 	
       
    49 
       
    50 =CALCULATE LADCP SHEAR PROFILE=
       
    51 
       
    52 The following simple example shows how to create separate shear profiles
       
    53 from an upward- and a downward-looking ADCP, as well as a BT-referenced
       
    54 velocity profile near the seabed:
       
    55 
       
    56 Input files:
       
    57 	001DL000.000	downlooker ADCP file
       
    58 	001UL000.000	uplooker ADCP file
       
    59 	001.cnv		CTD file
       
    60 
       
    61 LADCPproc -p 001DL.sh -b 001.BT 001DL000.000 001.cnv > /dev/null
       
    62 	- this example creates two files, 001DL.sh (shear profiles) and
       
    63 	  001.BT (bottom-track data)
       
    64 	- the default output (STDOUT) from [LADCPproc] is a list of
       
    65 	  valid shear samples, which is ignored (sent to /dev/null) in
       
    66 	  this example
       
    67 	- it is recommended that the diagnostic output (STDERR) is
       
    68 	  captured in a log file; refer to the manual of your shell on
       
    69 	  how to accomplish this
       
    70 
       
    71 LADCPproc -p 001UL.sh 001UL000.000 001.cnv > /dev/null
       
    72 	- this example creates one file, 001UL.sh (shear profiles)
       
    73 
       
    74 
       
    75 In this simple example, processing is carried out with standard
       
    76 parameters. Some of the important parameters can be modified with
       
    77 [LADCPproc] options, which are listed when [LADCPproc] is ran without
       
    78 input parameters. The following are the most important [LADCPproc]
       
    79 options:
       
    80 	-d 		generate diagnostic output (recommended)
       
    81 	-r		use RDI BT data instead of echo amplitudes to find
       
    82 	        	seabed and determine CTD velocity
       
    83 	-o <dz>		output grid resolution (defaults to 5m)
       
    84 	-p <shearprof>	generate shear profile output
       
    85 	-b <btm_track>	generate BT output
       
    86 	-s <setup_file>	read additional non-default processing parameters
       
    87 			from <setup_file>
       
    88 
       
    89 However, there are many more processing parameters than can be modified
       
    90 with options --- a full list with comments can be found in
       
    91 [LADCPproc.defaults]. To change any of the default parameter values,
       
    92 create a perl-file with variable assignments (see [LADCPproc.defaults]
       
    93 for syntax) and use the -s <setup_file> option in [LADCPproc].
       
    94 
       
    95 
       
    96 =CALCULATE LADCP VELOCITY PROFILE=
       
    97 
       
    98 Given the output from the above steps, different full-depth velocity
       
    99 profiles can be produced as follows:
       
   100 
       
   101 LADCPintsh 001DL.sh > 001DL.bc
       
   102 	- this creates baroclinic (zero vertical mean) velocity profile
       
   103 	  from the DL shear data
       
   104 
       
   105 LADCPintsh -r 001.BT 001DL.sh > 001DL.vel
       
   106 	- this creates a BT-referenced absolute velocity profile from
       
   107 	  the DL shear data
       
   108 
       
   109 LADCPintsh -r 001.BT -u 001UL.sh 001DL.sh > 001.vel
       
   110 	- this creates a BT-referenced absolute velocity profile from
       
   111 	  the combined DL/UL shear data
       
   112 
       
   113 It is also possible to use SADCP data to reference the velocity
       
   114 profiles, although it is up to the user to create an input data file
       
   115 in one of the supported formats. Note that it is *not* possible to use
       
   116 multiple simultaneous referencing constraints with [LADCPintsh]. 
       
   117 
       
   118 The following are common [LADCPintsh] options:
       
   119 	-u		use uplooker shear (in addition to downlooker,
       
   120 			which is always used)
       
   121 	-r <file>	use reference-velocity data to reference baroclinic
       
   122 			velocity profiles; the following file formats
       
   123 			are supported 1) bottom-track output produced by
       
   124 			the -b option of [LADCPproc], 2) bottom-track
       
   125 			output produced by the LDEO processing software
       
   126 			(.bot files). SADCP data can be used, too, but
       
   127 			they have to be supplied in one of the two
       
   128 			supported file formats.
       
   129 	-n <samp>	set minimum number of shear samples to use
       
   130 	-m <samp>	set minimum BT samples to use
       
   131 
       
   132 
       
   133 =QUALITY CHECKS=
       
   134 
       
   135 After processing, the quality of the resulting profiles must be
       
   136 assessed. The following steps are recommended:
       
   137 
       
   138 1) Compare the down- and up-cast profiles of velocity. Vertical
       
   139 velocity is particularly useful in this context as problematic casts
       
   140 often show a striking "X" pattern.
       
   141 
       
   142 2) Inspect the standard deviation profiles of the binned shear and
       
   143 determine (by comparison with similar data) whether the standard
       
   144 deviations have the correct magnitude.
       
   145 
       
   146 3) Calculate and compare independent solutions from the uplooker and
       
   147 downlooker data. This will only validate the baroclinic velocities (i.e.
       
   148 the vertical shear).
       
   149 
       
   150 4) Compare to velocity profiles calculated with different software (e.g.
       
   151 with the LDEO_IX velocity inversion code).