LWplot_corr
changeset 6 4d48ffde2471
parent 5 509cc9966b68
child 7 e4f531c66f2c
child 8 efdd251dd5b7
--- a/LWplot_corr	Thu Oct 27 10:29:18 2011 -0400
+++ b/LWplot_corr	Fri Mar 09 08:34:30 2012 +0000
@@ -2,14 +2,15 @@
 #======================================================================
 #					 L W P L O T _ C O R R 
 #                    doc: Sat Oct 15 13:42:50 2011
-#                    dlm: Mon Oct 24 10:29:12 2011
+#                    dlm: Mon Nov 14 14:58:54 2011
 #                    (c) 2011 A.M. Thurnherr
-#                    uE-Info: 13 0 NIL 0 0 72 2 2 4 NIL ofnI
+#                    uE-Info: 75 26 NIL 0 0 72 2 2 4 NIL ofnI
 #======================================================================
 
 # HISTORY:
 #	Oct 19, 2011: - adapted from [LWplot_residuals]
 #   Oct 24, 2011: - LADCP_bin_length => ADCP_bin_length
+#	Nov 14, 2011: - added file layout check
 
 #--------------------------------------------------
 # Usage
@@ -18,6 +19,7 @@
 USAGE="Usage: $0 <eps-file> [in-file]"
 	[ $# -eq 2 ] && exec <"$2" "$0" "$1"
 	[ $# -ne 1 ] && { echo $USAGE >&2; exit 1; }
+eps_file="$1"
 
 #--------------------------------------------------
 # Read Header Data
@@ -37,6 +39,12 @@
 	[ -z "$fields" ] && fields=`expr -- "$line" : '#ANTS#FIELDS# \(.*\)' | sed -e s/{//g -e s/}//g`
 done
 
+set -- $fields
+[ "$1" = ensemble -a "$4" = depth -a "${14}" = correlation ] || {
+	echo "$0: file layout error" >&2
+	exit 1
+}
+
 #--------------------------------------------------
 # Plot Data
 #--------------------------------------------------
@@ -64,10 +72,10 @@
 ens_width=`echo "scale=5;10/($max_ens-$min_ens+1)"|bc`
 bin_length=`echo "scale=5;10*$ADCP_bin_length/($max_depth-$min_depth+$ADCP_bin_length)"|bc`
 
-awk "{print \$1, \$4, \$12, $ens_width, $bin_length}" $TMPFILE \
-	| psxy -P -K $R $X $C -Sr > "$1"
+awk "{print \$1, \$4, \$14, $ens_width, $bin_length}" $TMPFILE \
+	| psxy -P -K $R $X $C -Sr > "$eps_file"
 	
-echo 0.02 0.02 12 0 0 TL $out_basename $run_label | pstext -O -K $U $X >> "$1"
+echo 0.02 0.02 12 0 0 TL $out_basename $run_label | pstext -O -K $U $X >> "$eps_file"
 
 if [ 0 -eq `echo "($max_depth-$min_depth)>1000"|bc` ]
 then
@@ -78,9 +86,9 @@
 	ens_tics=f500a1000
 fi
 
-cat `which LWplot_corr | sed 's@LWplot_corr$@corr_scale.skel@'` >> "$1"
+cat `which LWplot_corr | sed 's@LWplot_corr$@corr_scale.skel@'` >> "$eps_file"
 
-psbasemap -O $R $X -B$ens_tics:"Ensemble":/$depth_tics:"Depth [m]":WeSn >> "$1"
+psbasemap -O $R $X -B$ens_tics:"Ensemble":/$depth_tics:"Depth [m]":WeSn >> "$eps_file"
 
 rm $TMPFILE
-fixbb "$1"
+fixbb "$eps_file"