ERF_ApplySurfaceTreatment_BulkCoeff.cpp File Reference

#include <AMReX_MultiFab.H>
#include <ERF_SrcHeaders.H>
#include <AMReX_ParmParse.H>

Include dependency graph for ERF_ApplySurfaceTreatment_BulkCoeff.cpp:

Functions
void	ApplySurfaceTreatment_BulkCoeff_Mom (const Box &tbx, const Box &tby, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< const Real > &rho_u, const Array4< const Real > &rho_v, const Array4< const Real > &cons_state, const Array4< const Real > &z_phys_nd, const Array4< const Real > &surface_state_arr)
void	ApplySurfaceTreatment_BulkCoeff_CC (const Box &bx, const Array4< Real > &cell_rhs, const Array4< const Real > &cons_state, const Array4< const Real > &z_phys_cc, const Array4< const Real > &surface_state_arr)

Function Documentation

ApplySurfaceTreatment_BulkCoeff_CC()

void ApplySurfaceTreatment_BulkCoeff_CC	(	const Box &	bx,
		const Array4< Real > &	cell_rhs,
		const Array4< const Real > &	cons_state,
		const Array4< const Real > &	z_phys_cc,
		const Array4< const Real > &	surface_state_arr
	)

{
     ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
        if(k == 0) {
            Real dz = z_phys_cc(i,j,1)-z_phys_cc(i,j,0);
            Real ls_mask = surface_state_arr(i,j,0);
            Real Ch = 0.0015*(1.0-ls_mask);
            Real Ce = 0.0015*(1.0-ls_mask);
 
            Real rho = cons_state(i,j,k, Rho_comp);
            Real rhotheta = cons_state(i,j,k, RhoTheta_comp);
            Real rhoqv = cons_state(i,j,k, RhoQ1_comp);
            Real theta = rhotheta/rho;
            Real qv = rhoqv/rho;
            Real temp = getTgivenRandRTh(rho, rhotheta, qv);
            Real uvel = cons_state(i,j,k,1)/cons_state(i,j,k,0);
            Real vvel = cons_state(i,j,k,2)/cons_state(i,j,k,0);
            Real velmag = std::sqrt(uvel*uvel + vvel*vvel);
            Real dT = max(0.0, 301.0 - temp);
            Real dq = max(0.0, 0.024 - qv);
            cell_rhs(i, j, k, RhoTheta_comp) += (theta/(1005.0*temp))*rho*Ch*velmag*dT/dz;
            cell_rhs(i, j, k, RhoQ1_comp) += rho*Ce*velmag*dq/dz;
        }
    });
}

Referenced by make_sources().

Here is the call graph for this function:

Here is the caller graph for this function:

ApplySurfaceTreatment_BulkCoeff_Mom()

void ApplySurfaceTreatment_BulkCoeff_Mom	(	const Box &	tbx,
		const Box &	tby,
		const Array4< Real > &	rho_u_rhs,
		const Array4< Real > &	rho_v_rhs,
		const Array4< const Real > &	rho_u,
		const Array4< const Real > &	rho_v,
		const Array4< const Real > &	cons_state,
		const Array4< const Real > &	z_phys_nd,
		const Array4< const Real > &	surface_state_arr
	)

{
    int ndrag = 1;
    Real Cd_sea = 0.001;
    Real Cd_land = 0.01;
 
    ParallelFor(tbx, [=] AMREX_GPU_DEVICE(int i, int j, int k)
    {
        if(k <= ndrag) {
            Real dz = z_phys_nd(i,j,1) - z_phys_nd(i,j,0);
            Real ls_mask = (surface_state_arr(i-1,j,0)+surface_state_arr(i,j,0))/2.0;
            Real weight = 1.0 - Real(k)/ndrag;  // linear decrease
            Real fac = k==0? 1.0 : 0.0;
            Real Cd = fac*Cd_sea*(1.0-ls_mask) + Cd_land*ls_mask;
            Real rho_for_u = (cons_state(i-1,j,k,0)+cons_state(i,j,k,0))/2.0;
            Real rho_for_v = (cons_state(i,j-1,k,0)+cons_state(i,j,k,0))/2.0;
            Real uvel = rho_u(i,j,k)/rho_for_u;
            Real vvel = rho_v(i,j,k)/rho_for_v;
            Real velmag = std::sqrt(uvel*uvel + vvel*vvel);
            rho_u_rhs(i, j, k) += -1.0*weight*Cd*velmag*rho_for_u*uvel/dz;
        }
    });
 
 
    ParallelFor(tby, [=] AMREX_GPU_DEVICE(int i, int j, int k)
    {
       if(k <= ndrag) {
            Real dz = z_phys_nd(i,j,1) - z_phys_nd(i,j,0);
            Real ls_mask = (surface_state_arr(i,j-1,0)+surface_state_arr(i,j,0))/2.0;
            Real fac = k==0? 1.0 : 0.0;
            Real Cd = fac*Cd_sea*(1.0-ls_mask) + Cd_land*ls_mask;
            Real weight = 1.0 - Real(k)/ndrag;  // linear decrease
            Real rho_for_u = (cons_state(i-1,j,k,0)+cons_state(i,j,k,0))/2.0;
            Real rho_for_v = (cons_state(i,j-1,k,0)+cons_state(i,j,k,0))/2.0;
            Real uvel = rho_u(i,j,k)/rho_for_u;
            Real vvel = rho_v(i,j,k)/rho_for_v;
            Real velmag = std::sqrt(uvel*uvel + vvel*vvel);
            rho_v_rhs(i, j, k) += -1.0*weight*Cd*velmag*rho_for_v*vvel/dz;
        }
    });
}

Referenced by make_mom_sources().

Here is the caller graph for this function: