plot_wprof.pl
changeset 47 2ccb81b7cea5
parent 46 cc6c4309828a
child 49 5006e9158207
equal deleted inserted replaced
46:cc6c4309828a 47:2ccb81b7cea5
     1 #======================================================================
     1 #======================================================================
     2 #                    P L O T _ W P R O F . P L 
     2 #                    P L O T _ W P R O F . P L 
     3 #                    doc: Sun Jul 26 11:08:50 2015
     3 #                    doc: Sun Jul 26 11:08:50 2015
     4 #                    dlm: Tue May 24 22:31:14 2016
     4 #                    dlm: Thu May 26 11:29:32 2016
     5 #                    (c) 2015 A.M. Thurnherr
     5 #                    (c) 2015 A.M. Thurnherr
     6 #                    uE-Info: 19 51 NIL 0 0 72 2 2 4 NIL ofnI
     6 #                    uE-Info: 20 0 NIL 0 0 72 2 2 4 NIL ofnI
     7 #======================================================================
     7 #======================================================================
     8 
     8 
     9 # HISTORY:
     9 # HISTORY:
    10 #	Jul 26, 2015: - created from LWplot_prof_2beam
    10 #	Jul 26, 2015: - created from LWplot_prof_2beam
    11 #	Jul 30, 2015: - moved main label outside plot area
    11 #	Jul 30, 2015: - moved main label outside plot area
    15 #	May 18, 2016: - replaced rms tilt by mean tilt with traffic background
    15 #	May 18, 2016: - replaced rms tilt by mean tilt with traffic background
    16 #				  - added plot_wprof_ymin tweakable
    16 #				  - added plot_wprof_ymin tweakable
    17 #	May 24, 2016: - BUG: ymin did not work for nsamp
    17 #	May 24, 2016: - BUG: ymin did not work for nsamp
    18 #				  - fixed for partial-depth profiles
    18 #				  - fixed for partial-depth profiles
    19 #				  - suppress plotting of nsamp == 0
    19 #				  - suppress plotting of nsamp == 0
       
    20 #	May 26, 2016: - added instrument coord system to plot labels
    20 
    21 
    21 # Tweakables:
    22 # Tweakables:
    22 #
    23 #
    23 # $plot_wprof_xmin = -0.27;
    24 # $plot_wprof_xmin = -0.27;
    24 # $plot_wprof_ymin = 3600;
    25 # $plot_wprof_ymin = 3600;
   137 			print(GMT "0.02 0.10 BT\n");
   138 			print(GMT "0.02 0.10 BT\n");
   138 	}
   139 	}
   139 
   140 
   140 	GMT_pstext('-F+f9,Helvetica,CornFlowerBlue+jTL -N');
   141 	GMT_pstext('-F+f9,Helvetica,CornFlowerBlue+jTL -N');
   141 		printf(GMT "0.64 1.020 $LADCP{BEAM_FREQUENCY}kHz $LADCP{INSTRUMENT_TYPE} $P{ADCP_orientation}\n");
   142 		printf(GMT "0.64 1.020 $LADCP{BEAM_FREQUENCY}kHz $LADCP{INSTRUMENT_TYPE} $P{ADCP_orientation}\n");
   142 		printf(GMT "0.64 1.055 bin setup\n		0.77 1.055 : %.1fm/%1.fm/%1.fm\n",
   143 #		printf(GMT "0.64 1.055 %s\n		0.77 1.055 : %.1fm/%1.fm/%1.fm\n",
       
   144 		printf(GMT "0.64 1.055 %s [%.1fm/%1.fm/%1.fm]\n",
       
   145 			$LADCP{BEAM_COORDINATES} ? 'beam vels' : 'Earth vels',
   143 			$LADCP{BLANKING_DISTANCE},$LADCP{TRANSMITTED_PULSE_LENGTH},$LADCP{BIN_LENGTH});
   146 			$LADCP{BLANKING_DISTANCE},$LADCP{TRANSMITTED_PULSE_LENGTH},$LADCP{BIN_LENGTH});
   144 		print(GMT "0.64 1.090 mean tilt\n 		0.77 1.096 :\n");
   147 		print(GMT "0.64 1.090 mean tilt\n 		0.77 1.096 :\n");
   145 		print(GMT "0.64 1.130 rms a\@-pkg\@-\n	0.77 1.1315 :\n");
   148 		print(GMT "0.64 1.130 rms a\@-pkg\@-\n	0.77 1.1315 :\n");
   146 
   149 
   147 	if ($P{dc_mean_tilt} < 4) {
   150 	if ($P{dc_mean_tilt} < 4) {