Adding Spice Object makes position vector incorrect

[natasha-proteus]

I have some code that prints the position vector of a satellite, based on the scenarioBasicOrbit.py . I need to add a spice object because I am going to be getting eclipse times of the satellite. My code without the spice object is below:

    taskName = "unitTask"               # arbitrary name (don't change)
    processname = "TestProcess"         # arbitrary name (don't change)

    epoch = startTime

    # Create a sim module as an empty container
    scenarioSim = SimulationBaseClass.SimBaseClass()
    
    scenarioSim.SetProgressBar(True)

    # create the dynamics task and specify the integration update time
    dynProcess = scenarioSim.CreateNewProcess(processname)
    simulationTimeStep = macros.sec2nano(10.) # update process rate update time
    dynProcess.addTask(scenarioSim.CreateNewTask(taskName, simulationTimeStep))

    gravFactory = simIncludeGravBody.gravBodyFactory()
    earth = gravFactory.createEarth()
    earth.isCentralBody = True          # ensure this is the central gravitational body
    earth_mu = earth.mu

    # Create a spacecraft
    scObject = spacecraft.Spacecraft()

    scObject.ModelTag = name

    # add spacecraft object to the simulation process
    scenarioSim.AddModelToTask(taskName, scObject)

    #   setup orbit using orbitalMotion library
   oe = orbitalMotion.ClassicElements()
   oe.a = 7000. * 1000 
   oe.e = 0.0001
   oe.i = 33.3 * macros.D2R
   oe.Omega = 48.2 * macros.D2R
   oe.omega = 347.8 * macros.D2R
   oe.f = 85.3 * macros.D2R
   rN, vN = orbitalMotion.elem2rv(mu, oe)
   oe = orbitalMotion.rv2elem(mu, rN, vN) 

   scObject.hub.r_CN_NInit = rN
   scObject.hub.v_CN_NInit = vN

    gravFactory.addBodiesTo(scObject)

    # ****** where spice stuff goes *****

    # create a logging task object of the spacecraft output message at the desired down sampling ratio
    dataRec = scObject.scStateOutMsg.recorder()
    scenarioSim.AddModelToTask(taskName, dataRec)

    scenarioSim.InitializeSimulation()

    # Set the simulation time    
    scenarioSim.ConfigureStopTime(macros.hour2nano(48)) # how long after start time to end sim

    # Begin the simulation time run set above
    scenarioSim.ExecuteSimulation()

    posData = dataRec.r_BN_N
    print(posData[0])

When running this, what's printed is: [-2816801.60102349 5248174.84691615 3677157.2646773 ]

However, when I replace that line that says
# ****** where spice stuff goes *****
with
spiceObject = gravFactory.createSpiceInterface(time=epoch) scenarioSim.AddModelToTask(taskName, spiceObject)

what's printed is [ 6.48448729e+10 -1.25834874e+11 -5.45126724e+10] which is incorrect. How do I fix this?

This is the plot of positions over time for one period without the spice object

This is the plot of positions over time for one period with the spice object

[salomonLydonBusquets]

I had a similar problem recently. I solved it by adding spiceObject.zeroBase = 'Earth' and spiceObject.referenceBase = 'J2000' before adding the spiceObject to the task.