LWplot_prof_2beam
changeset 18 8818acdcd587
parent 14 fea65697bc7b
child 20 d6dc9c9da138
equal deleted inserted replaced
17:fc83e436a800 18:8818acdcd587
     1 #!/bin/ksh
     1 #!/bin/ksh
     2 #======================================================================
     2 #======================================================================
     3 #					 L W P L O T _ P R O F _ 2 B E A M 
     3 #					 L W P L O T _ P R O F _ 2 B E A M 
     4 #                    doc: Fri Oct 14 09:42:36 2011
     4 #                    doc: Fri Oct 14 09:42:36 2011
     5 #                    dlm: Thu Nov 21 10:22:15 2013
     5 #                    dlm: Mon Nov  3 22:01:45 2014
     6 #                    (c) 2011 A.M. Thurnherr
     6 #                    (c) 2011 A.M. Thurnherr
     7 #                    uE-Info: 82 0 NIL 0 0 72 2 2 4 NIL ofnI
     7 #                    uE-Info: 13 59 NIL 0 0 72 2 2 4 NIL ofnI
     8 #======================================================================
     8 #======================================================================
     9 
     9 
    10 # HISTORY:
    10 # HISTORY:
    11 #	May 15, 2013: - created from [LWplot_prof]
    11 #	May 15, 2013: - created from [LWplot_prof]
    12 #	Oct 30, 2103: - got rid of non-portable echo -e
    12 #	Oct 30, 2103: - got rid of non-portable echo -e
       
    13 #	Nov  3, 2014: - adapted to updated layout of .prof file
    13 
    14 
    14 # NOTES:
    15 # NOTES:
    15 #	- this version plots 2-beam solutions instead of final w
    16 #	- this version plots 2-beam solutions instead of final w
    16 
    17 
    17 #--------------------------------------------------
    18 #--------------------------------------------------
    39 done
    40 done
    40 
    41 
    41 set -- $fields
    42 set -- $fields
    42 [ "$1" = depth -a "$7" = dc_w12 -a "$8" = dc_w34 -a "$5" = dc_w.mad -a "$6" = dc_w.nsamp -a \
    43 [ "$1" = depth -a "$7" = dc_w12 -a "$8" = dc_w34 -a "$5" = dc_w.mad -a "$6" = dc_w.nsamp -a \
    43   "${14}" = uc_w12 -a "${15}" = uc_w34 -a "${12}" = uc_w.mad -a "${13}" = uc_w.nsamp -a \
    44   "${14}" = uc_w12 -a "${15}" = uc_w34 -a "${12}" = uc_w.mad -a "${13}" = uc_w.nsamp -a \
    44   "${20}" = BT_w  -a "${21}" = BT_w.mad -a "${22}" = BT_w.nsamp ] || {
    45   "${16}" = BT_w  -a "${17}" = BT_w.mad -a "${18}" = BT_w.nsamp ] || {
    45 		echo "$0: file layout error" >&2
    46 		echo "$0: file layout error" >&2
    46 		exit 1
    47 		exit 1
    47 }
    48 }
    48 
    49 
    49 #--------------------------------------------------
    50 #--------------------------------------------------
    83 # VERTICAL VELOCITIES (2-BEAM SOLUTIONS)
    84 # VERTICAL VELOCITIES (2-BEAM SOLUTIONS)
    84 awk '{print $7, $1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,coral,6_2:0 >> "$eps_file"
    85 awk '{print $7, $1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,coral,6_2:0 >> "$eps_file"
    85 awk '{print $8, $1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,coral,4_6:0 >> "$eps_file"
    86 awk '{print $8, $1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,coral,4_6:0 >> "$eps_file"
    86 awk '{print $14,$1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,SeaGreen,6_2:0 >> "$eps_file"
    87 awk '{print $14,$1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,SeaGreen,6_2:0 >> "$eps_file"
    87 awk '{print $15,$1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,SeaGreen,4_6:0 >> "$eps_file"
    88 awk '{print $15,$1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,SeaGreen,4_6:0 >> "$eps_file"
    88 awk '{print $20,$1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,black >> "$eps_file"
    89 awk '{print $16,$1}' $TMPFILE | psxy -O -K -N -Mn $R $X -W4,black >> "$eps_file"
    89 	
    90 	
    90 # MEAN ABSOLUTE DEVIATIONS (COMBINED SOLUTION)
    91 # MEAN ABSOLUTE DEVIATIONS (COMBINED SOLUTION)
    91 awk '{print  $5,$1, $4}' $TMPFILE | grep -vi nan | psxy -O -K -N  $R $X -Sc0.1c -Gcoral >> "$eps_file"
    92 awk '{print  $5,$1, $4}' $TMPFILE | grep -vi nan | psxy -O -K -N  $R $X -Sc0.1c -Gcoral >> "$eps_file"
    92 awk '{print $12,$1,$11}' $TMPFILE | grep -vi nan | psxy -O -K -N  $R $X -Sc0.1c -GSeaGreen >> "$eps_file"
    93 awk '{print $12,$1,$11}' $TMPFILE | grep -vi nan | psxy -O -K -N  $R $X -Sc0.1c -GSeaGreen >> "$eps_file"
    93 awk '{print $21,$1,$20}' $TMPFILE | grep -vi nan | psxy -O -K -N  $R $X -Sc0.1c -Gblack >> "$eps_file"
    94 awk '{print $17,$1,$20}' $TMPFILE | grep -vi nan | psxy -O -K -N  $R $X -Sc0.1c -Gblack >> "$eps_file"
    94 
    95 
    95 # NUMBER OF SAMPLES (COMBINED SOLUTION)
    96 # NUMBER OF SAMPLES (COMBINED SOLUTION)
    96 awk '{print  $6,$1, $4}' $TMPFILE | sed '/nan/s/.*/nan/' | psxy -O -K -N -Mn $R2 $X -W1/coral >> "$eps_file"
    97 awk '{print  $6,$1, $4}' $TMPFILE | sed '/nan/s/.*/nan/' | psxy -O -K -N -Mn $R2 $X -W1/coral >> "$eps_file"
    97 awk '{print $13,$1,$11}' $TMPFILE | sed '/nan/s/.*/nan/' | psxy -O -K -N -Mn $R2 $X -W1/SeaGreen >> "$eps_file"
    98 awk '{print $13,$1,$11}' $TMPFILE | sed '/nan/s/.*/nan/' | psxy -O -K -N -Mn $R2 $X -W1/SeaGreen >> "$eps_file"
    98 awk '{print $22,$1,$20}' $TMPFILE | sed '/nan/s/.*/nan/' | psxy -O -K -N -Mn $R2 $X -W1/black >> "$eps_file"
    99 awk '{print $18,$1,$20}' $TMPFILE | sed '/nan/s/.*/nan/' | psxy -O -K -N -Mn $R2 $X -W1/black >> "$eps_file"
    99 
   100 
   100 # LABELS
   101 # LABELS
   101 echo 0.02 0.02 12 0 0 TL $out_basename $run_label | pstext -O -K $U $X >> "$eps_file"
   102 echo 0.02 0.02 12 0 0 TL $out_basename $run_label | pstext -O -K $U $X >> "$eps_file"
   102 echo 0.6 0.98 12 0 0 BR m.a.d. | pstext -O -K $U $X >> "$eps_file"
   103 echo 0.6 0.98 12 0 0 BR m.a.d. | pstext -O -K $U $X >> "$eps_file"
   103 
   104