FW: error bar confusion

Dear community,

problem solved in the earlier parts of the script, so sorry for the noise and the incomplete email.

Best regards,

From: [email protected]
To: [email protected]
Subject: error bar confusion
Date: Fri, 20 Dec 2013 14:06:55 +0000




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]