raster_grid.cpp

Go to the documentation of this file.

 
/* ===============================================================================================================================
   This file is part of CoastalME, the Coastal Modelling Environment.
 
   CoastalME is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version.
 
   This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
===============================================================================================================================*/
#include <cstdio>
 
#include <new>
using std::bad_alloc;
 
#include "cme.h"
#include "cell.h"
#include "simulation.h"
#include "raster_grid.h"
 
CGeomRasterGrid* CGeomCell::m_pGrid = NULL; // Initialise m_pGrid, the static member of CGeomCell

CGeomRasterGrid::CGeomRasterGrid(CSimulation* pSimIn)
    : m_dD50Fine(0),
      m_dD50Sand(0),
      m_dD50Coarse(0),
      m_pSim(pSimIn),
      m_Cell(NULL)
{
}

CGeomRasterGrid::~CGeomRasterGrid(void)
{
   int const nXMax = m_pSim->nGetGridXMax();
 
   // Free the m_Cell memory
   for (int nX = 0; nX < nXMax; nX++)
      delete[] m_Cell[nX];
 
   delete[] m_Cell;
}

CSimulation* CGeomRasterGrid::pGetSim(void)
{
   return m_pSim;
}
 
// CGeomCell* CGeomRasterGrid::pGetCell(int const nX, int const nY)
// {
// return &m_Cell[nX][nY];
// }

int CGeomRasterGrid::nCreateGrid(void)
{
   // Create the 2D CGeomCell array (this is faster than using 2D STL vectors)
   int const nXMax = m_pSim->nGetGridXMax();
   int const nYMax = m_pSim->nGetGridYMax();
 
   // TODO 038 Do better error handling if insufficient memory
   try
   {
      m_Cell = new CGeomCell*[nXMax];
 
      for (int nX = 0; nX < nXMax; nX++)
         m_Cell[nX] = new CGeomCell[nYMax];
   }
 
   catch (bad_alloc&)
   {
      // Uh-oh, not enough memory
      return RTN_ERR_MEMALLOC;
   }
 
   // Initialize the CGeomCell shared pointer to the CGeomRasterGrid object
   CGeomCell::m_pGrid = this;
 
   return RTN_OK;
}