2017-12-27 32 views
1

나는 Southamerica의 북쪽 부분에 대한 몇 가지 FOCAL MECHANISMS (FM)을 GMT에 만들었고이를 맵에 넣어야합니다.GMT (Generic Mapping Tools)에서 FOCAL MECHANISMS를 수행하려면 어떻게해야합니까?

내가 85 개 지진과 문서를 가지고있는 FM 들어

이 코드입니다 :

psmeca eventos2.dat -R-80/-66/-5/14 -JM6i -Sm0.4/0/0.05u -P -V > coltopomap.ps 
GMT psconvert coltopomap.ps -Tf 

그리고 이것은 내가

으로 시도했습니다

GMT grdcut topo15.grd -R-80/-66/-5/14 -Gtopo15_cut_col.grd 
GMT grdgradient topo15_cut_col.grd -Gtopo15_cut_shaded.grd -A310 -Ne0.6 
GMT grdimage topo15_cut_col.grd -JM6i -R-80/-66/-5/14 -Cglobe -P -Ba --MAP_FRAME_TYPE=plain --MAP_GRID_CROSS_SIZE_PRIMARY=0p -Itopo15_cut_shaded.grd -K > coltopomap.ps 
GMT pscoast -R-80/-66/-5/14 -J -Df -N1/1,black -W1/0.5,black -Lf-68/-4/-68/-4/200k+l -Tdg-78/12+w0.3i+f3+lW,E,S,N --FONT_ANNOT=10p --FONT_LABEL=10p --FONT_TITLE=10p -O >> coltopomap.ps 

가 그들과 합류하기 위해지도입니다

GMT grdcut topo15.grd -R-80/-66/-5/14 -Gtopo15_cut_col.grd 
GMT grdgradient topo15_cut_col.grd -Gtopo15_cut_shaded.grd -A310 -Ne0.6 
GMT grdimage topo15_cut_col.grd -JM6i -R-80/-66/-5/14 -Cglobe -P -Ba --MAP_FRAME_TYPE=plain --MAP_GRID_CROSS_SIZE_PRIMARY=0p -Itopo15_cut_shaded.grd -K > coltopomap.ps 
GMT pscoast -R-80/-66/-5/14 -J -Df -N1/1,black -W1/0.5,black -Lf-68/-4/-68/-4/200k+l -Tdg-78/12+w0.3i+f3+lW,E,S,N --FONT_ANNOT=10p --FONT_LABEL=10p --FONT_TITLE=10p -O >> coltopomap.ps 
GMT psmeca eventos2.dat -R-80/-66/-5/14 -JM6i -Sm0.4/0/0.05u -P -V >> coltopomap.ps 
GMT psconvert coltopomap.ps -Tf 

하지만 내가 얻을 수있는 유일한 정보는 다음 이미지에서 볼 수있는 것처럼 한 페이지에있는지도와 다른 페이지에있는 FM입니다.

https://i.stack.imgur.com/bAvx7.png

답변

0

당신이 당신의 GMT pscoast ... 명령에서 -K 플래그를 누락 보인다. GMT documentation에서

: 다층에게 플롯

을 구축 할 때

마지막 플롯 명령을 제외한 모든 필수