minor bugfix

Author: sajjadazimi

src/cache/precipitation_precomputed_quantities.jl

@@ -115,20 +115,13 @@ function set_precipitation_velocities!(
 
     # compute the precipitation terminal velocity [m/s]
     # TODO sedimentation of snow is based on the 1M scheme
-    @. ᶜwNᵣ = getindex(CM2.rain_terminal_velocity(
+    @. (ᶜwNᵣ, ᶜwᵣ) = CM2.rain_terminal_velocity(
         cm2p.sb,
         cm2p.tv,
         q_rai,
         Y.c.ρ,
         Y.c.N_rai,
-    ), 1)
-    @. ᶜwᵣ = getindex(CM2.rain_terminal_velocity(
-        cm2p.sb,
-        cm2p.tv,
-        q_rai,
-        Y.c.ρ,
-        Y.c.N_rai,
-    ), 2)
+    )
     @. ᶜwNₛ = CM1.terminal_velocity(
         cm1p.ps,
         cm1p.tv.snow,

src/parameterized_tendencies/microphysics/cloud_condensate.jl

@@ -55,6 +55,6 @@ function cloud_condensate_tendency!(
     @. Yₜ.c.ρq_liq += Y.c.ρ * cloud_sources(cmc.liquid, thp, ᶜts, q_rai, dt)
     @. Yₜ.c.ρq_ice += Y.c.ρ * cloud_sources(cmc.ice, thp, ᶜts, q_sno, dt)
 
-    @. Yₜ.c.N_liq += aerosol_activation_sources(cmc.liquid, thp, ᶜts, Y.c.ρ, q_rai, Y.c.N_ice, cmc.N_cloud_liquid_droplets, dt)
+    @. Yₜ.c.N_liq += aerosol_activation_sources(cmc.liquid, thp, ᶜts, Y.c.ρ, q_rai, Y.c.N_liq, cmc.N_cloud_liquid_droplets, dt)
     @. Yₜ.c.N_ice += aerosol_activation_sources(cmc.ice, thp, ᶜts, Y.c.ρ, q_sno, Y.c.N_ice, cmc.N_cloud_liquid_droplets, dt)
 end