Attempting to use TauP on a Windows system

I say attempting because the boffins that wrote this have no appreciation of the real world as amply demonstrated by the fact that they consider a GUI would be a useful addition (just a little!) and that it is a problem as browsers don’t support Java 1.1 (No the world moved on and Java is now 7.21 but never mind).

In attempting to run this I unpacked it and tried to run taup_time.bat. It failed miserably because the program was not written to take account of modern folder naming. These days we actually are allowed to have spaces in the names. Cute isn’t it?

I moved the TauP-2.1.1 folder and it’s contents to the root of a drive and tried again.

Please bear in mind here that this is a command line package and thus you need to be in at the command line and NOT try running this from the Windows environment. In XP this means going to the Run part of the Start Bar

taup01

and entering cmd and pressing return.

taup02

Once in cmd you must then navigate to the TauP system folder. If you have it placed on a drive other then C: (as I have) you must first change to that drive by entering the drive letter and a colon and pressing return.

The you change directories to the system files for Taup. You do this using the CD command

CD TauP-2.1.1\bin

You should then be at the command prompt – >

E:\TauP-2.1.1\bin>taup_time.bat
Enter:
h for new depth
r to recalculate
p to append phases,
c to clear phases
l to list phases
s for new station lat lon
e for new event lat lon
a for new azimuth
b for new back azimuth
m for new model or
q to quit.

Enter Distance or Option [hrpclseabmq]:

Note that the > cursor is replaced by a colon, so we are not in the program.

After a bit of experimenting I discovered that you enter the option and press return and then enter the data. I entered the depth, station lat/lon and event lat/lon

Here is what the cmd window looked like just before I hit return to start it.

E:\>cd taup-2.1.1\bin

E:\TauP-2.1.1\bin>taup_time.bat
Enter:
h for new depth
r to recalculate
p to append phases,
c to clear phases
l to list phases
s for new station lat lon
e for new event lat lon
a for new azimuth
b for new back azimuth
m for new model or
q to quit.

Enter Distance or Option [hrpclseabmq]: h
Enter Depth: 20
Enter Distance or Option [hrpclseabmq]: s
Enter station lat and lon: 74.69 -94.9
Enter Distance or Option [hrpclseabmq]: e
Enter event lat and lon: 25.79 121.08
Enter Distance or Option [hrpclseabmq]: l
16 phases.
p,s,P,S,Pn,Sn,PcP,ScS,Pdiff,Sdiff,PKP,SKS,PKiKP,SKiKS,PKIKP,SKIKS
Enter Distance or Option [hrpclseabmq]:

The earthquake event was the Magnit 6.2 in Taiwan on 02 Jun 2013 and the station was Resolute in Canada. CN.RES

Note that I used the letter L to see what phases were in there. I pressed return to set it going.

After 20 minutes waiting I gave up and hit Ctrl+C to stop it. Lo and behold the results popped up

Distance   Depth   Phase   Travel    Ray Param  Takeoff  Incident  Purist    Purist
  (deg)     (km)   Name    Time (s)  p (s/deg)   (deg)    (deg)   Distance   Name
-----------------------------------------------------------------------------------
   76.86    20.0   P        710.59     5.636     19.30    17.10    76.86   = P

   76.86    20.0   PcP      721.46     4.326     14.70    13.04    76.86   = PcP

   76.86    20.0   PKiKP   1053.00     1.506      5.07     4.51    76.86   = PKiKP
   76.86    20.0   S       1297.82    10.900     21.64    19.23    76.86   = S

   76.86    20.0   SKS     1320.70     7.223     14.14    12.61    76.86   = SKS

   76.86    20.0   ScS     1325.91     8.080     15.86    14.13    76.86   = ScS

   76.86    20.0   SKiKS   1480.80     1.620      3.14     2.81    76.86   = SKiKS

Problem therefore is how long do you leave it?

I looked at the rudimentary GUI version and strongly suggest you don’t go there as you can only enter the number of degrees between the station and the event so unless you want to measure that – inaccurately as it transpired – on Google Earth, you are stuck.

This is NOT in any way shape or form a replacement for the USGS TT calculator. Many users have no idea how to use the command line, and I am not even sure if you can get to it in Windows 8.

My final test was to use the command line to send to a file as that way one might be able to determine when it had finished.

Using the same data I entered the command line so:

taup_time.bat -h 20 -sta 74.69 -94.9 -evt 23.79 121.08 -o e:\taiwan-tt.txt

In a fraction of a second the prompt retuned and the file had been created

Model: iasp91
Distance   Depth   Phase   Travel    Ray Param  Takeoff  Incident  Purist    Purist
  (deg)     (km)   Name    Time (s)  p (s/deg)   (deg)    (deg)   Distance   Name 
-----------------------------------------------------------------------------------
   78.84    20.0   P        721.58     5.488     18.77    16.63    78.84   = P    
   78.84    20.0   PcP      730.02     4.348     14.77    13.11    78.84   = PcP  
   78.84    20.0   PKiKP   1056.01     1.534      5.16     4.59    78.84   = PKiKP
   78.84    20.0   S       1319.09    10.656     21.13    18.78    78.84   = S    
   78.84    20.0   SKS     1334.81     7.063     13.82    12.32    78.84   = SKS  
   78.84    20.0   ScS     1341.91     8.125     15.95    14.21    78.84   = ScS  
   78.84    20.0   SKiKS   1484.03     1.647      3.19     2.85    78.84   = SKiKS

One thing in the manual amused me

taup.latlon.precision precision for latitude and longitude output, the default is 2 decimal digits. Note that this is precision, not accuracy. Just because you get more digits doesn’t imply that they have any meaning.

(Emphasis by me – I do realise this is about output)

No, but inputting more than 2 decimal digits sure improves the accuracy and those extra digits have meaning whatever comes out of the sausage machine at the other end.

Don’t forget as always this is all theoretical.

Advertisements

About PuterMan

A retired programmer.
This entry was posted in Reports and tagged . Bookmark the permalink.

One Response to Attempting to use TauP on a Windows system

  1. highace says:

    I’m pleased you figured it out, but that stuff is beyond me. Lets hope they bring the TT Calculator back at some stage.

Leave a Reply (All Comments are moderated by a human)

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s