1:split
#一个逗号代表一个空格
ls -1 *.SAC|gawk '
{x=$1; split(x,aa,"\."); print
"mv",$1,aa[1]"."aa[2]"."aa[3]"."aa[8] "." aa[9] "." aa[10]}'
|sh2:lh dist
#获取某个头文件值输入至某文件
ls *.?HZ |gawk '{print
"rh",$1;print "lh dist";}'|sac2000 |gawk -f
~/bin/extSacInfo.awk |sort -n --key=1 >
dist.info
3:if
#if双重否定
cat dist.info |gawk '{if ($2 <=
100&&$2 -n) print
$1,$2}'>100.list
4:plot sac files
#参数的意义参考psbasemap help , pssac2 help
psbasemap -R-5/100/0/380 -JX6/9 -B5f5:"time
(sec)":/20f10:"epicenter (deg)":WeNs -K -P >
z.eqr.epi.ps
ls -1 *.eqr | /opt/pssac2/pssac2 -R -JX -M0.5
-G0/0/255/0.001/-5/30 -Edt1 -O -K -P >> z.eqr.epi.ps
5:plot earthquakes all over the
earth
pscoast -JM6 -R70/140/0/55 -G220 -S0/0/200 -W1 -Di -N1 -V -P -K
-B5f1/2f5:.earthquake: -Y2 -X1 >station.ps
cat location2.dat | gawk '{print $1,$2}' | psxy -JM -R -Sc0.1
-G225/0/0 -V -O -K
-P>>station.ps
pstext station1.dat -Jm0.1i -R70/140/0/55 -G200/100/0 -K -X-0.4
-Y-0.7 -P -O >>station.ps
pstext station2.dat -Jm0.1i -R70/140/0/55 -G200/100/0 -K -X0.2
-Y-0.1 -P -O >>station.ps
format of stations files
97.50 31.00 10 3 0 1 CAD
90.76 37.86 10 3 0 1 HTG
95.00 35.00 10 3 0 1 GOM
92.25 32.25 10 3 0 1 NAQ
format of location files
81.44 30.59
74.65 35.52
74.60 35.42
74.51 35.41
74.49 35.30
77.39 39.48
77.43 39.36
77.34 39.52
77.22 39.43
80.65 35.80
77.58 39.41
78.44 36.79
83.67 30.69
83.63 30.71
81.15 31.02
经度 纬度
6:文件中的数据加减
cat ptime.info | gawk '{x=$3-$2;if(x -n) print x }'
>p.info
后面的if语句是保证结果非零!
7:合并两个文件
>paste 1.info
2.info >3.info