TT

来源：互联网发布：c#面试算法题编辑：程序博客网时间：2024/05/01 00:29

#!/bin/shcat stat.nau.d | gawk '{print $3,$2,$1}' > stat.naucat stat.au.d | gawk '{print $3,$2,$1}' > stat.aupsfile=fig/map.psR=110/160/-45/-10J=B135/-25/-36/-18/5ievlo=131.9755evla=-25.9665# set GMT defaultgmtset FONT_ANNOT_PRIMARY 12p# plot the basemap pscoast -J$J -R$R -K -B10/5 -I1 -W1/1p -N1/red -N2 \-A500 -Gwhite -Slightblue -Cblue \-FAU.NSW,AU.QLD,AU.SA,AU.TAS,AU.VIC,AU.WA,AU.ACT,AU.NT \-P -Xc -Yc --FORMAT_GEO_MAP=dddF > $psfile# event location or hypocentrepsxy -J$J -R$R -O -K -Sa0.15i -Gyellow << EOF >> $psfile$evlo $evlaEOFwhile read linedosta=`echo $line | awk '{print $3}'`stla=`echo $line | awk '{print $2}'`stlo=`echo $line | awk '{print $1}'`gmt project -C$evlo/$evla -E$stlo/$stla -G10 -Q > great_circle.xypgmt psxy -J$J -R$R -K -O -W0.8,blue -Gblue great_circle.xyp >> $psfilestlo=`echo $line | awk '{print $1+2}'`stla=`echo $line | awk '{print $2+1}'`echo $stlo $stla $sta | pstext -J$J -R$R -K -O -F+f8,blue >> $psfiledone < stat.auwhile read linedosta=`echo $line | awk '{print $3}'`stla=`echo $line | awk '{print $2}'`stlo=`echo $line | awk '{print $1}'`gmt project -C$evlo/$evla -E$stlo/$stla -G10 -Q > great_circle.xypgmt psxy -J$J -R$R -K -O -W0.8,red -Gblue great_circle.xyp >> $psfilestlo=`echo $line | awk '{print $1+2}'`stla=`echo $line | awk '{print $2+1}'`echo $stlo $stla $sta | pstext -J$J -R$R -K -O -F+f8,red >> $psfiledone < stat.nau# plot stationpsxy -J$J -R$R -O -K -Sd0.1i -Gblue stat.au >> $psfilepsxy -J$J -R$R -O -Sd0.1i -Gred stat.nau >> $psfile# convert ps file to jpegfor file in `ls fig/*.ps`dops2raster -A -P $filedone

0 0