Re: usage of -i option with psxy

Hello,

I try again to use -i option in place of AWK command, with GMT51.0 psvelo command, in order to draw GPS vectors (data file attached).

With the awk command :

awk '{print($1,$2,$3,$4,$5,$6,$7,"")}' JGR2007-Gan_GPS_Tibet.dat |\
gmt psvelo $region $projection -A$vecsize -Se0.05/0.95/12p -Gfirebrick -Wthick,black -L -K -O -V >> $carte

vectors are drawn correctly, without text specifications listed in the 7th column of the data file. 

So, I suppose that the simpler command :

gmt psvelo JGR2007-Gan_GPS_Tibet.dat $region $projection -A$vecsize -Se0.05/0.95/12p -i0,1,2,3,4,5,6 -Gfirebrick -Wthick,black -L -K -O -V >> $carte

should give the same output. But now, text data appear on the map, as if the -i option was not taken into account.

I don’t see where’s my mistake.
Thanks for help.
Cyril



Le 7 déc. 2013 à 00:13, Paul Wessel <[email protected]> a écrit :

> You are off by one column.  In GMT the first column is 0. So -i0,1,5s0.025
> 
> -p
> 
> On Dec 6, 2013, at 12:45 PM, Cyril Langlois <[email protected]> wrote:
> 
>> Hello,
>> 
>> I’ve tried to map with GMT5.1.0 the seamounts database Global_Seamounts_JGR2001_revised.txt published by Paul Wessel in 2001, with seamount radius as symbol size.
>> 
>> So, using awk, I write :
>> 
>> awk '{if((substr($1,1,1)!="#")&&(substr($1,1,1)!=">")&&($1!=""))print($1,$2,0.0025*$6)}' Global_Seamounts_JGR2001_revised.txt  |\
>> gmt psxy -R -J -Scc -Wthinner,black -Gfirebrick -O -V >> $map
>> 
>> It works (attached map below). 
>> 
>> Reading the GMT5.1.0 manual, I believed I could do the same without awk, with the psxy -i option, choosing longitude (1), latitude (2) and seamount size (6th) columns of the file and resizing seamount radius. So, I write :
>> 
>> gmt psxy -R -J -Scc -i1,2,6s0.025 -hi  -Wthinner,black-Gfirebrick -O -V >> $map
>> 
>> but it doesn’t work : data symbols are not plotted with the expected size nor at the good location.
>> 
>> Where am I wrong ?
>> Thanks for help,
>> 
>> Cyril

Cyril LANGLOIS
[email protected]
Tél : 06 89 05 68 22

http://www.normalesup.org/~clanglois/Index.html
§ ---------------------- §
"Les hommes ne souffrent que d'une seule maladie : l'ignorance. A ce mal, il n'y a qu'un seul remède: l'instruction. Mais il ne s'agit pas de prendre ce remède à doses homéopathiques. Il faut le prendre par seaux entiers, et par barriques de quarante seaux." Dmitri Pissariev

Dr. Cyril LANGLOIS
Professeur agrégé - enseignant de Sciences de la Terre 
[email protected]
Tél. portable : 06 89 05 68 22
Tél. bureau : 04 72 72 85 28
§ --------------------------- §
Préparation à l'Agrégation SV-STU
Ecole Normale Supérieure de Lyon -  Site Monod
46, allée d'Italie
69364 Lyon cedex 07
§ --------------------------- §


To unsubscribe, send the message "signoff gmt-help" to [email protected]