Index: trunk/BNC/src/bncbiassinex.cpp
===================================================================
--- trunk/BNC/src/bncbiassinex.cpp	(revision 9653)
+++ trunk/BNC/src/bncbiassinex.cpp	(revision 9653)
@@ -0,0 +1,164 @@
+/* -------------------------------------------------------------------------
+ * BNC
+ * -------------------------------------------------------------------------
+ *
+ * Class:      bncBiasSinex
+ *
+ * Purpose:    writes SINEX Bias files
+ *
+ * Author:     A. Stuerze
+ *
+ * Created:    01-Mar-2022
+ *
+ * Changes:
+ *
+ * -----------------------------------------------------------------------*/
+
+
+
+#include "bncbiassinex.h"
+
+#include <math.h>
+#include <iomanip>
+
+
+using namespace BNC_PPP;
+using namespace std;
+
+
+// Constructor
+////////////////////////////////////////////////////////////////////////////
+bncBiasSinex::bncBiasSinex(const QString& sklFileName, const QString& intr,
+    int sampl)
+  : bncoutf(sklFileName, intr, sampl) {
+  _sampl =  sampl;
+
+  if (!_sampl) {
+    _sampl = 5;
+  }cout << _sampl << endl;
+
+  _agency = agencyFromFileName();
+}
+
+
+// Destructor
+////////////////////////////////////////////////////////////////////////////
+bncBiasSinex::~bncBiasSinex() {
+  bncoutf::closeFile();
+}
+
+
+t_irc bncBiasSinex::write(int GPSweek, double GPSweeks, const QString& prn, const QString& obsCode, double bias) {
+
+  if (reopen(GPSweek, GPSweeks) == success) {
+
+      QDateTime datTimStart = dateAndTimeFromGPSweek(GPSweek, GPSweeks);
+
+      int daysec    = int(fmod(GPSweeks, 86400.0));
+      int dayOfYear = datTimStart.date().dayOfYear();
+      QString yyyy  = datTimStart.toString("yyyy");
+      QString timeStrStart  = QString("       %1:%2:%3").arg(yyyy)
+                                           .arg(dayOfYear, 3, 10, QLatin1Char('0'))
+                                           .arg(daysec   , 5, 10, QLatin1Char('0'));
+      QString timeStrEnd = QString(" %1:%2:%3").arg(yyyy)
+                                               .arg(dayOfYear, 3, 10, QLatin1Char('0'))
+                                               .arg(daysec+_sampl , 5, 10, QLatin1Char('0'));
+
+      _out << " OSB       " << prn.toLatin1().data()
+           << "           " << obsCode.toLatin1().data()
+           << timeStrStart.toLatin1().data()
+           << timeStrEnd.toLatin1().data()
+           << " ns   "
+           << setw(21) << setprecision(4) << (bias / t_CST::c) * 1.0e+9
+           << endl;
+
+    return success;
+  }
+  else {
+    cout << "return failure;" << endl;
+    return failure;
+  }
+  return success;
+}
+
+void bncBiasSinex::writeHeader(const QDateTime& datTim) {
+  int    GPSWeek;
+  double GPSWeeks;
+  bncSettings settings;
+  GPSweekFromDateAndTime(datTim, GPSWeek, GPSWeeks);
+  int daysec    = int(fmod(GPSWeeks, 86400.0));
+  int dayOfYear = datTim.date().dayOfYear();
+  QString yy    = datTim.toString("yy");
+  QString creationTime = QString("%1:%2:%3").arg(yy)
+                                            .arg(dayOfYear, 3, 10, QLatin1Char('0'))
+                                            .arg(daysec   , 5, 10, QLatin1Char('0'));
+  QString startTime = creationTime;
+  QString intStr = settings.value("uploadIntr").toString();
+  int intr, indHlp = 0;
+  if      ((indHlp = intStr.indexOf("min")) != -1) {
+    intr = intStr.left(indHlp-1).toInt();
+    intr *= 60;
+  }
+  else if ((indHlp = intStr.indexOf("hour")) != -1) {
+    intr = intStr.left(indHlp-1).toInt();
+    intr *= 3600;
+  }
+  else if ((indHlp = intStr.indexOf("day")) != -1) {
+    intr = intStr.left(indHlp-1).toInt();
+    intr *= 86400;
+  }
+
+    int nominalStartSec = daysec - (int(fmod(double(daysec), double(intr))));
+  int nominalEndSec = nominalStartSec  + intr - _sampl;
+  QString endTime = QString("%1:%2:%3").arg(yy)
+                                       .arg(dayOfYear     , 3, 10, QLatin1Char('0'))
+                                       .arg(nominalEndSec , 5, 10, QLatin1Char('0'));
+  int numEpochs = ((nominalEndSec - daysec) / _sampl) +1;
+  QString epo  = QString("%1").arg(numEpochs, 5, 10, QLatin1Char('0'));
+
+  _out << "%=BIA 1.00 " << _agency.toStdString() << " "
+       << creationTime.toStdString() << " " << _agency.toStdString() << " "
+       << startTime.toStdString()    << " " << endTime.toStdString() << " A "
+       << epo.toStdString() << endl;
+
+  _out << "+FILE/REFERENCE" << endl;
+  _out << "*INFO_TYPE_________ INFO________________________________________________________" << endl;
+  _out << " DESCRIPTION       " << " GNSS Satellite Code and Phase Biases from SSR stream  " << endl;
+  _out << " INPUT             " << " SSR satellite biases provided by " << _agency.toStdString() <<endl;
+  _out << " OUTPUT            " << " Absolute (observation-specific) bias parameters" << endl;
+  _out << " SOFTWARE          " << " " << BNCPGMNAME << endl;
+  _out << "-FILE/REFERENCE" << endl << endl;
+
+  _out << "+BIAS/DESCRIPTION" << endl;
+  _out << "*KEYWORD________________________________ VALUE(S)_______________________________"  << endl;
+  _out << " PARAMETER_SAMPLING                     " << " " << fixed   << setw(12) << _sampl  << endl;
+  _out << " DETERMINATION_METHOD                   " << " AC SPECIFIC" << endl;
+  _out << " BIAS_MODE                              " << " ABSOLUTE" << endl;
+  _out << " TIME_SYSTEM                            " << " G" << endl;
+  _out << "-BIAS/DESCRIPTION" << endl << endl;
+
+  _out << "+BIAS/SOLUTION" << endl;
+  _out << "*BIAS SVN_ PRN STATION__ OBS1 OBS2 BIAS_START____ BIAS_END______ UNIT __ESTIMATED_VALUE____ _STD_DEV___ __ESTIMATED_SLOPE____ _STD_DEV__" << endl;
+  return;
+}
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Index: trunk/BNC/src/bncbiassinex.h
===================================================================
--- trunk/BNC/src/bncbiassinex.h	(revision 9653)
+++ trunk/BNC/src/bncbiassinex.h	(revision 9653)
@@ -0,0 +1,39 @@
+/*
+ * bncBiasSinex.h
+ *
+ *  Created on: Feb 18, 2022
+ *      Author: stuerze
+ */
+
+#ifndef SRC_BNCBIASSINEX_H_
+#define SRC_BNCBIASSINEX_H_
+
+#include <fstream>
+#include <newmat.h>
+#include <iostream>
+
+#include "bncoutf.h"
+#include "bncversion.h"
+#include "pppOptions.h"
+#include "bncsettings.h"
+#include "bncantex.h"
+
+using namespace BNC_PPP;
+
+
+class bncBiasSinex : public bncoutf {
+ public:
+  bncBiasSinex(const QString& sklFileName, const QString& intr,
+              int sampl);
+  virtual ~bncBiasSinex();
+  virtual t_irc write(int GPSweek, double GPSweeks, const QString& prn,
+                      const QString& obsCode, double bias);
+
+ private:
+  virtual void writeHeader(const QDateTime& datTim);
+  int     _sampl;
+  QString _agency;
+};
+
+#endif /* SRC_BNCBIASSINEX_H_ */
+