debug flag fixed

PrincetonUniversity · Jan 5, 2024 · e999640 · e999640
1 parent 67f7012
commit e999640
Show file tree

Hide file tree

Showing 4 changed files with 10 additions and 14 deletions.
diff --git a/README.md b/README.md
@@ -1,4 +1,4 @@
-# Tristan v2.6
+# Tristan v2.7
 
 [![DOI](https://zenodo.org/badge/234551890.svg)](https://zenodo.org/badge/latestdoi/234551890)
 
@@ -120,6 +120,8 @@ or in the VSCode environment (see the extension list in the `.vscode/settings.js
 __@TODO__
 
 ## Latest Releases
+* `v2.6.1` __Aug 2023__
+  * Fixed a buggy ordering of synchrotron cooling term and particle coordinate update (very minor correction)
 * `v2.6` __Jan 2023__
   * Minor cleanup + bugfixes
 * `v2.6r1` __Nov 2022__

diff --git a/configure.py b/configure.py
@@ -267,7 +267,7 @@
     makefile_options[
         "WARNING_FLAGS"
     ] += "-Wall -Wextra -Wconversion -pedantic -Wno-compare-reals -Wno-unused-dummy-argument "
-    if int(args["d0ebug"]) >= 0:
+    if int(args["debug"]) >= 0:
         makefile_options["PREPROCESSOR_FLAGS"] += "-DDEBUG "
     if int(args["debug"]) >= 1:
         if args["flags"] == "intel":

diff --git a/src/algorithms/fieldsolvers/absorb_ampere.F08 b/src/algorithms/fieldsolvers/absorb_ampere.F08
@@ -37,22 +37,19 @@
         lam1 = (1.0 + lam) / (1.0 - lam)
         lam2 = 1.0 / (1.0 - lam)
         ex(i, j, k) = lam1 * ex(i, j, k) + lam2 * const * &
-                      (-bz(i, jm1, k) + bz(i, j, k)) - &
-                      2.0 * lam * lam2 * ex_target
+                      (-bz(i, jm1, k) + bz(i, j, k))
 
         lam = 0.5 * lambdaAbsorb(xg, yg + 0.5, zg)
         lam1 = (1.0 + lam) / (1.0 - lam)
         lam2 = 1.0 / (1.0 - lam)
         ey(i, j, k) = lam1 * ey(i, j, k) + lam2 * const * &
-                      (bz(im1, j, k) - bz(i, j, k)) - &
-                      2.0 * lam * lam2 * ey_target
+                      (bz(im1, j, k) - bz(i, j, k))
 
         lam = 0.5 * lambdaAbsorb(xg, yg, zg)
         lam1 = (1.0 + lam) / (1.0 - lam)
         lam2 = 1.0 / (1.0 - lam)
         ez(i, j, k) = lam1 * ez(i, j, k) + lam2 * const * &
-                      (bx(i, jm1, k) - bx(i, j, k) - by(im1, j, k) + by(i, j, k)) - &
-                      2.0 * lam * lam2 * ez_target
+                      (bx(i, jm1, k) - bx(i, j, k) - by(im1, j, k) + by(i, j, k))
       end do
     end do
 #elif defined(threeD)

diff --git a/src/algorithms/fieldsolvers/absorb_faraday.F08 b/src/algorithms/fieldsolvers/absorb_faraday.F08
@@ -37,22 +37,19 @@
         lam1 = (1.0 + lam) / (1.0 - lam)
         lam2 = 1.0 / (1.0 - lam)
         bx(i, j, k) = lam1 * bx(i, j, k) + lam2 * const * &
-                      (-ez(i, jp1, k) + ez(i, j, k)) - &
-                      2.0 * lam * lam2 * bx_target
+                      (-ez(i, jp1, k) + ez(i, j, k))
 
         lam = 0.25 * lambdaAbsorb(xg + 0.5, yg, zg)
         lam1 = (1.0 + lam) / (1.0 - lam)
         lam2 = 1.0 / (1.0 - lam)
         by(i, j, k) = lam1 * by(i, j, k) + lam2 * const * &
-                      (ez(ip1, j, k) - ez(i, j, k)) - &
-                      2.0 * lam * lam2 * by_target
+                      (ez(ip1, j, k) - ez(i, j, k))
 
         lam = 0.25 * lambdaAbsorb(xg + 0.5, yg + 0.5, zg)
         lam1 = (1.0 + lam) / (1.0 - lam)
         lam2 = 1.0 / (1.0 - lam)
         bz(i, j, k) = lam1 * bz(i, j, k) + lam2 * const * &
-                      (ex(i, jp1, k) - ex(i, j, k) - ey(ip1, j, k) + ey(i, j, k)) - &
-                      2.0 * lam * lam2 * bz_target
+                      (ex(i, jp1, k) - ex(i, j, k) - ey(ip1, j, k) + ey(i, j, k))
       end do
     end do
 #elif defined(threeD)