Polarization.cpp

Go to the documentation of this file.
00001 #include <iostream>
00002 #include <fstream>
00003 #include <iomanip>
00004 #include "CMtuFormat.h"
00005 using namespace std;
00006 
00007 int main()
00008 {
00009         string infilename;
00010         CMtuFormat TsData;
00011         ofstream hpol, epol, exout, eyout, hxout,hyout;
00012         double phi;
00013         const int npoints = 480000;
00014         
00015         cout << "Ts File: ";
00016         cin >> infilename;
00017         TsData.GetData(infilename.c_str());
00018         hpol.open((infilename+"_hpol.dat").c_str());
00019         epol.open((infilename+"_epol.dat").c_str());
00020         exout.open((infilename+"_ex.dat").c_str());
00021         eyout.open((infilename+"_ey.dat").c_str());
00022         hxout.open((infilename+"_hx.dat").c_str());
00023         hyout.open((infilename+"_hy.dat").c_str());
00024         for (int i = 0; i < npoints; ++i)
00025         {
00026                 phi = atan(TsData.Ey.data.at(i)/TsData.Ex.data.at(i));
00027                 epol << setprecision(15) << setw(18) << TsData.t.at(i)  - TsData.t.at(0)<< " " << phi << endl;
00028                 phi = atan(TsData.Hy.data.at(i)/TsData.Hx.data.at(i));
00029                 hpol <<  setprecision(15) << setw(18) << TsData.t.at(i) - TsData.t.at(0)<< " " << phi << endl;
00030                 exout <<  setprecision(15) << setw(18) << TsData.t.at(i) - TsData.t.at(0) << " " << TsData.Ex.data.at(i) << endl;
00031                 eyout <<  setprecision(15) << setw(18) << TsData.t.at(i) - TsData.t.at(0) << " " << TsData.Ey.data.at(i) << endl;
00032                 hxout <<  setprecision(15) << setw(18) << TsData.t.at(i) - TsData.t.at(0) << " " << TsData.Hx.data.at(i) << endl;
00033                 hyout <<  setprecision(15) << setw(18) << TsData.t.at(i) - TsData.t.at(0) << " " << TsData.Hy.data.at(i) << endl;
00034         }
00035         epol.close();
00036         hpol.close();
00037         exout.close();
00038         eyout.close();
00039         hxout.close();
00040         hyout.close();
00041 }