REMORA_advance_2d_onestep.cpp

Go to the documentation of this file.

#include <REMORA.H>
 
using namespace amrex;
 
/**
 * @param[in] lev            level of refinement
 * @param[in] dt_lev         baroclinic time step at level
 * @param[in] dtfast_lev     barotropic time step at level
 * @param[in] my_iif         current barotropic time step
 * @param[in] nfast_counter  how many bartropicc time steps to do per baroclinic step
 */
void REMORA::advance_2d_onestep (int lev, Real /*dt_lev*/, Real dtfast_lev, int my_iif, int nfast_counter)
{
    bool first_2d_step=(my_iif==0);
 
    //Predictor
    bool predictor_2d_step=true;
    int next_indx1 = 0;
    advance_2d(lev,
               vec_rhoS[lev].get() , vec_rhoA[lev].get(),
               vec_ru2d[lev].get()   , vec_rv2d[lev].get(),
               vec_rufrc[lev].get(), vec_rvfrc[lev].get(),
               vec_Zt_avg1[lev].get(),
               vec_DU_avg1[lev], vec_DU_avg2[lev],
               vec_DV_avg1[lev], vec_DV_avg2[lev],
               vec_rubar[lev], vec_rvbar[lev], vec_rzeta[lev],
                vec_ubar[lev],  vec_vbar[lev],  vec_zeta[lev].get(),
               vec_hOfTheConfusingName[lev].get(),
               vec_pm[lev].get(), vec_pn[lev].get(),
               vec_fcor[lev].get(),
               vec_visc2_p[lev].get(), vec_visc2_r[lev].get(),
               vec_mskr[lev].get(), vec_msku[lev].get(), vec_mskv[lev].get(),
               vec_mskp[lev].get(),
               dtfast_lev, predictor_2d_step, first_2d_step, my_iif, next_indx1);
 
    //Corrector. Skip it on last fast step
    predictor_2d_step=false;
    if (my_iif < nfast_counter - 1) {
 
        // These are needed to pass ctests
 
        advance_2d(lev,
                   vec_rhoS[lev].get(), vec_rhoA[lev].get(),
                   vec_ru2d[lev].get(), vec_rv2d[lev].get(),
                   vec_rufrc[lev].get(), vec_rvfrc[lev].get(),
                   vec_Zt_avg1[lev].get(),
                   vec_DU_avg1[lev], vec_DU_avg2[lev],
                   vec_DV_avg1[lev], vec_DV_avg2[lev],
                   vec_rubar[lev], vec_rvbar[lev], vec_rzeta[lev],
                    vec_ubar[lev],  vec_vbar[lev],  vec_zeta[lev].get(),
                   vec_hOfTheConfusingName[lev].get(),
                   vec_pm[lev].get(), vec_pn[lev].get(),
                    vec_fcor[lev].get(),
                   vec_visc2_p[lev].get(), vec_visc2_r[lev].get(),
                   vec_mskr[lev].get(), vec_msku[lev].get(), vec_mskv[lev].get(),
                   vec_mskp[lev].get(),
                   dtfast_lev, predictor_2d_step, first_2d_step, my_iif, next_indx1);
    }
}