multiply by rho to make it correct for AUX

zhichen3 · Jun 5, 2024 · 3b29ec9 · 3b29ec9
1 parent e21b85e
commit 3b29ec9
Showing 1 changed file with 1 addition and 10 deletions.
diff --git a/Exec/science/subchandra/problem_initialize_state_data.H b/Exec/science/subchandra/problem_initialize_state_data.H
@@ -59,19 +59,10 @@ void problem_initialize_state_data (int i, int j, int k,
     set_aux_comp_from_X(burn_state);
 
     for (int n = 0; n < NumAux; n++) {
-        state(i,j,k,UFX+n) = burn_state.aux[n];
+        state(i,j,k,UFX+n) = state(i,j,k,URHO) * burn_state.aux[n];
     }
 #endif
 
-
-    // burn_t burn_state;
-
-    // burn_state.rho = state(i,j,k,URHO);
-    // burn_state.T = state(i,j,k,UTEMP);
-    // for (int n = 0; n < NumSpec; n++) {
-    //     burn_state.xn[n] = state(i,j,k,UFS+n);
-    // }
-
     eos(eos_input_rt, burn_state);
 
     state(i,j,k,UEINT) = state(i,j,k,URHO) * burn_state.e;