1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70

	#!/bin/bash
	#gmtset BASEMAP_TYPE fancy
	gmtset ANOT_FONT_SIZE 8p
	gmtset LABEL_FONT_SIZE 8p
	gmtset BASEMAP_TYPE plain
	PS=slip.ps
	
	R=172/175.5/-43.5/-40.5

	# slip.txt为滑移量文件，分别为经度 维度 滑移量
	gawk '{print $2,$1,$3}' slip1.txt | surface -Gslip1.grd -I0.01/0.01 -R  -V 
	gawk '{print $2,$1,$3}' slip2.txt | surface -Gslip2.grd -I0.01/0.01 -R  -V
	gawk '{print $2,$1,$3}' slip3.txt | surface -Gslip3.grd -I0.01/0.01 -R  -V
	gawk '{print $2,$1,$3}' slip4.txt | surface -Gslip4.grd -I0.01/0.01 -R  -V
	gawk '{print $2,$1,$3}' slip5.txt | surface -Gslip5.grd -I0.01/0.01 -R  -V
	
	grd2cpt slip2.grd -Crainbow -S0/6/0.1 -Z > slip.cpt

	#pre.txt为圈定区域位置的文件，分别为滑移面四个顶点经维度组成的闭合区域 
	grdmask pre1.txt -Gfactor1.nc -R -I0.01/0.01 -Nnan/1/1
	grdmask pre2.txt -Gfactor2.nc -R -I0.01/0.01 -Nnan/1/1
	grdmask pre3.txt -Gfactor3.nc -R -I0.01/0.01 -Nnan/1/1
	grdmask pre4.txt -Gfactor4.nc -R -I0.01/0.01 -Nnan/1/1
	grdmask pre5.txt -Gfactor5.nc -R -I0.01/0.01 -Nnan/1/1
	
	grdmath factor1.nc slip1.grd MUL = test1.grd
	grdmath factor2.nc slip2.grd MUL = test2.grd
	grdmath factor3.nc slip3.grd MUL = test3.grd
	grdmath factor4.nc slip4.grd MUL = test4.grd
	grdmath factor5.nc slip5.grd MUL = test5.grd
	
	psbasemap -JX5id -R$R -B1f0.5/1f0.5WSen -K > $PS
	
	grdimage newzealand.nc -R -J -Cchina1.cpt  -K -O >> $PS
	
	psclip pre1.txt -R -J -K -O  -P >> $PS
	grdimage test1.grd -J -R  -Cslip.cpt -K -O >> $PS
	psclip -C -O -K -P >> $PS
	
	psclip pre2.txt -R -J -K -O  -P >> $PS
	grdimage test2.grd -J -R  -Cslip.cpt -K -O >> $PS
	psclip -C -O -K -P >> $PS
	
	psclip pre3.txt -R -J -K -O  -P >> $PS
	grdimage test3.grd -J -R  -Cslip.cpt -K -O >> $PS
	psclip -C -O -K -P >> $PS
	
	psclip pre4.txt -R -J -K -O  -P >> $PS
	grdimage test4.grd -J -R  -Cslip.cpt -K -O >> $PS
	psclip -C -O -K -P >> $PS
	
	psclip pre5.txt -R -J -K -O  -P >> $PS
	grdimage test5.grd -J -R  -Cslip.cpt -K -O >> $PS
	psclip -C -O -K -P >> $PS
	
	psxy line1.txt  -J -R -O -K -Wblack -M'999' >> $PS
	psxy line2.txt  -J -R -O -K -Wgrey -M'999' >> $PS
	psxy line3.txt  -J -R -O -K -Wgrey -M'999' >> $PS
	psxy line4.txt  -J -R -O -K -W1.4/black -M'999' >> $PS
	psxy line5.txt  -J -R -O -K -W1.2/black -M'999' >> $PS
	
	pscoast -R -J  -O -K -Dh -A300 -I1 -W0.25p -N1/1p  >> $PS
	
	psscale  -Cslip.cpt -D3.4i/0.5i/6.0/0.3h -B1/:"Slip/m": -K -O >> $PS
	
	psxy -J -R faultnzland.dat  -W1.5/139/71/38 -M -O -K >> $PS
	
	echo  173.054 -42.737 5 219 38 128 2.6 173.077 -42.757  | psmeca -R$R -J -Sa1/8p/1c -Gred -Ewhite  -W1.2/red -L1  -C1.5/red  -O  >> $PS
	
	ps2raster $PS -A -P -Tg