USAGE:  

        1  splist-plot [options] <file.dat> [<file.dat>...]

           eventually with accompanying <file.cfd> after each <file.dat>
           Plots a sequence of spectra (auto-assigned colours)
 
        2  splist-plot [options] + [splist-options] \
                       ++ <file.sp> <TRG1> [<TRG2>]
 plots one spectrum with TRG2 responsible for error bars 

 N.B. There is a plotting section in cal-seis
 The primary purpose of the present script is to
 plot a spectrum with confidence limits
 or plot a sequence of spectra that have been created with
 sasm03 or sasm06.
 
 OPTIONS:
 
  -T               - show commands, don't execute
 
  -L               - plot lines (default is dots)
 
  -yll #y          - crop the lower y-range at #y (default = -1e10)
  -yr #ylo,#yhi    - crop y-axis  (default = -1e10,+1e10)
 
  -x ext           - the file extensions in alt. 1 (default is `dat´ )
      
  -logf            - logarithmic frequency axis
       
 {-pha | -adm [unit] | -coh |-pwr [unit] }  
                   - selects appropriate y-axis title.
                     First three chars after `-' are significant.
              unit - is added to the axis label, `[ ]´ will be provided.
                     For -pha it's always `[deg]´  and for -coh `[-]´ 
                     -gai      is synonymous with -adm
                     -pow -psp  "      "      "   -pwr

  -Freq unit       - is added to the freq-axis title, `[ ]´ will be provided. 
 
  -o <opt>       -2- Send option <opt> to axlim, accumulative.
                     E.g. -o -T0/1  to widen the data range with
                     one tick interval at the upper end
                     
  -O odir        -2- the output directory for splist
 
  -P png-dir       - the directory for resulting png-files     [~/www/4me]

  -Leg leg-file  -2- A file for a legend. Default is no legend.
 
  -E[H] [<trg>]  -2- for error bars, import the confidence
                     data tagged <trg>
                     (you can specify the <trg> alternately at
                     the end of the command line as <trg2>)
                     -EH KOC  - error bars for coherence
                     -E       - for admittance or phase
 
  -EPW -         -2- for error bars on power (the limits
                     are computed in splist); therefore a second
                     tag <trg2> must not be specified.
                   !  The `-' sign after -EPW is required !
 
  -C <opt> + -EC -2- for error bars computed from statistic info
                     specify the option that will generate it
                     (if -EC won't do)     (does anyone understand this??"
                                           The -EPW option is easier to understand.
                                           Better perhaps to have special options
                                           for specific spectal properties.)
 
 LEGEND:
    This is optional in case 2 and skipped in case 1. A file must be provided containing legend information in the order of the arguments.
    Example:
T : Cross-spectral analysis SG-Q330HR / GCF-Acceleration
1 : 3.3
2 : 3.6
3 : 3.7
4 : 3.8
5 : 3.9
6 : 4.0
7 : 4.1
8 : 4.2
    Advantage may be made of file names' parts, like in this example:
/bin/cp -f q-seis-scg.leg tmp/q-seis-scg.leg'
set gainfiles = ( `ls -t tmp/q-gain-seis-scg-*.dat | head -8 | sort` )
ls  $gainfiles | awk '{sub(/.*-/,"",$0);sub(/\.dat/,"",$0);print NR,":",$0}' >> tmp/q-seis-scg.leg
splist-plot -gain dB -Freq Hz -logf -Leg tmp/q-seis-scg.leg -P ~/www/4me/seis-scg-181010/ -yr -35/15 -o -X1,1 -L $gainfiles
   where the template q-seis-scg.leg  contains the title.
.

 EXAMPLE  1:
 
   splist -q -f0 -cpd -EPW -AdB o/g090615-OPNEND-1h-ochy.ra.sp PSP > ! o/urtap-ochy-ra.spf
   splist -q -f0 -cpd -EPW -AdB o/eterna_residuals.sp PSP > ! o/eterna_residuals.spf
   splist-plot -x spf -L -pwr -logf o/urtap-ochy-ra.spf o/eterna_residuals.spf
 
 
 EXAMPLES 2:
 
   splist-plot -gain -P ~/www/4me/cal-seis$SUBDIR -O tmp -E \
               + -Hz -f0 -adB \
               ++ o-110910/zadm-0.01745.sp CAD ADC
 
   splist-plot -pha  -P ~/www/4me/cal-seis$SUBDIR -O tmp -E \
               + -Hz -f0 -Pd \
               ++ o-110910/zpha-0.01745.sp CPH APH
 
   splist-plot -coherence -o -b0 -O tmp -EH KOC \
               + -Hz -f0 ++ o-110910/cohc-0.01745.sp KOH
 
   splist-plot -pow -logf -EPW - + -f0 -Hz -AdB ++ o/psp-hf.sp PSY
 
 Recently:
 
   splist-plot -gain -P ~/www/4me/seis-scg-181010/ -yr -15/15 -O tmp -E \
               + -Hz -f0 -adB ++ o/g-gcf-1s.csp CAD ADC

   where sasm06-seis-scg-1Hz.ins provided the csp-file.