Merge pull request igfuw#457 from pdziekan/blk1m_dotth_fix

pdziekan · web-flow · commit 414591e85a3c · 2025-01-07T13:55:37.000+01:00
blk1m: dot_th depending on rr_to_rv, not full dot_rv
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;
       }
     }
 

Original file line number	Diff line number	Diff line change
`@@ -123,7 +123,7 @@ namespace libcloudphxx`
`123`	`123`
`124`	`124`	`dot_rv += rr_to_rv;`
`125`	`125`	`dot_rr -= rr_to_rv;`
`126`		`- dot_th -= const_cp::l_v(T) / (moist_air::c_pd<real_t>() * theta_std::exner(p)) * dot_rv / si::kelvins;`
	`126`	`+ dot_th -= const_cp::l_v(T) / (moist_air::c_pd<real_t>() * theta_std::exner(p)) * rr_to_rv / si::kelvins;`
`127`	`127`	`}`
`128`	`128`	`}`
`129`	`129`