equal
deleted
inserted
replaced
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: Sat May 23 11:11:42 2020 |
4 # dlm: Tue Mar 23 08:28:12 2021 |
5 # (c) 2015 A.M. Thurnherr |
5 # (c) 2015 A.M. Thurnherr |
6 # uE-Info: 89 108 NIL 0 0 72 0 2 4 NIL ofnI |
6 # uE-Info: 25 39 NIL 0 0 72 0 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 |
20 # May 26, 2016: - added instrument coord system to plot labels |
20 # May 26, 2016: - added instrument coord system to plot labels |
21 # Mar 20, 2018: - BUG: units of vertical package acceleration were wrong |
21 # Mar 20, 2018: - BUG: units of vertical package acceleration were wrong |
22 # - added blue background for likely in-ice package accelerations |
22 # - added blue background for likely in-ice package accelerations |
23 # May 16, 2020: - added residual profile data to background |
23 # May 16, 2020: - added residual profile data to background |
24 # May 23, 2020: - BUG: windows without samples made program bomb |
24 # May 23, 2020: - BUG: windows without samples made program bomb |
|
25 # Mar 23, 2021: - BUG: instrument frequency was rounded to 100kHz |
25 |
26 |
26 # Tweakables: |
27 # Tweakables: |
27 # |
28 # |
28 # $plot_wprof_xmin = -0.27; |
29 # $plot_wprof_xmin = -0.27; |
29 # $plot_wprof_ymin = 3600; |
30 # $plot_wprof_ymin = 3600; |
229 print(GMT "0.02 0.98 b.track\n"); |
230 print(GMT "0.02 0.98 b.track\n"); |
230 } |
231 } |
231 |
232 |
232 GMT_pstext('-F+f9,Helvetica,CornFlowerBlue+jTL -N'); |
233 GMT_pstext('-F+f9,Helvetica,CornFlowerBlue+jTL -N'); |
233 printf(GMT "0.64 1.020 %d kHz $LADCP{INSTRUMENT_TYPE} $P{ADCP_orientation}\n", |
234 printf(GMT "0.64 1.020 %d kHz $LADCP{INSTRUMENT_TYPE} $P{ADCP_orientation}\n", |
234 round($LADCP{BEAM_FREQUENCY},100)); |
235 round($LADCP{BEAM_FREQUENCY},50)); |
235 printf(GMT "0.64 1.055 %s [%.1fm/%1.fm/%1.fm]\n", |
236 printf(GMT "0.64 1.055 %s [%.1fm/%1.fm/%1.fm]\n", |
236 $LADCP{BEAM_COORDINATES} ? 'beam vels' : 'Earth vels', |
237 $LADCP{BEAM_COORDINATES} ? 'beam vels' : 'Earth vels', |
237 $LADCP{BLANKING_DISTANCE},$LADCP{TRANSMITTED_PULSE_LENGTH},$LADCP{BIN_LENGTH}); |
238 $LADCP{BLANKING_DISTANCE},$LADCP{TRANSMITTED_PULSE_LENGTH},$LADCP{BIN_LENGTH}); |
238 print(GMT "0.64 1.090 mean tilt\n 0.77 1.096 :\n"); |
239 print(GMT "0.64 1.090 mean tilt\n 0.77 1.096 :\n"); |
239 print(GMT "0.64 1.130 rms a\@-pkg\@-\n 0.77 1.1315 :\n"); |
240 print(GMT "0.64 1.130 rms a\@-pkg\@-\n 0.77 1.1315 :\n"); |