gmtset FRAME_WIDTH 0.1c
gmtset DOTS_PR_INCH 600
gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
makecpt -Crelief -T-8000/8000/100 >1.cpt
grdimage ETOPO1_Bed_g.g98 -R65/145/0/60 -Jm0.2c -B10g0WSne -C1.cpt -P -K -V  >atml.ps
pscoast -R -Jm -B -W1 -I1  -Df -A500 -Cblue -Slightblue  -V -K -O >>atml.ps
psvelo atml.out -R -W1/255/97/0  -Se1.0/0.3/8 -B -Jm -A0.01/0.12/0.03 -P -K  -L -V -O -D2 >>atml.ps
psvelo arrow.txt -V  -Se1.0/0.3/8 -A0.01/0.12/0.03 -W1/255/97/0 -R -Jm -P -U"Made by L.Xing"  -O -K -L>>atml.ps
psbasemap -R -Jm -O -P -Lf75/5/0/2000+lKM -B -K -V >>atml.ps
psxy ChinaBoundary.dat -R -J -B -W0.1p,red -mNan -O -P -V >>atml.ps