Skip to content

Commit

Permalink
Fix minor typos found with debug compile. (#1610)
Browse files Browse the repository at this point in the history
  • Loading branch information
AMLattanzi authored May 9, 2024
1 parent 8915c90 commit a893072
Show file tree
Hide file tree
Showing 3 changed files with 15 additions and 15 deletions.
12 changes: 6 additions & 6 deletions Source/BoundaryConditions/ABLMost.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -131,24 +131,24 @@ ABLMost::impose_most_bcs (const int& lev,
if (flux_type == FluxCalcType::MOENG) {
moeng_flux flux_comp(klo);
compute_most_bcs(lev, mfs,
xzmom_flux, xzmom_flux,
yzmom_flux, yzmom_flux,
xzmom_flux, zxmom_flux,
yzmom_flux, zymom_flux,
heat_flux,
eddyDiffs,
z_phys, m_geom[lev].CellSize(2), flux_comp);
} else if (flux_type == FluxCalcType::DONELAN) {
donelan_flux flux_comp(klo);
compute_most_bcs(lev, mfs,
xzmom_flux, xzmom_flux,
yzmom_flux, yzmom_flux,
xzmom_flux, zxmom_flux,
yzmom_flux, zymom_flux,
heat_flux,
eddyDiffs,
z_phys, m_geom[lev].CellSize(2), flux_comp);
} else {
custom_flux flux_comp(klo);
compute_most_bcs(lev, mfs,
xzmom_flux, xzmom_flux,
yzmom_flux, yzmom_flux,
xzmom_flux, zxmom_flux,
yzmom_flux, zymom_flux,
heat_flux,
eddyDiffs,
z_phys, m_geom[lev].CellSize(2), flux_comp);
Expand Down
6 changes: 3 additions & 3 deletions Source/BoundaryConditions/MOSTStress.H
Original file line number Diff line number Diff line change
Expand Up @@ -600,7 +600,7 @@ struct moeng_flux
const int& n,
const int& icomp,
const amrex::Real& dz,
const amrex::Real& dz1,
const amrex::Real& /*dz1*/,
const bool& exp_most,
const amrex::Array4<const amrex::Real>& eta_arr,
const amrex::Array4<const amrex::Real>& cons_arr,
Expand Down Expand Up @@ -890,7 +890,7 @@ struct donelan_flux
const int& n,
const int& icomp,
const amrex::Real& dz,
const amrex::Real& dz1,
const amrex::Real& /*dz1*/,
const bool& exp_most,
const amrex::Array4<const amrex::Real>& eta_arr,
const amrex::Array4<const amrex::Real>& cons_arr,
Expand Down Expand Up @@ -1161,7 +1161,7 @@ struct custom_flux
const int& n,
const int& icomp,
const amrex::Real& dz,
const amrex::Real& dz1,
const amrex::Real& /*dz1*/,
const bool& exp_most,
const amrex::Array4<const amrex::Real>& eta_arr,
const amrex::Array4<const amrex::Real>& cons_arr,
Expand Down
12 changes: 6 additions & 6 deletions Source/Microphysics/SAM/Init_SAM.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -218,8 +218,8 @@ void SAM::Compute_Coefficients ()
{
Real Prefactor;
Real pratio = sqrt(1.29 / rho1d_t(k));
Real rrr1 = 393.0/(tabs1d_t(k)+120.0)*std::pow((tabs1d_t(k)/273.0),1.5);
Real rrr2 = std::pow((tabs1d_t(k)/273.0),1.94)*(1000.0/pres1d_t(k));
//Real rrr1 = 393.0/(tabs1d_t(k)+120.0)*std::pow((tabs1d_t(k)/273.0),1.5);
//Real rrr2 = std::pow((tabs1d_t(k)/273.0),1.94)*(1000.0/pres1d_t(k));
Real estw = 100.0*erf_esatw(tabs1d_t(k));
Real esti = 100.0*erf_esati(tabs1d_t(k));

Expand All @@ -235,8 +235,8 @@ void SAM::Compute_Coefficients ()
coef2 = R_v * R_d / (diffelq * esti);
Prefactor = 2.0 * PI * nzeros / (rho1d_t(k) * (coef1 + coef2));
Prefactor *= (2.0/PI); // Shape factor snow
evapr1_t(k) = Prefactor * 0.65 * sqrt(rho1d_t(k) / (PI * rhos * nzeros));
evapr2_t(k) = Prefactor * 0.44 * sqrt(a_snow * rho1d_t(k) / muelq) * gams2
evaps1_t(k) = Prefactor * 0.65 * sqrt(rho1d_t(k) / (PI * rhos * nzeros));
evaps2_t(k) = Prefactor * 0.44 * sqrt(a_snow * rho1d_t(k) / muelq) * gams2
* sqrt(pratio) * pow(rho1d_t(k) / (PI * rhos * nzeros) , ((5.0+b_snow)/8.0));

// accretion by graupel:
Expand All @@ -249,8 +249,8 @@ void SAM::Compute_Coefficients ()
coef1 = (lsub/(tabs1d_t(k)*R_v)-1.0)*lsub/(therco*tabs1d_t(k));
coef2 = R_v * R_d / (diffelq * esti);
Prefactor = 2.0 * PI * nzerog / (rho1d_t(k) * (coef1 + coef2)); // Shape factor for graupel is 1
evapr1_t(k) = Prefactor * 0.78 * sqrt(rho1d_t(k) / (PI * rhog * nzerog));
evapr2_t(k) = Prefactor * 0.31 * sqrt(a_grau * rho1d_t(k) / muelq) * gamg2
evapg1_t(k) = Prefactor * 0.78 * sqrt(rho1d_t(k) / (PI * rhog * nzerog));
evapg2_t(k) = Prefactor * 0.31 * sqrt(a_grau * rho1d_t(k) / muelq) * gamg2
* sqrt(pratio) * pow(rho1d_t(k) / (PI * rhog * nzerog) , ((5.0+b_grau)/8.0));

// accretion by rain:
Expand Down

0 comments on commit a893072

Please sign in to comment.