test_gridded_fields.cc

Go to the documentation of this file.

#include "gridded_fields.h"
#include "matpack.h"

using std::cout;
using std::endl;

int main(void) {
  // Creating two gridded fields
  GriddedField1 gfone("I'm a GriddedField1");
  GriddedField2 gftwo;

  gftwo.set_name("I'm a GriddedField2");

  // Initializing the grids
  Vector gfonegrid(1, 5, 1);     // gfonegrid = [1,2,3,4,5]
  gfone.set_grid(0, gfonegrid);  // Set grid for the vector elements.

  Vector gftwogrid0(1, 5, 1);  // gftwogrid0 = [1,2,3,4,5]
  ArrayOfString gftwogrid1{"Chan1", "Chan2", "Chan3"};

  gftwo.set_grid(0, gftwogrid0);  // Set grid for the matrix rows.
  gftwo.set_grid(1, gftwogrid1);  // Set grid for the matrix columns.

  gfone.set_grid_name(0, "Pressure");

  gftwo.set_grid_name(0, "Pressure");
  gftwo.set_grid_name(1, "Channel");

  // Initializing the data
  Vector avector(1, 4, 0.5);  // avector = [1,1.5,2,2.5]

  gfone.data = avector;

  cout << gfone;

  Matrix amatrix(5, 3, 4.);  // amatrix = [[4,4,4],[4,4,4],...]

  gftwo.data = amatrix;

  // Consistency check

  if (!gfone.checksize())
    cout << gfone.get_name() << ": Sizes of grid and data don't match" << endl;

  // This should fail!
  if (!gftwo.checksize())
    cout << gftwo.get_name() << ": Sizes of grids and data don't match" << endl;

  // Output

  cout << "GriddedField1: " << gfone << endl;
  cout << "GriddedField2: " << gftwo << endl;

  return 0;
}