Changeset 263 for SH_SHM/trunk
- Timestamp:
- 08/20/2010 09:29:33 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
SH_SHM/trunk/util/locate_by_ttfit.csh
r46 r263 4 4 # =================== 5 5 # 6 # $Revision: 31 $, $Date: 2008-03-03 20:47:24 +0100 (Mo, 03 MÀr 2008) $6 # $Revision: 136 $, $Date: 2010-08-20 09:29:00 +0200 (Fr, 20 Aug 2010) $ 7 7 # 8 8 # Locate event via travel time fit. … … 13 13 exit 14 14 endif 15 16 # check for GMT command access. If GMT is defined, use it as prefix for 17 # all commands. 18 set testGMT=`which GMT` 19 set GMT="" 20 [ -x "$testGMT" ] && set GMT=$testGMT 15 21 16 22 #set echo … … 53 59 chdir $SH_SCRATCH 54 60 # set units for plot, want to change this to cm 55 gmtdefaults -Du >.gmtdefaults61 $GMT gmtdefaults -Du >.gmtdefaults 56 62 57 63 # delete previously existing result file … … 238 244 if ("$method" == "grid") then 239 245 grep -v '^\!' $fitxyz | grep " $mindep depth" | \ 240 psxy -: -R$lolon/$hilon/$lolat/$hilat \246 $GMT psxy -: -R$lolon/$hilon/$lolat/$hilat \ 241 247 $fmt -Ss$symsiz -C$cpt -U \ 242 248 -B${lab}:longitude:/${lab}:latitude::."$areaname":8 \ … … 244 250 else 245 251 grep -v '^\!' $fitxyz | \ 246 psxy -: -R$lolon/$hilon/$lolat/$hilat \252 $GMT psxy -: -R$lolon/$hilon/$lolat/$hilat \ 247 253 $fmt -Ss$symsiz -C$cpt -U \ 248 254 -B${lab}:longitude:/${lab}:latitude::."$areaname":8 \ … … 259 265 260 266 # plot coastlines and borders 261 pscoast -O -R$xlolon/$xhilon/$lolat/$hilat $fmt -A10 \267 $GMT pscoast -O -R$xlolon/$xhilon/$lolat/$hilat $fmt -A10 \ 262 268 -Df -N1/4/0/0/0/dashed -W1/0/0/0/solid -Ia/0.1/100 -K >>$psfile 263 269 264 270 # plot cities 265 271 awk '{print $1,$2}' $SH_UTIL/cities.dat \ 266 | psxy -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt \272 | $GMT psxy -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt \ 267 273 -W5/$citycol -G$citycol -Sc0.1 >>$psfile 268 274 awk '{print $1,$2,"14","0","4","1","\ ",$3}' $SH_UTIL/cities.dat \ 269 | pstext -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt -G$citycol >>$psfile275 | $GMT pstext -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt -G$citycol >>$psfile 270 276 271 277 # integrate text into psfile … … 279 285 if ($lcnt == 1) set cont="" 280 286 echo "6.0 $ypos.0 12 0.0 8 1 $line" | \ 281 pstext -O $cont -JX10/8 -R0/10/0/100 >>$psfile287 $GMT pstext -O $cont -JX10/8 -R0/10/0/100 >>$psfile 282 288 @ lcnt = $lcnt - 1 283 289 @ ypos = $ypos + 2
Note: See TracChangeset
for help on using the changeset viewer.