From 15ddf4c158a140604537962171a86cde9c072e2d Mon Sep 17 00:00:00 2001 From: Jay Kalinani Date: Wed, 4 Sep 2024 20:22:30 -0400 Subject: [PATCH] AsterSeeds: fix zlocal --- AsterSeeds/src/magTOV.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/AsterSeeds/src/magTOV.cxx b/AsterSeeds/src/magTOV.cxx index f2d0be89..d67d7a1f 100644 --- a/AsterSeeds/src/magTOV.cxx +++ b/AsterSeeds/src/magTOV.cxx @@ -25,7 +25,6 @@ extern "C" void AsterSeeds_InitializeCenteredAvec_TOV(CCTK_ARGUMENTS) { [=] CCTK_HOST(const Loop::PointDesc &p) CCTK_ATTRIBUTE_ALWAYS_INLINE { CCTK_REAL x_local = p.x - dipole_x[0]; CCTK_REAL y_local = p.y - dipole_y[0]; - CCTK_REAL z_local = p.z - dipole_z[0]; CCTK_REAL Pcut = press_max * press_cut; CCTK_REAL Pdiff = std::max(press(p.I) - Pcut, 0.0); CCTK_REAL Aphi_local = Ab * pow(Pdiff, Avec_kappa); @@ -42,6 +41,7 @@ extern "C" void AsterSeeds_InitializeCenteredAvec_TOV(CCTK_ARGUMENTS) { [=] CCTK_HOST(const Loop::PointDesc &p) CCTK_ATTRIBUTE_ALWAYS_INLINE { CCTK_REAL x_local = p.x - dipole_x[0]; CCTK_REAL y_local = p.y - dipole_y[0]; + CCTK_REAL z_local = p.z - dipole_z[0]; CCTK_REAL cylrad2 = x_local * x_local + y_local * y_local; CCTK_REAL rsph = sqrt(x_local * x_local + y_local * y_local + z_local * z_local);