Context Navigation

source: ntrip/trunk/BNC/src/upload/bncrtnetuploadcaster.cpp@ 8542

Last change on this file since 8542 was 8542, checked in by stuerze, 5 years ago
minor changes to consider a satellite clock drift that is introduced via RTNET format
File size: 99.4 KB

Line
1	/* -------------------------------------------------------------------------
2	* BKG NTRIP Server
3	* -------------------------------------------------------------------------
4	*
5	* Class: bncRtnetUploadCaster
6	*
7	* Purpose: Connection to NTRIP Caster
8	*
9	* Author: L. Mervart
10	*
11	* Created: 29-Mar-2011
12	*
13	* Changes:
14	*
15	* -----------------------------------------------------------------------*/
16
17	#include <math.h>
18	#include "bncrtnetuploadcaster.h"
19	#include "bncsettings.h"
20	#include "bncephuser.h"
21	#include "bncclockrinex.h"
22	#include "bncsp3.h"
23	#include "gnss.h"
24
25	using namespace std;
26
27	// Constructor
28	////////////////////////////////////////////////////////////////////////////
29	bncRtnetUploadCaster::bncRtnetUploadCaster(const QString& mountpoint,
30	const QString& outHost, int outPort,
31	const QString& ntripVersion,
32	const QString& userName, const QString& password,
33	const QString& crdTrafo, bool CoM, const QString& sp3FileName,
34	const QString& rnxFileName, int PID, int SID, int IOD, int iRow) :
35	bncUploadCaster(mountpoint, outHost, outPort, ntripVersion, userName, password, iRow, 0) {
36
37
38	if (!mountpoint.isEmpty()) {
39	_casterID += mountpoint;
40	}
41	if (!outHost.isEmpty()) {
42	_casterID += " " + outHost;
43	if (outPort) {
44	_casterID += ":" + QString("%1").arg(outPort, 10);
45	}
46	}
47	if (!crdTrafo.isEmpty()) {
48	_casterID += " " + crdTrafo;
49	}
50	if (!sp3FileName.isEmpty()) {
51	_casterID += " " + sp3FileName;
52	}
53	if (!rnxFileName.isEmpty()) {
54	_casterID += " " + rnxFileName;
55	}
56
57	_crdTrafo = crdTrafo;
58	_CoM = CoM;
59	_PID = PID;
60	_SID = SID;
61	_IOD = IOD;
62
63	// Member that receives the ephemeris
64	// ----------------------------------
65	_ephUser = new bncEphUser(true);
66
67	bncSettings settings;
68	QString intr = settings.value("uploadIntr").toString();
69	QStringList hlp = settings.value("cmbStreams").toStringList();
70	_samplRtcmEphCorr = settings.value("uploadSamplRtcmEphCorr").toDouble();
71	if (hlp.size() > 1) { // combination stream upload
72	_samplRtcmClkCorr = settings.value("cmbSampl").toInt();
73	}
74	else { // single stream upload or sp3 file generation
75	_samplRtcmClkCorr = 5; // default
76	}
77	int samplClkRnx = settings.value("uploadSamplClkRnx").toInt();
78	int samplSp3 = settings.value("uploadSamplSp3").toInt() * 60;
79
80	if (_samplRtcmEphCorr == 0.0) {
81	_usedEph = 0;
82	}
83	else {
84	_usedEph = new QMap<QString, const t_eph*>;
85	}
86
87	// RINEX writer
88	// ------------
89	if (!rnxFileName.isEmpty()) {
90	_rnx = new bncClockRinex(rnxFileName, intr, samplClkRnx);
91	}
92	else {
93	_rnx = 0;
94	}
95
96	// SP3 writer
97	// ----------
98	if (!sp3FileName.isEmpty()) {
99	_sp3 = new bncSP3(sp3FileName, intr, samplSp3);
100	}
101	else {
102	_sp3 = 0;
103	}
104
105	// Set Transformation Parameters
106	// -----------------------------
107	// Transformation Parameters from ITRF2014 to ETRF2000
108	if (_crdTrafo == "ETRF2000") {
109	_dx = 0.0537;
110	_dy = 0.0512;
111	_dz = -0.0551;
112	_dxr = 0.0001;
113	_dyr = 0.0001;
114	_dzr = -0.0019;
115	_ox = -0.000891;
116	_oy = -0.005390;
117	_oz = 0.008712;
118	_oxr = -0.000081;
119	_oyr = -0.000490;
120	_ozr = +0.000792;
121	_sc = 1.02;
122	_scr = 0.11;
123	_t0 = 2000.0;
124	}
125	// Transformation Parameters from ITRF2014 to GDA2020 (Ryan Ruddick, GA)
126	else if (_crdTrafo == "GDA2020") {
127	_dx = 0.0;
128	_dy = 0.0;
129	_dz = 0.0;
130	_dxr = 0.0;
131	_dyr = 0.0;
132	_dzr = 0.0;
133	_ox = 0.0;
134	_oy = 0.0;
135	_oz = 0.0;
136	_oxr = 0.00150379;
137	_oyr = 0.00118346;
138	_ozr = 0.00120716;
139	_sc = 0.0;
140	_scr = 0.0;
141	_t0 = 2020.0;
142	}
143	// Transformation Parameters from IGb14 to SIRGAS2000 (Sonia Costa, BRA)
144	else if (_crdTrafo == "SIRGAS2000") {
145	_dx = 0.0026;
146	_dy = 0.0018;
147	_dz = -0.0061;
148	_dxr = 0.0000;
149	_dyr = 0.0000;
150	_dzr = 0.0000;
151	_ox = 0.000170;
152	_oy = -0.000030;
153	_oz = 0.000070;
154	_oxr = 0.000000;
155	_oyr = 0.000000;
156	_ozr = 0.000000;
157	_sc = -1.000;
158	_scr = 0.000;
159	_t0 = 2000.4;
160	}
161	// Transformation Parameters from ITRF2014 to DREF91
162	else if (_crdTrafo == "DREF91") {
163	_dx = 0.0537;
164	_dy = 0.0512;
165	_dz = -0.0551;
166	_dxr = 0.0001;
167	_dyr = 0.0001;
168	_dzr = -0.0019;
169	_ox = -0.000233;
170	_oy = -0.005598;
171	_oz = 0.009467;
172	_oxr = -0.000081;
173	_oyr = -0.000490;
174	_ozr = 0.000792;
175	_sc = 1.02;
176	_scr = 0.11;
177	_t0 = 2000.0;
178	}
179	else if (_crdTrafo == "Custom") {
180	_dx = settings.value("trafo_dx").toDouble();
181	_dy = settings.value("trafo_dy").toDouble();
182	_dz = settings.value("trafo_dz").toDouble();
183	_dxr = settings.value("trafo_dxr").toDouble();
184	_dyr = settings.value("trafo_dyr").toDouble();
185	_dzr = settings.value("trafo_dzr").toDouble();
186	_ox = settings.value("trafo_ox").toDouble();
187	_oy = settings.value("trafo_oy").toDouble();
188	_oz = settings.value("trafo_oz").toDouble();
189	_oxr = settings.value("trafo_oxr").toDouble();
190	_oyr = settings.value("trafo_oyr").toDouble();
191	_ozr = settings.value("trafo_ozr").toDouble();
192	_sc = settings.value("trafo_sc").toDouble();
193	_scr = settings.value("trafo_scr").toDouble();
194	_t0 = settings.value("trafo_t0").toDouble();
195	}
196	}
197
198	// Destructor
199	////////////////////////////////////////////////////////////////////////////
200	bncRtnetUploadCaster::~bncRtnetUploadCaster() {
201	if (isRunning()) {
202	wait();
203	}
204	delete _rnx;
205	delete _sp3;
206	delete _ephUser;
207	delete _usedEph;
208	}
209
210	//
211	////////////////////////////////////////////////////////////////////////////
212	void bncRtnetUploadCaster::decodeRtnetStream(char* buffer, int bufLen) {
213
214	QMutexLocker locker(&_mutex);
215
216	// Append to internal buffer
217	// -------------------------
218	_rtnetStreamBuffer.append(QByteArray(buffer, bufLen));
219
220	// Select buffer part that contains last epoch
221	// -------------------------------------------
222	QStringList lines;
223	int iEpoBeg = _rtnetStreamBuffer.lastIndexOf('*'); // begin of last epoch
224	if (iEpoBeg == -1) {
225	_rtnetStreamBuffer.clear();
226	return;
227	}
228	int iEpoBegEarlier = _rtnetStreamBuffer.indexOf('*');
229	if (iEpoBegEarlier != -1 && iEpoBegEarlier < iEpoBeg) { // are there two epoch lines in buffer?
230	_rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoBegEarlier);
231	}
232	else {
233	_rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoBeg);
234	}
235
236	int iEpoEnd = _rtnetStreamBuffer.lastIndexOf("EOE"); // end of last epoch
237	if (iEpoEnd == -1) {
238	return;
239	}
240	else {
241	lines = _rtnetStreamBuffer.left(iEpoEnd).split('\n',
242	QString::SkipEmptyParts);
243	_rtnetStreamBuffer = _rtnetStreamBuffer.mid(iEpoEnd + 3);
244	}
245
246	if (lines.size() < 2) {
247	return;
248	}
249
250	// Read first line (with epoch time)
251	// ---------------------------------
252	QTextStream in(lines[0].toLatin1());
253	QString hlp;
254	int year, month, day, hour, min;
255	double sec;
256	in >> hlp >> year >> month >> day >> hour >> min >> sec;
257	bncTime epoTime;
258	epoTime.set(year, month, day, hour, min, sec);
259
260	emit(newMessage(
261	"bncRtnetUploadCaster: decode " + QByteArray(epoTime.datestr().c_str())
262	+ " " + QByteArray(epoTime.timestr().c_str()) + " "
263	+ _casterID.toLatin1(), false));
264
265	struct ClockOrbit co;
266	memset(&co, 0, sizeof(co));
267	co.EpochTime[CLOCKORBIT_SATGPS] = static_cast<int>(epoTime.gpssec());
268	double gt = epoTime.gpssec() + 3 * 3600 - gnumleap(year, month, day);
269	co.EpochTime[CLOCKORBIT_SATGLONASS] = static_cast<int>(fmod(gt, 86400.0));
270	co.EpochTime[CLOCKORBIT_SATGALILEO] = static_cast<int>(epoTime.gpssec());
271	co.EpochTime[CLOCKORBIT_SATQZSS] = static_cast<int>(epoTime.gpssec());
272	co.EpochTime[CLOCKORBIT_SATSBAS] = static_cast<int>(epoTime.gpssec());
273	co.EpochTime[CLOCKORBIT_SATBDS] = static_cast<int>(epoTime.bdssec());
274	co.Supplied[COBOFS_CLOCK] = 1;
275	co.Supplied[COBOFS_ORBIT] = 1;
276	co.SatRefDatum = DATUM_ITRF;
277	co.SSRIOD = _IOD;
278	co.SSRProviderID = _PID; // 256 .. BKG, 257 ... EUREF
279	co.SSRSolutionID = _SID;
280
281	struct CodeBias bias;
282	memset(&bias, 0, sizeof(bias));
283	bias.EpochTime[CLOCKORBIT_SATGPS] = co.EpochTime[CLOCKORBIT_SATGPS];
284	bias.EpochTime[CLOCKORBIT_SATGLONASS] = co.EpochTime[CLOCKORBIT_SATGLONASS];
285	bias.EpochTime[CLOCKORBIT_SATGALILEO] = co.EpochTime[CLOCKORBIT_SATGALILEO];
286	bias.EpochTime[CLOCKORBIT_SATQZSS] = co.EpochTime[CLOCKORBIT_SATQZSS];
287	bias.EpochTime[CLOCKORBIT_SATSBAS] = co.EpochTime[CLOCKORBIT_SATSBAS];
288	bias.EpochTime[CLOCKORBIT_SATBDS] = co.EpochTime[CLOCKORBIT_SATBDS];
289	bias.SSRIOD = _IOD;
290	bias.SSRProviderID = _PID;
291	bias.SSRSolutionID = _SID;
292
293	struct PhaseBias phasebias;
294	memset(&phasebias, 0, sizeof(phasebias));
295	unsigned int dispersiveBiasConsistenyIndicator = 0;
296	unsigned int mwConsistencyIndicator = 0;
297	phasebias.EpochTime[CLOCKORBIT_SATGPS] = co.EpochTime[CLOCKORBIT_SATGPS];
298	phasebias.EpochTime[CLOCKORBIT_SATGLONASS] = co.EpochTime[CLOCKORBIT_SATGLONASS];
299	phasebias.EpochTime[CLOCKORBIT_SATGALILEO] = co.EpochTime[CLOCKORBIT_SATGALILEO];
300	phasebias.EpochTime[CLOCKORBIT_SATQZSS] = co.EpochTime[CLOCKORBIT_SATQZSS];
301	phasebias.EpochTime[CLOCKORBIT_SATSBAS] = co.EpochTime[CLOCKORBIT_SATSBAS];
302	phasebias.EpochTime[CLOCKORBIT_SATBDS] = co.EpochTime[CLOCKORBIT_SATBDS];
303	phasebias.SSRIOD = _IOD;
304	phasebias.SSRProviderID = _PID;
305	phasebias.SSRSolutionID = _SID;
306
307	struct VTEC vtec;
308	memset(&vtec, 0, sizeof(vtec));
309	vtec.EpochTime = static_cast<int>(epoTime.gpssec());
310	vtec.SSRIOD = _IOD;
311	vtec.SSRProviderID = _PID;
312	vtec.SSRSolutionID = _SID;
313
314	// Default Update Interval
315	// -----------------------
316	int clkUpdInd = 2; // 5 sec
317	int ephUpdInd = clkUpdInd; // default
318
319	if (_samplRtcmClkCorr > 5.0 && _samplRtcmEphCorr <= 5.0) { // combined orb and clock
320	ephUpdInd = determineUpdateInd(_samplRtcmClkCorr);
321	}
322	if (_samplRtcmClkCorr > 5.0) {
323	clkUpdInd = determineUpdateInd(_samplRtcmClkCorr);
324	}
325	if (_samplRtcmEphCorr > 5.0) {
326	ephUpdInd = determineUpdateInd(_samplRtcmEphCorr);
327	}
328
329	co.UpdateInterval = clkUpdInd;
330	bias.UpdateInterval = clkUpdInd;
331	phasebias.UpdateInterval = clkUpdInd;
332
333	for (int ii = 1; ii < lines.size(); ii++) {
334	QString key; // prn or key VTEC, IND (phase bias indicators)
335	double rtnUra; // [m]
336	ColumnVector rtnAPC; // [m, m, m]
337	ColumnVector rtnVel; // [m/s, m/s, m/s]]
338	ColumnVector rtnCoM; // [m, m, m]
339	ColumnVector rtnClk; // [m, m/s, m/s²]
340	ColumnVector rtnClkSig; // [m, m/s, m/s²]
341	t_prn prn;
342
343	QTextStream in(lines[ii].toLatin1());
344
345	in >> key;
346
347	// non-satellite specific parameters
348	if (key.contains("IND", Qt::CaseSensitive)) {
349	in >> dispersiveBiasConsistenyIndicator >> mwConsistencyIndicator;
350	continue;
351	}
352	// non-satellite specific parameters
353	if (key.contains("VTEC", Qt::CaseSensitive)) {
354	double ui;
355	in >> ui >> vtec.NumLayers;
356	vtec.UpdateInterval = (unsigned int) determineUpdateInd(ui);
357	for (unsigned ll = 0; ll < vtec.NumLayers; ll++) {
358	int dummy;
359	in >> dummy >> vtec.Layers[ll].Degree >> vtec.Layers[ll].Order
360	>> vtec.Layers[ll].Height;
361	for (unsigned iDeg = 0; iDeg <= vtec.Layers[ll].Degree; iDeg++) {
362	for (unsigned iOrd = 0; iOrd <= vtec.Layers[ll].Order; iOrd++) {
363	in >> vtec.Layers[ll].Cosinus[iDeg][iOrd];
364	}
365	}
366	for (unsigned iDeg = 0; iDeg <= vtec.Layers[ll].Degree; iDeg++) {
367	for (unsigned iOrd = 0; iOrd <= vtec.Layers[ll].Order; iOrd++) {
368	in >> vtec.Layers[ll].Sinus[iDeg][iOrd];
369	}
370	}
371	}
372	continue;
373	}
374	// satellite specific parameters
375	char sys = key.mid(0, 1).at(0).toLatin1();
376	int number = key.mid(1, 2).toInt();
377	int flags = 0;
378	if (sys == 'E') { // I/NAV
379	flags = 1;
380	}
381	prn.set(sys, number, flags);
382	QString prnInternalStr = QString::fromStdString(prn.toInternalString());
383	QString prnStr = QString::fromStdString(prn.toString());
384
385	const t_eph* ephLast = _ephUser->ephLast(prnInternalStr);
386	const t_eph* ephPrev = _ephUser->ephPrev(prnInternalStr);
387	const t_eph* eph = ephLast;
388	if (eph) {
389
390	// Use previous ephemeris if the last one is too recent
391	// ----------------------------------------------------
392	const int MINAGE = 60; // seconds
393	if (ephPrev && eph->receptDateTime().isValid()
394	&& eph->receptDateTime().secsTo(currentDateAndTimeGPS()) < MINAGE) {
395	eph = ephPrev;
396	}
397
398	// Make sure the clock messages refer to same IOD as orbit messages
399	// ----------------------------------------------------------------
400	if (_usedEph) {
401	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
402	(*_usedEph)[prnInternalStr] = eph;
403	}
404	else {
405	eph = 0;
406	if (_usedEph->contains(prnInternalStr)) {
407	const t_eph* usedEph = _usedEph->value(prnInternalStr);
408	if (usedEph == ephLast) {
409	eph = ephLast;
410	}
411	else if (usedEph == ephPrev) {
412	eph = ephPrev;
413	}
414	}
415	}
416	}
417	}
418
419	if (eph) {
420
421	QMap<QString, double> codeBiases;
422	QList<phaseBiasSignal> phaseBiasList;
423	phaseBiasesSat pbSat;
424
425	while (true) {
426	QString key;
427	int numVal = 0;
428	in >> key;
429	if (in.status() != QTextStream::Ok) {
430	break;
431	}
432	if (key == "APC") {
433	in >> numVal;
434	rtnAPC.ReSize(3);
435	for (int ii = 0; ii < numVal; ii++) {
436	in >> rtnAPC[ii];
437	}
438	}
439	else if (key == "Ura") {
440	in >> numVal;
441	if (numVal == 1)
442	in >> rtnUra;
443	}
444	else if (key == "Clk") {
445	rtnClk.ReSize(3);
446	for (int ii = 0; ii < 3; ii++) {
447	rtnClk[ii] = 0.0;
448	}
449	in >> numVal;
450	for (int ii = 0; ii < numVal; ii++) {
451	in >> rtnClk[ii];
452	}
453	}
454	else if (key == "ClkSig") {
455	rtnClkSig.ReSize(3);
456	for (int ii = 0; ii < 3; ii++) {
457	rtnClkSig[ii] = 0.0;
458	}
459	in >> numVal;
460	for (int ii = 0; ii < numVal; ii++) {
461	in >> rtnClkSig[ii];
462	}
463	}
464	else if (key == "Vel") {
465	rtnVel.ReSize(3);
466	in >> numVal;
467	for (int ii = 0; ii < numVal; ii++) {
468	in >> rtnVel[ii];
469	}
470	}
471	else if (key == "CoM") {
472	rtnCoM.ReSize(3);
473	in >> numVal;
474	for (int ii = 0; ii < numVal; ii++) {
475	in >> rtnCoM[ii];
476	}
477	}
478	else if (key == "CodeBias") {
479	in >> numVal;
480	for (int ii = 0; ii < numVal; ii++) {
481	QString type;
482	double value;
483	in >> type >> value;
484	codeBiases[type] = value;
485	}
486	}
487	else if (key == "YawAngle") {
488	in >> numVal >> pbSat.yawAngle;
489	if (pbSat.yawAngle < 0.0) {
490	pbSat.yawAngle += (2*M_PI);
491	}
492	else if (pbSat.yawAngle > 2*M_PI) {
493	pbSat.yawAngle -= (2*M_PI);
494	}
495	}
496	else if (key == "YawRate") {
497	in >> numVal >> pbSat.yawRate;
498	}
499	else if (key == "PhaseBias") {
500	in >> numVal;
501	for (int ii = 0; ii < numVal; ii++) {
502	phaseBiasSignal pb;
503	in >> pb.type >> pb.bias >> pb.integerIndicator
504	>> pb.wlIndicator >> pb.discontinuityCounter;
505	phaseBiasList.append(pb);
506	}
507	}
508	else {
509	emit(newMessage(" RTNET format error: "
510	+ lines[ii].toLatin1(), false));
511	qDebug() << "bncRtnetUploadCaster: decode " + QByteArray(epoTime.datestr().c_str())
512	+ " " + QByteArray(epoTime.timestr().c_str()) + " "
513	+ _casterID.toLatin1();
514	qDebug() << " RTNET format error: " + lines[ii].toLatin1();
515	}
516	}
517
518
519	struct ClockOrbit::SatData* sd = 0;
520	if (prn.system() == 'G') {
521	sd = co.Sat + co.NumberOfSat[CLOCKORBIT_SATGPS];
522	++co.NumberOfSat[CLOCKORBIT_SATGPS];
523	}
524	else if (prn.system() == 'R') {
525	sd = co.Sat + CLOCKORBIT_NUMGPS + co.NumberOfSat[CLOCKORBIT_SATGLONASS];
526	++co.NumberOfSat[CLOCKORBIT_SATGLONASS];
527	}
528	else if (prn.system() == 'E') {
529	sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
530	+ co.NumberOfSat[CLOCKORBIT_SATGALILEO];
531	++co.NumberOfSat[CLOCKORBIT_SATGALILEO];
532	}
533	else if (prn.system() == 'J') {
534	sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
535	+ CLOCKORBIT_NUMGALILEO + co.NumberOfSat[CLOCKORBIT_SATQZSS];
536	++co.NumberOfSat[CLOCKORBIT_SATQZSS];
537	}
538	else if (prn.system() == 'S') {
539	sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
540	+ CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
541	+ co.NumberOfSat[CLOCKORBIT_SATSBAS];
542	++co.NumberOfSat[CLOCKORBIT_SATSBAS];
543	}
544	else if (prn.system() == 'C') {
545	sd = co.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
546	+ CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS + CLOCKORBIT_NUMSBAS
547	+ co.NumberOfSat[CLOCKORBIT_SATBDS];
548	++co.NumberOfSat[CLOCKORBIT_SATBDS];
549	}
550	if (sd) {
551	QString outLine;
552	t_irc irc = processSatellite(eph, epoTime.gpsw(), epoTime.gpssec(), prnStr, rtnAPC,
553	rtnUra, rtnClk, rtnVel, rtnCoM, rtnClkSig, sd, outLine);
554	if (irc != success) {
555	continue;
556	}
557	}
558
559	// Code Biases
560	// -----------
561	struct CodeBias::BiasSat* biasSat = 0;
562	if (!codeBiases.isEmpty()) {
563	if (prn.system() == 'G') {
564	biasSat = bias.Sat + bias.NumberOfSat[CLOCKORBIT_SATGPS];
565	++bias.NumberOfSat[CLOCKORBIT_SATGPS];
566	}
567	else if (prn.system() == 'R') {
568	biasSat = bias.Sat + CLOCKORBIT_NUMGPS
569	+ bias.NumberOfSat[CLOCKORBIT_SATGLONASS];
570	++bias.NumberOfSat[CLOCKORBIT_SATGLONASS];
571	}
572	else if (prn.system() == 'E') {
573	biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
574	+ bias.NumberOfSat[CLOCKORBIT_SATGALILEO];
575	++bias.NumberOfSat[CLOCKORBIT_SATGALILEO];
576	}
577	else if (prn.system() == 'J') {
578	biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
579	+ CLOCKORBIT_NUMGALILEO + bias.NumberOfSat[CLOCKORBIT_SATQZSS];
580	++bias.NumberOfSat[CLOCKORBIT_SATQZSS];
581	}
582	else if (prn.system() == 'S') {
583	biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
584	+ CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS
585	+ bias.NumberOfSat[CLOCKORBIT_SATSBAS];
586	++bias.NumberOfSat[CLOCKORBIT_SATSBAS];
587	}
588	else if (prn.system() == 'C') {
589	biasSat = bias.Sat + CLOCKORBIT_NUMGPS + CLOCKORBIT_NUMGLONASS
590	+ CLOCKORBIT_NUMGALILEO + CLOCKORBIT_NUMQZSS + CLOCKORBIT_NUMSBAS
591	+ bias.NumberOfSat[CLOCKORBIT_SATBDS];
592	++bias.NumberOfSat[CLOCKORBIT_SATBDS];
593	}
594	}
595
596	if (biasSat) {
597	biasSat->ID = prn.number();
598	biasSat->NumberOfCodeBiases = 0;
599	if (prn.system() == 'G') {
600	QMapIterator<QString, double> it(codeBiases);
601	while (it.hasNext()) {
602	it.next();
603	if (it.key() == "1C") {
604	int ii = biasSat->NumberOfCodeBiases;
605	if (ii >= CLOCKORBIT_NUMBIAS)
606	break;
607	biasSat->NumberOfCodeBiases += 1;
608	biasSat->Biases[ii].Type = CODETYPEGPS_L1_CA;
609	biasSat->Biases[ii].Bias = it.value();
610	}
611	else if (it.key() == "1S") {
612	int ii = biasSat->NumberOfCodeBiases;
613	if (ii >= CLOCKORBIT_NUMBIAS)
614	break;
615	biasSat->NumberOfCodeBiases += 1;
616	biasSat->Biases[ii].Type = CODETYPEGPS_L1C_D;
617	biasSat->Biases[ii].Bias = it.value();
618	}
619	else if (it.key() == "1L") {
620	int ii = biasSat->NumberOfCodeBiases;
621	if (ii >= CLOCKORBIT_NUMBIAS)
622	break;
623	biasSat->NumberOfCodeBiases += 1;
624	biasSat->Biases[ii].Type = CODETYPEGPS_L1C_P;
625	biasSat->Biases[ii].Bias = it.value();
626	}
627	else if (it.key() == "1X") {
628	int ii = biasSat->NumberOfCodeBiases;
629	if (ii >= CLOCKORBIT_NUMBIAS)
630	break;
631	biasSat->NumberOfCodeBiases += 1;
632	biasSat->Biases[ii].Type = CODETYPEGPS_L1C_DP;
633	biasSat->Biases[ii].Bias = it.value();
634	}
635	else if (it.key() == "1P") {
636	int ii = biasSat->NumberOfCodeBiases;
637	if (ii >= CLOCKORBIT_NUMBIAS)
638	break;
639	biasSat->NumberOfCodeBiases += 1;
640	biasSat->Biases[ii].Type = CODETYPEGPS_L1_P;
641	biasSat->Biases[ii].Bias = it.value();
642	}
643	else if (it.key() == "1W") {
644	int ii = biasSat->NumberOfCodeBiases;
645	if (ii >= CLOCKORBIT_NUMBIAS)
646	break;
647	biasSat->NumberOfCodeBiases += 1;
648	biasSat->Biases[ii].Type = CODETYPEGPS_L1_Z;
649	biasSat->Biases[ii].Bias = it.value();
650	}
651	else if (it.key() == "2C") {
652	int ii = biasSat->NumberOfCodeBiases;
653	if (ii >= CLOCKORBIT_NUMBIAS)
654	break;
655	biasSat->NumberOfCodeBiases += 1;
656	biasSat->Biases[ii].Type = CODETYPEGPS_L2_CA;
657	biasSat->Biases[ii].Bias = it.value();
658	}
659	else if (it.key() == "2D") {
660	int ii = biasSat->NumberOfCodeBiases;
661	if (ii >= CLOCKORBIT_NUMBIAS)
662	break;
663	biasSat->NumberOfCodeBiases += 1;
664	biasSat->Biases[ii].Type = CODETYPEGPS_SEMI_CODELESS;
665	biasSat->Biases[ii].Bias = it.value();
666	}
667	else if (it.key() == "2S") {
668	int ii = biasSat->NumberOfCodeBiases;
669	if (ii >= CLOCKORBIT_NUMBIAS)
670	break;
671	biasSat->NumberOfCodeBiases += 1;
672	biasSat->Biases[ii].Type = CODETYPEGPS_L2C_M;
673	biasSat->Biases[ii].Bias = it.value();
674	}
675	else if (it.key() == "2L") {
676	int ii = biasSat->NumberOfCodeBiases;
677	if (ii >= CLOCKORBIT_NUMBIAS)
678	break;
679	biasSat->NumberOfCodeBiases += 1;
680	biasSat->Biases[ii].Type = CODETYPEGPS_L2C_L;
681	biasSat->Biases[ii].Bias = it.value();
682	}
683	else if (it.key() == "2X") {
684	int ii = biasSat->NumberOfCodeBiases;
685	if (ii >= CLOCKORBIT_NUMBIAS)
686	break;
687	biasSat->NumberOfCodeBiases += 1;
688	biasSat->Biases[ii].Type = CODETYPEGPS_L2C_ML;
689	biasSat->Biases[ii].Bias = it.value();
690	}
691	else if (it.key() == "2P") {
692	int ii = biasSat->NumberOfCodeBiases;
693	if (ii >= CLOCKORBIT_NUMBIAS)
694	break;
695	biasSat->NumberOfCodeBiases += 1;
696	biasSat->Biases[ii].Type = CODETYPEGPS_L2_P;
697	biasSat->Biases[ii].Bias = it.value();
698	}
699	else if (it.key() == "2W") {
700	int ii = biasSat->NumberOfCodeBiases;
701	if (ii >= CLOCKORBIT_NUMBIAS)
702	break;
703	biasSat->NumberOfCodeBiases += 1;
704	biasSat->Biases[ii].Type = CODETYPEGPS_L2_Z;
705	biasSat->Biases[ii].Bias = it.value();
706	}
707	else if (it.key() == "5I") {
708	int ii = biasSat->NumberOfCodeBiases;
709	if (ii >= CLOCKORBIT_NUMBIAS)
710	break;
711	biasSat->NumberOfCodeBiases += 1;
712	biasSat->Biases[ii].Type = CODETYPEGPS_L5_I;
713	biasSat->Biases[ii].Bias = it.value();
714	}
715	else if (it.key() == "5Q") {
716	int ii = biasSat->NumberOfCodeBiases;
717	if (ii >= CLOCKORBIT_NUMBIAS)
718	break;
719	biasSat->NumberOfCodeBiases += 1;
720	biasSat->Biases[ii].Type = CODETYPEGPS_L5_Q;
721	biasSat->Biases[ii].Bias = it.value();
722	}
723	else if (it.key() == "5X") {
724	int ii = biasSat->NumberOfCodeBiases;
725	if (ii >= CLOCKORBIT_NUMBIAS)
726	break;
727	biasSat->NumberOfCodeBiases += 1;
728	biasSat->Biases[ii].Type = CODETYPEGPS_L5_IQ;
729	biasSat->Biases[ii].Bias = it.value();
730	}
731	}
732	}
733	else if (prn.system() == 'R') {
734	QMapIterator<QString, double> it(codeBiases);
735	while (it.hasNext()) {
736	it.next();
737	if (it.key() == "1C") {
738	int ii = biasSat->NumberOfCodeBiases;
739	if (ii >= CLOCKORBIT_NUMBIAS)
740	break;
741	biasSat->NumberOfCodeBiases += 1;
742	biasSat->Biases[ii].Type = CODETYPEGLONASS_L1_CA;
743	biasSat->Biases[ii].Bias = it.value();
744	}
745	else if (it.key() == "1P") {
746	int ii = biasSat->NumberOfCodeBiases;
747	if (ii >= CLOCKORBIT_NUMBIAS)
748	break;
749	biasSat->NumberOfCodeBiases += 1;
750	biasSat->Biases[ii].Type = CODETYPEGLONASS_L1_P;
751	biasSat->Biases[ii].Bias = it.value();
752	}
753	else if (it.key() == "2C") {
754	int ii = biasSat->NumberOfCodeBiases;
755	if (ii >= CLOCKORBIT_NUMBIAS)
756	break;
757	biasSat->NumberOfCodeBiases += 1;
758	biasSat->Biases[ii].Type = CODETYPEGLONASS_L2_CA;
759	biasSat->Biases[ii].Bias = it.value();
760	}
761	else if (it.key() == "2P") {
762	int ii = biasSat->NumberOfCodeBiases;
763	if (ii >= CLOCKORBIT_NUMBIAS)
764	break;
765	biasSat->NumberOfCodeBiases += 1;
766	biasSat->Biases[ii].Type = CODETYPEGLONASS_L2_P;
767	biasSat->Biases[ii].Bias = it.value();
768	}
769	}
770	}
771	else if (prn.system() == 'E') {
772	QMapIterator<QString, double> it(codeBiases);
773	while (it.hasNext()) {
774	it.next();
775	if (it.key() == "1A") {
776	int ii = biasSat->NumberOfCodeBiases;
777	if (ii >= CLOCKORBIT_NUMBIAS)
778	break;
779	biasSat->NumberOfCodeBiases += 1;
780	biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_A;
781	biasSat->Biases[ii].Bias = it.value();
782	}
783	else if (it.key() == "1B") {
784	int ii = biasSat->NumberOfCodeBiases;
785	if (ii >= CLOCKORBIT_NUMBIAS)
786	break;
787	biasSat->NumberOfCodeBiases += 1;
788	biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_B;
789	biasSat->Biases[ii].Bias = it.value();
790	}
791	else if (it.key() == "1C") {
792	int ii = biasSat->NumberOfCodeBiases;
793	if (ii >= CLOCKORBIT_NUMBIAS)
794	break;
795	biasSat->NumberOfCodeBiases += 1;
796	biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_C;
797	biasSat->Biases[ii].Bias = it.value();
798	}
799	else if (it.key() == "1X") {
800	int ii = biasSat->NumberOfCodeBiases;
801	if (ii >= CLOCKORBIT_NUMBIAS)
802	break;
803	biasSat->NumberOfCodeBiases += 1;
804	biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_BC;
805	biasSat->Biases[ii].Bias = it.value();
806	}
807	else if (it.key() == "1Z") {
808	int ii = biasSat->NumberOfCodeBiases;
809	if (ii >= CLOCKORBIT_NUMBIAS)
810	break;
811	biasSat->NumberOfCodeBiases += 1;
812	biasSat->Biases[ii].Type = CODETYPEGALILEO_E1_ABC;
813	biasSat->Biases[ii].Bias = it.value();
814	}
815	else if (it.key() == "5I") {
816	int ii = biasSat->NumberOfCodeBiases;
817	if (ii >= CLOCKORBIT_NUMBIAS)
818	break;
819	biasSat->NumberOfCodeBiases += 1;
820	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_I;
821	biasSat->Biases[ii].Bias = it.value();
822	}
823	else if (it.key() == "5Q") {
824	int ii = biasSat->NumberOfCodeBiases;
825	if (ii >= CLOCKORBIT_NUMBIAS)
826	break;
827	biasSat->NumberOfCodeBiases += 1;
828	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_Q;
829	biasSat->Biases[ii].Bias = it.value();
830	}
831	else if (it.key() == "5X") {
832	int ii = biasSat->NumberOfCodeBiases;
833	if (ii >= CLOCKORBIT_NUMBIAS)
834	break;
835	biasSat->NumberOfCodeBiases += 1;
836	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_IQ;
837	biasSat->Biases[ii].Bias = it.value();
838	}
839	else if (it.key() == "7I") {
840	int ii = biasSat->NumberOfCodeBiases;
841	if (ii >= CLOCKORBIT_NUMBIAS)
842	break;
843	biasSat->NumberOfCodeBiases += 1;
844	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_I;
845	biasSat->Biases[ii].Bias = it.value();
846	}
847	else if (it.key() == "7Q") {
848	int ii = biasSat->NumberOfCodeBiases;
849	if (ii >= CLOCKORBIT_NUMBIAS)
850	break;
851	biasSat->NumberOfCodeBiases += 1;
852	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_Q;
853	biasSat->Biases[ii].Bias = it.value();
854	}
855	else if (it.key() == "7X") {
856	int ii = biasSat->NumberOfCodeBiases;
857	if (ii >= CLOCKORBIT_NUMBIAS)
858	break;
859	biasSat->NumberOfCodeBiases += 1;
860	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_IQ;
861	biasSat->Biases[ii].Bias = it.value();
862	}
863	else if (it.key() == "8I") {
864	int ii = biasSat->NumberOfCodeBiases;
865	if (ii >= CLOCKORBIT_NUMBIAS)
866	break;
867	biasSat->NumberOfCodeBiases += 1;
868	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_I;
869	biasSat->Biases[ii].Bias = it.value();
870	}
871	else if (it.key() == "8Q") {
872	int ii = biasSat->NumberOfCodeBiases;
873	if (ii >= CLOCKORBIT_NUMBIAS)
874	break;
875	biasSat->NumberOfCodeBiases += 1;
876	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_Q;
877	biasSat->Biases[ii].Bias = it.value();
878	}
879	else if (it.key() == "8X") {
880	int ii = biasSat->NumberOfCodeBiases;
881	if (ii >= CLOCKORBIT_NUMBIAS)
882	break;
883	biasSat->NumberOfCodeBiases += 1;
884	biasSat->Biases[ii].Type = CODETYPEGALILEO_E5_IQ;
885	biasSat->Biases[ii].Bias = it.value();
886	}
887	else if (it.key() == "6A") {
888	int ii = biasSat->NumberOfCodeBiases;
889	if (ii >= CLOCKORBIT_NUMBIAS)
890	break;
891	biasSat->NumberOfCodeBiases += 1;
892	biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_A;
893	biasSat->Biases[ii].Bias = it.value();
894	}
895	else if (it.key() == "6B") {
896	int ii = biasSat->NumberOfCodeBiases;
897	if (ii >= CLOCKORBIT_NUMBIAS)
898	break;
899	biasSat->NumberOfCodeBiases += 1;
900	biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_B;
901	biasSat->Biases[ii].Bias = it.value();
902	}
903	else if (it.key() == "6C") {
904	int ii = biasSat->NumberOfCodeBiases;
905	if (ii >= CLOCKORBIT_NUMBIAS)
906	break;
907	biasSat->NumberOfCodeBiases += 1;
908	biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_C;
909	biasSat->Biases[ii].Bias = it.value();
910	}
911	else if (it.key() == "6X") {
912	int ii = biasSat->NumberOfCodeBiases;
913	if (ii >= CLOCKORBIT_NUMBIAS)
914	break;
915	biasSat->NumberOfCodeBiases += 1;
916	biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_BC;
917	biasSat->Biases[ii].Bias = it.value();
918	}
919	else if (it.key() == "6Z") {
920	int ii = biasSat->NumberOfCodeBiases;
921	if (ii >= CLOCKORBIT_NUMBIAS)
922	break;
923	biasSat->NumberOfCodeBiases += 1;
924	biasSat->Biases[ii].Type = CODETYPEGALILEO_E6_ABC;
925	biasSat->Biases[ii].Bias = it.value();
926	}
927	}
928	}
929	else if (prn.system() == 'J') {
930	QMapIterator<QString, double> it(codeBiases);
931	while (it.hasNext()) {
932	it.next();
933	if (it.key() == "1C") {
934	int ii = biasSat->NumberOfCodeBiases;
935	if (ii >= CLOCKORBIT_NUMBIAS)
936	break;
937	biasSat->NumberOfCodeBiases += 1;
938	biasSat->Biases[ii].Type = CODETYPEQZSS_L1_CA;
939	biasSat->Biases[ii].Bias = it.value();
940	}
941	else if (it.key() == "1S") {
942	int ii = biasSat->NumberOfCodeBiases;
943	if (ii >= CLOCKORBIT_NUMBIAS)
944	break;
945	biasSat->NumberOfCodeBiases += 1;
946	biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_D;
947	biasSat->Biases[ii].Bias = it.value();
948	}
949	else if (it.key() == "1L") {
950	int ii = biasSat->NumberOfCodeBiases;
951	if (ii >= CLOCKORBIT_NUMBIAS)
952	break;
953	biasSat->NumberOfCodeBiases += 1;
954	biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_P;
955	biasSat->Biases[ii].Bias = it.value();
956	}
957	else if (it.key() == "1X") {
958	int ii = biasSat->NumberOfCodeBiases;
959	if (ii >= CLOCKORBIT_NUMBIAS)
960	break;
961	biasSat->NumberOfCodeBiases += 1;
962	biasSat->Biases[ii].Type = CODETYPEQZSS_L1C_DP;
963	biasSat->Biases[ii].Bias = it.value();
964	}
965	else if (it.key() == "2S") {
966	int ii = biasSat->NumberOfCodeBiases;
967	if (ii >= CLOCKORBIT_NUMBIAS)
968	break;
969	biasSat->NumberOfCodeBiases += 1;
970	biasSat->Biases[ii].Type = CODETYPEQZSS_L2C_M;
971	biasSat->Biases[ii].Bias = it.value();
972	}
973	else if (it.key() == "2L") {
974	int ii = biasSat->NumberOfCodeBiases;
975	if (ii >= CLOCKORBIT_NUMBIAS)
976	break;
977	biasSat->NumberOfCodeBiases += 1;
978	biasSat->Biases[ii].Type = CODETYPEQZSS_L2C_L;
979	biasSat->Biases[ii].Bias = it.value();
980	}
981	else if (it.key() == "2X") {
982	int ii = biasSat->NumberOfCodeBiases;
983	if (ii >= CLOCKORBIT_NUMBIAS)
984	break;
985	biasSat->NumberOfCodeBiases += 1;
986	biasSat->Biases[ii].Type = CODETYPEQZSS_L2C_ML;
987	biasSat->Biases[ii].Bias = it.value();
988	}
989	else if (it.key() == "5I") {
990	int ii = biasSat->NumberOfCodeBiases;
991	if (ii >= CLOCKORBIT_NUMBIAS)
992	break;
993	biasSat->NumberOfCodeBiases += 1;
994	biasSat->Biases[ii].Type = CODETYPEQZSS_L5_I;
995	biasSat->Biases[ii].Bias = it.value();
996	}
997	else if (it.key() == "5Q") {
998	int ii = biasSat->NumberOfCodeBiases;
999	if (ii >= CLOCKORBIT_NUMBIAS)
1000	break;
1001	biasSat->NumberOfCodeBiases += 1;
1002	biasSat->Biases[ii].Type = CODETYPEQZSS_L5_Q;
1003	biasSat->Biases[ii].Bias = it.value();
1004	}
1005	else if (it.key() == "5X") {
1006	int ii = biasSat->NumberOfCodeBiases;
1007	if (ii >= CLOCKORBIT_NUMBIAS)
1008	break;
1009	biasSat->NumberOfCodeBiases += 1;
1010	biasSat->Biases[ii].Type = CODETYPEQZSS_L5_IQ;
1011	biasSat->Biases[ii].Bias = it.value();
1012	}
1013	else if (it.key() == "6S") {
1014	int ii = biasSat->NumberOfCodeBiases;
1015	if (ii >= CLOCKORBIT_NUMBIAS)
1016	break;
1017	biasSat->NumberOfCodeBiases += 1;
1018	biasSat->Biases[ii].Type = CODETYPEQZSS_L6_D;
1019	biasSat->Biases[ii].Bias = it.value();
1020	}
1021	else if (it.key() == "6L") {
1022	int ii = biasSat->NumberOfCodeBiases;
1023	if (ii >= CLOCKORBIT_NUMBIAS)
1024	break;
1025	biasSat->NumberOfCodeBiases += 1;
1026	biasSat->Biases[ii].Type = CODETYPEQZSS_L6_P;
1027	biasSat->Biases[ii].Bias = it.value();
1028	}
1029	else if (it.key() == "6X") {
1030	int ii = biasSat->NumberOfCodeBiases;
1031	if (ii >= CLOCKORBIT_NUMBIAS)
1032	break;
1033	biasSat->NumberOfCodeBiases += 1;
1034	biasSat->Biases[ii].Type = CODETYPEQZSS_L6_DP;
1035	biasSat->Biases[ii].Bias = it.value();
1036	}
1037	}
1038	}
1039	else if (prn.system() == 'S') {
1040	QMapIterator<QString, double> it(codeBiases);
1041	while (it.hasNext()) {
1042	it.next();
1043	if (it.key() == "1C") {
1044	int ii = biasSat->NumberOfCodeBiases;
1045	if (ii >= CLOCKORBIT_NUMBIAS)
1046	break;
1047	biasSat->NumberOfCodeBiases += 1;
1048	biasSat->Biases[ii].Type = CODETYPE_SBAS_L1_CA;
1049	biasSat->Biases[ii].Bias = it.value();
1050	}
1051	else if (it.key() == "5I") {
1052	int ii = biasSat->NumberOfCodeBiases;
1053	if (ii >= CLOCKORBIT_NUMBIAS)
1054	break;
1055	biasSat->NumberOfCodeBiases += 1;
1056	biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_I;
1057	biasSat->Biases[ii].Bias = it.value();
1058	}
1059	else if (it.key() == "5Q") {
1060	int ii = biasSat->NumberOfCodeBiases;
1061	if (ii >= CLOCKORBIT_NUMBIAS)
1062	break;
1063	biasSat->NumberOfCodeBiases += 1;
1064	biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_Q;
1065	biasSat->Biases[ii].Bias = it.value();
1066	}
1067	else if (it.key() == "5X") {
1068	int ii = biasSat->NumberOfCodeBiases;
1069	if (ii >= CLOCKORBIT_NUMBIAS)
1070	break;
1071	biasSat->NumberOfCodeBiases += 1;
1072	biasSat->Biases[ii].Type = CODETYPE_SBAS_L5_IQ;
1073	biasSat->Biases[ii].Bias = it.value();
1074	}
1075	}
1076	}
1077	else if (prn.system() == 'C') {
1078	QMapIterator<QString, double> it(codeBiases);
1079	while (it.hasNext()) {
1080	it.next();
1081	if (it.key() == "2I") {
1082	int ii = biasSat->NumberOfCodeBiases;
1083	if (ii >= CLOCKORBIT_NUMBIAS)
1084	break;
1085	biasSat->NumberOfCodeBiases += 1;
1086	biasSat->Biases[ii].Type = CODETYPE_BDS_B1_I;
1087	biasSat->Biases[ii].Bias = it.value();
1088	}
1089	else if (it.key() == "2Q") {
1090	int ii = biasSat->NumberOfCodeBiases;
1091	if (ii >= CLOCKORBIT_NUMBIAS)
1092	break;
1093	biasSat->NumberOfCodeBiases += 1;
1094	biasSat->Biases[ii].Type = CODETYPE_BDS_B1_Q;
1095	biasSat->Biases[ii].Bias = it.value();
1096	}
1097	else if (it.key() == "2X") {
1098	int ii = biasSat->NumberOfCodeBiases;
1099	if (ii >= CLOCKORBIT_NUMBIAS)
1100	break;
1101	biasSat->NumberOfCodeBiases += 1;
1102	biasSat->Biases[ii].Type = CODETYPE_BDS_B1_IQ;
1103	biasSat->Biases[ii].Bias = it.value();
1104	}
1105	else if (it.key() == "6I") {
1106	int ii = biasSat->NumberOfCodeBiases;
1107	if (ii >= CLOCKORBIT_NUMBIAS)
1108	break;
1109	biasSat->NumberOfCodeBiases += 1;
1110	biasSat->Biases[ii].Type = CODETYPE_BDS_B3_I;
1111	biasSat->Biases[ii].Bias = it.value();
1112	}
1113	else if (it.key() == "6Q") {
1114	int ii = biasSat->NumberOfCodeBiases;
1115	if (ii >= CLOCKORBIT_NUMBIAS)
1116	break;
1117	biasSat->NumberOfCodeBiases += 1;
1118	biasSat->Biases[ii].Type = CODETYPE_BDS_B3_Q;
1119	biasSat->Biases[ii].Bias = it.value();
1120	}
1121	else if (it.key() == "6X") {
1122	int ii = biasSat->NumberOfCodeBiases;
1123	if (ii >= CLOCKORBIT_NUMBIAS)
1124	break;
1125	biasSat->NumberOfCodeBiases += 1;
1126	biasSat->Biases[ii].Type = CODETYPE_BDS_B3_IQ;
1127	biasSat->Biases[ii].Bias = it.value();
1128	}
1129	else if (it.key() == "7I") {
1130	int ii = biasSat->NumberOfCodeBiases;
1131	if (ii >= CLOCKORBIT_NUMBIAS)
1132	break;
1133	biasSat->NumberOfCodeBiases += 1;
1134	biasSat->Biases[ii].Type = CODETYPE_BDS_B2_I;
1135	biasSat->Biases[ii].Bias = it.value();
1136	}
1137	else if (it.key() == "7Q") {
1138	int ii = biasSat->NumberOfCodeBiases;
1139	if (ii >= CLOCKORBIT_NUMBIAS)
1140	break;
1141	biasSat->NumberOfCodeBiases += 1;
1142	biasSat->Biases[ii].Type = CODETYPE_BDS_B2_Q;
1143	biasSat->Biases[ii].Bias = it.value();
1144	}
1145	else if (it.key() == "7X") {
1146	int ii = biasSat->NumberOfCodeBiases;
1147	if (ii >= CLOCKORBIT_NUMBIAS)
1148	break;
1149	biasSat->NumberOfCodeBiases += 1;
1150	biasSat->Biases[ii].Type = CODETYPE_BDS_B2_IQ;
1151	biasSat->Biases[ii].Bias = it.value();
1152	}
1153	}
1154	}
1155	}
1156	// Phase Biases
1157	// ------------
1158	struct PhaseBias::PhaseBiasSat* phasebiasSat = 0;
1159	if (!phaseBiasList.isEmpty()) {
1160	if (prn.system() == 'G') {
1161	phasebiasSat = phasebias.Sat
1162	+ phasebias.NumberOfSat[CLOCKORBIT_SATGPS];
1163	++phasebias.NumberOfSat[CLOCKORBIT_SATGPS];
1164	}
1165	else if (prn.system() == 'R') {
1166	phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
1167	+ phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS];
1168	++phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS];
1169	}
1170	else if (prn.system() == 'E') {
1171	phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
1172	+ CLOCKORBIT_NUMGLONASS
1173	+ phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO];
1174	++phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO];
1175	}
1176	else if (prn.system() == 'J') {
1177	phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
1178	+ CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
1179	+ phasebias.NumberOfSat[CLOCKORBIT_SATQZSS];
1180	++phasebias.NumberOfSat[CLOCKORBIT_SATQZSS];
1181	}
1182	else if (prn.system() == 'S') {
1183	phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
1184	+ CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
1185	+ CLOCKORBIT_NUMQZSS + phasebias.NumberOfSat[CLOCKORBIT_SATSBAS];
1186	++phasebias.NumberOfSat[CLOCKORBIT_SATSBAS];
1187	}
1188	else if (prn.system() == 'C') {
1189	phasebiasSat = phasebias.Sat + CLOCKORBIT_NUMGPS
1190	+ CLOCKORBIT_NUMGLONASS + CLOCKORBIT_NUMGALILEO
1191	+ CLOCKORBIT_NUMQZSS + CLOCKORBIT_NUMSBAS
1192	+ phasebias.NumberOfSat[CLOCKORBIT_SATBDS];
1193	++phasebias.NumberOfSat[CLOCKORBIT_SATBDS];
1194	}
1195	}
1196
1197	if (phasebiasSat) {
1198	phasebias.DispersiveBiasConsistencyIndicator = dispersiveBiasConsistenyIndicator;
1199	phasebias.MWConsistencyIndicator = mwConsistencyIndicator;
1200	phasebiasSat->ID = prn.number();
1201	phasebiasSat->NumberOfPhaseBiases = 0;
1202	phasebiasSat->YawAngle = pbSat.yawAngle;
1203	phasebiasSat->YawRate = pbSat.yawRate;
1204	if (prn.system() == 'G') {
1205	QListIterator<phaseBiasSignal> it(phaseBiasList);
1206	while (it.hasNext()) {
1207	const phaseBiasSignal &pbSig = it.next();
1208	if (pbSig.type == "1C") {
1209	int ii = phasebiasSat->NumberOfPhaseBiases;
1210	if (ii >= CLOCKORBIT_NUMBIAS)
1211	break;
1212	phasebiasSat->NumberOfPhaseBiases += 1;
1213	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_CA;
1214	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1215	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1216	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1217	pbSig.wlIndicator;
1218	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1219	pbSig.discontinuityCounter;
1220	}
1221	else if (pbSig.type == "1S") {
1222	int ii = phasebiasSat->NumberOfPhaseBiases;
1223	if (ii >= CLOCKORBIT_NUMBIAS)
1224	break;
1225	phasebiasSat->NumberOfPhaseBiases += 1;
1226	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1C_D;
1227	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1228	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1229	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1230	pbSig.wlIndicator;
1231	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1232	pbSig.discontinuityCounter;
1233	}
1234	else if (pbSig.type == "1L") {
1235	int ii = phasebiasSat->NumberOfPhaseBiases;
1236	if (ii >= CLOCKORBIT_NUMBIAS)
1237	break;
1238	phasebiasSat->NumberOfPhaseBiases += 1;
1239	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1C_P;
1240	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1241	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1242	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1243	pbSig.wlIndicator;
1244	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1245	pbSig.discontinuityCounter;
1246	}
1247	else if (pbSig.type == "1X") {
1248	int ii = phasebiasSat->NumberOfPhaseBiases;
1249	if (ii >= CLOCKORBIT_NUMBIAS)
1250	break;
1251	phasebiasSat->NumberOfPhaseBiases += 1;
1252	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1C_P;
1253	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1254	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1255	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1256	pbSig.wlIndicator;
1257	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1258	pbSig.discontinuityCounter;
1259	}
1260	else if (pbSig.type == "1P") {
1261	int ii = phasebiasSat->NumberOfPhaseBiases;
1262	if (ii >= CLOCKORBIT_NUMBIAS)
1263	break;
1264	phasebiasSat->NumberOfPhaseBiases += 1;
1265	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_P;
1266	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1267	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1268	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1269	pbSig.wlIndicator;
1270	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1271	pbSig.discontinuityCounter;
1272	}
1273	else if (pbSig.type == "1W") {
1274	int ii = phasebiasSat->NumberOfPhaseBiases;
1275	if (ii >= CLOCKORBIT_NUMBIAS)
1276	break;
1277	phasebiasSat->NumberOfPhaseBiases += 1;
1278	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L1_Z;
1279	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1280	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1281	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1282	pbSig.wlIndicator;
1283	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1284	pbSig.discontinuityCounter;
1285	}
1286	else if (pbSig.type == "2C") {
1287	int ii = phasebiasSat->NumberOfPhaseBiases;
1288	if (ii >= CLOCKORBIT_NUMBIAS)
1289	break;
1290	phasebiasSat->NumberOfPhaseBiases += 1;
1291	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_CA;
1292	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1293	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1294	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1295	pbSig.wlIndicator;
1296	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1297	pbSig.discontinuityCounter;
1298	}
1299	else if (pbSig.type == "2D") {
1300	int ii = phasebiasSat->NumberOfPhaseBiases;
1301	if (ii >= CLOCKORBIT_NUMBIAS)
1302	break;
1303	phasebiasSat->NumberOfPhaseBiases += 1;
1304	phasebiasSat->Biases[ii].Type = CODETYPEGPS_SEMI_CODELESS;
1305	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1306	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1307	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1308	pbSig.wlIndicator;
1309	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1310	pbSig.discontinuityCounter;
1311	}
1312	else if (pbSig.type == "2S") {
1313	int ii = phasebiasSat->NumberOfPhaseBiases;
1314	if (ii >= CLOCKORBIT_NUMBIAS)
1315	break;
1316	phasebiasSat->NumberOfPhaseBiases += 1;
1317	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2C_M;
1318	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1319	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1320	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1321	pbSig.wlIndicator;
1322	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1323	pbSig.discontinuityCounter;
1324	}
1325	else if (pbSig.type == "2L") {
1326	int ii = phasebiasSat->NumberOfPhaseBiases;
1327	if (ii >= CLOCKORBIT_NUMBIAS)
1328	break;
1329	phasebiasSat->NumberOfPhaseBiases += 1;
1330	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2C_L;
1331	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1332	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1333	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1334	pbSig.wlIndicator;
1335	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1336	pbSig.discontinuityCounter;
1337	}
1338	else if (pbSig.type == "2X") {
1339	int ii = phasebiasSat->NumberOfPhaseBiases;
1340	if (ii >= CLOCKORBIT_NUMBIAS)
1341	break;
1342	phasebiasSat->NumberOfPhaseBiases += 1;
1343	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2C_ML;
1344	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1345	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1346	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1347	pbSig.wlIndicator;
1348	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1349	pbSig.discontinuityCounter;
1350	}
1351	else if (pbSig.type == "2P") {
1352	int ii = phasebiasSat->NumberOfPhaseBiases;
1353	if (ii >= CLOCKORBIT_NUMBIAS)
1354	break;
1355	phasebiasSat->NumberOfPhaseBiases += 1;
1356	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_P;
1357	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1358	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1359	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1360	pbSig.wlIndicator;
1361	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1362	pbSig.discontinuityCounter;
1363	}
1364	else if (pbSig.type == "2W") {
1365	int ii = phasebiasSat->NumberOfPhaseBiases;
1366	if (ii >= CLOCKORBIT_NUMBIAS)
1367	break;
1368	phasebiasSat->NumberOfPhaseBiases += 1;
1369	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L2_Z;
1370	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1371	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1372	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1373	pbSig.wlIndicator;
1374	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1375	pbSig.discontinuityCounter;
1376	}
1377	else if (pbSig.type == "5I") {
1378	int ii = phasebiasSat->NumberOfPhaseBiases;
1379	if (ii >= CLOCKORBIT_NUMBIAS)
1380	break;
1381	phasebiasSat->NumberOfPhaseBiases += 1;
1382	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_I;
1383	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1384	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1385	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1386	pbSig.wlIndicator;
1387	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1388	pbSig.discontinuityCounter;
1389	}
1390	else if (pbSig.type == "5Q") {
1391	int ii = phasebiasSat->NumberOfPhaseBiases;
1392	if (ii >= CLOCKORBIT_NUMBIAS)
1393	break;
1394	phasebiasSat->NumberOfPhaseBiases += 1;
1395	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_Q;
1396	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1397	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1398	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1399	pbSig.wlIndicator;
1400	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1401	pbSig.discontinuityCounter;
1402	}
1403	else if (pbSig.type == "5X") {
1404	int ii = phasebiasSat->NumberOfPhaseBiases;
1405	if (ii >= CLOCKORBIT_NUMBIAS)
1406	break;
1407	phasebiasSat->NumberOfPhaseBiases += 1;
1408	phasebiasSat->Biases[ii].Type = CODETYPEGPS_L5_IQ;
1409	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1410	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1411	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1412	pbSig.wlIndicator;
1413	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1414	pbSig.discontinuityCounter;
1415	}
1416	}
1417	}
1418	if (prn.system() == 'R') {
1419	QListIterator<phaseBiasSignal> it(phaseBiasList);
1420	while (it.hasNext()) {
1421	const phaseBiasSignal &pbSig = it.next();
1422	if (pbSig.type == "1C") {
1423	int ii = phasebiasSat->NumberOfPhaseBiases;
1424	if (ii >= CLOCKORBIT_NUMBIAS)
1425	break;
1426	phasebiasSat->NumberOfPhaseBiases += 1;
1427	phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1_CA;
1428	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1429	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1430	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1431	pbSig.wlIndicator;
1432	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1433	pbSig.discontinuityCounter;
1434	}
1435	else if (pbSig.type == "1P") {
1436	int ii = phasebiasSat->NumberOfPhaseBiases;
1437	if (ii >= CLOCKORBIT_NUMBIAS)
1438	break;
1439	phasebiasSat->NumberOfPhaseBiases += 1;
1440	phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L1_P;
1441	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1442	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1443	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1444	pbSig.wlIndicator;
1445	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1446	pbSig.discontinuityCounter;
1447	}
1448	else if (pbSig.type == "2C") {
1449	int ii = phasebiasSat->NumberOfPhaseBiases;
1450	if (ii >= CLOCKORBIT_NUMBIAS)
1451	break;
1452	phasebiasSat->NumberOfPhaseBiases += 1;
1453	phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2_CA;
1454	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1455	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1456	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1457	pbSig.wlIndicator;
1458	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1459	pbSig.discontinuityCounter;
1460	}
1461	else if (pbSig.type == "2P") {
1462	int ii = phasebiasSat->NumberOfPhaseBiases;
1463	if (ii >= CLOCKORBIT_NUMBIAS)
1464	break;
1465	phasebiasSat->NumberOfPhaseBiases += 1;
1466	phasebiasSat->Biases[ii].Type = CODETYPEGLONASS_L2_P;
1467	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1468	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1469	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1470	pbSig.wlIndicator;
1471	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1472	pbSig.discontinuityCounter;
1473	}
1474	}
1475	}
1476	if (prn.system() == 'E') {
1477	QListIterator<phaseBiasSignal> it(phaseBiasList);
1478	while (it.hasNext()) {
1479	const phaseBiasSignal &pbSig = it.next();
1480	if (pbSig.type == "1A") {
1481	int ii = phasebiasSat->NumberOfPhaseBiases;
1482	if (ii >= CLOCKORBIT_NUMBIAS)
1483	break;
1484	phasebiasSat->NumberOfPhaseBiases += 1;
1485	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_A;
1486	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1487	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1488	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1489	pbSig.wlIndicator;
1490	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1491	pbSig.discontinuityCounter;
1492	}
1493	else if (pbSig.type == "1B") {
1494	int ii = phasebiasSat->NumberOfPhaseBiases;
1495	if (ii >= CLOCKORBIT_NUMBIAS)
1496	break;
1497	phasebiasSat->NumberOfPhaseBiases += 1;
1498	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_B;
1499	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1500	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1501	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1502	pbSig.wlIndicator;
1503	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1504	pbSig.discontinuityCounter;
1505	}
1506	else if (pbSig.type == "1C") {
1507	int ii = phasebiasSat->NumberOfPhaseBiases;
1508	if (ii >= CLOCKORBIT_NUMBIAS)
1509	break;
1510	phasebiasSat->NumberOfPhaseBiases += 1;
1511	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_C;
1512	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1513	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1514	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1515	pbSig.wlIndicator;
1516	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1517	pbSig.discontinuityCounter;
1518	}
1519	else if (pbSig.type == "1X") {
1520	int ii = phasebiasSat->NumberOfPhaseBiases;
1521	if (ii >= CLOCKORBIT_NUMBIAS)
1522	break;
1523	phasebiasSat->NumberOfPhaseBiases += 1;
1524	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_BC;
1525	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1526	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1527	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1528	pbSig.wlIndicator;
1529	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1530	pbSig.discontinuityCounter;
1531	}
1532	else if (pbSig.type == "1Z") {
1533	int ii = phasebiasSat->NumberOfPhaseBiases;
1534	if (ii >= CLOCKORBIT_NUMBIAS)
1535	break;
1536	phasebiasSat->NumberOfPhaseBiases += 1;
1537	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E1_ABC;
1538	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1539	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1540	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1541	pbSig.wlIndicator;
1542	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1543	pbSig.discontinuityCounter;
1544	}
1545	else if (pbSig.type == "5I") {
1546	int ii = phasebiasSat->NumberOfPhaseBiases;
1547	if (ii >= CLOCKORBIT_NUMBIAS)
1548	break;
1549	phasebiasSat->NumberOfPhaseBiases += 1;
1550	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_I;
1551	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1552	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1553	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1554	pbSig.wlIndicator;
1555	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1556	pbSig.discontinuityCounter;
1557	}
1558	else if (pbSig.type == "5Q") {
1559	int ii = phasebiasSat->NumberOfPhaseBiases;
1560	if (ii >= CLOCKORBIT_NUMBIAS)
1561	break;
1562	phasebiasSat->NumberOfPhaseBiases += 1;
1563	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_Q;
1564	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1565	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1566	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1567	pbSig.wlIndicator;
1568	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1569	pbSig.discontinuityCounter;
1570	}
1571	else if (pbSig.type == "5X") {
1572	int ii = phasebiasSat->NumberOfPhaseBiases;
1573	if (ii >= CLOCKORBIT_NUMBIAS)
1574	break;
1575	phasebiasSat->NumberOfPhaseBiases += 1;
1576	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5A_IQ;
1577	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1578	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1579	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1580	pbSig.wlIndicator;
1581	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1582	pbSig.discontinuityCounter;
1583	}
1584	else if (pbSig.type == "7I") {
1585	int ii = phasebiasSat->NumberOfPhaseBiases;
1586	if (ii >= CLOCKORBIT_NUMBIAS)
1587	break;
1588	phasebiasSat->NumberOfPhaseBiases += 1;
1589	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_I;
1590	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1591	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1592	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1593	pbSig.wlIndicator;
1594	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1595	pbSig.discontinuityCounter;
1596	}
1597	else if (pbSig.type == "7Q") {
1598	int ii = phasebiasSat->NumberOfPhaseBiases;
1599	if (ii >= CLOCKORBIT_NUMBIAS)
1600	break;
1601	phasebiasSat->NumberOfPhaseBiases += 1;
1602	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_Q;
1603	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1604	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1605	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1606	pbSig.wlIndicator;
1607	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1608	pbSig.discontinuityCounter;
1609	}
1610	else if (pbSig.type == "7X") {
1611	int ii = phasebiasSat->NumberOfPhaseBiases;
1612	if (ii >= CLOCKORBIT_NUMBIAS)
1613	break;
1614	phasebiasSat->NumberOfPhaseBiases += 1;
1615	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5B_IQ;
1616	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1617	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1618	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1619	pbSig.wlIndicator;
1620	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1621	pbSig.discontinuityCounter;
1622	}
1623	else if (pbSig.type == "8I") {
1624	int ii = phasebiasSat->NumberOfPhaseBiases;
1625	if (ii >= CLOCKORBIT_NUMBIAS)
1626	break;
1627	phasebiasSat->NumberOfPhaseBiases += 1;
1628	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_I;
1629	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1630	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1631	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1632	pbSig.wlIndicator;
1633	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1634	pbSig.discontinuityCounter;
1635	}
1636	else if (pbSig.type == "8Q") {
1637	int ii = phasebiasSat->NumberOfPhaseBiases;
1638	if (ii >= CLOCKORBIT_NUMBIAS)
1639	break;
1640	phasebiasSat->NumberOfPhaseBiases += 1;
1641	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_Q;
1642	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1643	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1644	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1645	pbSig.wlIndicator;
1646	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1647	pbSig.discontinuityCounter;
1648	}
1649	else if (pbSig.type == "8X") {
1650	int ii = phasebiasSat->NumberOfPhaseBiases;
1651	if (ii >= CLOCKORBIT_NUMBIAS)
1652	break;
1653	phasebiasSat->NumberOfPhaseBiases += 1;
1654	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E5_IQ;
1655	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1656	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1657	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1658	pbSig.wlIndicator;
1659	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1660	pbSig.discontinuityCounter;
1661	}
1662	else if (pbSig.type == "6A") {
1663	int ii = phasebiasSat->NumberOfPhaseBiases;
1664	if (ii >= CLOCKORBIT_NUMBIAS)
1665	break;
1666	phasebiasSat->NumberOfPhaseBiases += 1;
1667	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_A;
1668	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1669	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1670	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1671	pbSig.wlIndicator;
1672	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1673	pbSig.discontinuityCounter;
1674	}
1675	else if (pbSig.type == "6B") {
1676	int ii = phasebiasSat->NumberOfPhaseBiases;
1677	if (ii >= CLOCKORBIT_NUMBIAS)
1678	break;
1679	phasebiasSat->NumberOfPhaseBiases += 1;
1680	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_B;
1681	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1682	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1683	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1684	pbSig.wlIndicator;
1685	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1686	pbSig.discontinuityCounter;
1687	}
1688	else if (pbSig.type == "6C") {
1689	int ii = phasebiasSat->NumberOfPhaseBiases;
1690	if (ii >= CLOCKORBIT_NUMBIAS)
1691	break;
1692	phasebiasSat->NumberOfPhaseBiases += 1;
1693	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_C;
1694	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1695	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1696	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1697	pbSig.wlIndicator;
1698	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1699	pbSig.discontinuityCounter;
1700	}
1701	else if (pbSig.type == "6X") {
1702	int ii = phasebiasSat->NumberOfPhaseBiases;
1703	if (ii >= CLOCKORBIT_NUMBIAS)
1704	break;
1705	phasebiasSat->NumberOfPhaseBiases += 1;
1706	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_BC;
1707	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1708	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1709	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1710	pbSig.wlIndicator;
1711	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1712	pbSig.discontinuityCounter;
1713	}
1714	else if (pbSig.type == "6Z") {
1715	int ii = phasebiasSat->NumberOfPhaseBiases;
1716	if (ii >= CLOCKORBIT_NUMBIAS)
1717	break;
1718	phasebiasSat->NumberOfPhaseBiases += 1;
1719	phasebiasSat->Biases[ii].Type = CODETYPEGALILEO_E6_ABC;
1720	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1721	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1722	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1723	pbSig.wlIndicator;
1724	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1725	pbSig.discontinuityCounter;
1726	}
1727	}
1728	}
1729	if (prn.system() == 'J') {
1730	QListIterator<phaseBiasSignal> it(phaseBiasList);
1731	while (it.hasNext()) {
1732	const phaseBiasSignal &pbSig = it.next();
1733	if (pbSig.type == "1C") {
1734	int ii = phasebiasSat->NumberOfPhaseBiases;
1735	if (ii >= CLOCKORBIT_NUMBIAS)
1736	break;
1737	phasebiasSat->NumberOfPhaseBiases += 1;
1738	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1_CA;
1739	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1740	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1741	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1742	pbSig.wlIndicator;
1743	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1744	pbSig.discontinuityCounter;
1745	}
1746	else if (pbSig.type == "1S") {
1747	int ii = phasebiasSat->NumberOfPhaseBiases;
1748	if (ii >= CLOCKORBIT_NUMBIAS)
1749	break;
1750	phasebiasSat->NumberOfPhaseBiases += 1;
1751	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_D;
1752	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1753	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1754	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1755	pbSig.wlIndicator;
1756	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1757	pbSig.discontinuityCounter;
1758	}
1759	else if (pbSig.type == "1L") {
1760	int ii = phasebiasSat->NumberOfPhaseBiases;
1761	if (ii >= CLOCKORBIT_NUMBIAS)
1762	break;
1763	phasebiasSat->NumberOfPhaseBiases += 1;
1764	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_P;
1765	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1766	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1767	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1768	pbSig.wlIndicator;
1769	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1770	pbSig.discontinuityCounter;
1771	}
1772	else if (pbSig.type == "1X") {
1773	int ii = phasebiasSat->NumberOfPhaseBiases;
1774	if (ii >= CLOCKORBIT_NUMBIAS)
1775	break;
1776	phasebiasSat->NumberOfPhaseBiases += 1;
1777	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L1C_DP;
1778	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1779	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1780	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1781	pbSig.wlIndicator;
1782	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1783	pbSig.discontinuityCounter;
1784	}
1785	else if (pbSig.type == "2S") {
1786	int ii = phasebiasSat->NumberOfPhaseBiases;
1787	if (ii >= CLOCKORBIT_NUMBIAS)
1788	break;
1789	phasebiasSat->NumberOfPhaseBiases += 1;
1790	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2C_M;
1791	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1792	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1793	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1794	pbSig.wlIndicator;
1795	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1796	pbSig.discontinuityCounter;
1797	}
1798	else if (pbSig.type == "2L") {
1799	int ii = phasebiasSat->NumberOfPhaseBiases;
1800	if (ii >= CLOCKORBIT_NUMBIAS)
1801	break;
1802	phasebiasSat->NumberOfPhaseBiases += 1;
1803	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2C_L;
1804	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1805	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1806	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1807	pbSig.wlIndicator;
1808	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1809	pbSig.discontinuityCounter;
1810	}
1811	else if (pbSig.type == "2X") {
1812	int ii = phasebiasSat->NumberOfPhaseBiases;
1813	if (ii >= CLOCKORBIT_NUMBIAS)
1814	break;
1815	phasebiasSat->NumberOfPhaseBiases += 1;
1816	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L2C_ML;
1817	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1818	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1819	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1820	pbSig.wlIndicator;
1821	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1822	pbSig.discontinuityCounter;
1823	}
1824	else if (pbSig.type == "5I") {
1825	int ii = phasebiasSat->NumberOfPhaseBiases;
1826	if (ii >= CLOCKORBIT_NUMBIAS)
1827	break;
1828	phasebiasSat->NumberOfPhaseBiases += 1;
1829	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_I;
1830	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1831	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1832	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1833	pbSig.wlIndicator;
1834	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1835	pbSig.discontinuityCounter;
1836	}
1837	else if (pbSig.type == "5Q") {
1838	int ii = phasebiasSat->NumberOfPhaseBiases;
1839	if (ii >= CLOCKORBIT_NUMBIAS)
1840	break;
1841	phasebiasSat->NumberOfPhaseBiases += 1;
1842	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_Q;
1843	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1844	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1845	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1846	pbSig.wlIndicator;
1847	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1848	pbSig.discontinuityCounter;
1849	}
1850	else if (pbSig.type == "5X") {
1851	int ii = phasebiasSat->NumberOfPhaseBiases;
1852	if (ii >= CLOCKORBIT_NUMBIAS)
1853	break;
1854	phasebiasSat->NumberOfPhaseBiases += 1;
1855	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L5_IQ;
1856	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1857	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1858	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1859	pbSig.wlIndicator;
1860	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1861	pbSig.discontinuityCounter;
1862	}
1863	else if (pbSig.type == "6S") {
1864	int ii = phasebiasSat->NumberOfPhaseBiases;
1865	if (ii >= CLOCKORBIT_NUMBIAS)
1866	break;
1867	phasebiasSat->NumberOfPhaseBiases += 1;
1868	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_D;
1869	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1870	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1871	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1872	pbSig.wlIndicator;
1873	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1874	pbSig.discontinuityCounter;
1875	}
1876	else if (pbSig.type == "6L") {
1877	int ii = phasebiasSat->NumberOfPhaseBiases;
1878	if (ii >= CLOCKORBIT_NUMBIAS)
1879	break;
1880	phasebiasSat->NumberOfPhaseBiases += 1;
1881	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_P;
1882	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1883	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1884	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1885	pbSig.wlIndicator;
1886	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1887	pbSig.discontinuityCounter;
1888	}
1889	else if (pbSig.type == "6X") {
1890	int ii = phasebiasSat->NumberOfPhaseBiases;
1891	if (ii >= CLOCKORBIT_NUMBIAS)
1892	break;
1893	phasebiasSat->NumberOfPhaseBiases += 1;
1894	phasebiasSat->Biases[ii].Type = CODETYPEQZSS_L6_DP;
1895	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1896	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1897	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1898	pbSig.wlIndicator;
1899	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1900	pbSig.discontinuityCounter;
1901	}
1902	}
1903	}
1904	if (prn.system() == 'S') {
1905	QListIterator<phaseBiasSignal> it(phaseBiasList);
1906	while (it.hasNext()) {
1907	const phaseBiasSignal &pbSig = it.next();
1908	if (pbSig.type == "1C") {
1909	int ii = phasebiasSat->NumberOfPhaseBiases;
1910	if (ii >= CLOCKORBIT_NUMBIAS)
1911	break;
1912	phasebiasSat->NumberOfPhaseBiases += 1;
1913	phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L1_CA;
1914	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1915	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1916	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1917	pbSig.wlIndicator;
1918	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1919	pbSig.discontinuityCounter;
1920	}
1921	else if (pbSig.type == "5I") {
1922	int ii = phasebiasSat->NumberOfPhaseBiases;
1923	if (ii >= CLOCKORBIT_NUMBIAS)
1924	break;
1925	phasebiasSat->NumberOfPhaseBiases += 1;
1926	phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_I;
1927	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1928	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1929	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1930	pbSig.wlIndicator;
1931	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1932	pbSig.discontinuityCounter;
1933	}
1934	else if (pbSig.type == "5Q") {
1935	int ii = phasebiasSat->NumberOfPhaseBiases;
1936	if (ii >= CLOCKORBIT_NUMBIAS)
1937	break;
1938	phasebiasSat->NumberOfPhaseBiases += 1;
1939	phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_Q;
1940	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1941	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1942	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1943	pbSig.wlIndicator;
1944	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1945	pbSig.discontinuityCounter;
1946	}
1947	else if (pbSig.type == "5X") {
1948	int ii = phasebiasSat->NumberOfPhaseBiases;
1949	if (ii >= CLOCKORBIT_NUMBIAS)
1950	break;
1951	phasebiasSat->NumberOfPhaseBiases += 1;
1952	phasebiasSat->Biases[ii].Type = CODETYPE_SBAS_L5_IQ;
1953	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1954	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1955	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1956	pbSig.wlIndicator;
1957	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1958	pbSig.discontinuityCounter;
1959	}
1960	}
1961	}
1962	if (prn.system() == 'C') {
1963	QListIterator<phaseBiasSignal> it(phaseBiasList);
1964	while (it.hasNext()) {
1965	const phaseBiasSignal &pbSig = it.next();
1966	if (pbSig.type == "2I") {
1967	int ii = phasebiasSat->NumberOfPhaseBiases;
1968	if (ii >= CLOCKORBIT_NUMBIAS)
1969	break;
1970	phasebiasSat->NumberOfPhaseBiases += 1;
1971	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_I;
1972	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1973	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1974	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1975	pbSig.wlIndicator;
1976	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1977	pbSig.discontinuityCounter;
1978	}
1979	else if (pbSig.type == "2Q") {
1980	int ii = phasebiasSat->NumberOfPhaseBiases;
1981	if (ii >= CLOCKORBIT_NUMBIAS)
1982	break;
1983	phasebiasSat->NumberOfPhaseBiases += 1;
1984	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_Q;
1985	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1986	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
1987	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
1988	pbSig.wlIndicator;
1989	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
1990	pbSig.discontinuityCounter;
1991	}
1992	else if (pbSig.type == "2X") {
1993	int ii = phasebiasSat->NumberOfPhaseBiases;
1994	if (ii >= CLOCKORBIT_NUMBIAS)
1995	break;
1996	phasebiasSat->NumberOfPhaseBiases += 1;
1997	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B1_IQ;
1998	phasebiasSat->Biases[ii].Bias = pbSig.bias;
1999	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2000	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2001	pbSig.wlIndicator;
2002	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2003	pbSig.discontinuityCounter;
2004	}
2005	else if (pbSig.type == "6I") {
2006	int ii = phasebiasSat->NumberOfPhaseBiases;
2007	if (ii >= CLOCKORBIT_NUMBIAS)
2008	break;
2009	phasebiasSat->NumberOfPhaseBiases += 1;
2010	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_I;
2011	phasebiasSat->Biases[ii].Bias = pbSig.bias;
2012	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2013	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2014	pbSig.wlIndicator;
2015	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2016	pbSig.discontinuityCounter;
2017	}
2018	else if (pbSig.type == "6Q") {
2019	int ii = phasebiasSat->NumberOfPhaseBiases;
2020	if (ii >= CLOCKORBIT_NUMBIAS)
2021	break;
2022	phasebiasSat->NumberOfPhaseBiases += 1;
2023	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_Q;
2024	phasebiasSat->Biases[ii].Bias = pbSig.bias;
2025	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2026	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2027	pbSig.wlIndicator;
2028	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2029	pbSig.discontinuityCounter;
2030	}
2031	else if (pbSig.type == "6X") {
2032	int ii = phasebiasSat->NumberOfPhaseBiases;
2033	if (ii >= CLOCKORBIT_NUMBIAS)
2034	break;
2035	phasebiasSat->NumberOfPhaseBiases += 1;
2036	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B3_IQ;
2037	phasebiasSat->Biases[ii].Bias = pbSig.bias;
2038	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2039	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2040	pbSig.wlIndicator;
2041	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2042	pbSig.discontinuityCounter;
2043	}
2044	else if (pbSig.type == "7I") {
2045	int ii = phasebiasSat->NumberOfPhaseBiases;
2046	if (ii >= CLOCKORBIT_NUMBIAS)
2047	break;
2048	phasebiasSat->NumberOfPhaseBiases += 1;
2049	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_I;
2050	phasebiasSat->Biases[ii].Bias = pbSig.bias;
2051	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2052	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2053	pbSig.wlIndicator;
2054	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2055	pbSig.discontinuityCounter;
2056	}
2057	else if (pbSig.type == "7Q") {
2058	int ii = phasebiasSat->NumberOfPhaseBiases;
2059	if (ii >= CLOCKORBIT_NUMBIAS)
2060	break;
2061	phasebiasSat->NumberOfPhaseBiases += 1;
2062	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_Q;
2063	phasebiasSat->Biases[ii].Bias = pbSig.bias;
2064	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2065	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2066	pbSig.wlIndicator;
2067	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2068	pbSig.discontinuityCounter;
2069	}
2070	else if (pbSig.type == "7X") {
2071	int ii = phasebiasSat->NumberOfPhaseBiases;
2072	if (ii >= CLOCKORBIT_NUMBIAS)
2073	break;
2074	phasebiasSat->NumberOfPhaseBiases += 1;
2075	phasebiasSat->Biases[ii].Type = CODETYPE_BDS_B2_IQ;
2076	phasebiasSat->Biases[ii].Bias = pbSig.bias;
2077	phasebiasSat->Biases[ii].SignalIntegerIndicator = pbSig.integerIndicator;
2078	phasebiasSat->Biases[ii].SignalsWideLaneIntegerIndicator =
2079	pbSig.wlIndicator;
2080	phasebiasSat->Biases[ii].SignalDiscontinuityCounter =
2081	pbSig.discontinuityCounter;
2082	}
2083	}
2084	}
2085	}
2086	}
2087	}
2088
2089	QByteArray hlpBufferCo;
2090
2091	// Orbit and Clock Corrections together
2092	// ------------------------------------
2093	if (_samplRtcmEphCorr == 0.0) {
2094	if (co.NumberOfSat[CLOCKORBIT_SATGPS] > 0
2095	\|\| co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0
2096	\|\| co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0
2097	\|\| co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0
2098	\|\| co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0
2099	\|\| co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
2100	char obuffer[CLOCKORBIT_BUFFERSIZE];
2101	int len = MakeClockOrbit(&co, COTYPE_AUTO, 0, obuffer, sizeof(obuffer));
2102	if (len > 0) {
2103	hlpBufferCo = QByteArray(obuffer, len);
2104	}
2105	}
2106	}
2107
2108	// Orbit and Clock Corrections separately
2109	// --------------------------------------
2110	else {
2111	if (co.NumberOfSat[CLOCKORBIT_SATGPS] > 0) {
2112	char obuffer[CLOCKORBIT_BUFFERSIZE];
2113	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2114	co.UpdateInterval = ephUpdInd;
2115	int len1 = MakeClockOrbit(&co, COTYPE_GPSORBIT, 1, obuffer,
2116	sizeof(obuffer));
2117	co.UpdateInterval = clkUpdInd;
2118	if (len1 > 0) {
2119	hlpBufferCo += QByteArray(obuffer, len1);
2120	}
2121	}
2122	int mmsg = (co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0) ? 1 : 0;
2123	int len2 = MakeClockOrbit(&co, COTYPE_GPSCLOCK, mmsg, obuffer,
2124	sizeof(obuffer));
2125	if (len2 > 0) {
2126	hlpBufferCo += QByteArray(obuffer, len2);
2127	}
2128	}
2129	if (co.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0) {
2130	char obuffer[CLOCKORBIT_BUFFERSIZE];
2131	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2132	co.UpdateInterval = ephUpdInd;
2133	int len1 = MakeClockOrbit(&co, COTYPE_GLONASSORBIT, 1, obuffer,
2134	sizeof(obuffer));
2135	co.UpdateInterval = clkUpdInd;
2136	if (len1 > 0) {
2137	hlpBufferCo += QByteArray(obuffer, len1);
2138	}
2139	}
2140	int mmsg = (co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0) ? 1 : 0;
2141	int len2 = MakeClockOrbit(&co, COTYPE_GLONASSCLOCK, mmsg, obuffer,
2142	sizeof(obuffer));
2143	if (len2 > 0) {
2144	hlpBufferCo += QByteArray(obuffer, len2);
2145	}
2146	}
2147	if (co.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0) {
2148	char obuffer[CLOCKORBIT_BUFFERSIZE];
2149	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2150	co.UpdateInterval = ephUpdInd;
2151	int len1 = MakeClockOrbit(&co, COTYPE_GALILEOORBIT, 1, obuffer,
2152	sizeof(obuffer));
2153	co.UpdateInterval = clkUpdInd;
2154	if (len1 > 0) {
2155	hlpBufferCo += QByteArray(obuffer, len1);
2156	}
2157	}
2158	int mmsg = (co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0) ? 1 : 0;
2159	int len2 = MakeClockOrbit(&co, COTYPE_GALILEOCLOCK, mmsg, obuffer,
2160	sizeof(obuffer));
2161	if (len2 > 0) {
2162	hlpBufferCo += QByteArray(obuffer, len2);
2163	}
2164	}
2165	if (co.NumberOfSat[CLOCKORBIT_SATQZSS] > 0) {
2166	char obuffer[CLOCKORBIT_BUFFERSIZE];
2167	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2168	co.UpdateInterval = ephUpdInd;
2169	int len1 = MakeClockOrbit(&co, COTYPE_QZSSORBIT, 1, obuffer,
2170	sizeof(obuffer));
2171	co.UpdateInterval = clkUpdInd;
2172	if (len1 > 0) {
2173	hlpBufferCo += QByteArray(obuffer, len1);
2174	}
2175	}
2176	int mmsg = (co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0) ? 1 : 0;
2177	int len2 = MakeClockOrbit(&co, COTYPE_QZSSCLOCK, mmsg, obuffer,
2178	sizeof(obuffer));
2179	if (len2 > 0) {
2180	hlpBufferCo += QByteArray(obuffer, len2);
2181	}
2182	}
2183	if (co.NumberOfSat[CLOCKORBIT_SATSBAS] > 0) {
2184	char obuffer[CLOCKORBIT_BUFFERSIZE];
2185	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2186	co.UpdateInterval = ephUpdInd;
2187	int len1 = MakeClockOrbit(&co, COTYPE_SBASORBIT, 1, obuffer,
2188	sizeof(obuffer));
2189	co.UpdateInterval = clkUpdInd;
2190	if (len1 > 0) {
2191	hlpBufferCo += QByteArray(obuffer, len1);
2192	}
2193	}
2194	int mmsg = (co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) ? 1 : 0;
2195	int len2 = MakeClockOrbit(&co, COTYPE_SBASCLOCK, mmsg, obuffer,
2196	sizeof(obuffer));
2197	if (len2 > 0) {
2198	hlpBufferCo += QByteArray(obuffer, len2);
2199	}
2200	}
2201	if (co.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
2202	char obuffer[CLOCKORBIT_BUFFERSIZE];
2203	if (fmod(epoTime.gpssec(), _samplRtcmEphCorr) == 0.0) {
2204	co.UpdateInterval = ephUpdInd;
2205	int len1 = MakeClockOrbit(&co, COTYPE_BDSORBIT, 1, obuffer,
2206	sizeof(obuffer));
2207	co.UpdateInterval = clkUpdInd;
2208	if (len1 > 0) {
2209	hlpBufferCo += QByteArray(obuffer, len1);
2210	}
2211	}
2212	int len2 = MakeClockOrbit(&co, COTYPE_BDSCLOCK, 0, obuffer,
2213	sizeof(obuffer));
2214	if (len2 > 0) {
2215	hlpBufferCo += QByteArray(obuffer, len2);
2216	}
2217	}
2218	}
2219
2220	// Code Biases
2221	// -----------
2222	QByteArray hlpBufferBias;
2223	if (bias.NumberOfSat[CLOCKORBIT_SATGPS] > 0
2224	\|\| bias.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0
2225	\|\| bias.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0
2226	\|\| bias.NumberOfSat[CLOCKORBIT_SATQZSS] > 0
2227	\|\| bias.NumberOfSat[CLOCKORBIT_SATSBAS] > 0
2228	\|\| bias.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
2229	char obuffer[CLOCKORBIT_BUFFERSIZE];
2230	int len = MakeCodeBias(&bias, BTYPE_AUTO, 0, obuffer, sizeof(obuffer));
2231	if (len > 0) {
2232	hlpBufferBias = QByteArray(obuffer, len);
2233	}
2234	}
2235
2236	// Phase Biases
2237	// ------------
2238	QByteArray hlpBufferPhaseBias;
2239	if (phasebias.NumberOfSat[CLOCKORBIT_SATGPS] > 0
2240	\|\| phasebias.NumberOfSat[CLOCKORBIT_SATGLONASS] > 0
2241	\|\| phasebias.NumberOfSat[CLOCKORBIT_SATGALILEO] > 0
2242	\|\| phasebias.NumberOfSat[CLOCKORBIT_SATQZSS] > 0
2243	\|\| phasebias.NumberOfSat[CLOCKORBIT_SATSBAS] > 0
2244	\|\| phasebias.NumberOfSat[CLOCKORBIT_SATBDS] > 0) {
2245	char obuffer[CLOCKORBIT_BUFFERSIZE];
2246	int len = MakePhaseBias(&phasebias, PBTYPE_AUTO, 0, obuffer,
2247	sizeof(obuffer));
2248	if (len > 0) {
2249	hlpBufferPhaseBias = QByteArray(obuffer, len);
2250	}
2251	}
2252
2253	// VTEC
2254	// ----
2255	QByteArray hlpBufferVtec;
2256	if (vtec.NumLayers > 0) {
2257	char obuffer[CLOCKORBIT_BUFFERSIZE];
2258	int len = MakeVTEC(&vtec, 0, obuffer, sizeof(obuffer));
2259	if (len > 0) {
2260	hlpBufferVtec = QByteArray(obuffer, len);
2261	}
2262	}
2263
2264	_outBuffer += hlpBufferCo + hlpBufferBias + hlpBufferPhaseBias
2265	+ hlpBufferVtec;
2266	}
2267
2268	//
2269	////////////////////////////////////////////////////////////////////////////
2270	t_irc bncRtnetUploadCaster::processSatellite(const t_eph* eph, int GPSweek,
2271	double GPSweeks, const QString& prn, const ColumnVector& rtnAPC,
2272	double rtnUra, const ColumnVector& rtnClk, const ColumnVector& rtnVel,
2273	const ColumnVector& rtnCoM, const ColumnVector& rtnClkSig,
2274	struct ClockOrbit::SatData* sd, QString& outLine) {
2275
2276	// Broadcast Position and Velocity
2277	// -------------------------------
2278	ColumnVector xB(6);
2279	ColumnVector vB(3);
2280	t_irc irc = eph->getCrd(bncTime(GPSweek, GPSweeks), xB, vB, false);
2281
2282	if (irc != success) {
2283	return irc;
2284	}
2285
2286	// Precise Position
2287	// ----------------
2288	ColumnVector xP = _CoM ? rtnCoM : rtnAPC;
2289
2290	if (xP.size() == 0) {
2291	return failure;
2292	}
2293
2294	double dc = 0.0;
2295	if (_crdTrafo != "IGS14") {
2296	crdTrafo(GPSweek, xP, dc);
2297	}
2298
2299	// Difference in xyz
2300	// -----------------
2301	ColumnVector dx = xB.Rows(1, 3) - xP;
2302	ColumnVector dv = vB - rtnVel;
2303
2304	// Difference in RSW
2305	// -----------------
2306	ColumnVector rsw(3);
2307	XYZ_to_RSW(xB.Rows(1, 3), vB, dx, rsw);
2308
2309	ColumnVector dotRsw(3);
2310	XYZ_to_RSW(xB.Rows(1, 3), vB, dv, dotRsw);
2311
2312	// Clock Correction
2313	// ----------------
2314	double dClkA0 = rtnClk(1) - (xB(4) - dc) * t_CST::c;
2315	double dClkA1 = rtnClk(2) - xB(5) * t_CST::c;
2316	double dClkA2 = rtnClk(3) - xB(6) * t_CST::c;
2317
2318	if (sd) {
2319	sd->ID = prn.mid(1).toInt();
2320	sd->IOD = eph->IOD();
2321	sd->Clock.DeltaA0 = dClkA0;
2322	sd->Clock.DeltaA1 = dClkA1;
2323	sd->Clock.DeltaA2 = dClkA2;
2324	sd->UserRangeAccuracy = rtnUra;
2325	sd->Orbit.DeltaRadial = rsw(1);
2326	sd->Orbit.DeltaAlongTrack = rsw(2);
2327	sd->Orbit.DeltaCrossTrack = rsw(3);
2328	sd->Orbit.DotDeltaRadial = dotRsw(1);
2329	sd->Orbit.DotDeltaAlongTrack = dotRsw(2);
2330	sd->Orbit.DotDeltaCrossTrack = dotRsw(3);
2331	}
2332
2333	outLine.sprintf("%d %.1f %s %u %10.3f %8.3f %8.3f %8.3f %8.3f %8.3f\n", GPSweek,
2334	GPSweeks, eph->prn().toString().c_str(), eph->IOD(), dClkA0, dClkA1, dClkA2,
2335	rsw(1), rsw(2), rsw(3));
2336
2337	// RTNET full clock for RINEX and SP3 file
2338	// ---------------------------------------
2339	double relativity = -2.0 * DotProduct(xP, rtnVel) / t_CST::c;
2340	double clkRnx = (rtnClk[0] - relativity) / t_CST::c; // [s]
2341	double clkRnxRate = rtnClk[1] / t_CST::c; // [s/s = -]
2342	double clkRnxAcc = rtnClk[2] / t_CST::c; // [s/s² ) -/s]
2343
2344	if (_rnx) {
2345	double clkRnxSig, clkRnxRateSig, clkRnxAccSig;
2346	int s = rtnClkSig.size();
2347	switch (s) {
2348	case 1:
2349	clkRnxSig = rtnClkSig[0] / t_CST::c; // [s]
2350	clkRnxRateSig = 0.0; // [s/s = -]
2351	clkRnxAccSig = 0.0; // [s/s² ) -/s]
2352	break;
2353	case 2:
2354	clkRnxSig = rtnClkSig[0] / t_CST::c; // [s]
2355	clkRnxRateSig = rtnClkSig[1] / t_CST::c; // [s/s = -]
2356	clkRnxAccSig = 0.0; // [s/s² ) -/s]
2357	break;
2358	case 3:
2359	clkRnxSig = rtnClkSig[0] / t_CST::c; // [s]
2360	clkRnxRateSig = rtnClkSig[1] / t_CST::c; // [s/s = -]
2361	clkRnxAccSig = rtnClkSig[2] / t_CST::c; // [s/s² ) -/s]
2362	break;
2363	}
2364	_rnx->write(GPSweek, GPSweeks, prn, clkRnx, clkRnxRate, clkRnxAcc,
2365	clkRnxSig, clkRnxRateSig, clkRnxAccSig);
2366	}
2367	if (_sp3) {
2368	_sp3->write(GPSweek, GPSweeks, prn, rtnCoM, clkRnx, rtnVel, clkRnxRate);
2369	}
2370	return success;
2371	}
2372
2373	// Transform Coordinates
2374	////////////////////////////////////////////////////////////////////////////
2375	void bncRtnetUploadCaster::crdTrafo(int GPSWeek, ColumnVector& xyz,
2376	double& dc) {
2377
2378	// Current epoch minus 2000.0 in years
2379	// ------------------------------------
2380	double dt = (GPSWeek - (1042.0 + 6.0 / 7.0)) / 365.2422 * 7.0 + 2000.0 - _t0;
2381
2382	ColumnVector dx(3);
2383
2384	dx(1) = _dx + dt * _dxr;
2385	dx(2) = _dy + dt * _dyr;
2386	dx(3) = _dz + dt * _dzr;
2387
2388	static const double arcSec = 180.0 * 3600.0 / M_PI;
2389
2390	double ox = (_ox + dt * _oxr) / arcSec;
2391	double oy = (_oy + dt * _oyr) / arcSec;
2392	double oz = (_oz + dt * _ozr) / arcSec;
2393
2394	double sc = 1.0 + _sc * 1e-9 + dt * _scr * 1e-9;
2395
2396	// Specify approximate center of area
2397	// ----------------------------------
2398	ColumnVector meanSta(3);
2399
2400	if (_crdTrafo == "ETRF2000") {
2401	meanSta(1) = 3661090.0;
2402	meanSta(2) = 845230.0;
2403	meanSta(3) = 5136850.0;
2404	}
2405	else if (_crdTrafo == "GDA2020") {
2406	meanSta(1) = -4052050.0;
2407	meanSta(2) = 4212840.0;
2408	meanSta(3) = -2545110.0;
2409	}
2410	else if (_crdTrafo == "SIRGAS2000") {
2411	meanSta(1) = 3740860.0;
2412	meanSta(2) = -4964290.0;
2413	meanSta(3) = -1425420.0;
2414	}
2415	else if (_crdTrafo == "DREF91") {
2416	meanSta(1) = 3959579.0;
2417	meanSta(2) = 721719.0;
2418	meanSta(3) = 4931539.0;
2419	}
2420	else if (_crdTrafo == "Custom") {
2421	meanSta(1) = 0.0; // TODO
2422	meanSta(2) = 0.0; // TODO
2423	meanSta(3) = 0.0; // TODO
2424	}
2425
2426	// Clock correction proportional to topocentric distance to satellites
2427	// -------------------------------------------------------------------
2428	double rho = (xyz - meanSta).norm_Frobenius();
2429	dc = rho * (sc - 1.0) / sc / t_CST::c;
2430
2431	Matrix rMat(3, 3);
2432	rMat(1, 1) = 1.0;
2433	rMat(1, 2) = -oz;
2434	rMat(1, 3) = oy;
2435	rMat(2, 1) = oz;
2436	rMat(2, 2) = 1.0;
2437	rMat(2, 3) = -ox;
2438	rMat(3, 1) = -oy;
2439	rMat(3, 2) = ox;
2440	rMat(3, 3) = 1.0;
2441
2442	xyz = sc * rMat * xyz + dx;
2443	}
2444
2445	int bncRtnetUploadCaster::determineUpdateInd(double samplingRate) {
2446
2447	if (samplingRate == 10.0) {
2448	return 3;
2449	}
2450	else if (samplingRate == 15.0) {
2451	return 4;
2452	}
2453	else if (samplingRate == 30.0) {
2454	return 5;
2455	}
2456	else if (samplingRate == 60.0) {
2457	return 6;
2458	}
2459	else if (samplingRate == 120.0) {
2460	return 7;
2461	}
2462	else if (samplingRate == 240.0) {
2463	return 8;
2464	}
2465	else if (samplingRate == 300.0) {
2466	return 9;
2467	}
2468	else if (samplingRate == 600.0) {
2469	return 10;
2470	}
2471	else if (samplingRate == 900.0) {
2472	return 11;
2473	}
2474	else if (samplingRate == 1800.0) {
2475	return 12;
2476	}
2477	else if (samplingRate == 3600.0) {
2478	return 13;
2479	}
2480	else if (samplingRate == 7200.0) {
2481	return 14;
2482	}
2483	else if (samplingRate == 10800.0) {
2484	return 15;
2485	}
2486	return 2; // default
2487	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: