diff --git a/src/moments_eq_rhs.F90 b/src/moments_eq_rhs.F90
index 1914e3852ba05e2b4214554bafdfa4f6894ca067..6d54d19486719c32c41fb16a37b339a60653233a 100644
--- a/src/moments_eq_rhs.F90
+++ b/src/moments_eq_rhs.F90
@@ -123,7 +123,7 @@ SUBROUTINE moments_eq_rhs_e
           !! Adding non linearity
           IF ( NON_LIN ) THEN
             moments_rhs_e(ip,ij,ikx,iky,iz,updatetlevel) = &
-              moments_rhs_e(ip,ij,ikx,iky,iz,updatetlevel) + Sepj(ip,ij,ikx,iky,iz)
+              moments_rhs_e(ip,ij,ikx,iky,iz,updatetlevel) - Sepj(ip,ij,ikx,iky,iz)
           ENDIF
 
           END DO kyloope
@@ -262,7 +262,7 @@ SUBROUTINE moments_eq_rhs_i
           !! Adding non linearity
           IF ( NON_LIN ) THEN
            moments_rhs_i(ip,ij,ikx,iky,iz,updatetlevel) = &
-             moments_rhs_i(ip,ij,ikx,iky,iz,updatetlevel) + Sipj(ip,ij,ikx,iky,iz)
+             moments_rhs_i(ip,ij,ikx,iky,iz,updatetlevel) - Sipj(ip,ij,ikx,iky,iz)
           ENDIF
 
           END DO kyloopi