Function WriteFab

Tools/Postprocessing/C_Src/PlotfileToMatLab.cpp:93–127

}

static
void
WriteFab (std::ostream&         os,
          const FArrayBox& fab,
          const char*      name)

{
  int nx = fab.box().length(0);
  int ny = fab.box().length(1);
  int dim = BL_SPACEDIM;


#if (BL_SPACEDIM == 2)
  for (int i = 0; i < nx; i++)
    for (int j = 0; j < ny; j++)
    {
      int index = j*nx + i;
      const Real * ptr = fab.dataPtr();
      os.write((char*)(ptr+index),sizeof(Real));
    }

#elif (BL_SPACEDIM == 3)
  int nz = fab.box().length(2);
  os.write((char*)&nz,sizeof(int));

  for (int k = 0; k < nz; k++)
    for (int i = 0; i < nx; i++)
      for (int j = 0; j < ny; j++)
      {
        int index = k*(nx*ny) + j*nx + i;
        const Real * ptr = fab.dataPtr();
        os.write((char*)(ptr+index),sizeof(Real));
      }
#endif
}

static
void