Changeset 860 in ntrip

Timestamp:

Apr 26, 2008, 11:31:00 AM (16 years ago)

Author:

mervart

Message:

* empty log message *

Location:

trunk/BNS

Files:

• bnsoutf.h (modified) (2 diffs)
• bnssp3.cpp (modified) (1 diff)
• bnsutils.cpp (modified) (1 diff)
• bnsutils.h (modified) (1 diff)

trunk/BNS/bnsoutf.h

@@ -21 +21 @@
   virtual void closeFile();
   std::ofstream _out;
+  int           _sampl;
 
  private:
@@ -28 +29 @@
   void resolveFileName(int GPSweek, const QDateTime& datTim);
 
-  int           _sampl;
   bool          _headerWritten;
   QDateTime     _nextCloseEpoch;

trunk/BNS/bnssp3.cpp

@@ -79 +79 @@
 ////////////////////////////////////////////////////////////////////////////
 void bnsSP3::writeHeader(const QDateTime& datTim) {
-  _out << "#cP2007  7  1  0  0  0.00000000      96 ORBIT IGS05 HLM  IGS"
-       << "## 1434      0.00000000   900.00000000 54282 0.0000000000000"
+
+  int    GPSWeek;
+  double GPSWeeks;
+  GPSweekFromDateAndTime(datTim, GPSWeek, GPSWeeks);
+
+  double sec = fmod(GPSWeeks, 60.0);
+
+  int    mjd;
+  double dayfrac;
+  mjdFromDateAndTime(datTim, mjd, dayfrac);
+
+  _out << "#aP" << datTim.toString("yyyy MM dd hh mm").toAscii().data() 
+       << setw(12) << setprecision(8) << sec
+       << "      96 ORBIT IGS05 HLM  IGS" << endl;
+
+  _out << "## " 
+       << setw(4)  << GPSWeek
+       << setw(16) << setprecision(8) << GPSWeeks
+       << setw(15) << setprecision(8) << double(_sampl)
+       << setw(6)  << mjd
+       << setw(16) << setprecision(13) << dayfrac << endl;
 
   _out << "+   32   G01G02G03G04G05G06G07G08G09G10G11G12G13G14G15G16G17\n"

trunk/BNS/bnsutils.cpp

@@ -92 +92 @@
         currTime.second()                          +
         currTime.msec()                   / 1000.0;
+}
+
+// 
+////////////////////////////////////////////////////////////////////////////
+void mjdFromDateAndTime(const QDateTime& dateTime, int& mjd, double& dayfrac) {
+
+  double mjddec = dateTime.date().toJulianDay() - 2400000.5 + 
+                   (dateTime.time().hour() +
+                    (dateTime.time().minute() +
+                     (dateTime.time().second() + 
+                      dateTime.time().msec() / 1000.0) / 60.0) / 60.0) / 24.0;
+
+  mjd     = int(mjddec);
+  dayfrac = mjddec - mjd;
 }
 

trunk/BNS/bnsutils.h

@@ -19 +19 @@
                             int& GPSWeek, double& GPSWeeks);
 
+void mjdFromDateAndTime(const QDateTime& dateTime, int& mjd, double& dayfrac);
+
 void currentGPSWeeks(int& week, double& sec);