diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90 index 6f478b7c6e1bdd1f1f558d707b6498efe17cec94..72b71046f30035b4117cec67bd8d33f35731b812 100644 --- a/src/geometry_mod.F90 +++ b/src/geometry_mod.F90 @@ -258,7 +258,7 @@ CONTAINS hatB_NL(iz,eo) = 1._dp ! Derivative of the magnetic field strenght - gradxB(iz,eo) = 0._dp ! Gene put a factor hatB^2 or 1/hatR^2 in this + gradxB(iz,eo) = -1._dp ! Gene put a factor hatB^2 or 1/hatR^2 in this gradyB(iz,eo) = 0._dp gradzB(iz,eo) = 0._dp ! Gene put a factor hatB or 1/hatR in this @@ -267,7 +267,7 @@ CONTAINS ky = kyarray(iky) DO ikx= ikxs, ikxe kx = kxarray(ikx) - Ckxky(iky, ikx, iz,eo) = -ky * hatB(iz,eo) ! .. multiply by hatB to cancel the 1/ hatB factor in moments_eqs_rhs.f90 routine + Ckxky(iky, ikx, iz,eo) = -ky ENDDO ENDDO ! coefficient in the front of parallel derivative