plot_mean_residuals.pl
author Andreas Thurnherr <ant@ldeo.columbia.edu>
Sat, 10 Apr 2021 06:00:45 -0400
changeset 55 2d8e1139acd5
parent 54 828e5466391b
permissions -rw-r--r--
adapted to Nortek files
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     1
#======================================================================
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     2
#                    P L O T _ M E A N _ R E S I D U A L S . P L 
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     3
#                    doc: Tue Jul 28 13:21:09 2015
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
     4
#                    dlm: Fri May 15 19:06:51 2020
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     5
#                    (c) 2015 A.M. Thurnherr
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
     6
#                    uE-Info: 122 47 NIL 0 0 72 2 2 4 NIL ofnI
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     7
#======================================================================
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     8
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     9
# HISTORY:
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    10
#	Jul 28, 2015: - created from [LWplot_BR]
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    11
#	Jul 29, 2015: - finished
31
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    12
#	Jul 30, 2015: - added bin_tics
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    13
#				  - added outGrid_* support
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    14
#	Jan 22, 2015: - many changes
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    15
#				  - added quality assessment label
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    16
#	Jan 25, 2016: - added return on no data
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
    17
#	Mar 16, 2016: - adapted to gmt5
45
6d49c7420a6c stable V1.3beta, just before implementing bin interpolation
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 41
diff changeset
    18
#   May 18, 2016: - added version
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    19
#	May 15, 2020: - adapted to bin-residuals separate per beam pair
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    20
#				  - added orange range
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    21
#				  - slightly relaxed green range
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    22
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    23
require "$ANTS/libGMT.pl";
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    24
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    25
sub plot_mean_residuals($)
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    26
{
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    27
	my($pfn) = @_;
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    28
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    29
	return unless ($P{BR_max_bin});
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    30
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    31
	my($xmin) = -0.05;
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    32
	my($xmax) =  0.05;
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    33
	my($ymin) =  0.5;
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    34
	my($ymax) = $P{BR_max_bin} + 0.5;
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    35
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    36
	my($R) = "-R$xmin/$xmax/$ymin/$ymax";
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    37
	GMT_begin($pfn,'-JX10/-10',$R,'-P');
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    38
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    39
	if ($outGrid_firstBin>$LADCP_firstBin || $outGrid_lastBin<$LADCP_lastBin) {		# mark used bins
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
    40
		GMT_psxy('-G200 -L');
31
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    41
		printf(GMT ">\n%f %f\n%f %f\n%f %f\n%f %f\n",
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    42
			$xmin,$LADCP_firstBin-0.5,
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    43
			$xmax,$LADCP_firstBin-0.5,
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    44
			$xmax,$outGrid_firstBin-0.5,
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    45
			$xmin,$outGrid_firstBin-0.5)
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    46
				if ($outGrid_firstBin>$LADCP_firstBin);
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    47
		printf(GMT ">\n%f %f\n%f %f\n%f %f\n%f %f\n",
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    48
			$xmin,$LADCP_lastBin+0.5,
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    49
			$xmax,$LADCP_lastBin+0.5,
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    50
			$xmax,$outGrid_lastBin+0.5,
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    51
			$xmin,$outGrid_lastBin+0.5)
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    52
				if ($outGrid_lastBin<$LADCP_lastBin);
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    53
	}
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
    54
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    55
	GMT_psxy('-W1');																# plot zero line
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    56
	printf(GMT "0 $ymin\n0 $ymax\n");
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    57
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    58
	GMT_psxy('-W2,coral,8_2:0');
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    59
		for (my($bin)=$LADCP_firstBin; $bin<@dc_bres12; $bin++) {			
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    60
			next if ($bin+1<$outGrid_firstBin || $bin+1>$outGrid_lastBin);
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    61
			next unless ($dc_bres12_nsamp[$bin] >= $dc_bres12_max_nsamp/3);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    62
			printf(GMT "%f %d\n",$dc_avg_bres12[$bin],$bin+1);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    63
        }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    64
	GMT_psxy('-W2,coral,1_1:0');
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    65
		for (my($bin)=$LADCP_firstBin; $bin<@dc_bres34; $bin++) {			
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    66
			next if ($bin+1<$outGrid_firstBin || $bin+1>$outGrid_lastBin);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    67
			next unless ($dc_bres34_nsamp[$bin] >= $dc_bres34_max_nsamp/3);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    68
			printf(GMT "%f %d\n",$dc_avg_bres34[$bin],$bin+1);
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    69
        }
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
    70
	GMT_psxy('-Ex0.2/2,coral');
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    71
		for (my($bin)=$LADCP_firstBin-1; $bin<@dc_bres12; $bin++) {
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    72
			printf(GMT "%f %d %f\n",
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    73
							$dc_avg_bres12[$bin],
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    74
							$bin+1,
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    75
							($dc_bres12_nsamp[$bin] > 1) ?
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    76
								$dc_sig_bres12[$bin]/sqrt($dc_bres12_nsamp[$bin]-1) : 0);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    77
		}
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    78
		for (my($bin)=$LADCP_firstBin-1; $bin<@dc_bres34; $bin++) {
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    79
			printf(GMT "%f %d %f\n",
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    80
							$dc_avg_bres34[$bin],
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    81
							$bin+1,
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    82
							($dc_bres34_nsamp[$bin] > 1) ?
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    83
								$dc_sig_bres34[$bin]/sqrt($dc_bres34_nsamp[$bin]-1) : 0);
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    84
        }
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    85
	GMT_psxy('-W2,SeaGreen,8_2:0');
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    86
		for (my($bin)=$LADCP_firstBin; $bin<@uc_bres12; $bin++) {			
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
    87
			next if ($bin+1<$outGrid_firstBin || $bin+1>$outGrid_lastBin);
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    88
			next unless ($uc_bres12_nsamp[$bin] >= $uc_bres12_max_nsamp/3);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    89
			printf(GMT "%f %d\n",$uc_avg_bres12[$bin],$bin+1);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    90
        }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    91
	GMT_psxy('-W2,SeaGreen,1_1:0');
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    92
		for (my($bin)=$LADCP_firstBin; $bin<@uc_bres34; $bin++) {			
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    93
			next if ($bin+1<$outGrid_firstBin || $bin+1>$outGrid_lastBin);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    94
			next unless ($uc_bres34_nsamp[$bin] >= $uc_bres34_max_nsamp/3);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    95
			printf(GMT "%f %d\n",$uc_avg_bres34[$bin],$bin+1);
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    96
        }
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
    97
	GMT_psxy('-Ex0.2/2,SeaGreen');
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
    98
		for (my($bin)=$LADCP_firstBin-1; $bin<@uc_bres12; $bin++) {
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    99
			printf(GMT "%f %d %f\n",
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   100
							$uc_avg_bres12[$bin],
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   101
							$bin+1,
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   102
							($uc_bres12_nsamp[$bin] > 1) ?
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   103
								$uc_sig_bres12[$bin]/sqrt($uc_bres12_nsamp[$bin]-1) : 0);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   104
		}
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   105
		for (my($bin)=$LADCP_firstBin-1; $bin<@uc_bres34; $bin++) {
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   106
			printf(GMT "%f %d %f\n",
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   107
							$uc_avg_bres34[$bin],
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   108
							$bin+1,
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   109
							($uc_bres34_nsamp[$bin] > 1) ?
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   110
								$uc_sig_bres34[$bin]/sqrt($uc_bres34_nsamp[$bin]-1) : 0);
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   111
        }
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   112
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   113
	GMT_unitcoords();																	# LABELS
45
6d49c7420a6c stable V1.3beta, just before implementing bin interpolation
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 41
diff changeset
   114
	GMT_pstext('-F+f9,Helvetica,orange+jTR -N -Gwhite');
6d49c7420a6c stable V1.3beta, just before implementing bin interpolation
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 41
diff changeset
   115
        print(GMT "0.99 0.01 V$VERSION\n");
6d49c7420a6c stable V1.3beta, just before implementing bin interpolation
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 41
diff changeset
   116
        
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   117
	GMT_pstext('-F+f14,Helvetica,blue+jBL -N');											# profile id
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   118
		print(GMT "0.0 -0.03 $P{out_basename} $P{run_label}\n");
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
   119
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   120
	GMT_pstext('-F+f12,Helvetica-Bold,black+jBL -Gwhite'); 							# rms residuals
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   121
		print(GMT "0.1 0.88 beams <1,2>\n");
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   122
	GMT_pstext('-F+f12,Helvetica-Bold,black+jBL -Gwhite');
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   123
		print(GMT "0.7 0.88 beams <3,4>\n");
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   124
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   125
	GMT_pstext('-F+f12,Helvetica,coral+jBL');											# rms residuals
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   126
		print(GMT "0.01 0.93 dc\n");
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   127
	if ($dc_bres12_rms >= 0.005) { 		GMT_pstext('-F+f12,Helvetica,white+jBL -Gred'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   128
	elsif ($dc_bres12_rms >= 0.003) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gorange'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   129
	elsif ($dc_bres12_rms >= 0.0015) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gyellow'); }
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   130
	else {								GMT_pstext('-F+f12,Helvetica,black+jBL -Ggreen'); }
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   131
		printf(GMT "0.10 0.93 %.1f mm/s rms\n",1000*$dc_bres12_rms);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   132
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   133
	if ($dc_bres34_rms >= 0.005) { 		GMT_pstext('-F+f12,Helvetica,white+jBL -Gred'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   134
	elsif ($dc_bres34_rms >= 0.003) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gorange'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   135
	elsif ($dc_bres34_rms >= 0.0015) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gyellow'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   136
	else {								GMT_pstext('-F+f12,Helvetica,black+jBL -Ggreen'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   137
		printf(GMT "0.70 0.93 %.1f mm/s rms\n",1000*$dc_bres34_rms);
34
e550db661c17 pre-Tampa
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 31
diff changeset
   138
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   139
	GMT_pstext('-F+f12,Helvetica,SeaGreen+jBL');
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   140
		print(GMT "0.01 0.98 uc\n");
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   141
	if ($uc_bres12_rms >= 0.005) { 		GMT_pstext('-F+f12,Helvetica,white+jBL -Gred'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   142
	elsif ($uc_bres12_rms >= 0.003) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gorange'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   143
	elsif ($uc_bres12_rms >= 0.0015) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gyellow'); }
41
6bddb82924e3 V1.2beta5
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 34
diff changeset
   144
	else {								GMT_pstext('-F+f12,Helvetica,black+jBL -Ggreen'); }
54
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   145
		printf(GMT "0.10 0.98 %.1f mm/s rms\n",1000*$uc_bres12_rms);
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   146
	if ($uc_bres34_rms >= 0.005) { 		GMT_pstext('-F+f12,Helvetica,white+jBL -Gred'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   147
	elsif ($uc_bres34_rms >= 0.003) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gorange'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   148
	elsif ($uc_bres34_rms >= 0.0015) { 	GMT_pstext('-F+f12,Helvetica,black+jBL -Gyellow'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   149
	else {								GMT_pstext('-F+f12,Helvetica,black+jBL -Ggreen'); }
828e5466391b before A20
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 45
diff changeset
   150
		printf(GMT "0.70 0.98 %.1f mm/s rms\n",1000*$uc_bres34_rms);
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   151
31
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
   152
	my($bin_tics) = ($ymax <= 20) ? 'f1a1' : 'f1a2';
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   153
	GMT_setR($R);																		# FINISH PLOT
31
d0ae3cb99021 DoMORE-2
A.M. Thurnherr <athurnherr@yahoo.com>
parents: 30
diff changeset
   154
	GMT_end("-Bf0.005a0.02:'Residual Vertical Velocity [m/s]':/$bin_tics:'Bin [#]':WeSn");
30
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   155
}
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   156
7fb67e771d85 LWplots expunged
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
   157
1; # return true on require