Program Hypotable: iaspei-tauThe original 1991 release of the iaspei-tau package and its 1992 update are hard-wired to handle only the iasp91 velocity structure, which it did through polynomial fits within designated layers for the velocities and density. In 1991 I changed subroutine emdlv.f so that it could read in a v-z model from a file. (In 1992 I replaced the linear interpolation in this subroutine with cubic-spline interpolation.) File ./build-tables/iasp91.mod is the iasp91 model in the format expected by emdlv. The February 1996 version of the package added a few core phases and included a variant of my emdlv.f for reading velocity models from files. (It also includes model and table files for model ak135.) For more information about the iaspei-tau package, look at ./iaspei-tau/README_tau, written (mostly) by Ray Buland (USGS). The 1996 release for Unix and earlier releases for other operating systems can be gotten by anonymous ftp from the USGS. (The 1996 release used here is in the subdirectory experimental.) Program Hypotable requires three input files: the ascii text velocity model file (xxx.mod), and two binary files produced by the iaspei-tau package (xxx.hed and xxx.tbl). If one simply wants to use the iasp91 model and iaspei-tau-package format tables, one can just use the files for the iasp91 model which come with this distribution. If one wants to use a different velocity structure, one follows the following procedure: First one creates an ascii model file in the same format as iasp91.mod. One then runs two program from the iaspei-tau package: first remodl and then setbrn. All other input comes from prompts. The output from this two-program sequence is two binary files: here IASP91.hed and IASP91.tbl. The code is optimized for the iasp91 model, so this procedure may not work for models which differ significantly from that model. (Using rsetbrn a stripped-down version of setbrn with only the main P and S branches, may help.) Cases for which this procedure has been applied by the authors have not used core phases, so the models differed from iasp91 only in the (two-layer) crust and the uppermost mantle. Linking program Hypotable requires subroutines included in libiasp.a (which is produced by running the Makefile in ./iaspei-tau). As written, it requires that the the .tbl, .mod, and .hed files are all in the directory from which the program is being run. The files (with their most recent modification dates) in directory ./iaspei-tau are as follows:
Modif. date name comment Jun 22 1992 emdlv.f jas/vt version of the iasp91 subroutine Apr 29 2002 hypotable.f source code in f77 Jun 5 1992 libsun.f part of the iasp91 package Oct 28 1998 libtau.f part of the iasp91 package Jun 5 1992 limits.inc part of the iasp91 package Apr 21 2001 Makefile unix makefile for libiasp.a and all executables Jun 17 1995 README_tau Ray Buland's readme file for the iasp91 package Jan 14 2000 remodl.f source file in iasp91 package Jan 14 2000 rsetbrn.f source for stripped-down setbern (jas/vt) Jan 14 2000 setbrn.f source file in iasp91 package Jun 22 1993 tau.man man page realted to iasp91 package Apr 26 2002 ttimes.f source file for my version of ttimes Jun 5 1992 ttlim.inc part of iasp91 package |