error bar confusion
"Gery ." <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
Dear community,
this seems to be an old, discussed (solved perhaps) question, but some results are new for me. After using psxy to plot asymmetric error bars, just one bar is different in each case, ie., the last one (first place in the images), here my 2 cases:
output: case1.jpg
code:
echo "49 -0.155 1.12 49 -0.155 -1.43" | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
echo "37 -1.215 -0.89 37 -1.215 -1.54" | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
echo "23 -1.955 -1.43 23 -1.955 -2.48" | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
echo "16 -1.205 -0.31 16 -1.205 -2.1" | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
echo "12 2.08 4.58 12 2.08 -0.42" | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
echo "5 5 11 5 5 -1" | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
output: case2.jpg
code:
echo "49 -0.155 1.12 49 -0.155 -1.43" | psxy -J -R -EY5p -Gblack -Ss10p -K -O >> $ps
echo "37 -1.215 -0.89 37 -1.215 -1.54" | psxy -J -R -EY5p -Gblack -Ss10p -K -O >> $ps
echo "23 -1.955 -1.43 23 -1.955 -2.48" | psxy -J -R -EY5p -Gblack -Ss10p -K -O >> $ps
echo "16 -1.205 -0.31 16 -1.205 -2.1" | psxy -J -R -EY5p -Gblack -Ss10p -K -O >> $ps
echo "12 2.08 4.58 12 2.08 -0.42" | psxy -J -R -EY5p -Gblack -Ss10p -K -O >> $ps
echo "5 5 11 5 5 -1" | psxy -J -R -EY5p -Gblack -Ss10p -K -O >> $ps
for both cases:
JXscale="-JX20/15"
Rregion="-R0/90/-15/15"
As you can see, the only difference is "y" and "Y".
The problem is that I have a loop that goes over several rows, so the above is just an example. This is the raw format of one row:
MN,5,11,-1
and this is what does the job:
awk '{ if ( $1 == "MN" ) print $2, ($3+$4)/2, $3, $2, ($3+$4)/2, $4 }' FS="," myfile | psxy -J -R -Ey5p -Gblack -Ss10p -K -O >> $ps
so, what am I doing wrong? how can I avoid such difference? my GMT version is
To unsubscribe, send the message "signoff gmt-help" to [email protected]
case1.jpg
(image/jpeg, 289.6 KB) - not displayed
case2.jpg
(image/jpeg, 275.2 KB) - not displayed