keepcommon.cpp

Go to the documentation of this file.

#include <iostream>
#include <vector>
#include "MTStationList.h"
#include "Util.h"
using namespace std;
using namespace gplib;

int main()
  {
    try
      {
        MTStationList MTSites;
        std::string infilename = AskFilename("Station Filename: ");

        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 (FatalException &e)
      {
        cerr << e.what() << endl;
        return -1;
      }
  }