Changeset 263 for SH_SHM/trunk

Timestamp:

08/20/2010 09:29:33 AM (13 years ago)

Author:

marcus

Message:

r136 | walther | 2010-08-20 09:29:00 +0200 (Fr, 20 Aug 2010) | 1 line

Adopted script to use GMT prefix if available.

File:

• SH_SHM/trunk/util/locate_by_ttfit.csh (modified) (7 diffs)

SH_SHM/trunk/util/locate_by_ttfit.csh

@@ -4 +4 @@
 #      ===================
 #
-# $Revision: 31 $, $Date: 2008-03-03 20:47:24 +0100 (Mo, 03 MÀr 2008) $
+# $Revision: 136 $, $Date: 2010-08-20 09:29:00 +0200 (Fr, 20 Aug 2010) $
 #
 # Locate event via travel time fit.
@@ -13 +13 @@
         exit
 endif
+
+# check for GMT command access. If GMT is defined, use it as prefix for
+# all commands.
+set testGMT=`which GMT`
+set GMT=""
+[ -x "$testGMT" ] && set GMT=$testGMT
 
 #set echo
@@ -53 +59 @@
 chdir $SH_SCRATCH
 # set units for plot, want to change this to cm
-gmtdefaults -Du >.gmtdefaults
+$GMT gmtdefaults -Du >.gmtdefaults
 
 # delete previously existing result file
@@ -238 +244 @@
 if  ("$method" == "grid")  then
         grep -v '^\!' $fitxyz | grep " $mindep depth" | \
-                psxy -: -R$lolon/$hilon/$lolat/$hilat \
+                $GMT psxy -: -R$lolon/$hilon/$lolat/$hilat \
                 $fmt -Ss$symsiz -C$cpt -U \
                 -B${lab}:longitude:/${lab}:latitude::."$areaname":8 \
@@ -244 +250 @@
 else
         grep -v '^\!' $fitxyz | \
-                psxy -: -R$lolon/$hilon/$lolat/$hilat \
+                $GMT psxy -: -R$lolon/$hilon/$lolat/$hilat \
                 $fmt -Ss$symsiz -C$cpt -U \
                 -B${lab}:longitude:/${lab}:latitude::."$areaname":8 \
@@ -259 +265 @@
 
 # plot coastlines and borders
-pscoast -O -R$xlolon/$xhilon/$lolat/$hilat $fmt -A10 \
+$GMT pscoast -O -R$xlolon/$xhilon/$lolat/$hilat $fmt -A10 \
         -Df -N1/4/0/0/0/dashed -W1/0/0/0/solid -Ia/0.1/100 -K >>$psfile
 
 # plot cities
 awk '{print $1,$2}' $SH_UTIL/cities.dat \
-        | psxy -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt \
+        | $GMT psxy -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt \
         -W5/$citycol -G$citycol -Sc0.1 >>$psfile
 awk '{print $1,$2,"14","0","4","1","\ ",$3}' $SH_UTIL/cities.dat \
-        | pstext  -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt -G$citycol >>$psfile
+        | $GMT pstext  -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt -G$citycol >>$psfile
 
 # integrate text into psfile
@@ -279 +285 @@
         if  ($lcnt == 1)  set cont=""
         echo "6.0 $ypos.0 12 0.0 8 1  $line" | \
-                pstext -O $cont -JX10/8 -R0/10/0/100 >>$psfile
+                $GMT pstext -O $cont -JX10/8 -R0/10/0/100 >>$psfile
         @ lcnt = $lcnt - 1
         @ ypos = $ypos + 2