Context Navigation

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

Last change on this file since 8710 was 8710, checked in by stuerze, 5 years ago
minor changes
File size: 105.4 KB

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

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

Download in other formats: