blk1m: dot_th depending on rr_to_rv, not full dot_rv

igfuw · Jan 7, 2025 · 63579bb · 63579bb
1 parent e3e0dd1
commit 63579bb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/include/libcloudph++/blk_1m/rhs_cellwise.hpp b/include/libcloudph++/blk_1m/rhs_cellwise.hpp
@@ -123,7 +123,7 @@ namespace libcloudphxx
 
         dot_rv += rr_to_rv;
         dot_rr -= rr_to_rv;
-        dot_th -= const_cp::l_v(T) / (moist_air::c_pd<real_t>() * theta_std::exner(p)) * dot_rv / si::kelvins;
+        dot_th -= const_cp::l_v(T) / (moist_air::c_pd<real_t>() * theta_std::exner(p)) * rr_to_rv / si::kelvins;
       }
     }