320 const int i = team.league_rank();
322 const Real inv_qc_relvar_max = 10;
323 const Real inv_qc_relvar_min = 0.001;
325 const int nlev_packs = ekat::npack<Spack>(
nlev);
326 Kokkos::parallel_for(Kokkos::TeamVectorRange(team, nlev_packs), [&] (
const Int& k)
332 qv(i,k) =
qw(i,k) -
qc(i,k);
339 const auto condition = (
qc(i,k) != 0 &&
qc2(i,k) != 0);
340 if (condition.any()) {
342 ekat::min(inv_qc_relvar_max,
343 ekat::max(inv_qc_relvar_min,
344 ekat::square(
qc(i,k))/
qc2(i,k))));
351 T_mid(i,k) = PF::calculate_temperature_from_dse(dse_ik,z_mid_ik,phis_i);
int Int
Definition: ERF_ShocInterface.H:17
amrex::Real Real
Definition: ERF_ShocInterface.H:16
typename SHF::Spack Spack
Definition: ERF_ShocInterface.H:40
view_1d_const surf_sens_flux
Definition: ERF_ShocInterface.H:379
view_1d ice_flux
Definition: ERF_ShocInterface.H:382
view_1d vapor_flux
Definition: ERF_ShocInterface.H:380
view_1d_const phis
Definition: ERF_ShocInterface.H:376
int nlev
Definition: ERF_ShocInterface.H:366
view_1d heat_flux
Definition: ERF_ShocInterface.H:383
view_2d cldfrac_liq
Definition: ERF_ShocInterface.H:372
view_2d_const qw
Definition: ERF_ShocInterface.H:369
bool compute_mass_and_energy_fluxes
Definition: ERF_ShocInterface.H:377
view_1d water_flux
Definition: ERF_ShocInterface.H:381
view_1d_const surf_evap
Definition: ERF_ShocInterface.H:378
view_2d_const qc_copy
Definition: ERF_ShocInterface.H:369
view_2d_const z_mid
Definition: ERF_ShocInterface.H:375
view_2d_const dse
Definition: ERF_ShocInterface.H:375
view_2d tke
Definition: ERF_ShocInterface.H:368
view_2d T_mid
Definition: ERF_ShocInterface.H:374
view_2d inv_qc_relvar
Definition: ERF_ShocInterface.H:373
view_2d qc
Definition: ERF_ShocInterface.H:368
view_2d qv
Definition: ERF_ShocInterface.H:368
view_2d_const qc2
Definition: ERF_ShocInterface.H:371
view_2d_const tke_copy
Definition: ERF_ShocInterface.H:369