Context Navigation

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

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

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

Download in other formats: