From 8d42c92af8eecb1f281ad5b112f93e884eb2d41d Mon Sep 17 00:00:00 2001 From: "Aaron M. Lattanzi" <103702284+AMLattanzi@users.noreply.github.com> Date: Fri, 2 Feb 2024 12:10:07 -0700 Subject: [PATCH] Fix wrf init for moisture vars. (#1421) --- Source/Initialization/ERF_init_from_wrfinput.cpp | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/Source/Initialization/ERF_init_from_wrfinput.cpp b/Source/Initialization/ERF_init_from_wrfinput.cpp index dd4eb777d..780d4633f 100644 --- a/Source/Initialization/ERF_init_from_wrfinput.cpp +++ b/Source/Initialization/ERF_init_from_wrfinput.cpp @@ -294,11 +294,13 @@ init_state_from_wrfinput (int lev, if (moisture_type != MoistureType::None) { state_fab.template copy(NC_QVAPOR_fab[idx], 0, RhoQ1_comp, 1); - state_fab.template plus(NC_QCLOUD_fab[idx], 0, RhoQ1_comp, 1); state_fab.template mult(NC_rho_fab[idx] , 0, RhoQ1_comp, 1); - state_fab.template copy(NC_QRAIN_fab[idx], 0, RhoQ2_comp, 1); - state_fab.template mult(NC_rho_fab[idx] , 0, RhoQ2_comp, 1); + state_fab.template copy(NC_QCLOUD_fab[idx], 0, RhoQ2_comp, 1); + state_fab.template mult(NC_rho_fab[idx] , 0, RhoQ2_comp, 1); + + state_fab.template copy(NC_QRAIN_fab[idx], 0, RhoQ3_comp, 1); + state_fab.template mult(NC_rho_fab[idx] , 0, RhoQ3_comp, 1); } } // idx }