Context Navigation

source: ntrip/branches/BNC_2.12/src/upload/bncrtnetuploadcaster.cpp@ 8992

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

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

Download in other formats: