From cba159e02fb6f301e7cb8327f394bb667d49f9c5 Mon Sep 17 00:00:00 2001
From: nickishch <nickishch@gmail.com>
Date: Thu, 3 Apr 2025 21:55:52 +0200
Subject: [PATCH] now OWAR ignores dedner cleaning, to avoid them counteracting
each other
---
src/mhd/DirectInduction/mhd_iact.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/mhd/DirectInduction/mhd_iact.h b/src/mhd/DirectInduction/mhd_iact.h
index 5f9f65cebe..3af83c99a3 100644
--- a/src/mhd/DirectInduction/mhd_iact.h
+++ b/src/mhd/DirectInduction/mhd_iact.h
@@ -658,8 +658,8 @@ __attribute__((always_inline)) INLINE static void runner_iact_mhd_force(
for (int i = 0; i < 3; i++) {
pi->mhd_data.Adv_B_source[i] += mj * dB_dt_pref_i * dB_dt_i[i];
pj->mhd_data.Adv_B_source[i] += mi * dB_dt_pref_j * dB_dt_j[i];
- pi->mhd_data.Adv_B_source[i] -= mj * grad_psi_i * dx[i];
- pj->mhd_data.Adv_B_source[i] -= mi * grad_psi_j * dx[i];
+ //pi->mhd_data.Adv_B_source[i] -= mj * grad_psi_i * dx[i];
+ //pj->mhd_data.Adv_B_source[i] -= mi * grad_psi_j * dx[i];
pi->mhd_data.Diff_B_source[i] += mj * dB_dt_pref_PR_i * wi_dr * dB[i];
pj->mhd_data.Diff_B_source[i] -= mi * dB_dt_pref_PR_j * wj_dr * dB[i];
pi->mhd_data.Diff_B_source[i] += mj * art_diff_pref_i * dB[i];
@@ -913,7 +913,7 @@ __attribute__((always_inline)) INLINE static void runner_iact_nonsym_mhd_force(
for (int i = 0; i < 3; i++) {
pi->mhd_data.Adv_B_source[i] += mj * dB_dt_pref_i * dB_dt_i[i];
- pi->mhd_data.Adv_B_source[i] -= mj * grad_psi_i * dx[i];
+ //pi->mhd_data.Adv_B_source[i] -= mj * grad_psi_i * dx[i];
pi->mhd_data.Diff_B_source[i] += mj * dB_dt_pref_PR * wi_dr * dB[i];
pi->mhd_data.Diff_B_source[i] += mj * art_diff_pref * dB[i];
pi->mhd_data.Delta_B[i] += mj * dB_dt_pref_Lap * wi_dr * dB[i];
--
GitLab