float **grid;

CGrid::CGrid (int numX, int numY)
{
grid=new float*[numX];
for (int i=0; i<numX; i++)
   grid[i]=new float[numY];
}

CGrid::~CGrid()
{
for ( int i = 0; i < numX; ++i )
   delete[] grid[i];
delete[] grid;
}

class CGrid   
{
public:
   //====== Конструктор экспорта
   CGrid (LPCTSTR fileName);
   //====== Конструктор для расчета
   CGrid (int b1, int b2);
   //Деструктор
   ~CGrid();

   float **grid;
   double minX, minY, xSize, ySize, minZ, maxZ;
   long numX, numY;

   float getX(int i) const;
   float getY(int j) const;
   void calculate(vector<CPoint3D> data, int approx, float n);

   float* CGrid::operator [] (int i)
   {
      assert (i >= 0 && i <= numX);
      return grid[i];
   }
private:
   float calc(float src, float a, int approx, float n);
   float distance(float x1, float y1, float x2, float y2);
   int   lags,
      num_close_wells;

};

CGrid::CGrid (int b1, int b2)
{
   lags=14;
   numX=b1;
   numY=b2;
   num_close_wells=10;
   grid=new float*[numX];
   for (int i=0; i<numX; i++)
      grid[i]=new float[numY];
}
CGrid::CGrid (LPCTSTR fileName)
{
   FILE *f;
   f = fopen(fileName, "rb");

   fseek(f, 5*4, SEEK_SET);
   fread(&numY, sizeof(long), 1, f);
   fread(&numX, sizeof(long), 1, f);
   fread(&minX, sizeof(double), 1, f);
   fread(&minY, sizeof(double), 1, f);
   fread(&xSize, sizeof(double), 1, f);
   fread(&ySize, sizeof(double), 1, f);
   fread(&minZ, sizeof(double), 1, f);
   fread(&maxZ, sizeof(double), 1, f);
   fseek(f, 24, SEEK_CUR);

   grid=new float*[numX];
   for (int i=0; i<numX; i++)
      grid[i]=new float[numY];

   double *tmp = new double[numX*numY];
   fread(tmp, sizeof(double), numX*numY, f);
   int k=0;
   for(int j=0; j<numY; j++)
      for(int i=0; i<numX; i++,k++)
         grid[i][j]=tmp[k];
   delete tmp;
   fclose(f);
}
CGrid::~CGrid()
{
   for ( int i = 0; i < numX; ++i )
      delete[] grid[i];
   delete[] grid;
}

class CGrid
{
private:
    float *grid_;
public:
    CGrid() { /*...*/ grid_ = new float[numX*numY]; }
    float *operator[] ( size_t i ) { return &grid_[i*numY]; }
};