/*=========================================================================

  Program:   Visualization Toolkit
  Module:    RGrid.cxx

  Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
  All rights reserved.
  See Copyright.txt or http://www.kitware.com/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notice for more information.

=========================================================================*/
// This example shows how to create a rectilinear grid.
//

#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkFloatArray.h"
#include "vtkRectilinearGrid.h"
#include "vtkRectilinearGridGeometryFilter.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkProperty.h"
#include "vtkCamera.h"

int main()
{
  int i;
  static double x[47]={
         -1.22396, -1.17188, -1.11979, -1.06771, -1.01562, -0.963542,
         -0.911458, -0.859375, -0.807292, -0.755208, -0.703125, -0.651042,
         -0.598958, -0.546875, -0.494792, -0.442708, -0.390625, -0.338542,
         -0.286458, -0.234375, -0.182292, -0.130209, -0.078125, -0.026042,
          0.0260415, 0.078125, 0.130208, 0.182291, 0.234375, 0.286458,
          0.338542, 0.390625, 0.442708, 0.494792, 0.546875, 0.598958,
          0.651042, 0.703125, 0.755208, 0.807292, 0.859375, 0.911458,
          0.963542, 1.01562, 1.06771, 1.11979, 1.17188};
  static double y[33]={-1.25, -1.17188, -1.09375, -1.01562, -0.9375, -0.859375,
         -0.78125, -0.703125, -0.625, -0.546875, -0.46875, -0.390625,
         -0.3125, -0.234375, -0.15625, -0.078125, 0, 0.078125,
          0.15625, 0.234375, 0.3125, 0.390625, 0.46875, 0.546875,
          0.625, 0.703125, 0.78125, 0.859375, 0.9375, 1.01562,
          1.09375, 1.17188, 1.25};
  static double z[44]={0, 0.1, 0.2, 0.3, 0.4, 0.5,
         0.6, 0.7, 0.75, 0.8, 0.9, 1,
         1.1, 1.2, 1.3, 1.4, 1.5, 1.6,
         1.7, 1.75, 1.8, 1.9, 2, 2.1,
         2.2, 2.3, 2.4, 2.5, 2.6, 2.7,
         2.75, 2.8, 2.9, 3, 3.1, 3.2,
         3.3, 3.4, 3.5, 3.6, 3.7, 3.75,
         3.8, 3.9};

  // Create a rectilinear grid by defining three arrays specifying the
  // coordinates in the x-y-z directions.
  vtkFloatArray *xCoords = vtkFloatArray::New();
  for (i=0; i<47; i++)
  {
    xCoords->InsertNextValue(x[i]);
  }
  vtkFloatArray *yCoords = vtkFloatArray::New();
  for (i=0; i<33; i++)
  {
    yCoords->InsertNextValue(y[i]);
  }

  vtkFloatArray *zCoords = vtkFloatArray::New();
  for (i=0; i<44; i++)
  {
    zCoords->InsertNextValue(z[i]);
  }

  // The coordinates are assigned to the rectilinear grid. Make sure that
  // the number of values in each of the XCoordinates, YCoordinates,
  // and ZCoordinates is equal to what is defined in SetDimensions().
  //
  vtkRectilinearGrid *rgrid = vtkRectilinearGrid::New();
    rgrid->SetDimensions(47,33,44);
    rgrid->SetXCoordinates(xCoords);
    rgrid->SetYCoordinates(yCoords);
    rgrid->SetZCoordinates(zCoords);

    // Extract a plane from the grid to see what we've got.
  vtkRectilinearGridGeometryFilter *plane = vtkRectilinearGridGeometryFilter::New();
    plane->SetInputData(rgrid);
    plane->SetExtent(0,46, 16,16, 0,43);

  vtkPolyDataMapper *rgridMapper = vtkPolyDataMapper::New();
      rgridMapper->SetInputConnection(plane->GetOutputPort());

  vtkActor *wireActor = vtkActor::New();
      wireActor->SetMapper(rgridMapper);
      wireActor->GetProperty()->SetRepresentationToWireframe();
      wireActor->GetProperty()->SetColor(0,0,0);

  // Create the usual rendering stuff.
  vtkRenderer *renderer = vtkRenderer::New();
  vtkRenderWindow *renWin = vtkRenderWindow::New();
    renWin->AddRenderer(renderer);
  vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
    iren->SetRenderWindow(renWin);

  renderer->AddActor(wireActor);
  renderer->SetBackground(1,1,1);
  renderer->ResetCamera();
  renderer->GetActiveCamera()->Elevation(60.0);
  renderer->GetActiveCamera()->Azimuth(30.0);
  renderer->GetActiveCamera()->Zoom(1.0);

  renWin->SetSize(300,300);

  // interact with data
  renWin->Render();
  iren->Start();

  // Clean up
  plane->Delete();
  renderer->Delete();
  renWin->Delete();
  iren->Delete();
  xCoords->Delete();
  yCoords->Delete();
  zCoords->Delete();
  rgrid->Delete();
  rgridMapper->Delete();
  wireActor->Delete();

  return 0;
}
