Method get_old_grav_vector

Source/gravity/Gravity.cpp:853–928

}

void
Gravity::get_old_grav_vector(int level, MultiFab& grav_vector, Real time)
{
    BL_PROFILE("Gravity::get_old_grav_vector()");

    int ng = grav_vector.nGrow();

    // Fill data from the level below if we're not doing a solve on this level.

    if (level > gravity::max_solve_level) {

        LevelData[level]->FillCoarsePatch(grav_vector,0,time,Gravity_Type,0,3,ng);

        return;

    }

    // Note that grav_vector coming into this routine always has three components.
    // So we'll define a temporary MultiFab with AMREX_SPACEDIM dimensions.
    // Then at the end we'll copy in all AMREX_SPACEDIM dimensions from this into
    // the outgoing grav_vector, leaving any higher dimensions unchanged.

    MultiFab grav(grids[level], dmap[level], AMREX_SPACEDIM, ng);
    grav.setVal(0.0,ng);

    const Geometry& geom = parent->Geom(level);

    if (gravity::gravity_type == "ConstantGrav") {

        if (AMREX_SPACEDIM == 2 && geom.Coord() == 2) {
            // 2D spherical r-theta, we want g in the radial direction
            grav.setVal(gravity::const_grav, 0, 1, ng);
        } else {
            // Set to constant value in the AMREX_SPACEDIM direction and zero in all others.
            grav.setVal(gravity::const_grav, AMREX_SPACEDIM-1, 1, ng);
        }

    } else if (gravity::gravity_type == "MonopoleGrav") {

       const Real prev_time = LevelData[level]->get_state_data(State_Type).prevTime();
       make_radial_gravity(level,prev_time,radial_grav_old[level]);
       interpolate_monopole_grav(level,radial_grav_old[level],grav);

    } else if (gravity::gravity_type == "PoissonGrav") {

       amrex::average_face_to_cellcenter(grav, amrex::GetVecOfConstPtrs(grad_phi_prev[level]), geom);
       grav.mult(-1.0, ng); // g = - grad(phi)

    } else {
       amrex::Abort("Unknown gravity_type in get_old_grav_vector");
    }

    // Do the copy to the output vector.

    for (int dir = 0; dir < 3; dir++) {
        if (dir < AMREX_SPACEDIM) {
            MultiFab::Copy(grav_vector, grav, dir, dir, 1, ng);
        } else {