Use lwp instead of q_liq for the cldn function, use log - log

Author: szy21

src/callbacks/callbacks.jl

@@ -302,6 +302,14 @@ NVTX.@annotate function rrtmgp_model_callback!(integrator)
                 end
             end
 
+            lwp_col = p.scratch.temp_field_level
+            ᶜliquid_water_mass_concentration =
+                @. lazy(cloud_liquid_water_content * Y.c.ρ)
+            Operators.column_integral_definite!(
+                lwp_col,
+                ᶜliquid_water_mass_concentration,
+            )
+
             @. ᶜreliq = ifelse(
                 cloud_liquid_water_content > FT(0),
                 CM.CloudDiagnostics.effective_radius_Liu_Hallet_97(
@@ -313,8 +321,7 @@ NVTX.@annotate function rrtmgp_model_callback!(integrator)
                         dust_aero_conc,
                         seasalt_aero_conc,
                         SO4_aero_conc,
-                        cloud_liquid_water_content /
-                        max(eps(FT), cloud_fraction),
+                        lwp_col,
                     ),
                     FT(0),
                     FT(0),

src/parameterized_tendencies/microphysics/microphysics_wrappers.jl

@@ -54,10 +54,10 @@ function ml_N_cloud_liquid_droplets(cmc, c_dust, c_seasalt, c_SO4, q_liq)
     FT = eltype(N₀)
     return N₀ * (
         FT(1) +
-        α_dust * log(max(c_dust, eps(FT)) / c₀_dust) +
-        α_seasalt * log(max(c_seasalt, eps(FT)) / c₀_seasalt) +
-        α_SO4 * log(max(c_SO4, eps(FT)) / c₀_SO4) +
-        α_q_liq * log(max(q_liq, eps(FT)) / q₀_liq)
+        α_dust * (log(max(c_dust, eps(FT))) - log(c₀_dust)) +
+        α_seasalt * (log(max(c_seasalt, eps(FT))) - log(c₀_seasalt)) +
+        α_SO4 * (log(max(c_SO4, eps(FT))) - log(c₀_SO4)) +
+        α_q_liq * (log(max(q_liq, eps(FT))) - log(q₀_liq))
     )
 end