diff --git a/README.md b/README.md
index 1f85d46..eec5d58 100644
--- 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
index c52f475..8017442 100644
--- 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
index a2ee4d0..46a01e9 100644
--- 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
index 7f77faa..dab6027 100644
--- 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)