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: Sun Apr 11 06:52:26 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: 62 52 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 |
|
12 # Apr 11, 2021: - halved x-axis range |
12 |
13 |
13 require "$ANTS/libGMT.pl"; |
14 require "$ANTS/libGMT.pl"; |
14 |
15 |
15 sub plotDC($$) |
16 sub plotDC($$) |
16 { |
17 { |
47 my($pfn) = @_; |
48 my($pfn) = @_; |
48 |
49 |
49 my($yellow_light) = 0.004; |
50 my($yellow_light) = 0.004; |
50 my($red_light) = 0.01; |
51 my($red_light) = 0.01; |
51 |
52 |
52 my($xmin) = -0.05; |
53 my($xmin) = -0.01; |
53 my($xmax) = 0.05; |
54 my($xmax) = 0.01; |
54 my($ymin) = round(antsParam('min_depth')-25,50); |
55 my($ymin) = round(antsParam('min_depth')-25,50); |
55 my($ymax) = ($P{water_depth} > 0) ? |
56 my($ymax) = ($P{water_depth} > 0) ? |
56 round($P{water_depth}+25,50) : |
57 round($P{water_depth}+25,50) : |
57 round($P{max_depth}+$P{ADCP_bin_length}+25,50); |
58 round($P{max_depth}+$P{ADCP_bin_length}+25,50); |
58 |
59 |
59 my($R) = "-R$xmin/$xmax/$ymin/$ymax"; |
60 my($R) = "-R$xmin/$xmax/$ymin/$ymax"; |
60 my($depth_tics) = ($ymax < 1000 ) ? 'f10a100g100' : 'f100a500g500'; |
61 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"); |
62 GMT_begin($pfn,'-JX10/-10',$R,"-P -Bf0.001a0.005g0.005:'Residual Vertical Velocity [m/s]':/$depth_tics:'Depth [m]':WeSn"); |
62 |
63 |
63 GMT_psxy('-W2,CornflowerBlue'); # zero line |
64 GMT_psxy('-W2,CornflowerBlue'); # zero line |
64 printf(GMT "0 $ymin\n0 $ymax\n"); |
65 printf(GMT "0 $ymin\n0 $ymax\n"); |
65 |
66 |
66 GMT_psxy('-W1,coral,8_2:0'); my($r12DC) = plotDC('MEAN_RESIDUAL12',$opt_k); # dc residual12 (pitch plane) |
67 GMT_psxy('-W1,coral,8_2:0'); my($r12DC) = plotDC('MEAN_RESIDUAL12',$opt_k); # dc residual12 (pitch plane) |