: # use perl
eval 'exec perl -S $0 "$@"'
if 0;

require "getopts.pl" ;

$usage = "Usage: $0 [-p phases] model [model ...]\n" ; 
if ( ! &Getopts('h:p:v') || @ARGV < 1)
    { die $usage ;  }

if ( $opt_p ) {
    $phases = $opt_p ; 
} else {
    $phases = "PP,Pb,PcP,Pg,Pn,Sb,ScP,Sg,Sn" ;
}

foreach $model ( @ARGV ) {
    $cmd = "remodl $model" ; 
    print STDERR "$cmd\n" ; 
    $result = system ( $cmd ) ; 

    if ( $result != 0 ) {
	warn ( "remodl failed\n" ) ; 
    } else {
	$cmd = "setbrn > /tmp/setbrn.out" ; 
	print STDERR "$cmd\n" ; 
	$result = system ( $cmd ) ; 

	if ( $result != 0 ) {
	    warn ( "setbrn failed\n" ) ; 
	} else {
	    $cmd = "taup2locsat -p firstP $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 

	    $cmd = "taup2locsat -p firstpP $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 

	    $cmd = "taup2locsat -p firstsP $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 

	    $cmd = "taup2locsat -p firstS $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 

	    $cmd = "taup2locsat -p firstpS $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 

	    $cmd = "taup2locsat -p firstsS $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 

	    $cmd = "taup2locsat -p $phases $model" ; 
	    print STDERR "$cmd\n" ; 
	    $result = system ( $cmd ) ; 
	    if ( $result != 0 ) {
		warn ( "$cmd failed\n" ) ; 
	    } 
	    
	    print STDERR "Saving standard Lg table\n" ;
	    $dsap = $ENV{'DSAP'} ;
	    $Lg = "$dsap/data/tables/dblocsat/$model.Lg" ; 
	    open ( LG, ">$Lg" ) ; 
	    while ( <DATA> ) {
		print LG ; 
	    }
	    close LG ;
	}
    }
}
__END__
n # Lg     travel-time (and amplitude) tables
 2    # number of depth samples
    0.00    5.00
 41    # number of distance samples
    0.00    1.00    2.00    3.00    4.00    5.00    6.00    7.00    8.00    9.00
   10.00   11.00   12.00   13.00   14.00   15.00   16.00   17.00   18.00   19.00
   20.00   21.00   22.00   23.00   24.00   25.00   26.00   27.00   28.00   29.00
   30.00   31.00   32.00   33.00   34.00   35.00   36.00   37.00   38.00   39.00
   40.00
# Travel-time/amplitude for z =    0.00
       0.000	Lg
      31.323	Lg
      62.645	Lg
      93.968	Lg
     125.290	Lg
     156.613	Lg
     187.935	Lg
     219.258	Lg
     250.580	Lg
     281.903	Lg
     313.225	Lg
     344.548	Lg
     375.870	Lg
     407.193	Lg
     438.516	Lg
     469.838	Lg
     501.161	Lg
     532.483	Lg
     563.806	Lg
     595.128	Lg
     626.451	Lg
     657.773	Lg
     689.096	Lg
     720.418	Lg
     751.741	Lg
     783.063	Lg
     814.386	Lg
     845.708	Lg
     877.031	Lg
     908.354	Lg
     939.676	Lg
     970.999	Lg
    1002.321	Lg
    1033.644	Lg
    1064.966	Lg
    1096.289	Lg
    1127.611	Lg
    1158.934	Lg
    1190.256	Lg
    1221.579	Lg
    1252.901	Lg
# Travel-time/amplitude for z =    5.00
       1.408	Lg
      31.354	Lg
      62.661	Lg
      93.978	Lg
     125.298	Lg
     156.619	Lg
     187.940	Lg
     219.262	Lg
     250.584	Lg
     281.906	Lg
     313.229	Lg
     344.551	Lg
     375.873	Lg
     407.195	Lg
     438.518	Lg
     469.840	Lg
     501.163	Lg
     532.485	Lg
     563.807	Lg
     595.130	Lg
     626.452	Lg
     657.775	Lg
     689.097	Lg
     720.420	Lg
     751.742	Lg
     783.065	Lg
     814.387	Lg
     845.710	Lg
     877.032	Lg
     908.355	Lg
     939.677	Lg
     971.000	Lg
    1002.322	Lg
    1033.645	Lg
    1064.967	Lg
    1096.290	Lg
    1127.612	Lg
    1158.935	Lg
    1190.257	Lg
    1221.580	Lg
    1252.902	Lg


# $Id: mktt.sh,v 1.2 1998/02/16 20:13:17 danq Exp $ 
