@Test
  public void testOnlyPre1972Data() throws OrekitException {

    Utils.setDataRoot("USNO");
    TimeScalesFactory.addUTCTAIOffsetsLoader(
        new TAIUTCDatFilesLoader("tai-utc-only-pre-1972-data.dat"));

    // linear models between 1961 and 1972
    checkOffset(1961, 1, 2, -(1.422818 + 1 * 0.001296)); // MJD 37300 +   1
    checkOffset(1961, 8, 2, -(1.372818 + 213 * 0.001296)); // MJD 37300 + 213
    checkOffset(1962, 1, 2, -(1.845858 + 1 * 0.0011232)); // MJD 37665 +   1
    checkOffset(1963, 11, 2, -(1.945858 + 670 * 0.0011232)); // MJD 37665 + 670
    checkOffset(1964, 1, 2, -(3.240130 - 365 * 0.001296)); // MJD 38761 - 365
    checkOffset(1964, 4, 2, -(3.340130 - 274 * 0.001296)); // MJD 38761 - 274
    checkOffset(1964, 9, 2, -(3.440130 - 121 * 0.001296)); // MJD 38761 - 121
    checkOffset(1965, 1, 2, -(3.540130 + 1 * 0.001296)); // MJD 38761 +   1
    checkOffset(1965, 3, 2, -(3.640130 + 60 * 0.001296)); // MJD 38761 +  60
    checkOffset(1965, 7, 2, -(3.740130 + 182 * 0.001296)); // MJD 38761 + 182
    checkOffset(1965, 9, 2, -(3.840130 + 244 * 0.001296)); // MJD 38761 + 244
    checkOffset(1966, 1, 2, -(4.313170 + 1 * 0.002592)); // MJD 39126 +   1
    checkOffset(1968, 2, 2, -(4.213170 + 762 * 0.002592)); // MJD 39126 + 762

    // last linear drift is not stopped as we miss the first constant offset in 1972
    checkOffset(1972, 3, 5, -(4.213170 + 2255 * 0.002592)); // MJD 39126 +  2255
    checkOffset(1972, 7, 14, -(4.213170 + 2386 * 0.002592)); // MJD 39126 +  2386
    checkOffset(1979, 12, 31, -(4.213170 + 5112 * 0.002592)); // MJD 39126 +  5112
    checkOffset(1980, 1, 22, -(4.213170 + 5134 * 0.002592)); // MJD 39126 +  5134
    checkOffset(2006, 7, 7, -(4.213170 + 14797 * 0.002592)); // MJD 39126 + 14797
    checkOffset(2010, 7, 7, -(4.213170 + 16258 * 0.002592)); // MJD 39126 + 16258
    checkOffset(2012, 7, 7, -(4.213170 + 16989 * 0.002592)); // MJD 39126 + 16989
    checkOffset(2015, 7, 7, -(4.213170 + 18084 * 0.002592)); // MJD 39126 + 18084
  }
  @Test
  public void testModifiedLinearData() throws OrekitException {

    Utils.setDataRoot("USNO");
    TimeScalesFactory.addUTCTAIOffsetsLoader(
        new TAIUTCDatFilesLoader("tai-utc-modified-linear.dat"));

    // linear models between 1961 and 1972
    checkOffset(1961, 1, 2, -(1.4000000 + 1 * 0.001000)); // MJD 37300 +   1
    checkOffset(1961, 8, 2, -(1.5000000 + 213 * 0.001100)); // MJD 37300 + 213
    checkOffset(1962, 1, 2, -(1.6000000 + 1 * 0.001200)); // MJD 37665 +   1
    checkOffset(1963, 11, 2, -(1.7000000 + 670 * 0.001300)); // MJD 37665 + 670
    checkOffset(1964, 1, 2, -(1.8000000 - 365 * 0.001400)); // MJD 38761 - 365
    checkOffset(1964, 4, 2, -(1.9000000 - 274 * 0.001500)); // MJD 38761 - 274
    checkOffset(1964, 9, 2, -(2.0000000 - 121 * 0.001600)); // MJD 38761 - 121
    checkOffset(1965, 1, 2, -(2.1000000 + 1 * 0.001700)); // MJD 38761 +   1
    checkOffset(1965, 3, 2, -(2.2000000 + 60 * 0.001800)); // MJD 38761 +  60
    checkOffset(1965, 7, 2, -(2.3000000 + 182 * 0.001900)); // MJD 38761 + 182
    checkOffset(1965, 9, 2, -(2.4000000 + 244 * 0.002000)); // MJD 38761 + 244
    checkOffset(1966, 1, 2, -(2.5000000 + 1 * 0.002100)); // MJD 39126 +   1
    checkOffset(1968, 2, 2, -(2.6000000 + 762 * 0.002200)); // MJD 39126 + 762

    // last linear drift is not stopped as we miss the first constant offset in 1972
    checkOffset(1972, 3, 5, -(2.6000000 + 2255 * 0.002200)); // MJD 39126 +  2255
    checkOffset(1972, 7, 14, -(2.6000000 + 2386 * 0.002200)); // MJD 39126 +  2386
    checkOffset(1979, 12, 31, -(2.6000000 + 5112 * 0.002200)); // MJD 39126 +  5112
    checkOffset(1980, 1, 22, -(2.6000000 + 5134 * 0.002200)); // MJD 39126 +  5134
    checkOffset(2006, 7, 7, -(2.6000000 + 14797 * 0.002200)); // MJD 39126 + 14797
    checkOffset(2010, 7, 7, -(2.6000000 + 16258 * 0.002200)); // MJD 39126 + 16258
    checkOffset(2012, 7, 7, -(2.6000000 + 16989 * 0.002200)); // MJD 39126 + 16989
    checkOffset(2015, 7, 7, -(2.6000000 + 18084 * 0.002200)); // MJD 39126 + 18084
  }
Esempio n. 3
0
 /**
  * load orekit data and gravity field.
  *
  * @throws Exception on error.
  */
 @BeforeClass
 public static void setUpBefore() throws Exception {
   Utils.setDataRoot("earth:geoid:regular-data");
   GravityFieldFactory.clearPotentialCoefficientsReaders();
   GravityFieldFactory.addPotentialCoefficientsReader(new EGMFormatReader("egm96", false));
   potential = GravityFieldFactory.getConstantNormalizedProvider(maxDegree, maxOrder);
 }
  @Test
  public void testRegularFile() throws OrekitException {

    Utils.setDataRoot("USNO");

    // we arbitrary put UTC == TAI before 1961-01-01
    checkOffset(1950, 1, 1, 0);

    // linear models between 1961 and 1972
    checkOffset(1961, 1, 2, -(1.422818 + 1 * 0.001296)); // MJD 37300 +   1
    checkOffset(1961, 8, 2, -(1.372818 + 213 * 0.001296)); // MJD 37300 + 213
    checkOffset(1962, 1, 2, -(1.845858 + 1 * 0.0011232)); // MJD 37665 +   1
    checkOffset(1963, 11, 2, -(1.945858 + 670 * 0.0011232)); // MJD 37665 + 670
    checkOffset(1964, 1, 2, -(3.240130 - 365 * 0.001296)); // MJD 38761 - 365
    checkOffset(1964, 4, 2, -(3.340130 - 274 * 0.001296)); // MJD 38761 - 274
    checkOffset(1964, 9, 2, -(3.440130 - 121 * 0.001296)); // MJD 38761 - 121
    checkOffset(1965, 1, 2, -(3.540130 + 1 * 0.001296)); // MJD 38761 +   1
    checkOffset(1965, 3, 2, -(3.640130 + 60 * 0.001296)); // MJD 38761 +  60
    checkOffset(1965, 7, 2, -(3.740130 + 182 * 0.001296)); // MJD 38761 + 182
    checkOffset(1965, 9, 2, -(3.840130 + 244 * 0.001296)); // MJD 38761 + 244
    checkOffset(1966, 1, 2, -(4.313170 + 1 * 0.002592)); // MJD 39126 +   1
    checkOffset(1968, 2, 2, -(4.213170 + 762 * 0.002592)); // MJD 39126 + 762

    // since 1972-01-01, offsets are only whole seconds
    checkOffset(1972, 3, 5, -10);
    checkOffset(1972, 7, 14, -11);
    checkOffset(1979, 12, 31, -18);
    checkOffset(1980, 1, 22, -19);
    checkOffset(2006, 7, 7, -33);
    checkOffset(2010, 7, 7, -34);
    checkOffset(2012, 7, 7, -35);
    checkOffset(2015, 7, 7, -36);
  }
  @Deprecated
  @Test
  public void testDeprecatedAPI() throws OrekitException {

    Utils.setDataRoot("USNO");
    TimeScalesFactory.addUTCTAILoader(
        new UTCTAILoader() {

          @Override
          public boolean stillAcceptsData() {
            return false;
          }

          @Override
          public void loadData(InputStream input, String name) {}

          @Override
          public SortedMap<DateComponents, Integer> loadTimeSteps() {
            SortedMap<DateComponents, Integer> map = new TreeMap<DateComponents, Integer>();
            map.put(new DateComponents(1972, 1, 1), 10);
            map.put(new DateComponents(1972, 7, 1), 11);
            return map;
          }

          @Override
          public String getSupportedNames() {
            return "non-existant-name";
          }
        });

    checkOffset(1972, 1, 2, -10.0);
    checkOffset(1972, 6, 30, -10.0);
    checkOffset(1972, 7, 2, -11.0);
    checkOffset(3000, 1, 1, -11.0);
  }
 private void checkException(String name, OrekitMessages message) {
   Utils.setDataRoot("USNO");
   TimeScalesFactory.addUTCTAIOffsetsLoader(new TAIUTCDatFilesLoader(name));
   try {
     TimeScalesFactory.getUTC();
     Assert.fail("an exception should have been thrown");
   } catch (OrekitException oe) {
     Assert.assertEquals(message, oe.getSpecifier());
   }
 }
Esempio n. 7
0
 @Before
 public void setUp() throws OrekitException {
   Utils.setDataRoot("regular-data");
   geoTLE =
       new TLE(
           "1 27508U 02040A   12021.25695307 -.00000113  00000-0  10000-3 0  7326",
           "2 27508   0.0571 356.7800 0005033 344.4621 218.7816  1.00271798 34501");
   leoTLE =
       new TLE(
           "1 31135U 07013A   11003.00000000  .00000816  00000+0  47577-4 0    11",
           "2 31135   2.4656 183.9084 0021119 236.4164  60.4567 15.10546832    15");
 }
 @Before
 public void setUp() throws OrekitException {
   Utils.setDataRoot("regular-data:potential/shm-format");
   GravityFieldFactory.addPotentialCoefficientsReader(
       new SHMFormatReader("^eigen_cg03c_coef$", false));
   mu = GravityFieldFactory.getUnnormalizedProvider(0, 0).getMu();
   final Vector3D position = new Vector3D(7.0e6, 1.0e6, 4.0e6);
   final Vector3D velocity = new Vector3D(-500.0, 8000.0, 1000.0);
   initDate = AbsoluteDate.J2000_EPOCH;
   final Orbit orbit =
       new EquinoctialOrbit(
           new PVCoordinates(position, velocity), FramesFactory.getEME2000(), initDate, mu);
   initialState = new SpacecraftState(orbit);
   double[][] tolerance = NumericalPropagator.tolerances(0.001, orbit, OrbitType.EQUINOCTIAL);
   AdaptiveStepsizeIntegrator integrator =
       new DormandPrince853Integrator(0.001, 200, tolerance[0], tolerance[1]);
   integrator.setInitialStepSize(60);
   propagator = new NumericalPropagator(integrator);
   propagator.setInitialState(initialState);
 }
  @Before
  public void setUp() throws OrekitException {

    Utils.setDataRoot("regular-data");
    earth =
        new OneAxisEllipsoid(
            Constants.WGS84_EARTH_EQUATORIAL_RADIUS,
            Constants.WGS84_EARTH_FLATTENING,
            FramesFactory.getITRF(IERSConventions.IERS_2010, true));

    gp = new GeodeticPoint(FastMath.toRadians(51.0), FastMath.toRadians(66.6), 300.0);
    final TimeScale utc = TimeScalesFactory.getUTC();
    final Vector3D position = new Vector3D(-6142438.668, 3492467.56, -25767.257);
    final Vector3D velocity = new Vector3D(505.848, 942.781, 7435.922);
    final AbsoluteDate date = new AbsoluteDate(2003, 9, 16, utc);
    orbit =
        new EquinoctialOrbit(
            new PVCoordinates(position, velocity),
            FramesFactory.getEME2000(),
            date,
            Constants.EIGEN5C_EARTH_MU);
  }
Esempio n. 10
0
  @Before
  public void setUp() {
    try {

      Utils.setDataRoot("regular-data");

      // Computation date
      date =
          new AbsoluteDate(
              new DateComponents(2008, 04, 07), TimeComponents.H00, TimeScalesFactory.getUTC());

      // Body mu
      mu = 3.9860047e14;

      // Reference frame = ITRF
      frameItrf = FramesFactory.getITRF(IERSConventions.IERS_2010, true);

      // Elliptic earth shape
      earthSpheric = new OneAxisEllipsoid(6378136.460, 0., frameItrf);

    } catch (OrekitException oe) {
      Assert.fail(oe.getMessage());
    }
  }
 @Before
 public void setUp() {
   Utils.setDataRoot("regular-data");
 }
 /** set Orekit data for useful debugging messages from dates. */
 @BeforeClass
 public static void setUpBefore() {
   Utils.setDataRoot("regular-data");
 }
Esempio n. 13
0
 @Before
 public void setUp() throws OrekitException {
   Utils.setDataRoot("regular-data:potential/shm-format");
   GravityFieldFactory.addPotentialCoefficientsReader(
       new SHMFormatReader("^eigen_cg03c_coef$", false));
 }