ERF_UpdateRhoQtSources_RICO.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);
    // RICO (van Zanten et al. 2011): linear in z up to 2980 m and constant above. Above the
    // cutoff the source keeps its value at the cutoff, so it is continuous there
    Real moisture_cutoff            = amrex::Real(2980.0); pp_prob.query("moisture_source_cutoff", moisture_cutoff);
 
    // Only apply moisture source below nominal inversion height
    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)
        {
            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) =  amrex::Real(5.226199353716133e-12)*z_cc + advection_moisture_rate;
                } else {
                    qsrc_arr(i, j, k) = amrex::Real(5.226199353716133e-12)*moisture_cutoff + advection_moisture_rate;
                }
            });
        }
    }