Function main

Diagnostics/Sedov/main.cpp:118–382

}

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

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

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

    // Input arguments
    std::string pltfile, slcfile;
    bool sphr = false;

    GetInputArgs (argc, argv, pltfile, slcfile, sphr);

    auto center = GetCenter(pltfile);
    double xctr = center[0];
    double yctr = center[1];
    double zctr = center[2];

    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 = pf.cellSize(fine_level);
    auto problo = pf.probLo();
    auto probhi = pf.probHi();
    int coord = pf.coordSys();

    // compute the size of the radially-binned array -- we'll do it to
    // the furtherest corner of the domain

#if (AMREX_SPACEDIM == 1)
    double maxdist = std::abs(probhi[0] - problo[0]);

#elif (AMREX_SPACEDIM == 2)
    double x_maxdist = amrex::max(std::abs(probhi[0] - xctr),
                                  std::abs(problo[0] - xctr));
    double y_maxdist = amrex::max(std::abs(probhi[1] - yctr),
                                  std::abs(problo[1] - yctr));
    double maxdist = std::sqrt(x_maxdist*x_maxdist +
                               y_maxdist*y_maxdist);

    // For spherical 2D, the max distance away from blast center
    // always happens at max(θ - θ_0) for 0 < θ < π.
    if (coord == 2) {
        Real theta_maxdist = amrex::max(std::abs(problo[1] - yctr),
                                        std::abs(probhi[1] - yctr));
        maxdist = probhi[0] * probhi[0] + xctr * xctr -
                  2.0_rt * probhi[0] * xctr * std::cos(theta_maxdist);
    }
#else
    double x_maxdist = amrex::max(std::abs(probhi[0] - xctr),
                                  std::abs(problo[0] - xctr));
    double y_maxdist = amrex::max(std::abs(probhi[1] - yctr),
                                  std::abs(problo[1] - yctr));