00001 #include "ResPkModel.h"
00002 #include <iostream>
00003 #include <string>
00004 #include "Util.h"
00005 using namespace std;
00006 int main()
00007 {
00008 const double refslow = 0.05;
00009 ResPkModel Model;
00010 string modelfilename = AskFilename("Model Filename: ");
00011 Model.GetData(modelfilename);
00012 double sourcedepth = 0.0;
00013 double recdepth = 0.0;
00014 cout << "Source Depth : ";
00015 cin >> sourcedepth;
00016 cout << "Receiver Depth: ";
00017 cin >> recdepth;
00018 double stime = Model.CalcTravelTime(SeismicModel::DirectS, sourcedepth,
00019 recdepth, refslow);
00020 double ptime = Model.CalcTravelTime(SeismicModel::DirectP, sourcedepth,
00021 recdepth, refslow);
00022 cout << "Direct S: " << stime << " s Direct P: " << ptime << " s" << endl;
00023 }
