ERF_UpdateRhoQtSources_SineMassFlux.H

Go to the documentation of this file.

    ParmParse pp_prob("prob");
 
    Real advection_moisture_rate    =   zero; pp_prob.query("advection_moisture_rate", advection_moisture_rate);
    Real moisture_cutoff            = amrex::Real(500.0); pp_prob.query("moisture_source_cutoff", moisture_cutoff);
    Real moisture_cutoff_transition =  amrex::Real(50.0); pp_prob.query("moisture_source_cutoff_transition", moisture_cutoff_transition);
 
    for (MFIter mfi(*qsrc, TilingIfNotGPU()); mfi.isValid(); ++mfi)
    {
        const auto &box = mfi.tilebox();
        const Array4<Real>& qsrc_arr = qsrc->array(mfi);
        if (box.length(0) != 1)
        {
            qsrc->setVal(0.0);
        } else {
            ParallelFor(box, [=] AMREX_GPU_DEVICE (int i, int j, int k) {
                const Real z_cc = d_zlevels_arr[k];
                if (z_cc < moisture_cutoff) {
                    qsrc_arr(i, j, k) = advection_moisture_rate;
                } else if (z_cc < moisture_cutoff+moisture_cutoff_transition) {
                    Real slope = -advection_moisture_rate / moisture_cutoff_transition;
                    qsrc_arr(i, j, k) = (z_cc-moisture_cutoff) * slope + advection_moisture_rate;
                } else {
                    qsrc_arr(i, j, k) = zero;
                }
            });
        }
    }