Changeset 263


Ignore:
Timestamp:
08/20/10 09:29:33 (9 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:
1 edited

Legend:

Unmodified
Added
Removed
  • SH_SHM/trunk/util/locate_by_ttfit.csh

    r46 r263  
    44#      =================== 
    55# 
    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) $ 
    77# 
    88# Locate event via travel time fit. 
     
    1313        exit 
    1414endif 
     15 
     16# check for GMT command access. If GMT is defined, use it as prefix for 
     17# all commands. 
     18set testGMT=`which GMT` 
     19set GMT="" 
     20[ -x "$testGMT" ] && set GMT=$testGMT 
    1521 
    1622#set echo 
     
    5359chdir $SH_SCRATCH 
    5460# set units for plot, want to change this to cm 
    55 gmtdefaults -Du >.gmtdefaults 
     61$GMT gmtdefaults -Du >.gmtdefaults 
    5662 
    5763# delete previously existing result file 
     
    238244if  ("$method" == "grid")  then 
    239245        grep -v '^\!' $fitxyz | grep " $mindep depth" | \ 
    240                 psxy -: -R$lolon/$hilon/$lolat/$hilat \ 
     246                $GMT psxy -: -R$lolon/$hilon/$lolat/$hilat \ 
    241247                $fmt -Ss$symsiz -C$cpt -U \ 
    242248                -B${lab}:longitude:/${lab}:latitude::."$areaname":8 \ 
     
    244250else 
    245251        grep -v '^\!' $fitxyz | \ 
    246                 psxy -: -R$lolon/$hilon/$lolat/$hilat \ 
     252                $GMT psxy -: -R$lolon/$hilon/$lolat/$hilat \ 
    247253                $fmt -Ss$symsiz -C$cpt -U \ 
    248254                -B${lab}:longitude:/${lab}:latitude::."$areaname":8 \ 
     
    259265 
    260266# 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 \ 
    262268        -Df -N1/4/0/0/0/dashed -W1/0/0/0/solid -Ia/0.1/100 -K >>$psfile 
    263269 
    264270# plot cities 
    265271awk '{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 \ 
    267273        -W5/$citycol -G$citycol -Sc0.1 >>$psfile 
    268274awk '{print $1,$2,"14","0","4","1","\ ",$3}' $SH_UTIL/cities.dat \ 
    269         | pstext  -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt -G$citycol >>$psfile 
     275        | $GMT pstext  -: -O -K -R$lolon/$hilon/$lolat/$hilat $fmt -G$citycol >>$psfile 
    270276 
    271277# integrate text into psfile 
     
    279285        if  ($lcnt == 1)  set cont="" 
    280286        echo "6.0 $ypos.0 12 0.0 8 1  $line" | \ 
    281                 pstext -O $cont -JX10/8 -R0/10/0/100 >>$psfile 
     287                $GMT pstext -O $cont -JX10/8 -R0/10/0/100 >>$psfile 
    282288        @ lcnt = $lcnt - 1 
    283289        @ ypos = $ypos + 2 
Note: See TracChangeset for help on using the changeset viewer.