1 #====================================================================== |
1 #====================================================================== |
2 # P L O T _ R E S I D U A L _ P R O F S . P L |
2 # P L O T _ R E S I D U A L _ P R O F S . P L |
3 # doc: Wed May 18 18:43:33 2016 |
3 # doc: Wed May 18 18:43:33 2016 |
4 # dlm: Tue May 24 22:02:28 2016 |
4 # dlm: Thu Jul 1 13:25:20 2021 |
5 # (c) 2016 A.M. Thurnherr |
5 # (c) 2016 A.M. Thurnherr |
6 # uE-Info: 77 0 NIL 0 0 72 2 2 4 NIL ofnI |
6 # uE-Info: 54 42 NIL 0 0 72 2 2 4 NIL ofnI |
7 #====================================================================== |
7 #====================================================================== |
8 |
8 |
9 # HISTORY: |
9 # HISTORY: |
10 # May 18, 2016: - created from [plot_mean_residuals.pl] |
10 # May 18, 2016: - created from [plot_mean_residuals.pl] |
11 # May 24, 2016: - improved |
11 # May 24, 2016: - improved |
49 my($yellow_light) = 0.004; |
49 my($yellow_light) = 0.004; |
50 my($red_light) = 0.01; |
50 my($red_light) = 0.01; |
51 |
51 |
52 my($xmin) = -0.05; |
52 my($xmin) = -0.05; |
53 my($xmax) = 0.05; |
53 my($xmax) = 0.05; |
54 my($ymin) = round(antsParam('min_depth')-25,50); |
54 my($ymin) = round(antsParam('depth.min')-25,50); |
55 my($ymax) = ($P{water_depth} > 0) ? |
55 my($ymax) = ($P{water_depth} > 0) ? |
56 round($P{water_depth}+25,50) : |
56 round($P{water_depth}+25,50) : |
57 round($P{max_depth}+$P{ADCP_bin_length}+25,50); |
57 round($P{'depth.max'}+$P{ADCP_bin_length}+25,50); |
58 |
58 |
59 my($R) = "-R$xmin/$xmax/$ymin/$ymax"; |
59 my($R) = "-R$xmin/$xmax/$ymin/$ymax"; |
60 my($depth_tics) = ($ymax < 1000 ) ? 'f10a100g100' : 'f100a500g500'; |
60 my($depth_tics) = ($ymax < 1000 ) ? 'f10a100g100' : 'f100a500g500'; |
61 GMT_begin($pfn,'-JX10/-10',$R,"-P -Bf0.005a0.02g0.01:'Residual Vertical Velocity [m/s]':/$depth_tics:'Depth [m]':WeSn"); |
61 GMT_begin($pfn,'-JX10/-10',$R,"-P -Bf0.005a0.02g0.01:'Residual Vertical Velocity [m/s]':/$depth_tics:'Depth [m]':WeSn"); |
62 |
62 |
73 print(GMT "0.99 0.01 V$VERSION\n"); |
73 print(GMT "0.99 0.01 V$VERSION\n"); |
74 GMT_pstext('-F+f14,Helvetica,blue+jBL -N'); # profile id |
74 GMT_pstext('-F+f14,Helvetica,blue+jBL -N'); # profile id |
75 print(GMT "0.0 -0.03 $P{out_basename} $P{run_label}\n"); |
75 print(GMT "0.0 -0.03 $P{out_basename} $P{run_label}\n"); |
76 |
76 |
77 GMT_pstext('-F+f12,Helvetica-Bold,black+jBL -Gwhite'); # rms residuals |
77 GMT_pstext('-F+f12,Helvetica-Bold,black+jBL -Gwhite'); # rms residuals |
78 print(GMT "0.01 0.89 beams <1,2>\n"); |
78 print(GMT "0.1 0.89 beams <1,2>\n"); |
|
79 GMT_pstext('-F+f12,Helvetica-Bold,black+jBR -Gwhite'); |
|
80 print(GMT "0.99 0.89 beams <3,4>\n"); |
|
81 GMT_pstext('-F+f12,Helvetica,coral+jBL -Gwhite'); # rms residuals |
|
82 print(GMT "0.01 0.93 dc\n"); |
|
83 GMT_pstext('-F+f12,Helvetica,SeaGreen+jBL -Gwhite'); |
|
84 print(GMT "0.01 0.98 uc\n"); |
|
85 |
79 if ($r12DC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBL -Gred'); } |
86 if ($r12DC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBL -Gred'); } |
80 elsif ($r12DC >= $yellow_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBL -Gyellow'); } |
87 elsif ($r12DC >= $yellow_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBL -Gyellow'); } |
81 else { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBL -Gwhite'); } |
88 else { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBL -Gwhite'); } |
82 printf(GMT "0.01 0.935 %.1f mm/s rms\n",1000*$r12DC); |
89 printf(GMT "0.1 0.935 %.1f mm/s rms\n",1000*$r12DC); |
83 if ($r12UC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBL -Gred'); } |
90 if ($r12UC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBL -Gred'); } |
84 elsif ($r12UC >= $yellow_light) { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBL -Gyellow'); } |
91 elsif ($r12UC >= $yellow_light) { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBL -Gyellow'); } |
85 else { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBL -Gwhite'); } |
92 else { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBL -Gwhite'); } |
86 printf(GMT "0.01 0.98 %.1f mm/s rms\n",1000*$r12UC); |
93 printf(GMT "0.1 0.98 %.1f mm/s rms\n",1000*$r12UC); |
87 GMT_pstext('-F+f12,Helvetica-Bold,black+jBR -Gwhite'); |
94 |
88 print(GMT "0.99 0.89 beams <3,4>\n"); |
|
89 if ($r34DC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBR -Gred'); } |
95 if ($r34DC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBR -Gred'); } |
90 elsif ($r34DC >= $yellow_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBR -Gyellow'); } |
96 elsif ($r34DC >= $yellow_light) { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBR -Gyellow'); } |
91 else { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBR -Gwhite'); } |
97 else { GMT_pstext('-F+f12,Helvetica-Bold,coral+jBR -Gwhite'); } |
92 printf(GMT "0.99 0.935 %.1f mm/s rms\n",1000*$r34DC); |
98 printf(GMT "0.99 0.935 %.1f mm/s rms\n",1000*$r34DC); |
93 if ($r34UC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBR -Gred'); } |
99 if ($r34UC >= $red_light) { GMT_pstext('-F+f12,Helvetica-Bold,SeaGreen+jBR -Gred'); } |