writevtk.cpp

Go to the documentation of this file.

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

using namespace std;

string version = "$Id: writevtk.cpp 1453 2007-11-13 13:40:38Z 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;
        const double kmlatitude = 111194;
        const double kmlongitude = 71474;       
        
        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 CFatalException &e)
        {
                cerr << e.what() << endl;
        }
}

---