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 diff --git a/src/cssrlib/pppssr.py b/src/cssrlib/pppssr.py index cee5c1f..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]]) @@ -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: @@ -858,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) @@ -930,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))