QZSSAlmanac.java

/* Copyright 2002-2026 CS GROUP
 * Licensed to CS GROUP (CS) under one or more
 * contributor license agreements.  See the NOTICE file distributed with
 * this work for additional information regarding copyright ownership.
 * CS licenses this file to You under the Apache License, Version 2.0
 * (the "License"); you may not use this file except in compliance with
 * the License.  You may obtain a copy of the License at
 *
 *   http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
package org.orekit.propagation.analytical.gnss.data;

import org.hipparchus.CalculusFieldElement;
import org.orekit.frames.Frame;
import org.orekit.orbits.FieldKeplerianOrbit;
import org.orekit.orbits.KeplerianOrbit;
import org.orekit.time.FieldGNSSDate;
import org.orekit.time.GNSSDate;
import org.orekit.time.TimeScales;

import java.util.function.DoubleFunction;

/**
 * This class holds a QZSS almanac as read from YUMA files.
 *
 * @author Bryan Cazabonne
 * @since 10.0
 *
 */
public class QZSSAlmanac extends GNSSOrbitalElements<QZSSAlmanac> {

    /** Source of the almanac. */
    private final String source;

    /** Health status. */
    private final int health;

    /**
     * Constructor.
     * @param timeScales known time scales
     * @param prn        PRN number of the satellite
     * @param toe        tim of ephemeris (<em>must</em> be consistent with {@code orbit})
     * @param orbit      Keplerian orbit in Earth-frozen frame
     * @param aDot       change rate in semi-major axis (m/s)
     * @param deltaN0    delta of satellite mean motion
     * @param deltaN0Dot change rate in Δn₀
     * @param iDot       inclination rate (rad/s)
     * @param omegaDot   rate of right ascension (rad/s)
     * @param cuc        amplitude of the cosine harmonic correction term to the argument of latitude
     * @param cus        amplitude of the sine harmonic correction term to the argument of latitude
     * @param crc        amplitude of the cosine harmonic correction term to the orbit radius
     * @param crs        amplitude of the sine harmonic correction term to the orbit radius
     * @param cic        amplitude of the cosine harmonic correction term to the inclination
     * @param cis        amplitude of the sine harmonic correction term to the inclination
     * @param af0        zero-th order clock correction (s)
     * @param af1        first order clock correction (s/s)
     * @param af2        second order clock correction (s/s²)
     * @param tgd        group delay differential TGD for L1-L2 correction
     * @param toc        time of clock
     * @param source     source of the almanac
     * @param health     health status
     */
    public QZSSAlmanac(final TimeScales timeScales, final int prn,
                       final GNSSDate toe, final KeplerianOrbit orbit,
                       final double aDot, final double deltaN0, final double deltaN0Dot,
                       final double iDot, final double omegaDot,
                       final double cuc, final double cus,
                       final double crc, final double crs,
                       final double cic, final double cis,
                       final double af0, final double af1, final double af2,
                       final double tgd, final GNSSDate toc, final String source, final int health) {
        super(GNSSConstants.QZSS_AV, GNSSConstants.QZSS_WEEK_NB, timeScales, null, prn,
              toe, orbit, aDot, deltaN0, deltaN0Dot, iDot, omegaDot,
              cuc, cus, crc, crs, cic, cis, af0, af1, af2, tgd, toc);
        this.source = source;
        this.health = health;
    }

    /** Constructor from field instance.
     * @param <T> type of the field elements
     * @param original regular field instance
     */
    public <T extends CalculusFieldElement<T>> QZSSAlmanac(final FieldQZSSAlmanac<T> original) {
        super(original);
        source = original.getSource();
        health = original.getHealth();
    }

    /** {@inheritDoc} */
    @Override
    public <T extends CalculusFieldElement<T>>
        FieldQZSSAlmanac<T> toField(final FieldKeplerianOrbit<T> orbit,
                                    final T[] nonKeplerian,
                                    final DoubleFunction<T> converter) {
        return new FieldQZSSAlmanac<>(getAngularVelocity(), getWeeksInCycle(), getTimeScales(),
                                      getType(), getPrn(),
                                      new FieldGNSSDate<>(orbit.getDate().getField(), getTimeOfEphemeris()),
                                      orbit, nonKeplerian,
                                      converter.apply(getTgd()),
                                      new FieldGNSSDate<>(orbit.getDate().getField(), getTimeOfClock()),
                                      getSource(), getHealth());
    }

    /** Get the source of this QZSS almanac.
     * @return the source of this QZSS almanac
     */
    public String getSource() {
        return source;
    }

    /** Get the Health status.
     * @return the Health status
     */
    public int getHealth() {
        return health;
    }

    /** {@inheritDoc} */
    @Override
    public QZSSAlmanacFactory baseFactory(final Frame inertial, final Frame bodyFixed) {
        return new QZSSAlmanacFactory(getTimeScales(), getTimeOfEphemeris().getSystem(), inertial, bodyFixed);
    }

}