dike.cpp

Go to the documentation of this file.

#include "CAnalyticDike.h"
#include "CTwoQuarter.h"
#include <iostream>
#include <fstream>

using namespace std;

CAnalyticDike   AnalyticDike;
CTwoQuarter             TwoQuarter;

int main (void)
{
        ofstream outfile;
        AnalyticDike.rho1 = 0.1;
        AnalyticDike.rho2 = 160;
        AnalyticDike.length = 1000;
        AnalyticDike.height = 1000;
        AnalyticDike.sub = 0;
        TwoQuarter.rho1 = 1;
        TwoQuarter.rho2 = 100;
        TwoQuarter.height = 1000;
        TwoQuarter.sub = 0;
        
        outfile.open("Zphase.out");
        for (int i = 0; i < 1000; i+=10)
        {
                AnalyticDike.CalcResponse(1.28e-4,i);
                outfile << i << " " << arg(AnalyticDike.Z1) *180/PI << " " << 
                        arg(AnalyticDike.Z2) *180/PI << endl;
                cout << "Z1: " << AnalyticDike.Z1 << " Phase: " << arg(AnalyticDike.Z1) *180/PI << endl;
                cout << "Z2: " << AnalyticDike.Z2 << " Phase: " << arg(AnalyticDike.Z2) *180/PI << endl;
        }
        /*for (int i = -1000; i < 1000; i+=50)
        {
                TwoQuarter.CalcResponse(1.28,i);
                outfile << i << " " << arg(TwoQuarter.Z1) *180/PI +180<< " " << 
                        arg(TwoQuarter.Z2) *180/PI +180<< endl;
                cout << "Z1: " << TwoQuarter.Z1 << " Phase: " << arg(TwoQuarter.Z1) *180/PI << endl;
                cout << "Z2: " << TwoQuarter.Z2 << " Phase: " << arg(TwoQuarter.Z2) *180/PI << endl;
        }*/
        outfile.close();
        return(0);
}

---