writevtk.cpp

Go to the documentation of this file.

#include <iostream>
#include <string>
#include "ThreeDMTModel.h"
#include "FatalException.h"
#include <fstream>

using namespace std;
using namespace gplib;

string version = "$Id: writevtk.cpp 1816 2009-09-07 11:28:35Z mmoorkamp $";

int main(void)
  {
    cout
        << " This is writevtk, convert Mackie or WinGlink model files to vtk format for display"
        << endl;
    cout
        << " Output format is legacy VTK format, ending '.vtk' is appended to output"
        << endl;
    cout << " This is Version: " << version << endl << endl;
    ThreeDMTModel Model;
    string infilename, outfilename;

    int mode;

    cout << "Modelfile: ";
    cin >> infilename;
    cout << "Type 1 for old Mackie, 2 for Winglink: ";
    cin >> mode;

    cout << "Outfile: ";
    cin >> outfilename;
    try
      {
        if (mode == 1)
          Model.ReadMackie(infilename);
        else
          Model.ReadWinGLink(infilename);
        Model.WriteVTK(outfilename + ".vtk");
        Model.WriteNetCDF(outfilename + ".nc");
      } catch (const FatalException &e)
      {
        cerr << e.what() << endl;
      }
  }