GLOFNavUT1TimeOffset.cpp

Go to the documentation of this file.

//==============================================================================
//
//  This file is part of GNSSTk, the ARL:UT GNSS Toolkit.
//
//  The GNSSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 3.0 of the License, or
//  any later version.
//
//  The GNSSTk is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GNSSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin.
//  Copyright 2004-2022, The Board of Regents of The University of Texas System
//
//==============================================================================
 
 
//==============================================================================
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin, under contract to an agency or agencies
//  within the U.S. Department of Defense. The U.S. Government retains all
//  rights to use, duplicate, distribute, disclose, or release this software.
//
//  Pursuant to DoD Directive 523024
//
//  DISTRIBUTION STATEMENT A: This software has been approved for public
//                            release, distribution is unlimited.
//
//==============================================================================
#include "GLOFNavUT1TimeOffset.hpp"
#include "TimeString.hpp"
 
using namespace std;
 
namespace gnsstk
{
   GLOFNavUT1TimeOffset ::
   GLOFNavUT1TimeOffset()
         : tauc(std::numeric_limits<double>::quiet_NaN()),
           B1(std::numeric_limits<double>::quiet_NaN()),
           B2(std::numeric_limits<double>::quiet_NaN()),
           KP(-1), NT(-1), NA(-1)
   {
      msgLenSec = 2.0;
      weekFmt = "";
   }
 
 
   void GLOFNavUT1TimeOffset ::
   dump(std::ostream& s, DumpDetail dl) const
   {
      const ios::fmtflags oldFlags = s.flags();
      if (dl == DumpDetail::Terse)
      {
         dumpTerse(s);
         return;
      }
      if (dl == DumpDetail::OneLine)
      {
         NavData::dump(s,dl);
         return;
      }
      if (dl == DumpDetail::Brief)
      {
         double offset;
         getOffset(TimeSystem::GLO, TimeSystem::UTC, timeStamp, offset);
            // brief just shows the offset as of the reference time.
         s << StringUtils::asString(TimeSystem::GLO) << "-"
           << StringUtils::asString(TimeSystem::UTC)
           << " offset = " << offset << endl;
         return;
      }
      s << "****************************************************************"
        << "************" << endl
        << "Time System Offset"
        << endl
        << endl
        << getSignalString() << endl
        << endl
        << "           TIMES OF INTEREST" << endl << endl
        << "              " << getDumpTimeHdr(dl) << endl
        << "Transmit:     " << getDumpTime(dl, timeStamp)
        << endl << endl
        << "           " << StringUtils::asString(TimeSystem::GLO) << " "
        << StringUtils::asString(TimeSystem::UTC) << " PARAMETERS" << endl
        << "Parameter              Value" << endl;
      s.setf(ios::scientific, ios::floatfield);
      s.setf(ios::right, ios::adjustfield);
      s.setf(ios::uppercase);
      s.precision(6);
      s.fill(' ');
      s << "tau_c         " << setw(14) << tauc    << " sec" << endl
        << "B1            " << setw(14) << B1      << " sec" << endl
        << "B2            " << setw(14) << B2      << " sec/sec" << endl
        << "KP              " << setw(12) << KP;
      if (KP==0)
      {
         s << " No leap second this quarter";
      }
      else if (KP==1)
      {
         s << " +1 leap second this quarter";
      }
      else if (KP==3)
      {
         s << " -1 leap second this quarter";
      }
      else
      {
         s << " Unknown meaning";
      }
      s << endl
        << "NT              " << setw(12) << NT << " days" << endl
        << "NA              " << setw(12) << NA << " days" << endl;
      s.flags(oldFlags);
   }
 
 
   void GLOFNavUT1TimeOffset ::
   dumpTerse(ostream& s) const
   {
      s << std::setw(4) << signal.sat.id
        << printTime(timeStamp,"%02m/%02d/%04Y %03j %02H:%02M:%02S") << "! ";
      s.setf(ios::scientific, ios::floatfield);
      s.setf(ios::right, ios::adjustfield);
      s.setf(ios::uppercase);
      s.precision(9);
      s.fill(' ');
      s << "B1:" << B1;
      s.precision(6);
      s << " s, B2:" << B2;
 
      s.precision(0);
      s << " s/s, KP:" << KP;
      s << " ";
   }
 
 
   bool GLOFNavUT1TimeOffset ::
   getOffset(TimeSystem fromSys, TimeSystem toSys,
             const CommonTime& when, double& offset)
      const
   {
      if (((fromSys == TimeSystem::GLO) && (toSys == TimeSystem::UTC)) ||
          ((fromSys == TimeSystem::UTC) && (toSys == TimeSystem::GLO)))
      {
         GNSSTK_ASSERT(when.getTimeSystem() == fromSys);
         long dt = NT-NA;
         double dUT1 = B1 + B2*dt;
            // T_GLO = T_UTC(SU) + 3 hours, or
            // T_UTC(SU) = T_GLO - 3 hours
            // UT1 = UTC(SU) - dUT1
         offset = -10800 + tauc + dUT1;
         if (fromSys == TimeSystem::UTC)
            offset = -offset;
         return true;
      }
      return false;
   }
 
 
   TimeCvtSet GLOFNavUT1TimeOffset ::
   getConversions() const
   {
      TimeCvtKey keyF(TimeSystem::GLO,TimeSystem::UTC);
      TimeCvtKey keyR(TimeSystem::UTC,TimeSystem::GLO);
      return TimeCvtSet({ keyF, keyR });
   }
 
 
   bool GLOFNavUT1TimeOffset ::
   validate() const
   {
      return true;
   }
}