Examples with DSSTForceModel org.orekit.propagation.semianalytical.dsst.forces.DSSTForceModel used on opensource projects

Example 16 with DSSTForceModel

use of org.orekit.propagation.semianalytical.dsst.forces.DSSTForceModel in project Orekit by CS-SI.

the class DSSTPropagatorTest method testGetInitialOsculatingState.

@Test
public void testGetInitialOsculatingState() throws IllegalArgumentException, OrekitException {
    final SpacecraftState initialState = getGEOState();
    // build integrator
    final double minStep = initialState.getKeplerianPeriod() * 0.1;
    final double maxStep = initialState.getKeplerianPeriod() * 10.0;
    final double[][] tol = DSSTPropagator.tolerances(0.1, initialState.getOrbit());
    AdaptiveStepsizeIntegrator integrator = new DormandPrince853Integrator(minStep, maxStep, tol[0], tol[1]);
    // build the propagator for the propagation of the mean elements
    DSSTPropagator prop = new DSSTPropagator(integrator, true);
    final UnnormalizedSphericalHarmonicsProvider provider = GravityFieldFactory.getUnnormalizedProvider(4, 0);
    final Frame earthFrame = CelestialBodyFactory.getEarth().getBodyOrientedFrame();
    DSSTForceModel zonal = new DSSTZonal(provider, 4, 3, 9);
    DSSTForceModel tesseral = new DSSTTesseral(earthFrame, Constants.WGS84_EARTH_ANGULAR_VELOCITY, provider, 4, 0, 4, 8, 4, 0, 2);
    prop.addForceModel(zonal);
    prop.addForceModel(tesseral);
    // Set the initial state as osculating
    prop.setInitialState(initialState, false);
    // Check the stored initial state is the osculating one
    Assert.assertEquals(initialState, prop.getInitialState());
    // Check that no propagation, i.e. propagation to the initial date, provides the initial
    // osculating state although the propagator is configured to propagate mean elements !!!
    Assert.assertEquals(initialState, prop.propagate(initialState.getDate()));
}