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 }