Function main

Tools/Plotfile/faverage.cpp:19–265

int main(int argc, char* argv[])
{

    amrex::Initialize(argc, argv, false);

    {

        // timer for profiling
        BL_PROFILE_VAR("main()", pmain);

        // Input arguments
        std::string varname{"density"};
        bool do_favre{false};

        const int narg = amrex::command_argument_count();

        int farg = 1;
        while (farg <= narg) {
            const std::string& name = amrex::get_command_argument(farg);
            if (name == "-v" || name == "--variable") {
                varname = amrex::get_command_argument(++farg);
            } else if (name == "-f" || name == "--favre") {
                do_favre = true;
            } else {
                break;
            }
            ++farg;
        }

        if (farg > narg) {
            PrintUsage();
            return -1;
        }

        const std::string& pltfile = amrex::get_command_argument(farg);
        std::string slcfile = pltfile + ".slice";

        PlotFileData pf(pltfile);

        int fine_level = pf.finestLevel();
        const int dim = pf.spaceDim();

        // get the index bounds and dx.
        Box domain = pf.probDomain(fine_level);
        auto dx_fine = pf.cellSize(fine_level);
        auto problo = pf.probLo();
        auto probhi = pf.probHi();

        // compute the size of the radially-binned array -- we'll take the
        // vertical direction to be the dimensionality

        int nbins = static_cast<int>(std::abs(probhi[dim-1] - problo[dim-1]) / dx_fine[dim-1]);

        // height coordinate
        Vector<Real> h(nbins);

        for (auto i = 0; i < nbins; i++) {
            h[i] = (i + 0.5) * dx_fine[dim-1];