equal
deleted
inserted
replaced
1 #====================================================================== |
1 #====================================================================== |
2 # P L O T _ C O R R E L A T I O N . P L |
2 # P L O T _ C O R R E L A T I O N . P L |
3 # doc: Tue Jul 28 13:21:09 2015 |
3 # doc: Tue Jul 28 13:21:09 2015 |
4 # dlm: Tue Jan 26 20:46:18 2016 |
4 # dlm: Wed Mar 16 16:26:09 2016 |
5 # (c) 2015 A.M. Thurnherr |
5 # (c) 2015 A.M. Thurnherr |
6 # uE-Info: 19 31 NIL 0 0 72 0 2 4 NIL ofnI |
6 # uE-Info: 67 37 NIL 0 0 72 0 2 4 NIL ofnI |
7 #====================================================================== |
7 #====================================================================== |
8 |
8 |
9 # HISTORY: |
9 # HISTORY: |
10 # Jul 28, 2015: - created from [LWplot_corr] |
10 # Jul 28, 2015: - created from [LWplot_corr] |
11 # Jan 26, 2016: - added return on no data to plot |
11 # Jan 26, 2016: - added return on no data to plot |
|
12 # Mar 16, 2016: - adapted to gmt5 |
12 |
13 |
13 require "$ANTS/libGMT.pl"; |
14 require "$ANTS/libGMT.pl"; |
14 |
15 |
15 sub plot_correlation($) |
16 sub plot_correlation($) |
16 { |
17 { |
51 GMT_psxy('-G204/153/102'); |
52 GMT_psxy('-G204/153/102'); |
52 print(GMT "$xmin $ymax\n$xmax $ymax\n$xmax $P{water_depth}\n $xmin $P{water_depth}\n"); |
53 print(GMT "$xmin $ymax\n$xmax $ymax\n$xmax $P{water_depth}\n $xmin $P{water_depth}\n"); |
53 } |
54 } |
54 |
55 |
55 GMT_unitcoords(); # LABELS |
56 GMT_unitcoords(); # LABELS |
56 GMT_pstext(-Gblue); |
57 GMT_pstext('-F+f14,Helvetica,blue+jTL -N'); |
57 print(GMT "0.02 0.98 12 0 0 BL $P{out_basename} $P{run_label}\n"); |
58 print(GMT "0.01 -0.06 $P{out_basename} $P{run_label}\n"); |
58 |
59 |
59 my($depth_tics) = ($ymax < 1000 ) ? 'f10a100' : 'f100a500'; # AXES |
60 my($depth_tics) = ($ymax < 1000 ) ? 'f10a100' : 'f100a500'; # AXES |
60 my($ens_tics) = ($ymax < 1000 ) ? 'f50a500' : 'f500a2000'; |
61 my($ens_tics) = ($ymax < 1000 ) ? 'f50a500' : 'f500a2000'; |
61 GMT_setR($R); |
62 GMT_setR($R); |
62 GMT_psbasemap("-B$ens_tics:'Ensemble [#]':/$depth_tics:'Depth [m]':WeSn"); |
63 GMT_psbasemap("-B$ens_tics:'Ensemble [#]':/$depth_tics:'Depth [m]':WeSn"); |
63 |
64 |
64 GMT_setAnnotFontSize(7); # SCALE BAR |
65 GMT_setAnnotFontSize(7); # SCALE BAR |
65 GMT_psscale("-E -D8/2/3/0.4 $C -B/:corr:"); |
66 # GMT_psscale("-E -D8/2/3/0.4 $C -B/:corr:"); |
|
67 GMT_psscale("-Dn0.85/0.1+w3/0.4+e $C -B/:corr:"); |
66 |
68 |
67 GMT_end(); # FINISH PLOT |
69 GMT_end(); # FINISH PLOT |
68 } |
70 } |
69 |
71 |
70 1; # return true on require |
72 1; # return true on require |