keepcommon.cpp

Go to the documentation of this file.

#include <iostream>
#include <vector>
#include "MTStationList.h"

using namespace std;
string infilename;

int main()
  {
    try
      {
        MTStationList MTSites;

        cout << "Station Filename: ";
        cin >> infilename;

        MTSites.GetData(infilename);
        const unsigned int nstats = MTSites.GetList().size();

        for (unsigned int i = 0; i < nstats; ++i)
          {
            trealdata freqs(MTSites.at(i).GetFrequencies());
            const int nfreq =MTSites.at(i).GetFrequencies().size();
            for (int j = 0; j < nfreq; ++j)
              {
                if ( find(MTSites.GetComFreqIndices().at(i).begin(),MTSites.GetComFreqIndices().at(i).end(),j)
                    == MTSites.GetComFreqIndices().at(i).end())
                freqs.at(j) = 0;
              }
            MTSites.at(i).SetFrequencies(freqs);
            MTSites.at(i).WriteAsMtt(MTSites.at(i).GetName());
          }
      }
    catch(CFatalException &e)
      {
        cerr << e.what() << endl;
        return -1;
      }
  }

---