Function main

Tools/Postprocessing/C_Src/IntegrateComp.cpp:103–204

};

int
main (int   argc,
      char* argv[])
{
  amrex::Initialize(argc,argv);
  {
    if (argc < 2)
      print_usage(argc,argv);

    ParmParse pp;

    const std::string farg = amrex::get_command_argument(1);
    if (farg == "-h" || farg == "--help")
      print_usage(argc,argv);

    // Create the AmrData object from a pltfile on disk
    std::string infile; pp.get("infile",infile);
    DataServices::SetBatchMode();
    Amrvis::FileType fileType(Amrvis::NEWPLT);
    DataServices dataServices(infile, fileType);
    if( ! dataServices.AmrDataOk()) {
      DataServices::Dispatch(DataServices::ExitRequest, NULL);
    }
    AmrData& amrData = dataServices.AmrDataRef();

    int nv = pp.countval("varNames");
    Vector<std::string> varNames(nv); pp.getarr("varNames",varNames,0,nv);

    // Make a data struct for just the variables needed
    AMReXDataHierarchy data(amrData,varNames);
    const AMReXMeshHierarchy& mesh = data.Mesh();



    // Compute the volume integrals
    const int nGrow = 0;
    const int finestLevel = mesh.FinestLevel();
    const int nLev = finestLevel + 1;
    const int nComp = varNames.size();

    Vector<Real> integrals(nComp,0); // Results, initialized to zero
    for (int lev=0; lev<nLev; ++lev) {
      const BoxArray ba = mesh.boxArray(lev);

      // Make boxes that are projection of finer ones (if exist)
      const BoxArray baf = lev < finestLevel
                                 ? BoxArray(mesh.boxArray(lev+1)).coarsen(mesh.RefRatio()[lev])
                                 : BoxArray();

      // Compute volume of a cell at this level
      Real vol=1.0;
      for (int d=0; d<AMREX_SPACEDIM; ++d) {
        vol *= mesh.ProbSize()[d] / mesh.ProbDomain()[lev].length(d);
      }

      // For each component listed...
      for (int n=0; n<nComp; ++n) {