From d5314490847eb6dd3a20abb1b3b9c43ab92f9393 Mon Sep 17 00:00:00 2001 From: AndreHauschild Date: Tue, 1 Oct 2024 12:19:42 +0200 Subject: [PATCH 1/3] Typo correction in comment --- src/cssrlib/ppp.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/cssrlib/ppp.py b/src/cssrlib/ppp.py index bb5e027..4e3cf50 100644 --- a/src/cssrlib/ppp.py +++ b/src/cssrlib/ppp.py @@ -210,7 +210,7 @@ def eci2ecef(tgps, erpv): dpsi, deps = nut_iau1980(t_, f) N = Rx(-eps-deps)@Rz(-dpsi)@Rx(eps) - # Greenwich aparent sidereal time [rad] + # Greenwich apparent sidereal time [rad] gmst = utc2gmst(tutc, erpv[2]) gast = gmst+dpsi*cos(eps) gast += (0.00264*sin(f[4])+0.000063*sin(2.0*f[4]))*gn.rCST.AS2R From 29c0a85d63c90dcda5ec2b623e2b3e3d3a06595d Mon Sep 17 00:00:00 2001 From: AndreHauschild Date: Tue, 1 Oct 2024 12:20:10 +0200 Subject: [PATCH 2/3] Use zero antenna corrections if nav.rcv_ant is None --- src/cssrlib/pppssr.py | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/src/cssrlib/pppssr.py b/src/cssrlib/pppssr.py index cee5c1f..5e184ce 100644 --- a/src/cssrlib/pppssr.py +++ b/src/cssrlib/pppssr.py @@ -659,8 +659,12 @@ def zdres(self, obs, cs, bsx, rs, vs, dts, rr, rtype=1): # Receiver/satellite antenna offset # - antrPR = antModelRx(self.nav, pos, e[i, :], sigsPR, rtype) - antrCP = antModelRx(self.nav, pos, e[i, :], sigsCP, rtype) + if self.nav.rcv_ant is None: + antrPR = [0.0 for _ in sigsPR] + antrCP = [0.0 for _ in sigsCP] + else: + antrPR = antModelRx(self.nav, pos, e[i, :], sigsPR, rtype) + antrCP = antModelRx(self.nav, pos, e[i, :], sigsCP, rtype) if self.nav.ephopt == 4: From d01d0d175b238bf24d0fa07ed6feaca4b896398d Mon Sep 17 00:00:00 2001 From: AndreHauschild Date: Tue, 22 Oct 2024 10:08:57 +0200 Subject: [PATCH 3/3] Fix typos --- src/cssrlib/pppssr.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/src/cssrlib/pppssr.py b/src/cssrlib/pppssr.py index 5e184ce..e475581 100644 --- a/src/cssrlib/pppssr.py +++ b/src/cssrlib/pppssr.py @@ -199,7 +199,7 @@ def initx(self, x0, v0, i): self.nav.P[j, i] = self.nav.P[i, j] = v0 if i == j else 0 def IB(self, s, f, na=3): - """ return index of phase ambguity """ + """ return index of phase ambiguity """ idx = na+uGNSS.MAXSAT*f+s-1 return idx @@ -335,7 +335,7 @@ def udstate(self, obs): if sys[i] == uGNSS.GLO: if sat[i] not in self.nav.glo_ch: - print("glonass channed not found: {:d}" + print("glonass channel not found: {:d}" .format(sat[i])) continue f1 = sig1.frequency(self.nav.glo_ch[sat[i]]) @@ -862,6 +862,7 @@ def sdres(self, obs, x, y, e, sat, el): self.IT(self.nav.na)]))) if self.nav.niono > 0: # iono is estimated + # SD ionosphere # idx_i = self.II(sat[i], self.nav.na) @@ -934,7 +935,7 @@ def sdres(self, obs, x, y, e, sat, el): nb[b] += 1 # counter for single-differences per signal nv += 1 # counter for single-difference observations - b += 1 # counter for signal (pseudrange+carrier-phase) + b += 1 # counter for signal (pseudorange+carrier-phase) v = np.resize(v, nv) H = np.resize(H, (nv, self.nav.nx))