some grdcontour bugs

[Alexander Kukharsky]

Dear colleagues,

I've found a non accurate (splitted) lines as a result of the grdcontour action.
(1) This wrong lines was drawn along the zero meridian when the pole is inside the boundaries (level 524). The correct line was drawn in case of North Pole was outside the boundaries with the same input file.
(2) Wrong indication of lows/highs (maybe only along the zero meridian too).

Please, see my added picture and explain how so solve this problems.

I also have some troubles with

(3) and with changing the color of ANNOTATION FONT (I see that -W parameter controls the color of local lows/highs and the color of the lines simultaneously).
(4) the elimination of the contour clipping at the boundaries

May be is it possible to control this parameters?

I am using GMT-4.5 on a 32-bit Windows system.

surface mean2.xyz -I1 -R -Gdata2.grd
grdcontour data2.grd -JS55/90/16c -R0/30/150/85r -P -C4 -Q100 -A+a0+s18+kBrown+n.1/.1 -T:LH -G5i/5i -S2  -Wa1.5p/Red -K -X0 -Y0>02.ps
pscoast -Dc -R0/30/150/85r -bo -B20g20/10g10 -JS55/90/16c -P -W0.5p  -O>>02.ps


Best Regards

Alex Kukharsky


To unsubscribe, send the message "signoff gmt-help" to listserv@...

[Paul Wessel-2]

Hi Alex-

A few issues:

1. The command you sent me does not reproduce you PNG.  Could you tell  
me which -R -J you used for that plot so I can reproduce.  Also your  
script passed -R to surface but specified no values.  I surmised from  
your grid it is -R0/360/0/90.
2. However, surface has NO knowledge of geographic boundary  
conditions, and furthermore cannot handle geographic poles.  Looking  
at your data I see it is already gridded at 2.5 degrees, so the  
solution is to just put it into grd format with xyz2grd mean2.xyz -
I2.5 -R0/360/0/90 -Gdata2.grd (this is just reformatting).
3. Once that is done the contours looks OK to me, but perhaps you can  
redo your exercise and send me a PNG and the exact commands used.

Once the lines look OK we can discuss if there are font/color issues.

Cheers,
Paul

On Nov 6, 2009, at 9:13 PM, Alexander Kukharsky wrote:

> Dear colleagues,
>
> I've found a non accurate (splitted) lines as a result of the  
> grdcontour action.
> (1) This wrong lines was drawn along the zero meridian when the pole  
> is inside the boundaries (level 524). The correct line was drawn in  
> case of North Pole was outside the boundaries with the same input  
> file.
> (2) Wrong indication of lows/highs (maybe only along the zero  
> meridian too).
>
> Please, see my added picture and explain how so solve this problems.
>
> I also have some troubles with
>
> (3) and with changing the color of ANNOTATION FONT (I see that -W  
> parameter controls the color of local lows/highs and the color of  
> the lines simultaneously).
> (4) the elimination of the contour clipping at the boundaries
>
> May be is it possible to control this parameters?
>
> I am using GMT-4.5 on a 32-bit Windows system.
>
> surface mean2.xyz -I1 -R -Gdata2.grd
> grdcontour data2.grd -JS55/90/16c -R0/30/150/85r -P -C4 -Q100 -A
> +a0+s18+kBrown+n.1/.1 -T:LH -G5i/5i -S2  -Wa1.5p/Red -K -X0 -Y0>02.ps
> pscoast -Dc -R0/30/150/85r -bo -B20g20/10g10 -JS55/90/16c -P -W0.5p  
> -O>>02.ps
>
>
> Best Regards
>
> Alex Kukharsky
>
>
> To unsubscribe, send the message "signoff gmt-help" to listserv@...
> <LOW-HIGH.ZIP><low-high.png>

To unsubscribe, send the message "signoff gmt-help" to listserv@...

[Alexander Kukharsky]

Hi, Paul

1. It was unnecessary to add the surface string in my script. The correct example should be produced using only grdcontour line (file data2.grd loaded in my archive). Sorry for my confused message.

Moreover, another image with script was posted for GMT mailing list {27 Oct 2009} without possibility of notification.
https://listserv.hawaii.edu/cgi-bin/wa?A2=ind0910&L=GMT-HELP&D=1&T=0&O=D&P=108109

2. This stripped lines was drawn in any cases of gridding. The lines with -I1 produced more smooth lines and needed detalization level.

Regards

Alex Kukharsky



-----Original Message-----

> Hi Alex-
>
> A few issues:
>
> 1. The command you sent me does not reproduce you PNG.  Could you tell  
> me which -R -J you used for that plot so I can reproduce.  Also your  
> script passed -R to surface but specified no values.  I surmised from  
> your grid it is -R0/360/0/90.
> 2. However, surface has NO knowledge of geographic boundary  
> conditions, and furthermore cannot handle geographic poles.  Looking  
> at your data I see it is already gridded at 2.5 degrees, so the  
> solution is to just put it into grd format with xyz2grd mean2.xyz -
> I2.5 -R0/360/0/90 -Gdata2.grd (this is just reformatting).
> 3. Once that is done the contours looks OK to me, but perhaps you can  
> redo your exercise and send me a PNG and the exact commands used.
>
> Once the lines look OK we can discuss if there are font/color issues.
>
> Cheers,
> Paul
>
> On Nov 6, 2009, at 9:13 PM, Alexander Kukharsky wrote:
>
> > Dear colleagues,
> >
> > I've found a non accurate (splitted) lines as a result of the  
> > grdcontour action.
> > (1) This wrong lines was drawn along the zero meridian when the pole  
> > is inside the boundaries (level 524). The correct line was drawn in  
> > case of North Pole was outside the boundaries with the same input  
> > file.
> > (2) Wrong indication of lows/highs (maybe only along the zero  
> > meridian too).
> >
> > Please, see my added picture and explain how so solve this problems.
> >
> > I also have some troubles with
> >
> > (3) and with changing the color of ANNOTATION FONT (I see that -W  
> > parameter controls the color of local lows/highs and the color of  
> > the lines simultaneously).
> > (4) the elimination of the contour clipping at the boundaries
> >
> > May be is it possible to control this parameters?
> >
> > I am using GMT-4.5 on a 32-bit Windows system.
> >
> > surface mean2.xyz -I1 -R -Gdata2.grd
> > grdcontour data2.grd -JS55/90/16c -R0/30/150/85r -P -C4 -Q100 -A
> > +a0+s18+kBrown+n.1/.1 -T:LH -G5i/5i -S2  -Wa1.5p/Red -K -X0 -Y0>02.ps
> > pscoast -Dc -R0/30/150/85r -bo -B20g20/10g10 -JS55/90/16c -P -W0.5p  
> > -O>>02.ps
> >
> >
> > Best Regards
> >
> > Alex Kukharsky
> >
> >
> > To unsubscribe, send the message "signoff gmt-help" to listserv@...
> > <LOW-HIGH.ZIP><low-high.png>
>
> To unsubscribe, send the message "signoff gmt-help" to listserv@...
>

To unsubscribe, send the message "signoff gmt-help" to listserv@...