Draw velocity using GMT
This paper introduces a method for velocity drawing by GMT .The main procedure used in this example is Psvelos. The velocity data were calculated by YUAN Yongdong using IGS data before and after the great earthquake happened in Japan last year.Bat code :set R=-R128/148/30/46Rem extract and combine DEM data over Japan area grdraster 32 -R128/140/40/46 -I0.5m -Gs1.ncgrdraster 33 -R128/140/30/40 -I0.5m -Gs2.ncgrdraster 34 -R140/148/30/40 -I0.5m -Gs3.ncgrdraster 35 -R140/148/40/46 -I0.5m -Gs4.ncgrdpaste s1.nc s2.nc -Gjapan1.ncgrdpaste s3.nc s4.nc -Gjapan2.ncgrdpaste japan1.nc japan2.nc -Gjapan.ncRem make a color palette makecpt -Ctopo -T-100/2500/200 -Z>japan.cptRem draw DEM using GTOPO DEM data grdimage -JM5i %R% -Cjapan.cpt japan.nc -P -K >japanVelo.pspscoast %R% -JM5i-Ba5g5/a5g5-Di -A5000 -Swhite -Na -P -V -K -O>>japanVelo.psRem you can give a legen
页:
[1]