Context Navigation

source: ntrip/trunk/BNC/src/rinex/rnxobsfile.cpp@ 6251

Last change on this file since 6251 was 6251, checked in by mervart, 10 years ago

File size: 34.5 KB

Rev	Line
[3716]	1	// Part of BNC, a utility for retrieving decoding and
	2	// converting GNSS data streams from NTRIP broadcasters.
	3	//
	4	// Copyright (C) 2007
	5	// German Federal Agency for Cartography and Geodesy (BKG)
	6	// http://www.bkg.bund.de
	7	// Czech Technical University Prague, Department of Geodesy
	8	// http://www.fsv.cvut.cz
	9	//
	10	// Email: euref-ip@bkg.bund.de
	11	//
	12	// This program is free software; you can redistribute it and/or
	13	// modify it under the terms of the GNU General Public License
	14	// as published by the Free Software Foundation, version 2.
	15	//
	16	// This program is distributed in the hope that it will be useful,
	17	// but WITHOUT ANY WARRANTY; without even the implied warranty of
	18	// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
	19	// GNU General Public License for more details.
	20	//
	21	// You should have received a copy of the GNU General Public License
	22	// along with this program; if not, write to the Free Software
	23	// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
	24
	25	/* -------------------------------------------------------------------------
	26	* BKG NTRIP Client
	27	* -------------------------------------------------------------------------
	28	*
	29	* Class: t_rnxObsFile
	30	*
	31	* Purpose: Reads RINEX Observation File
	32	*
	33	* Author: L. Mervart
	34	*
	35	* Created: 24-Jan-2012
	36	*
	37	* Changes:
	38	*
	39	* -----------------------------------------------------------------------*/
	40
	41	#include <iostream>
	42	#include <iomanip>
	43	#include <sstream>
[3722]	44	#include "rnxobsfile.h"
[3716]	45	#include "bncutils.h"
[5070]	46	#include "bnccore.h"
[6242]	47	#include "bncsettings.h"
[3716]	48
	49	using namespace std;
	50
	51	// Constructor
	52	////////////////////////////////////////////////////////////////////////////
[4480]	53	t_rnxObsHeader::t_rnxObsHeader() {
[6247]	54	_usedSystems = "GREJCS";
[3716]	55	_antNEU.ReSize(3); _antNEU = 0.0;
	56	_antXYZ.ReSize(3); _antXYZ = 0.0;
	57	_antBSG.ReSize(3); _antBSG = 0.0;
	58	_xyz.ReSize(3); _xyz = 0.0;
	59	_version = 0.0;
	60	_interval = 0.0;
[5742]	61	for (unsigned iPrn = 1; iPrn <= t_prn::MAXPRN_GPS; iPrn++) {
[3716]	62	_wlFactorsL1[iPrn] = 1;
	63	_wlFactorsL2[iPrn] = 1;
	64	}
	65	}
	66
	67	// Destructor
	68	////////////////////////////////////////////////////////////////////////////
[4480]	69	t_rnxObsHeader::~t_rnxObsHeader() {
[3716]	70	}
	71
	72	// Read Header
	73	////////////////////////////////////////////////////////////////////////////
[4480]	74	t_irc t_rnxObsHeader::read(QTextStream* stream, int maxLines) {
[4534]	75	_comments.clear();
[3716]	76	int numLines = 0;
[3718]	77	while ( stream->status() == QTextStream::Ok && !stream->atEnd() ) {
	78	QString line = stream->readLine(); ++ numLines;
	79	if (line.isEmpty()) {
[3716]	80	continue;
	81	}
[3718]	82	if (line.indexOf("END OF FILE") != -1) {
[3716]	83	break;
	84	}
[3718]	85	QString value = line.mid(0,60).trimmed();
	86	QString key = line.mid(60).trimmed();
[3716]	87	if (key == "END OF HEADER") {
	88	break;
	89	}
	90	else if (key == "RINEX VERSION / TYPE") {
[3718]	91	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	92	in >> _version;
	93	}
	94	else if (key == "MARKER NAME") {
	95	_markerName = value;
	96	}
[3930]	97	else if (key == "MARKER NUMBER") {
	98	_markerNumber = line.mid(0,20).trimmed();
	99	}
[3716]	100	else if (key == "ANT # / TYPE") {
[3930]	101	_antennaNumber = line.mid( 0,20).trimmed();
	102	_antennaName = line.mid(20,20).trimmed();
[3716]	103	}
[3930]	104	else if (key == "OBSERVER / AGENCY") {
	105	_observer = line.mid( 0,20).trimmed();
	106	_agency = line.mid(20,40).trimmed();
	107	}
	108	else if (key == "REC # / TYPE / VERS") {
	109	_receiverNumber = line.mid( 0,20).trimmed();
	110	_receiverType = line.mid(20,20).trimmed();
	111	_receiverVersion = line.mid(40,20).trimmed();
	112	}
[3716]	113	else if (key == "INTERVAL") {
[3718]	114	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	115	in >> _interval;
	116	}
[4494]	117	else if (key == "COMMENT") {
	118	_comments << line.mid(0,60).trimmed();
	119	}
[3716]	120	else if (key == "WAVELENGTH FACT L1/2") {
[3718]	121	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	122	int wlFactL1 = 0;
	123	int wlFactL2 = 0;
	124	int numSat = 0;
	125	in >> wlFactL1 >> wlFactL2 >> numSat;
	126	if (numSat == 0) {
[5742]	127	for (unsigned iPrn = 1; iPrn <= t_prn::MAXPRN_GPS; iPrn++) {
[3716]	128	_wlFactorsL1[iPrn] = wlFactL1;
	129	_wlFactorsL2[iPrn] = wlFactL2;
	130	}
	131	}
	132	else {
	133	for (int ii = 0; ii < numSat; ii++) {
[3718]	134	QString prn; in >> prn;
[3716]	135	if (prn[0] == 'G') {
	136	int iPrn;
	137	readInt(prn, 1, 2, iPrn);
	138	_wlFactorsL1[iPrn] = wlFactL1;
	139	_wlFactorsL2[iPrn] = wlFactL2;
	140	}
	141	}
	142	}
	143	}
	144	else if (key == "APPROX POSITION XYZ") {
[3718]	145	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	146	in >> _xyz[0] >> _xyz[1] >> _xyz[2];
	147	}
	148	else if (key == "ANTENNA: DELTA H/E/N") {
[3718]	149	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	150	in >> _antNEU[2] >> _antNEU[1] >> _antNEU[0];
	151	}
	152	else if (key == "ANTENNA: DELTA X/Y/Z") {
[3718]	153	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	154	in >> _antXYZ[0] >> _antXYZ[1] >> _antXYZ[2];
	155	}
	156	else if (key == "ANTENNA: B.SIGHT XYZ") {
[3718]	157	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3716]	158	in >> _antBSG[0] >> _antBSG[1] >> _antBSG[2];
	159	}
	160	else if (key == "# / TYPES OF OBSERV") {
[6223]	161	if (_version == 0.0) {
	162	_version = t_rnxObsHeader::defaultRnxObsVersion2;
	163	}
[3864]	164	QTextStream* in = new QTextStream(value.toAscii(), QIODevice::ReadOnly);
[3716]	165	int nTypes;
[3864]	166	*in >> nTypes;
[6245]	167	char sys0 = _usedSystems[0].toAscii();
[6125]	168	_obsTypes[sys0].clear();
[3716]	169	for (int ii = 0; ii < nTypes; ii++) {
[3864]	170	if (ii > 0 && ii % 9 == 0) {
	171	line = stream->readLine(); ++numLines;
	172	delete in;
[4152]	173	in = new QTextStream(line.left(60).toAscii(), QIODevice::ReadOnly);
[3864]	174	}
[3718]	175	QString hlp;
[3864]	176	*in >> hlp;
[6125]	177	_obsTypes[sys0].append(hlp);
[3716]	178	}
[6245]	179	for (int ii = 1; ii < _usedSystems.length(); ii++) {
	180	char sysI = _usedSystems[ii].toAscii();
[6125]	181	_obsTypes[sysI] = _obsTypes[sys0];
[6119]	182	}
[3716]	183	}
	184	else if (key == "SYS / # / OBS TYPES") {
[6223]	185	if (_version == 0.0) {
[6228]	186	_version = t_rnxObsHeader::defaultRnxObsVersion3;
[6223]	187	}
[3718]	188	QTextStream* in = new QTextStream(value.toAscii(), QIODevice::ReadOnly);
[3716]	189	char sys;
	190	int nTypes;
	191	*in >> sys >> nTypes;
[6119]	192	_obsTypes[sys].clear();
[3716]	193	for (int ii = 0; ii < nTypes; ii++) {
	194	if (ii > 0 && ii % 13 == 0) {
[3718]	195	line = stream->readLine(); ++numLines;
[3716]	196	delete in;
[3718]	197	in = new QTextStream(line.toAscii(), QIODevice::ReadOnly);
[3716]	198	}
[3718]	199	QString hlp;
[3716]	200	*in >> hlp;
[6119]	201	_obsTypes[sys].push_back(hlp);
[3716]	202	}
	203	delete in;
	204	}
[3837]	205	else if (key == "TIME OF FIRST OBS") {
	206	QTextStream in(value.toAscii(), QIODevice::ReadOnly);
	207	int year, month, day, hour, min;
	208	double sec;
	209	in >> year >> month >> day >> hour >> min >> sec;
	210	_startTime.set(year, month, day, hour, min, sec);
	211	}
[3716]	212	if (maxLines > 0 && numLines == maxLines) {
	213	break;
	214	}
	215	}
	216
[6245]	217	// Systems used
	218	// ------------
	219	_usedSystems.clear();
	220	QMapIterator<char, QStringList> it(_obsTypes);
	221	while (it.hasNext()) {
	222	it.next();
	223	_usedSystems += QChar(it.key());
	224	}
	225
[3716]	226	return success;
	227	}
	228
[6222]	229	// Set Default Header
	230	////////////////////////////////////////////////////////////////////////////
	231	void t_rnxObsHeader::setDefault(const QString& markerName, int version) {
	232
	233	_markerName = markerName;
	234
	235	if (version <= 2) {
	236	_version = t_rnxObsHeader::defaultRnxObsVersion2;
	237	}
	238	else {
	239	_version = t_rnxObsHeader::defaultRnxObsVersion3;
	240	}
	241
	242	_comments << "Default set of observation types used";
	243
	244	_obsTypes.clear();
	245	if (_version < 3.0) {
	246	_obsTypes['G'] << "C1" << "P1" << "L1" << "S1"
	247	<< "C2" << "P2" << "L2" << "S2";
	248	_obsTypes['R'] = _obsTypes['G'];
	249	_obsTypes['E'] = _obsTypes['G'];
	250	_obsTypes['J'] = _obsTypes['G'];
	251	_obsTypes['S'] = _obsTypes['G'];
	252	_obsTypes['C'] = _obsTypes['G'];
	253	}
	254	else {
	255	_obsTypes['G'] << "C1C" << "L1C" << "S1C"
	256	<< "C2W" << "L2W" << "S2W"
	257	<< "C5" << "L5" << "S5";
	258
	259	_obsTypes['J'] = _obsTypes['G'];
	260
	261	_obsTypes['R'] << "C1C" << "L1C" << "S1C"
	262	<< "C2P" << "L2P" << "S2P";
	263
	264	_obsTypes['E'] << "C1" << "L1" << "S1"
	265	<< "C5" << "L5" << "S5"
	266	<< "C7" << "L7" << "S7"
	267	<< "C8" << "L8" << "S8";
	268
	269	_obsTypes['S'] << "C1" << "L1" << "S1"
	270	<< "C5" << "L5" << "S5";
	271
	272	_obsTypes['C'] << "C1" << "L1" << "S1"
	273	<< "C6" << "L6" << "S6"
	274	<< "C7" << "L7" << "S7";
	275	}
	276	}
	277
	278	// Copy header
	279	////////////////////////////////////////////////////////////////////////////
	280	void t_rnxObsHeader::set(const t_rnxObsHeader& header, int version,
	281	const QStringList* useObsTypes) {
	282
	283	if (version <= 2) {
	284	_version = t_rnxObsHeader::defaultRnxObsVersion2;
	285	}
	286	else {
	287	_version = t_rnxObsHeader::defaultRnxObsVersion3;
	288	}
	289	_interval = header._interval;
	290	_antennaNumber = header._antennaNumber;
	291	_antennaName = header._antennaName;
	292	_markerName = header._markerName;
	293	_markerNumber = header._markerNumber;
	294	_antNEU = header._antNEU;
	295	_antXYZ = header._antXYZ;
	296	_antBSG = header._antBSG;
	297	_xyz = header._xyz;
	298	_observer = header._observer;
	299	_agency = header._agency;
	300	_receiverNumber = header._receiverNumber;
	301	_receiverType = header._receiverType;
	302	_receiverVersion = header._receiverVersion;
[6224]	303	_startTime = header._startTime;
	304	_comments = header._comments;
[6245]	305	_usedSystems = header._usedSystems;
[6222]	306	for (unsigned iPrn = 1; iPrn <= t_prn::MAXPRN_GPS; iPrn++) {
	307	_wlFactorsL1[iPrn] = header._wlFactorsL1[iPrn];
	308	_wlFactorsL2[iPrn] = header._wlFactorsL2[iPrn];
	309	}
	310
	311	// Set observation types
	312	// ---------------------
	313	_obsTypes.clear();
	314	if (!useObsTypes \|\| useObsTypes->size() == 0) {
	315	if (int(_version) == int(header._version)) {
	316	_obsTypes = header._obsTypes;
	317	}
	318	else {
[6250]	319	if (_version >= 3.0) {
	320	for (int iSys = 0; iSys < header.numSys(); iSys++) {
	321	char sys = header.system(iSys);
	322	for (int iType = 0; iType < header.nTypes(sys); iType++) {
	323	QString type = header.obsType(sys, iType, _version);
	324	if (!_obsTypes[sys].contains(type)) {
	325	_obsTypes[sys].push_back(type);
	326	}
[6249]	327	}
[6250]	328	}
	329	}
	330	else {
	331	for (int iSys = 0; iSys < header.numSys(); iSys++) {
	332	char sys = header.system(iSys);
	333	for (int iType = 0; iType < header.nTypes(sys); iType++) {
	334	QString type = header.obsType(sys, iType, _version);
[6249]	335	for (int jSys = 0; jSys < _usedSystems.length(); jSys++) {
	336	char thisSys = _usedSystems[jSys].toAscii();
	337	if (!_obsTypes[thisSys].contains(type)) {
	338	_obsTypes[thisSys].push_back(type);
	339	}
	340	}
	341	}
[6222]	342	}
	343	}
	344	}
	345	}
	346	else {
	347	for (int iType = 0; iType < useObsTypes->size(); iType++) {
	348	if (useObsTypes->at(iType).indexOf(":") != -1) {
	349	QStringList hlp = useObsTypes->at(iType).split(":", QString::SkipEmptyParts);
	350	if (hlp.size() == 2 && hlp[0].length() == 1) {
[6248]	351	if (_version >= 3.0) {
	352	char sys = hlp[0][0].toAscii();
	353	QString type = t_rnxObsFile::type2to3(sys, hlp[1]);
	354	if (!_obsTypes[sys].contains(type)) {
	355	_obsTypes[sys].push_back(type);
	356	}
	357	}
	358	else {
	359	for (int iSys = 0; iSys < _usedSystems.length(); iSys++) {
	360	char sys = _usedSystems[iSys].toAscii();
	361	QString type = t_rnxObsFile::type3to2(sys, hlp[1]);
	362	if (!_obsTypes[sys].contains(type)) {
	363	_obsTypes[sys].push_back(type);
	364	}
	365	}
	366	}
[6222]	367	}
	368	}
	369	else {
[6245]	370	for (int iSys = 0; iSys < _usedSystems.length(); iSys++) {
[6248]	371	char sys = _usedSystems[iSys].toAscii();
	372	QString type = _version >= 3.0 ? t_rnxObsFile::type2to3(sys, useObsTypes->at(iType)) :
	373	t_rnxObsFile::type3to2(sys, useObsTypes->at(iType));
	374	if (!_obsTypes[sys].contains(type)) {
	375	_obsTypes[sys].push_back(type);
	376	}
[6222]	377	}
	378	}
	379	}
[6245]	380	_usedSystems.clear();
[6244]	381	QMapIterator<char, QStringList> it(_obsTypes);
	382	while (it.hasNext()) {
	383	it.next();
[6245]	384	_usedSystems += QChar(it.key());
[6244]	385	}
[6222]	386	}
	387	}
	388
[4481]	389	// Write Header
	390	////////////////////////////////////////////////////////////////////////////
	391	void t_rnxObsHeader::write(QTextStream* stream,
	392	const QMap<QString, QString>* txtMap) const {
	393
[4493]	394	QStringList newComments;
[5068]	395	QString runBy = BNC_CORE->userName();
[4481]	396
	397	if (txtMap) {
	398	QMapIterator<QString, QString> it(*txtMap);
	399	while (it.hasNext()) {
	400	it.next();
	401	if (it.key() == "RUN BY") {
	402	runBy = it.value();
	403	}
	404	else if (it.key() == "COMMENT") {
[4493]	405	newComments = it.value().split("\\n", QString::SkipEmptyParts);
[4481]	406	}
	407	}
	408	}
	409
	410	*stream << QString("%1 Observation data Mixed")
	411	.arg(_version, 9, 'f', 2)
	412	.leftJustified(60)
	413	<< "RINEX VERSION / TYPE\n";
	414
	415	const QString fmtDate = (_version < 3.0) ? "dd-MMM-yy hh:mm"
	416	: "yyyyMMdd hhmmss UTC";
	417	*stream << QString("%1%2%3")
[5068]	418	.arg(BNC_CORE->pgmName(), -20)
[4481]	419	.arg(runBy.trimmed().left(20), -20)
	420	.arg(QDateTime::currentDateTime().toUTC().toString(fmtDate), -20)
	421	.leftJustified(60)
	422	<< "PGM / RUN BY / DATE\n";
	423
[4493]	424	QStringListIterator itCmnt(_comments + newComments);
[4481]	425	while (itCmnt.hasNext()) {
	426	*stream << itCmnt.next().trimmed().left(60).leftJustified(60) << "COMMENT\n";
	427	}
	428
	429	*stream << QString("%1")
	430	.arg(_markerName, -60)
	431	.leftJustified(60)
	432	<< "MARKER NAME\n";
	433
	434	if (!_markerNumber.isEmpty()) {
	435	*stream << QString("%1")
	436	.arg(_markerNumber, -20)
	437	.leftJustified(60)
	438	<< "MARKER NUMBER\n";
	439	}
	440
	441	*stream << QString("%1%2")
[4529]	442	.arg(_observer, -20)
	443	.arg(_agency, -40)
[4481]	444	.leftJustified(60)
	445	<< "OBSERVER / AGENCY\n";
	446
	447	*stream << QString("%1%2%3")
[4529]	448	.arg(_receiverNumber, -20)
	449	.arg(_receiverType, -20)
	450	.arg(_receiverVersion, -20)
[4481]	451	.leftJustified(60)
	452	<< "REC # / TYPE / VERS\n";
	453
	454	*stream << QString("%1%2")
[4529]	455	.arg(_antennaNumber, -20)
	456	.arg(_antennaName, -20)
[4481]	457	.leftJustified(60)
	458	<< "ANT # / TYPE\n";
	459
	460	*stream << QString("%1%2%3")
	461	.arg(_xyz(1), 14, 'f', 4)
	462	.arg(_xyz(2), 14, 'f', 4)
	463	.arg(_xyz(3), 14, 'f', 4)
	464	.leftJustified(60)
	465	<< "APPROX POSITION XYZ\n";
	466
	467	*stream << QString("%1%2%3")
	468	.arg(_antNEU(3), 14, 'f', 4)
	469	.arg(_antNEU(2), 14, 'f', 4)
	470	.arg(_antNEU(1), 14, 'f', 4)
	471	.leftJustified(60)
	472	<< "ANTENNA: DELTA H/E/N\n";
	473
	474	if (_version < 3.0) {
	475	int defaultWlFact1 = _wlFactorsL1[1];
	476	int defaultWlFact2 = _wlFactorsL2[1]; // TODO check all prns
	477	*stream << QString("%1%2")
	478	.arg(defaultWlFact1, 6)
	479	.arg(defaultWlFact2, 6)
	480	.leftJustified(60)
	481	<< "WAVELENGTH FACT L1/2\n";
	482	}
	483
	484	*stream << obsTypesStrings().join("");
	485
[4484]	486	if (_interval > 0) {
	487	*stream << QString("%1")
	488	.arg(_interval, 10, 'f', 3)
	489	.leftJustified(60)
	490	<< "INTERVAL\n";
	491	}
[4481]	492
	493	unsigned year, month, day, hour, min;
	494	double sec;
	495	_startTime.civil_date(year, month, day);
	496	_startTime.civil_time(hour, min, sec);
	497	*stream << QString("%1%2%3%4%5%6%7")
	498	.arg(year, 6)
	499	.arg(month, 6)
	500	.arg(day, 6)
	501	.arg(hour, 6)
	502	.arg(min, 6)
	503	.arg(sec, 13, 'f', 7)
	504	.arg("GPS", 8)
	505	.leftJustified(60)
	506	<< "TIME OF FIRST OBS\n";
	507
	508	*stream << QString()
	509	.leftJustified(60)
	510	<< "END OF HEADER\n";
	511	}
	512
[6119]	513	// Number of Different Systems
	514	////////////////////////////////////////////////////////////////////////////
	515	int t_rnxObsHeader::numSys() const {
	516	return _obsTypes.size();
	517	}
	518
[6130]	519	//
	520	////////////////////////////////////////////////////////////////////////////
	521	char t_rnxObsHeader::system(int iSys) const {
	522	int iSysLocal = -1;
[6230]	523	QMapIterator<char, QStringList> it(_obsTypes);
[6130]	524	while (it.hasNext()) {
	525	++iSysLocal;
	526	it.next();
	527	if (iSysLocal == iSys) {
	528	return it.key();
	529	}
	530	}
	531	return ' ';
	532	}
	533
[3716]	534	// Number of Observation Types (satellite-system specific)
	535	////////////////////////////////////////////////////////////////////////////
[4480]	536	int t_rnxObsHeader::nTypes(char sys) const {
[6119]	537	if (_obsTypes.contains(sys)) {
	538	return _obsTypes[sys].size();
[3716]	539	}
	540	else {
[6119]	541	return 0;
[3716]	542	}
	543	}
	544
	545	// Observation Type (satellite-system specific)
	546	////////////////////////////////////////////////////////////////////////////
[6126]	547	QString t_rnxObsHeader::obsType(char sys, int index, double version) const {
	548
	549	if (version == 0.0) {
	550	version = _version;
	551	}
[6119]	552	if (_obsTypes.contains(sys)) {
[6126]	553	QString origType = _obsTypes[sys].at(index);
	554	if (int(version) == int(_version)) {
	555	return origType;
	556	}
	557	else if (int(version) == 2) {
[6192]	558	return t_rnxObsFile::type3to2(sys, origType);
[6126]	559	}
	560	else if (int(version) == 3) {
[6192]	561	return t_rnxObsFile::type2to3(sys, origType);
[6126]	562	}
[3716]	563	}
[6126]	564	return "";
[3716]	565	}
	566
[4481]	567	// Write Observation Types
	568	////////////////////////////////////////////////////////////////////////////
	569	QStringList t_rnxObsHeader::obsTypesStrings() const {
	570
	571	QStringList strList;
	572
	573	if (_version < 3.0) {
[6245]	574	char sys0 = _usedSystems[0].toAscii();
[4481]	575	QString hlp;
[6125]	576	QTextStream(&hlp) << QString("%1").arg(_obsTypes[sys0].size(), 6);
	577	for (int ii = 0; ii < _obsTypes[sys0].size(); ii++) {
	578	QTextStream(&hlp) << QString("%1").arg(_obsTypes[sys0][ii], 6);
	579	if ((ii+1) % 9 == 0 \|\| ii == _obsTypes[sys0].size()-1) {
[4481]	580	strList.append(hlp.leftJustified(60) + "# / TYPES OF OBSERV\n");
	581	hlp = QString().leftJustified(6);
	582	}
	583	}
	584	}
	585	else {
[6219]	586	for (int iSys = 0; iSys < numSys(); iSys++) {
	587	char sys = system(iSys);
[4481]	588	QString hlp;
[6219]	589	QTextStream(&hlp) << QString("%1 %2").arg(sys).arg(nTypes(sys), 3);
	590	for (int iType = 0; iType < nTypes(sys); iType++) {
[6221]	591	QString type = obsType(sys, iType);
[6219]	592	QTextStream(&hlp) << QString(" %1").arg(type, -3);
	593	if ((iType+1) % 13 == 0 \|\| iType == nTypes(sys)-1) {
[4481]	594	strList.append(hlp.leftJustified(60) + "SYS / # / OBS TYPES\n");
	595	hlp = QString().leftJustified(6);
	596	}
	597	}
	598	}
	599	}
	600
	601	return strList;
	602	}
	603
[3716]	604	// Constructor
	605	////////////////////////////////////////////////////////////////////////////
[3843]	606	t_rnxObsFile::t_rnxObsFile(const QString& fileName, e_inpOut inpOut) {
	607	_inpOut = inpOut;
[3716]	608	_stream = 0;
	609	_flgPowerFail = false;
[3843]	610	if (_inpOut == input) {
	611	openRead(fileName);
	612	}
	613	else {
[3845]	614	openWrite(fileName);
[3843]	615	}
[3716]	616	}
	617
[3845]	618	// Open for input
[3716]	619	////////////////////////////////////////////////////////////////////////////
[3843]	620	void t_rnxObsFile::openRead(const QString& fileName) {
[3718]	621
	622	_fileName = fileName; expandEnvVar(_fileName);
	623	_file = new QFile(_fileName);
	624	_file->open(QIODevice::ReadOnly \| QIODevice::Text);
	625	_stream = new QTextStream();
	626	_stream->setDevice(_file);
	627
[3716]	628	_header.read(_stream);
	629
	630	// Guess Observation Interval
	631	// --------------------------
	632	if (_header._interval == 0.0) {
	633	bncTime ttPrev;
	634	for (int iEpo = 0; iEpo < 10; iEpo++) {
	635	const t_rnxEpo* rnxEpo = nextEpoch();
	636	if (!rnxEpo) {
[3718]	637	throw QString("t_rnxObsFile: not enough epochs");
[3716]	638	}
	639	if (iEpo > 0) {
	640	double dt = rnxEpo->tt - ttPrev;
	641	if (_header._interval == 0.0 \|\| dt < _header._interval) {
	642	_header._interval = dt;
	643	}
	644	}
	645	ttPrev = rnxEpo->tt;
	646	}
[3719]	647	_stream->seek(0);
[3716]	648	_header.read(_stream);
	649	}
[3837]	650
	651	// Time of first observation
	652	// -------------------------
	653	if (!_header._startTime.valid()) {
	654	const t_rnxEpo* rnxEpo = nextEpoch();
	655	if (!rnxEpo) {
	656	throw QString("t_rnxObsFile: not enough epochs");
	657	}
	658	_header._startTime = rnxEpo->tt;
	659	_stream->seek(0);
	660	_header.read(_stream);
	661	}
[3716]	662	}
	663
[3845]	664	// Open for output
	665	////////////////////////////////////////////////////////////////////////////
	666	void t_rnxObsFile::openWrite(const QString& fileName) {
	667
	668	_fileName = fileName; expandEnvVar(_fileName);
	669	_file = new QFile(_fileName);
	670	_file->open(QIODevice::WriteOnly \| QIODevice::Text);
	671	_stream = new QTextStream();
	672	_stream->setDevice(_file);
	673	}
	674
[3716]	675	// Destructor
	676	////////////////////////////////////////////////////////////////////////////
	677	t_rnxObsFile::~t_rnxObsFile() {
	678	close();
	679	}
	680
	681	// Close
	682	////////////////////////////////////////////////////////////////////////////
	683	void t_rnxObsFile::close() {
	684	delete _stream; _stream = 0;
[3718]	685	delete _file; _file = 0;
[3716]	686	}
	687
	688	// Handle Special Epoch Flag
	689	////////////////////////////////////////////////////////////////////////////
[4540]	690	void t_rnxObsFile::handleEpochFlag(int flag, const QString& line,
	691	bool& headerReRead) {
[3716]	692
[4540]	693	headerReRead = false;
	694
[3716]	695	// Power Failure
	696	// -------------
	697	if (flag == 1) {
	698	_flgPowerFail = true;
	699	}
	700
	701	// Start moving antenna
	702	// --------------------
	703	else if (flag == 2) {
	704	// no action
	705	}
	706
	707	// Re-Read Header
	708	// --------------
[5468]	709	else if (flag == 3 \|\| flag == 4 \|\| flag == 5) {
[3716]	710	int numLines = 0;
	711	if (version() < 3.0) {
	712	readInt(line, 29, 3, numLines);
	713	}
	714	else {
	715	readInt(line, 32, 3, numLines);
	716	}
[5468]	717	if (flag == 3 \|\| flag == 4) {
	718	_header.read(_stream, numLines);
	719	headerReRead = true;
	720	}
	721	else {
	722	for (int ii = 0; ii < numLines; ii++) {
	723	_stream->readLine();
	724	}
	725	}
[3716]	726	}
	727
	728	// Unhandled Flag
	729	// --------------
	730	else {
[3718]	731	throw QString("t_rnxObsFile: unhandled flag\n" + line);
[3716]	732	}
	733	}
	734
	735	// Retrieve single Epoch
	736	////////////////////////////////////////////////////////////////////////////
[3994]	737	t_rnxObsFile::t_rnxEpo* t_rnxObsFile::nextEpoch() {
[3716]	738	_currEpo.clear();
	739	if (version() < 3.0) {
	740	return nextEpochV2();
	741	}
	742	else {
	743	return nextEpochV3();
	744	}
	745	}
	746
	747	// Retrieve single Epoch (RINEX Version 3)
	748	////////////////////////////////////////////////////////////////////////////
[3994]	749	t_rnxObsFile::t_rnxEpo* t_rnxObsFile::nextEpochV3() {
[3716]	750
[3719]	751	while ( _stream->status() == QTextStream::Ok && !_stream->atEnd() ) {
[3716]	752
[3718]	753	QString line = _stream->readLine();
[3716]	754
[3718]	755	if (line.isEmpty()) {
[3716]	756	continue;
	757	}
	758
	759	int flag = 0;
	760	readInt(line, 31, 1, flag);
	761	if (flag > 0) {
[4540]	762	bool headerReRead = false;
	763	handleEpochFlag(flag, line, headerReRead);
	764	if (headerReRead) {
	765	continue;
	766	}
[3716]	767	}
	768
[3718]	769	QTextStream in(line.mid(1).toAscii(), QIODevice::ReadOnly);
[3716]	770
	771	// Epoch Time
	772	// ----------
	773	int year, month, day, hour, min;
	774	double sec;
	775	in >> year >> month >> day >> hour >> min >> sec;
[3720]	776	_currEpo.tt.set(year, month, day, hour, min, sec);
[3716]	777
	778	// Number of Satellites
	779	// --------------------
	780	int numSat;
	781	readInt(line, 32, 3, numSat);
	782
	783	_currEpo.rnxSat.resize(numSat);
	784
	785	// Observations
	786	// ------------
	787	for (int iSat = 0; iSat < numSat; iSat++) {
[3718]	788	line = _stream->readLine();
[6119]	789	t_prn prn; prn.set(line.left(3).toAscii().data());
	790	_currEpo.rnxSat[iSat].prn = prn;
	791	char sys = prn.system();
[3716]	792	for (int iType = 0; iType < _header.nTypes(sys); iType++) {
	793	int pos = 3 + 16*iType;
	794	double obsValue = 0.0;
	795	int lli = 0;
	796	int snr = 0;
	797	readDbl(line, pos, 14, obsValue);
	798	readInt(line, pos + 14, 1, lli);
	799	readInt(line, pos + 15, 1, snr);
	800	if (_flgPowerFail) {
	801	lli \|= 1;
	802	}
[6119]	803	QString type = obsType(sys, iType);
	804	_currEpo.rnxSat[iSat].obs[type].value = obsValue;
	805	_currEpo.rnxSat[iSat].obs[type].lli = lli;
	806	_currEpo.rnxSat[iSat].obs[type].snr = snr;
[3716]	807	}
	808	}
	809
	810	_flgPowerFail = false;
	811
	812	return &_currEpo;
	813	}
	814
	815	return 0;
	816	}
	817
	818	// Retrieve single Epoch (RINEX Version 2)
	819	////////////////////////////////////////////////////////////////////////////
[3994]	820	t_rnxObsFile::t_rnxEpo* t_rnxObsFile::nextEpochV2() {
[3716]	821
[3719]	822	while ( _stream->status() == QTextStream::Ok && !_stream->atEnd() ) {
[3716]	823
[3718]	824	QString line = _stream->readLine();
[3716]	825
[3718]	826	if (line.isEmpty()) {
[3716]	827	continue;
	828	}
	829
	830	int flag = 0;
	831	readInt(line, 28, 1, flag);
	832	if (flag > 0) {
[4540]	833	bool headerReRead = false;
	834	handleEpochFlag(flag, line, headerReRead);
	835	if (headerReRead) {
	836	continue;
	837	}
[3716]	838	}
	839
[3718]	840	QTextStream in(line.toAscii(), QIODevice::ReadOnly);
[3716]	841
	842	// Epoch Time
	843	// ----------
	844	int year, month, day, hour, min;
	845	double sec;
	846	in >> year >> month >> day >> hour >> min >> sec;
	847	if (year < 80) {
	848	year += 2000;
	849	}
	850	else if (year < 100) {
	851	year += 1900;
	852	}
[3720]	853	_currEpo.tt.set(year, month, day, hour, min, sec);
[3716]	854
	855	// Number of Satellites
	856	// --------------------
	857	int numSat;
	858	readInt(line, 29, 3, numSat);
	859
	860	_currEpo.rnxSat.resize(numSat);
	861
	862	// Read Satellite Numbers
	863	// ----------------------
	864	int pos = 32;
	865	for (int iSat = 0; iSat < numSat; iSat++) {
	866	if (iSat > 0 && iSat % 12 == 0) {
[3718]	867	line = _stream->readLine();
[3716]	868	pos = 32;
	869	}
	870
[6119]	871	char sys = line.toAscii()[pos];
	872	int satNum; readInt(line, pos + 1, 2, satNum);
	873	_currEpo.rnxSat[iSat].prn.set(sys, satNum);
	874
[3716]	875	pos += 3;
	876	}
	877
	878	// Read Observation Records
	879	// ------------------------
	880	for (int iSat = 0; iSat < numSat; iSat++) {
[6119]	881	char sys = _currEpo.rnxSat[iSat].prn.system();
[3718]	882	line = _stream->readLine();
[3716]	883	pos = 0;
[6119]	884	for (int iType = 0; iType < _header.nTypes(sys); iType++) {
[3716]	885	if (iType > 0 && iType % 5 == 0) {
[3718]	886	line = _stream->readLine();
[3716]	887	pos = 0;
	888	}
	889	double obsValue = 0.0;
	890	int lli = 0;
	891	int snr = 0;
	892	readDbl(line, pos, 14, obsValue);
	893	readInt(line, pos + 14, 1, lli);
	894	readInt(line, pos + 15, 1, snr);
	895
	896	if (_flgPowerFail) {
	897	lli \|= 1;
	898	}
	899
[6119]	900	QString type = obsType(sys, iType);
	901	_currEpo.rnxSat[iSat].obs[type].value = obsValue;
	902	_currEpo.rnxSat[iSat].obs[type].lli = lli;
	903	_currEpo.rnxSat[iSat].obs[type].snr = snr;
[3716]	904
	905	pos += 16;
	906	}
	907	}
	908
	909	_flgPowerFail = false;
	910
	911	return &_currEpo;
	912	}
	913
	914	return 0;
	915	}
	916
[3845]	917	// Write Data Epoch
	918	////////////////////////////////////////////////////////////////////////////
	919	void t_rnxObsFile::writeEpoch(const t_rnxEpo* epo) {
[6246]	920	if (epo == 0) {
	921	return;
	922	}
	923	t_rnxEpo epoLocal;
	924	epoLocal.tt = epo->tt;
	925	for (unsigned ii = 0; ii < epo->rnxSat.size(); ii++) {
	926	const t_rnxSat& rnxSat = epo->rnxSat[ii];
	927	if (_header._obsTypes[rnxSat.prn.system()].size() > 0) {
	928	epoLocal.rnxSat.push_back(rnxSat);
	929	}
	930	}
	931
[3866]	932	if (version() < 3.0) {
[6246]	933	return writeEpochV2(_stream, _header, &epoLocal);
[3866]	934	}
	935	else {
[6246]	936	return writeEpochV3(_stream, _header, &epoLocal);
[3866]	937	}
	938	}
[3855]	939
[3866]	940	// Write Data Epoch (RINEX Version 2)
	941	////////////////////////////////////////////////////////////////////////////
[6222]	942	void t_rnxObsFile::writeEpochV2(QTextStream* stream, const t_rnxObsHeader& header,
	943	const t_rnxEpo* epo) {
[3866]	944
[3855]	945	unsigned year, month, day, hour, min;
	946	double sec;
	947	epo->tt.civil_date(year, month, day);
	948	epo->tt.civil_time(hour, min, sec);
	949
	950	QString dateStr;
[3868]	951	QTextStream(&dateStr) << QString(" %1 %2 %3 %4 %5%6")
[3856]	952	.arg(int(fmod(year, 100)), 2, 10, QChar('0'))
[3868]	953	.arg(month, 2, 10, QChar('0'))
	954	.arg(day, 2, 10, QChar('0'))
	955	.arg(hour, 2, 10, QChar('0'))
	956	.arg(min, 2, 10, QChar('0'))
	957	.arg(sec, 11, 'f', 7);
[3855]	958
[3857]	959	int flag = 0;
[6222]	960	*stream << dateStr << QString("%1%2").arg(flag, 3).arg(epo->rnxSat.size(), 3);
[3858]	961	for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
[3869]	962	const t_rnxSat& rnxSat = epo->rnxSat[iSat];
[3858]	963	if (iSat > 0 && iSat % 12 == 0) {
[6222]	964	*stream << endl << QString().leftJustified(32);
[3857]	965	}
[6222]	966	*stream << rnxSat.prn.toString().c_str();
[3857]	967	}
[6222]	968	*stream << endl;
[3858]	969	for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
[3967]	970
[3858]	971	const t_rnxSat& rnxSat = epo->rnxSat[iSat];
[6119]	972	char sys = rnxSat.prn.system();
[3966]	973
[6222]	974	for (int iTypeV2 = 0; iTypeV2 < header.nTypes(sys); iTypeV2++) {
[6216]	975	if (iTypeV2 > 0 && iTypeV2 % 5 == 0) {
[6222]	976	*stream << endl;
[3967]	977	}
[6222]	978	QString typeV2 = header.obsType(sys, iTypeV2);
[6216]	979	bool found = false;
[6231]	980
[6241]	981	QString preferredAttrib = signalPriorities(sys);
[6231]	982	for (int iPref = 0; iPref < preferredAttrib.length(); iPref++) {
	983	QMapIterator<QString, t_rnxObs> itObs(rnxSat.obs);
	984	while (itObs.hasNext()) {
	985	itObs.next();
	986	const QString& type = itObs.key();
	987	const t_rnxObs& rnxObs = itObs.value();
[6232]	988	if ( preferredAttrib[iPref] == '?' \|\|
[6241]	989	(type.length() == 2 && preferredAttrib[iPref] == '_' ) \|\|
[6232]	990	(type.length() == 3 && preferredAttrib[iPref] == type[2]) ) {
	991	if (typeV2 == type3to2(sys, type)) {
	992	found = true;
	993	if (rnxObs.value == 0.0) {
	994	*stream << QString().leftJustified(16);
[6231]	995	}
	996	else {
[6232]	997	*stream << QString("%1").arg(rnxObs.value, 14, 'f', 3);
	998	if (rnxObs.lli != 0.0) {
	999	*stream << QString("%1").arg(rnxObs.lli,1);
	1000	}
	1001	else {
	1002	*stream << ' ';
	1003	}
	1004	if (rnxObs.snr != 0.0) {
	1005	*stream << QString("%1").arg(rnxObs.snr,1);
	1006	}
	1007	else {
	1008	*stream << ' ';
	1009	}
[6231]	1010	}
[6251]	1011	goto end_loop_iPref;
[6216]	1012	}
[6127]	1013	}
[6231]	1014	}
[6251]	1015	} end_loop_iPref:
[6216]	1016	if (!found) {
[6222]	1017	*stream << QString().leftJustified(16);
[6216]	1018	}
[3858]	1019	}
[6222]	1020	*stream << endl;
[3858]	1021	}
[3845]	1022	}
[3866]	1023
	1024	// Write Data Epoch (RINEX Version 3)
	1025	////////////////////////////////////////////////////////////////////////////
[6222]	1026	void t_rnxObsFile::writeEpochV3(QTextStream* stream, const t_rnxObsHeader& header,
	1027	const t_rnxEpo* epo) {
[3866]	1028
	1029	unsigned year, month, day, hour, min;
	1030	double sec;
	1031	epo->tt.civil_date(year, month, day);
	1032	epo->tt.civil_time(hour, min, sec);
	1033
[3867]	1034	QString dateStr;
[3868]	1035	QTextStream(&dateStr) << QString("> %1 %2 %3 %4 %5%6")
[3867]	1036	.arg(year, 4)
	1037	.arg(month, 2, 10, QChar('0'))
	1038	.arg(day, 2, 10, QChar('0'))
	1039	.arg(hour, 2, 10, QChar('0'))
	1040	.arg(min, 2, 10, QChar('0'))
	1041	.arg(sec, 11, 'f', 7);
[3866]	1042
[3867]	1043	int flag = 0;
[6222]	1044	*stream << dateStr << QString("%1%2\n").arg(flag, 3).arg(epo->rnxSat.size(), 3);
[3867]	1045
[3869]	1046	for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
	1047	const t_rnxSat& rnxSat = epo->rnxSat[iSat];
[6119]	1048	char sys = rnxSat.prn.system();
[3964]	1049
[6222]	1050	*stream << rnxSat.prn.toString().c_str();
[6220]	1051
[6222]	1052	for (int iTypeV3 = 0; iTypeV3 < header.nTypes(sys); iTypeV3++) {
	1053	QString typeV3 = header.obsType(sys, iTypeV3);
[6220]	1054	bool found = false;
	1055	QMapIterator<QString, t_rnxObs> itObs(rnxSat.obs);
	1056	while (itObs.hasNext()) {
	1057	itObs.next();
	1058	const QString& type = itObs.key();
	1059	const t_rnxObs& rnxObs = itObs.value();
	1060	if (typeV3 == type2to3(sys, type)) {
	1061	found = true;
	1062	if (rnxObs.value == 0.0) {
[6222]	1063	*stream << QString().leftJustified(16);
[6127]	1064	}
	1065	else {
[6222]	1066	*stream << QString("%1").arg(rnxObs.value, 14, 'f', 3);
[6220]	1067	if (rnxObs.lli != 0.0) {
[6222]	1068	*stream << QString("%1").arg(rnxObs.lli,1);
[6220]	1069	}
	1070	else {
[6222]	1071	*stream << ' ';
[6220]	1072	}
	1073	if (rnxObs.snr != 0.0) {
[6222]	1074	*stream << QString("%1").arg(rnxObs.snr,1);
[6220]	1075	}
	1076	else {
[6222]	1077	*stream << ' ';
[6220]	1078	}
[6127]	1079	}
[3979]	1080	}
[3869]	1081	}
[6220]	1082	if (!found) {
[6222]	1083	*stream << QString().leftJustified(16);
[6220]	1084	}
[3869]	1085	}
[6222]	1086	*stream << endl;
[3869]	1087	}
[3866]	1088	}
[3956]	1089
	1090	// Translate Observation Type v2 --> v3
	1091	////////////////////////////////////////////////////////////////////////////
[6193]	1092	QString t_rnxObsFile::type2to3(char sys, const QString& typeV2) {
[6147]	1093	if (typeV2 == "P1") {
[6193]	1094	return (sys == 'G') ? "C1W" : "C1P";
[3956]	1095	}
[6147]	1096	else if (typeV2 == "P2") {
[6193]	1097	return (sys == 'G') ? "C2W" : "C2P";
[3956]	1098	}
[6147]	1099	return typeV2;
[3956]	1100	}
[3962]	1101
	1102	// Translate Observation Type v3 --> v2
	1103	////////////////////////////////////////////////////////////////////////////
[6192]	1104	QString t_rnxObsFile::type3to2(char /* sys */, const QString& typeV3) {
[6193]	1105	if (typeV3 == "C1P" \|\| typeV3 == "C1W") {
[3962]	1106	return "P1";
	1107	}
[6193]	1108	else if (typeV3 == "C2P" \|\| typeV3 == "C2W") {
[3962]	1109	return "P2";
	1110	}
[6147]	1111	return typeV3.left(2);
[3962]	1112	}
[4053]	1113
[5883]	1114	// Set Observations from RINEX File
	1115	////////////////////////////////////////////////////////////////////////////
[6137]	1116	void t_rnxObsFile::setObsFromRnx(const t_rnxObsFile* rnxObsFile, const t_rnxObsFile::t_rnxEpo* epo,
	1117	const t_rnxObsFile::t_rnxSat& rnxSat, t_satObs& obs) {
[5883]	1118
[6137]	1119	obs._staID = rnxObsFile->markerName().toAscii().constData();
	1120	obs._prn = rnxSat.prn;
	1121	obs._time = epo->tt;
[5883]	1122
[6137]	1123	char sys = rnxSat.prn.system();
[5883]	1124
[6137]	1125	for (int iType = 0; iType < rnxObsFile->nTypes(sys); iType++) {
	1126	QString type = rnxObsFile->obsType(sys, iType);
[6119]	1127	if (rnxSat.obs.contains(type)) {
	1128	const t_rnxObs& rnxObs = rnxSat.obs[type];
[6188]	1129	if (rnxObs.value != 0.0) {
	1130	string type2ch(type.mid(1).toAscii().data());
	1131
	1132	t_frqObs* frqObs = 0;
	1133	for (unsigned iFrq = 0; iFrq < obs._obs.size(); iFrq++) {
	1134	if (obs._obs[iFrq]->_rnxType2ch == type2ch) {
	1135	frqObs = obs._obs[iFrq];
	1136	break;
	1137	}
	1138	}
	1139	if (frqObs == 0) {
	1140	frqObs = new t_frqObs;
	1141	frqObs->_rnxType2ch = type2ch;
	1142	obs._obs.push_back(frqObs);
	1143	}
	1144
	1145	switch( type.toAscii().data()[0] ) {
	1146	case 'C':
	1147	frqObs->_codeValid = true;
	1148	frqObs->_code = rnxObs.value;
[6137]	1149	break;
[6188]	1150	case 'L':
	1151	frqObs->_phaseValid = true;
	1152	frqObs->_phase = rnxObs.value;
	1153	frqObs->_slip = (rnxObs.lli & 1);
	1154	break;
	1155	case 'D':
	1156	frqObs->_dopplerValid = true;
	1157	frqObs->_doppler = rnxObs.value;
	1158	break;
	1159	case 'S':
	1160	frqObs->_snrValid = true;
	1161	frqObs->_snr = rnxObs.value;
	1162	break;
[6137]	1163	}
[6119]	1164	}
[5883]	1165	}
[6119]	1166	}
	1167	}
	1168
[6241]	1169	// Tracking Mode Priorities
	1170	////////////////////////////////////////////////////////////////////////////
	1171	QString t_rnxObsFile::signalPriorities(char sys) {
[6242]	1172
	1173	bncSettings settings;
	1174	QStringList priorList = settings.value("rnxV2Priority").toString().split(" ", QString::SkipEmptyParts);
[6243]	1175	if (priorList.empty()) {
	1176	priorList << "CWPX_?";
	1177	}
[6242]	1178
	1179	QString result;
	1180	for (int ii = 0; ii < priorList.size(); ii++) {
	1181	if (priorList[ii].indexOf(":") != -1) {
	1182	QStringList hlp = priorList[ii].split(":", QString::SkipEmptyParts);
	1183	if (hlp.size() == 2 && hlp[0].length() == 1 && hlp[0][0] == sys) {
	1184	result = hlp[1];
	1185	break;
	1186	}
	1187	}
	1188	else {
	1189	result = priorList[ii];
	1190	}
	1191	}
	1192
	1193	return result;
[6241]	1194	}

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

Download in other formats: