Re: grdvector with GMT5

Thanks Cyril, for exposing more issues. It was a mix of bugs and "features" with the new geo-vectors.

The weirdness had to do with the vector heads being longer than the length of the vector.  That is, the chosen vector head length of 0.25cm, when converted to km on the ground, in many places exceeded the length of the data vector projected to km on the ground.  This would bring up the issue of using the +n<norm> modifier to -Q, but it had bugs as well.  These are now all fixed I hope.  Here is the modified behavior:

1. If the vector head length exceed 90% of the vector length the head is skipped and a warning is issued
2. If +n<norm> is used with -Q then four things can happen:
	a) vector length >= norm and vector head length >= 90% of the vector length: Normal plotting
	b) vector length >= norm and vector head length < 90% of the vector length: Head is shrunk down to satisfy criteria; warning issued
	c) vector length <norm:
		Vector pen and head shrinks down as desired. Then
		d) if shrunk head length still >= 90% of the vector length: head shrinks down further, vector pen stays the same.; warning issued

I have added a new test to our suite for this behavior. PDF is attached and panels with blue line has +n such that the left side experiences shrinking.
The other panel just skips head when vector is too short.

In r12520.

Cheers, Paul

On Nov 15, 2013, at 6:36 AM, Langlois <[email protected]> wrote:

> Hello,
> 
> I’m getting trouble with grdvector use with GMT5.1.0 :
> I have two scripts with the same grdvector command to draw Nuvel-1 velocity vectors, one for a regional map, the other for a global map (attached). But the output maps look completely differents (see attached images) and I don’t undersand why.
> 
> My grdvector command is, for both map :
>  
> # add vector field
> vecsize=0.25+e+p0.25p
> velscale=i80
> 
> #Nuvel-1A
> gmt grdvector $rasterpath/plates/nuvel1.vy.1.5.grd $rasterpath/plates/nuvel1.vx.1.5.grd \
>  	 -T $region $projection -I30 -S$velscale -Q$vecsize \
>  	 -Wthick,black -Gorange --MAP_VECTOR_SHAPE=0.5 -O -K -V >> $carte
> 
> where "$rasterpath/plates/nuvel1.vy.1.5.grd » is the Nuvel-1 velocity grid (x and y components). Only the -I option are changed (-I30 for the global map, -I5 for the regional one).
> 
> How can the same script command gives so different results ?? Both scripts are also attached below.
> Thanks for help,
> 
> Cyril
> 
> 
> To unsubscribe, send the message "signoff gmt-help" to [email protected] <GMT5_Nuvel_Eu-fixe-JM.sh>
> To unsubscribe, send the message "signoff gmt-help" to [email protected]<GMT5_Tibet_GPS-Nuvel-Eurasie-2013-3.sh>
> To unsubscribe, send the message "signoff gmt-help" to [email protected]<GMT_Tibet_GPS_Nuvel_2013-simple.pdf>
> 
> To unsubscribe, send the message "signoff gmt-help" to [email protected]<GMT5_Nuvel-plaques-Eu-Mercator.pdf>
> 
> 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]


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