rename molmass_ratio to Rv_over_Rd

Author: szy21

src/initial_conditions/initial_conditions.jl

@@ -686,7 +686,7 @@ function moist_baroclinic_wave_values(z, ϕ, λ, params, perturb, deep_atmospher
     T = T_v / (1 + ε * q_tot) # This is the formula used in the paper.
 
     # This is the actual formula, which would be consistent with TD:
-    # T = T_v * (1 + q_tot) / (1 + q_tot * CAP.molmass_ratio(params))
+    # T = T_v * (1 + q_tot) / (1 + q_tot * CAP.Rv_over_Rd(params))
 
     return (; T, p, q_tot, u, v)
 end
@@ -1221,16 +1221,16 @@ function (initial_condition::TRMM_LBA)(params)
     # in Pressel et al., 2015). Note that the measured profiles are different from the
     # ones required for hydrostatic balance.
     # TODO: Move this to APL.
-    molmass_ratio = TD.Parameters.molmass_ratio(thermo_params)
+    Rv_over_Rd = TD.Parameters.Rv_over_Rd(thermo_params)
     measured_p = APL.TRMM_LBA_p(FT)
     measured_RH = APL.TRMM_LBA_RH(FT)
     measured_z_values = APL.TRMM_LBA_z(FT)
     measured_q_tot_values = map(measured_z_values) do z
         p_v_sat = TD.saturation_vapor_pressure(thermo_params, T(z), TD.Liquid())
         denominator =
             measured_p(z) - p_v_sat +
-            (1 / molmass_ratio) * p_v_sat * measured_RH(z) / 100
-        q_v_sat = p_v_sat * (1 / molmass_ratio) / denominator
+            (1 / Rv_over_Rd) * p_v_sat * measured_RH(z) / 100
+        q_v_sat = p_v_sat * (1 / Rv_over_Rd) / denominator
         return q_v_sat * measured_RH(z) / 100
     end
     q_tot = Intp.extrapolate(

src/parameters/Parameters.jl

@@ -137,7 +137,7 @@ for var in fieldnames(TD.Parameters.ThermodynamicsParameters)
     @eval $var(ps::ACAP) = TD.Parameters.$var(thermodynamics_params(ps))
 end
 # Thermodynamics derived parameters
-for var in [:molmass_ratio, :R_d, :R_v, :e_int_v0, :cp_d, :cv_v, :cv_l, :cv_d]
+for var in [:Rv_over_Rd, :R_d, :R_v, :e_int_v0, :cp_d, :cv_v, :cv_l, :cv_d]
     @eval $var(ps::ACAP) = TD.Parameters.$var(thermodynamics_params(ps))
 end
 

src/prognostic_equations/eddy_diffusion_closures.jl

@@ -120,7 +120,7 @@ function buoyancy_gradients(
     FT = eltype(bg_model)
 
     g = TDP.grav(thermo_params)
-    molmass_ratio = TDP.molmass_ratio(thermo_params)
+    Rv_over_Rd = TDP.Rv_over_Rd(thermo_params)
     R_d = TDP.R_d(thermo_params)
     R_v = TDP.R_v(thermo_params)
 
@@ -159,7 +159,7 @@ function buoyancy_gradients(
         qv_sat = get_qv_sat(thermo_params, bg_model)
         ∂b∂θli_sat = (
             ∂b∂θv *
-            (1 + molmass_ratio * (1 + lh / R_v / t_sat) * qv_sat - qt_sat) /
+            (1 + Rv_over_Rd * (1 + lh / R_v / t_sat) * qv_sat - qt_sat) /
             (1 + lh * lh / cp_m / R_v / t_sat / t_sat * qv_sat)
         )
         ∂b∂qt_sat =