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
and entering cmd and pressing return.
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
You should then be at the command prompt – >
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.