diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90
index 9774219befe420dcb28745f7a809346ca48f8491..4c90e1e520794f3e5278eeaab0703fedace44824 100644
--- a/src/geometry_mod.F90
+++ b/src/geometry_mod.F90
@@ -167,7 +167,7 @@ CONTAINS
 
         ! Nonlinear term prefactor
         ! (according to my derivations, there should be a metric dependent factor in front of the Poisson bracket)
-        hatB_NL(iz,eo) = 1._dp !Jacobian(iz,eo)*(gxx(iz,eo)*gyy(iz,eo) - gxy(iz,eo)**2)/hatB(iz,eo)
+        hatB_NL(iz,eo) = (gxx(iz,eo)*gyy(iz,eo) - gxy(iz,eo)**2)
 
       ENDDO
     ENDDO