x2sys confusion
杨磊 <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <CADXaJRYX4o2ktR+At9JiBMTL96UFz-jDP7d7ccv8jZZUkuqMKQ@mail.gmail.com> |
Dear community,
I use the GMT x2sys to calculate crossover of some tracks, but error
happens.
here is the confusing result [image: 内嵌图片 1]
two crossover points are right, but others near the equator are wrong.
The data and def file were attached.
the codes are :
gmt x2sys_init CROSSHY2 -Dxys -V -I1 -Cg -Etxt -Gg -Ndk -Wd8 -F
gmt x2sys_cross -TCROSSHY2 srl_test_gdr005.txt -Qi -W5 -Ia -V2 >
COE_orig.txt
gmt psbasemap -Rd -JQ4i -P -K -Bxa60 -Bya10g10 -BWSen >trackshy2.ps
gawk "NR>4 {print $1,$2,$11}" COE_orig.txt > data.d
gmt psxy -R -J -O -K data.d -Sc0.2i >> trackshy2.ps
gmt psxy srl_test_gdr005.txt -Sc0.01c -R -J -P -O >>trackshy2.ps
So, where is the point ? the GMT is 5.1.0
To unsubscribe, send the message "signoff gmt-help" to [email protected]
trackshy2.png
(image/png, 14.6 KB) - not displayed
srl_test_gdr005.txt
(text/plain, 628.2 KB) - not displayed
xys.def
(application/octet-stream, 426 B) - not displayed