def _getCellAreas(self):
areas = numerix.ones((4, self.numberOfCells), 'd')
areas[0] = self.dx * self.getCellCenters()[0]
areas[1] = self.dy * (self.getCellCenters()[0] + self.dx / 2)
areas[2] = self.dx * self.getCellCenters()[0]
areas[3] = self.dy * (self.getCellCenters()[0] - self.dx / 2)
return areas
def getCellVolumes(self):
return numerix.ones(self.numberOfCells, 'd') * self.dx * self.dy * self.getCellCenters()[0]
