REMORA_Advance.cpp

Go to the documentation of this file.

#include <REMORA.H>
 
using namespace amrex;
 
#ifdef REMORA_USE_FUNWAVE_FORT
#include <REMORA_funwave_Fortran_Interface.H>
#endif
 
/**
 * @param[in] lev            level of refinement
 * @param[in] time           simulation time at start of step
 * @param[in] dt_lev         baroclinic time step at level
 * @param[in] iteration      iteration in subcycling, if using
 * @param[in] ncycle         total number of subcycles, if using
 */
 void
REMORA::Advance (int lev, Real time, Real dt_lev, int /*iteration*/, int /*ncycle*/)
{
    BL_PROFILE("REMORA::Advance()");
 
    setup_step(lev, time, dt_lev);
 
    int nfast_counter=nfast + 1;
 
    //***************************************************
    //Compute fast timestep from dt_lev and ratio
    //***************************************************
    Real dtfast_lev=dt_lev/Real(ndtfast);
 
    //***************************************************
    //Advance nfast_counter steps of the 2d integrator
    //***************************************************
    for (int my_iif = 0; my_iif < nfast_counter; my_iif++) {
        advance_2d_onestep(lev, dt_lev, dtfast_lev, my_iif, nfast_counter);
    }
 
#ifdef REMORA_USE_FUNWAVE_FORT
    MultiFab* mf_rhoS = vec_rhoS[lev].get();
    for ( MFIter mfi(*mf_rhoS, TilingIfNotGPU()); mfi.isValid(); ++mfi )
    {
        Box bx = mfi.validbox();
        int ims = bx.smallEnd(0);
        int jms = bx.smallEnd(1);
        int kms = bx.smallEnd(2);
        int ime = bx.bigEnd(0);
        int jme = bx.bigEnd(1);
        int kme = bx.bigEnd(2);
 
        Array4<Real> const& rho_salt = mf_rhoS->array(mfi);
 
        funwave_advance_c(rho_salt.dataPtr(), ims, ime, jms, jme, kms, kme);
    }
#endif
 
    //***************************************************
    //Advance one step of the 3d integrator
    //***************************************************
    advance_3d_ml(lev, dt_lev);
 
}