try again

AMReX-Astro · Apr 21, 2024 · 2d0d90f · 2d0d90f
1 parent 956f363
commit 2d0d90f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/Exec/science/wdmerger/problem_initialize_state_data.H b/Exec/science/wdmerger/problem_initialize_state_data.H
@@ -151,11 +151,11 @@ void problem_initialize_state_data (int i, int j, int k,
     if (problem::problem != 1) {
 
         if (P_star_test || S_star_test) {
-            if (rho_P > rho_S && problem::mass_P > 0.0_rt) {
+            if (rho_P > rho_S && problem::mass_P > 0.0_rt && dist_P < problem::radius_P) {
                 state(i,j,k,UMX) += problem::vel_P[0] * state(i,j,k,URHO);
                 state(i,j,k,UMY) += problem::vel_P[1] * state(i,j,k,URHO);
                 state(i,j,k,UMZ) += problem::vel_P[2] * state(i,j,k,URHO);
-            } else if (problem::mass_S > 0.0_rt) {
+            } else if (problem::mass_S > 0.0_rt && dist_S < problem::radius_S) {
                 state(i,j,k,UMX) += problem::vel_S[0] * state(i,j,k,URHO);
                 state(i,j,k,UMY) += problem::vel_S[1] * state(i,j,k,URHO);
                 state(i,j,k,UMZ) += problem::vel_S[2] * state(i,j,k,URHO);