dike.cpp

Go to the documentation of this file.
00001 #include "CAnalyticDike.h"
00002 #include "CTwoQuarter.h"
00003 #include <iostream>
00004 #include <fstream>
00005 
00006 using namespace std;
00007 
00008 CAnalyticDike   AnalyticDike;
00009 CTwoQuarter             TwoQuarter;
00010 
00011 int main (void)
00012 {
00013         ofstream outfile;
00014         AnalyticDike.rho1 = 0.1;
00015         AnalyticDike.rho2 = 160;
00016         AnalyticDike.length = 1000;
00017         AnalyticDike.height = 1000;
00018         AnalyticDike.sub = 0;
00019         TwoQuarter.rho1 = 1;
00020         TwoQuarter.rho2 = 100;
00021         TwoQuarter.height = 1000;
00022         TwoQuarter.sub = 0;
00023         
00024         outfile.open("Zphase.out");
00025         for (int i = 0; i < 1000; i+=10)
00026         {
00027                 AnalyticDike.CalcResponse(1.28e-4,i);
00028                 outfile << i << " " << arg(AnalyticDike.Z1) *180/PI << " " << 
00029                         arg(AnalyticDike.Z2) *180/PI << endl;
00030                 cout << "Z1: " << AnalyticDike.Z1 << " Phase: " << arg(AnalyticDike.Z1) *180/PI << endl;
00031                 cout << "Z2: " << AnalyticDike.Z2 << " Phase: " << arg(AnalyticDike.Z2) *180/PI << endl;
00032         }
00033         /*for (int i = -1000; i < 1000; i+=50)
00034         {
00035                 TwoQuarter.CalcResponse(1.28,i);
00036                 outfile << i << " " << arg(TwoQuarter.Z1) *180/PI +180<< " " << 
00037                         arg(TwoQuarter.Z2) *180/PI +180<< endl;
00038                 cout << "Z1: " << TwoQuarter.Z1 << " Phase: " << arg(TwoQuarter.Z1) *180/PI << endl;
00039                 cout << "Z2: " << TwoQuarter.Z2 << " Phase: " << arg(TwoQuarter.Z2) *180/PI << endl;
00040         }*/
00041         outfile.close();
00042         return(0);
00043 }