diff --git a/src/propagation/jetcoeffs.jl b/src/propagation/jetcoeffs.jl
index 11944f69..7fc1c275 100644
--- a/src/propagation/jetcoeffs.jl
+++ b/src/propagation/jetcoeffs.jl
@@ -138,7 +138,7 @@ function TaylorIntegration._allocate_jetcoeffs!(::Val{RNp1BP_pN_A_J23E_J2S_ng_ep
     accY = Taylor1(identity(constant_term(zero_q_1)), order)
     accZ = Taylor1(identity(constant_term(zero_q_1)), order)
     local M_ = Array{S}(undef, 3, 3, N)
-    local M_[:, :, ea] = t2c_jpl_de430(dsj2k)
+    local M_[:, :, ea] = t2c_jpl_de430(dsj2k) .+ zero_q_1
     dq[1] = Taylor1(identity(constant_term(q[4])), order)
     dq[2] = Taylor1(identity(constant_term(q[5])), order)
     dq[3] = Taylor1(identity(constant_term(q[6])), order)
@@ -810,7 +810,7 @@ function TaylorIntegration.jetcoeffs!(::Val{RNp1BP_pN_A_J23E_J2S_ng_eph_threads!
     accZ.coeffs[1] = identity(constant_term(zero_q_1))
     accZ.coeffs[2:order + 1] .= zero(accZ.coeffs[1])
     local M_ = Array{S}(undef, 3, 3, N)
-    local M_[:, :, ea] = t2c_jpl_de430(dsj2k)
+    local M_[:, :, ea] = t2c_jpl_de430(dsj2k) .+ zero_q_1
     (dq[1]).coeffs[1] = identity(constant_term(q[4]))
     (dq[1]).coeffs[2:order + 1] .= zero((dq[1]).coeffs[1])
     (dq[2]).coeffs[1] = identity(constant_term(q[5]))
@@ -1749,7 +1749,7 @@ function TaylorIntegration._allocate_jetcoeffs!(::Val{RNp1BP_pN_A_J23E_J2S_eph_t
     accY = Taylor1(identity(constant_term(zero_q_1)), order)
     accZ = Taylor1(identity(constant_term(zero_q_1)), order)
     local M_ = Array{S}(undef, 3, 3, N)
-    local M_[:, :, ea] = t2c_jpl_de430(dsj2k)
+    local M_[:, :, ea] = t2c_jpl_de430(dsj2k) .+ zero_q_1
     dq[1] = Taylor1(identity(constant_term(q[4])), order)
     dq[2] = Taylor1(identity(constant_term(q[5])), order)
     dq[3] = Taylor1(identity(constant_term(q[6])), order)
@@ -2324,7 +2324,7 @@ function TaylorIntegration.jetcoeffs!(::Val{RNp1BP_pN_A_J23E_J2S_eph_threads!},
     accZ.coeffs[1] = identity(constant_term(zero_q_1))
     accZ.coeffs[2:order + 1] .= zero(accZ.coeffs[1])
     local M_ = Array{S}(undef, 3, 3, N)
-    local M_[:, :, ea] = t2c_jpl_de430(dsj2k)
+    local M_[:, :, ea] = t2c_jpl_de430(dsj2k) .+ zero_q_1
     (dq[1]).coeffs[1] = identity(constant_term(q[4]))
     (dq[1]).coeffs[2:order + 1] .= zero((dq[1]).coeffs[1])
     (dq[2]).coeffs[1] = identity(constant_term(q[5]))