LWplot_TL
changeset 11 9e5eba6992f7
parent 9 6ea152571f00
child 13 2788bf1bf1de
equal deleted inserted replaced
10:7e76fdef8953 11:9e5eba6992f7
     1 #!/bin/ksh
     1 #!/bin/ksh
     2 #======================================================================
     2 #======================================================================
     3 #                    L W P L O T _ T L 
     3 #                    L W P L O T _ T L 
     4 #                    doc: Thu Oct 13 10:51:49 2011
     4 #                    doc: Thu Oct 13 10:51:49 2011
     5 #                    dlm: Wed Apr 18 10:33:35 2012
     5 #                    dlm: Tue Oct 16 13:02:45 2012
     6 #                    (c) 2011 A.M. Thurnherr
     6 #                    (c) 2011 A.M. Thurnherr
     7 #                    uE-Info: 18 0 NIL 0 0 72 0 2 4 NIL ofnI
     7 #                    uE-Info: 20 52 NIL 0 0 72 0 2 4 NIL ofnI
     8 #======================================================================
     8 #======================================================================
     9 
     9 
    10 # HISTORY:
    10 # HISTORY:
    11 #	Oct 13, 2011: - created
    11 #	Oct 13, 2011: - created
    12 #	Oct 14, 2011: - fiddled
    12 #	Oct 14, 2011: - fiddled
    13 #	Oct 17, 2011: - turned into filter
    13 #	Oct 17, 2011: - turned into filter
    14 #				  - added %run_label
    14 #				  - added %run_label
    15 #	Nov 14, 2011: - added file layout check
    15 #	Nov 14, 2011: - added file layout check
    16 #	Mar 12, 2012: - sh -> ksh as this does not work with /bin/sh on MacOSX 10.6.4
    16 #	Mar 12, 2012: - sh -> ksh as this does not work with /bin/sh on MacOSX 10.6.4
    17 #   Apr 12, 2012: - made re-entrant
    17 #   Apr 12, 2012: - made re-entrant
       
    18 #	Oct 12, 2012: - split dc/uc time lagging
       
    19 #				  - reduced number of xtic labels 
       
    20 #	Oct 16, 2012: - renamed elapsed to elapsed.LADCP
    18 
    21 
    19 USAGE="Usage: $0 <eps-file> [in-file]"
    22 USAGE="Usage: $0 <eps-file> [in-file]"
    20 	[ $# -eq 2 ] && exec <"$2" "$0" "$1"
    23 	[ $# -eq 2 ] && exec <"$2" "$0" "$1"
    21 	[ $# -ne 1 ] && { echo $USAGE >&2; exit 1; }
    24 	[ $# -ne 1 ] && { echo $USAGE >&2; exit 1; }
    22 eps_file="$1"
    25 eps_file="$1"
    27     [ ! -t 1 ] && echo "$line"
    30     [ ! -t 1 ] && echo "$line"
    28 	[ -z "$out_basename" ] && out_basename=`expr -- "$line" : '#ANTS#PARAMS#.*out_basename{\([^}]*\)}`
    31 	[ -z "$out_basename" ] && out_basename=`expr -- "$line" : '#ANTS#PARAMS#.*out_basename{\([^}]*\)}`
    29 	[ -z "$run_label" ] && run_label=`expr -- "$line" : '#ANTS#PARAMS#.*run_label{\([^}]*\)}`
    32 	[ -z "$run_label" ] && run_label=`expr -- "$line" : '#ANTS#PARAMS#.*run_label{\([^}]*\)}`
    30 	[ -z "$min_elapsed" ] && min_elapsed=`expr -- "$line" : '#ANTS#PARAMS#.*elapsed.min{\([^}]*\)}`
    33 	[ -z "$min_elapsed" ] && min_elapsed=`expr -- "$line" : '#ANTS#PARAMS#.*elapsed.min{\([^}]*\)}`
    31 	[ -z "$max_elapsed" ] && max_elapsed=`expr -- "$line" : '#ANTS#PARAMS#.*elapsed.max{\([^}]*\)}`
    34 	[ -z "$max_elapsed" ] && max_elapsed=`expr -- "$line" : '#ANTS#PARAMS#.*elapsed.max{\([^}]*\)}`
    32 	[ -z "$best_scan_offset" ] && best_scan_offset=`expr -- "$line" : '#ANTS#PARAMS#.*best_scan_offset{\([^}]*\)}`
    35 	[ -z "$bot_elapsed" ] && bot_elapsed=`expr -- "$line" : '#ANTS#PARAMS#.*elapsed.bot{\([^}]*\)}`
       
    36 	[ -z "$best_scan_offset_dc" ] && best_scan_offset_dc=`expr -- "$line" : '#ANTS#PARAMS#.*best_scan_offset\.dc{\([^}]*\)}`
       
    37 	[ -z "$best_scan_offset_uc" ] && best_scan_offset_uc=`expr -- "$line" : '#ANTS#PARAMS#.*best_scan_offset\.uc{\([^}]*\)}`
    33 	[ -z "$fields" ] && fields=`expr -- "$line" : '#ANTS#FIELDS# \(.*\)' | sed -e s/{//g -e s/}//g`
    38 	[ -z "$fields" ] && fields=`expr -- "$line" : '#ANTS#FIELDS# \(.*\)' | sed -e s/{//g -e s/}//g`
    34 done
    39 done
    35 
    40 
    36 set -- $fields
    41 set -- $fields
    37 [ "$1" = elapsed -a "$2" = scan_offset -a "$4" = downcast ] || {
    42 [ "$1" = elapsed.LADCP -a "$2" = scan_offset -a "$4" = downcast ] || {
    38 	echo "$0: file layout error ($1,$2,$4)" >&2
    43 	echo "$0: file layout error ($1,$2,$4)" >&2
    39 	exit 1
    44 	exit 1
    40 }
    45 }
    41 
    46 
    42 eps_file="$PWD/$eps_file"	# make outfile name absolute (hopefully, it is not already...)
    47 eps_file="$PWD/$eps_file"	# make outfile name absolute (hopefully, it is not already...)
    61 U=-R0/1/0/1
    66 U=-R0/1/0/1
    62 J=-JX10
    67 J=-JX10
    63 
    68 
    64 awk '{if ($4 == 1) { print $1/60,$2; }}' $TMPFILE \
    69 awk '{if ($4 == 1) { print $1/60,$2; }}' $TMPFILE \
    65 	| psxy -P -K $J $R \
    70 	| psxy -P -K $J $R \
    66 		   -Bf1a15:"Elapsed Time [min]":/f1a5:"Best Offset [scans]":WeSn \
    71 		   -Bf1a30:"Elapsed Time [min]":/f1a5:"Best Offset [scans]":WeSn \
    67 		   -Sc0.1 -Gcoral \
    72 		   -Sc0.1 -Gcoral \
    68 	> "$eps_file"		   
    73 	> "$eps_file"		   
    69 	
    74 	
    70 awk '{if ($4 == 0) { print $1/60,$2; }}' $TMPFILE \
    75 awk '{if ($4 == 0) { print $1/60,$2; }}' $TMPFILE \
    71 	| psxy -O -K $J $R -Sc0.1 -GSeaGreen \
    76 	| psxy -O -K $J $R -Sc0.1 -GSeaGreen \
    72 	>> "$eps_file"		   
    77 	>> "$eps_file"		   
    73 	
    78 	
    74 echo 0.02 0.98 12 0 0 TL $out_basename $run_label | pstext -O -K $U $J >> "$eps_file"
    79 echo 0.02 0.98 12 0 0 TL $out_basename $run_label | pstext -O -K $U $J >> "$eps_file"
    75 
    80 
    76 echo -e "`echo $min_elapsed/60|bc` $best_scan_offset\n`echo $max_elapsed/60+1|bc` $best_scan_offset" \
    81 echo -e "`echo $min_elapsed/60|bc` $best_scan_offset_dc\n`echo $bot_elapsed/60+1|bc` $best_scan_offset_dc" \
    77 	| psxy -O $R $J >> "$eps_file"
    82 	| psxy -O -K $R $J -W4/coral >> "$eps_file"
       
    83 echo -e "`echo $bot_elapsed/60|bc` $best_scan_offset_uc\n`echo $max_elapsed/60+1|bc` $best_scan_offset_uc" \
       
    84 	| psxy -O $R $J -W4/SeaGreen >> "$eps_file"
    78 
    85 
    79 rm $TMPFILE
    86 rm $TMPFILE
    80 fixbb "$eps_file"
    87 fixbb "$eps_file"
    81 
    88 
    82 cd "$PWD"
    89 cd "$PWD"