321 const int i = team.league_rank();
323 const Real inv_qc_relvar_max = 10;
324 const Real inv_qc_relvar_min = 0.001;
326 const int nlev_packs = ekat::npack<Spack>(
nlev);
327 Kokkos::parallel_for(Kokkos::TeamVectorRange(team, nlev_packs), [&] (
const Int& k)
333 qv(i,k) =
qw(i,k) -
qc(i,k);
340 const auto condition = (
qc(i,k) != 0 &&
qc2(i,k) != 0);
341 if (condition.any()) {
343 ekat::min(inv_qc_relvar_max,
344 ekat::max(inv_qc_relvar_min,
345 ekat::square(
qc(i,k))/
qc2(i,k))));
352 T_mid(i,k) = PF::calculate_temperature_from_dse(dse_ik,z_mid_ik,phis_i);
int Int
Definition: ERF_ShocInterface.H:20
amrex::Real Real
Definition: ERF_ShocInterface.H:19
typename SHF::Spack Spack
Definition: ERF_ShocInterface.H:43
view_1d_const surf_sens_flux
Definition: ERF_ShocInterface.H:380
view_1d ice_flux
Definition: ERF_ShocInterface.H:383
view_1d vapor_flux
Definition: ERF_ShocInterface.H:381
view_1d_const phis
Definition: ERF_ShocInterface.H:377
int nlev
Definition: ERF_ShocInterface.H:367
view_1d heat_flux
Definition: ERF_ShocInterface.H:384
view_2d cldfrac_liq
Definition: ERF_ShocInterface.H:373
view_2d_const qw
Definition: ERF_ShocInterface.H:370
bool compute_mass_and_energy_fluxes
Definition: ERF_ShocInterface.H:378
view_1d water_flux
Definition: ERF_ShocInterface.H:382
view_1d_const surf_evap
Definition: ERF_ShocInterface.H:379
view_2d_const qc_copy
Definition: ERF_ShocInterface.H:370
view_2d_const z_mid
Definition: ERF_ShocInterface.H:376
view_2d_const dse
Definition: ERF_ShocInterface.H:376
view_2d tke
Definition: ERF_ShocInterface.H:369
view_2d T_mid
Definition: ERF_ShocInterface.H:375
view_2d inv_qc_relvar
Definition: ERF_ShocInterface.H:374
view_2d qc
Definition: ERF_ShocInterface.H:369
view_2d qv
Definition: ERF_ShocInterface.H:369
view_2d_const qc2
Definition: ERF_ShocInterface.H:372
view_2d_const tke_copy
Definition: ERF_ShocInterface.H:370