cssrlib 1.2.1__tar.gz → 1.3.2__tar.gz

{cssrlib-1.2.1/src/cssrlib.egg-info → cssrlib-1.3.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cssrlib
-Version: 1.2.1
+Version: 1.3.2
 Summary: Python Toolkit for PPP/PPP-RTK positioning
 Home-page: https://github.com/hirokawa/cssrlib
 Author: Rui Hirokawa
@@ -16,6 +16,7 @@ Requires-Dist: notebook
 Requires-Dist: numpy
 Requires-Dist: matplotlib
 Requires-Dist: cartopy
+Requires-Dist: cryptography
 Dynamic: license-file
 
 CSSRlib - Toolkit for PPP-RTK/RTK in Python using Compact SSR

{cssrlib-1.2.1 → cssrlib-1.3.2}/pyproject.toml

@@ -7,7 +7,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "cssrlib"
-version = "1.2.1"
+version = "1.3.2"
 description = "Python Toolkit for PPP/PPP-RTK positioning"
 authors = [
     {name="Rui Hirokawa", email="rui.hirokawa@gmail.com"},
@@ -23,4 +23,5 @@ dependencies = [
     "numpy",
     "matplotlib",
     "cartopy",
+    "cryptography",
 ]

{cssrlib-1.2.1 → cssrlib-1.3.2}/setup.cfg

@@ -1,6 +1,6 @@
 [metadata]
 name = cssrlib
-version = 1.2.1
+version = 1.3.2
 author = Rui Hirokawa
 author_email = rui.hirokawa@gmail.com
 description = Python Toolkit for PPP/PPP-RTK positioning

{cssrlib-1.2.1 → cssrlib-1.3.2}/src/cssrlib/__init__.py

@@ -1,3 +1,5 @@
+"""CSSRLIB package for GNSS SSR and PPP/PPP-RTK processing."""
+
 
 from cssrlib.gnss import (gpst2time,epoch2time,ecef2pos,
                           prn2sat)

{cssrlib-1.2.1 → cssrlib-1.3.2}/src/cssrlib/cssr_bds.py

@@ -1,7 +1,6 @@
-"""
-BDS PPP correction data decoder
+"""BDS PPP correction data decoder.
 
-[1] BeiDou Navigation Satellite System Signal In Space Interface Control Document
+[1] BeiDou Navigation Satellite System Signal In Space ICD
 Precise Point Positioning Service Signal PPP-B2b (Version 1.0), 2020
 """
 
@@ -12,7 +11,10 @@ from cssrlib.gnss import bdt2time, bdt2gpst, uGNSS, uSIG, uTYP, rSigRnx
 
 
 class cssr_bds(cssr):
+    """Class to decode BDS PPP (B2b) Compact SSR messages."""
+
     def __init__(self, foutname=None):
+        """Initialize BDS PPP SSR decoder."""
         super().__init__(foutname)
         self.MAXNET = 1
         self.cssrmode = sCSSRTYPE.BDS_PPP
@@ -23,6 +25,7 @@ class cssr_bds(cssr):
         self.cb_scl = 0.017
 
     def ssig2rsig(self, sys: sGNSS, utyp: uTYP, ssig):
+        """Convert SSR signal index to RINEX signal code."""
         gps_tbl = {
             0: uSIG.L1C,
             1: uSIG.L1P,
@@ -72,13 +75,14 @@ class cssr_bds(cssr):
         return rSigRnx(sys, utyp, usig_tbl[ssig])
 
     def sval(self, u, n, scl):
-        """ calculate signed value based on n-bit int, lsb """
+        """Calculate signed value based on n-bit int, lsb."""
         invalid = -2**(n-1)
         dnu = -(2**(n-1)-1)  # this value seems to be invalid
         y = np.nan if u == invalid or u == dnu else u*scl
         return y
 
     def slot2prn(self, slot):
+        """Convert BDS B2b satellite slot index to system and PRN."""
         prn = 0
         sys = uGNSS.NONE
         if slot >= 1 and slot <= 63:
@@ -96,6 +100,7 @@ class cssr_bds(cssr):
         return sys, prn
 
     def decode_head(self, msg, i, st=-1):
+        """Decode BDS PPP message header."""
         self.tod, _, iodssr = bs.unpack_from('u17u4u2', msg, i)
         i += 23
 
@@ -107,10 +112,11 @@ class cssr_bds(cssr):
             if self.week >= 0:
                 self.time = bdt2gpst(bdt2time(self.week, self.tow))
 
-        head = {'uint': 0, 'mi': 0, 'iodssr': iodssr}
+        head = {'udi': 0, 'mi': 0, 'iodssr': iodssr}
         return head, i
 
     def add_gnss(self, mask, blen, gnss):
+        """Add decoded GNSS satellite mask entries."""
         prn, nsat = self.decode_mask(mask, blen)
         self.nsat_g[gnss] = nsat
         self.nsat_n += nsat
@@ -124,7 +130,7 @@ class cssr_bds(cssr):
             self.gnss_n.append(gnss)
 
     def decode_cssr_mask(self, msg, i):
-        """decode MT1 Mask message """
+        """Decode MT1 Mask message."""
         head, i = self.decode_head(msg, i, sCSSR.MASK)
 
         self.iodp = bs.unpack_from('u4', msg, i)[0]
@@ -179,6 +185,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_orb_sat(self, msg, i, inet, sat_n):
+        """Decode satellite orbit correction for BDS PPP."""
         slot, iodn, iodc = bs.unpack_from('u9u10u3', msg, i)
         i += 22
 
@@ -216,7 +223,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_orb(self, msg, i, inet=0):
-        """decode MT2 orbit + URA message """
+        """Decode MT2 orbit + URA message."""
         head, i = self.decode_head(msg, i)
         sat_n = np.array(self.sat_n)
 
@@ -233,7 +240,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_cbias(self, msg, i, inet=0):
-        """decode MT3 Code Bias Correction message """
+        """Decode MT3 Code Bias Correction message."""
         head, i = self.decode_head(msg, i)
         nsat = self.nsat_n
         self.flg_net = False
@@ -270,7 +277,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_clk_sat(self, msg, i, inet, sat):
-        """ decode clock correction for satellite """
+        """Decode clock correction for satellite."""
         iodc, dclk = bs.unpack_from('u3s{:d}'.format(self.dclk_blen), msg, i)
         i += 3+self.dclk_blen
 
@@ -285,7 +292,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_clk(self, msg, i, inet=0):
-        """decode MT4 Clock Correction message """
+        """Decode MT4 Clock Correction message."""
         head, i = self.decode_head(msg, i)
         if self.iodssr != head['iodssr']:
             return -1
@@ -311,7 +318,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_ura(self, msg, i):
-        """decode MT5 URA message """
+        """Decode MT5 URA message."""
         head, i = self.decode_head(msg, i)
         if self.iodssr != head['iodssr']:
             return -1
@@ -332,7 +339,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_comb1(self, msg, i, inet=0):
-        """decode MT6 combined message #1 """
+        """Decode MT6 combined message #1."""
         numc, numo = bs.unpack_from('u5u3', msg, i)
         sat_n = np.array(self.sat_n)
 
@@ -354,7 +361,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_comb2(self, msg, i, inet=0):
-        """decode MT7 combined message #2 """
+        """Decode MT7 combined message #2."""
         numc, numo = bs.unpack_from('u5u3', msg, i)
         sat_n = np.array(self.sat_n)
 
@@ -377,6 +384,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr(self, msg, i=0):
+        """Decode BDS PPP B2b message."""
         mt = bs.unpack_from('u6', msg, i)[0]
         i += 6
 

{cssrlib-1.2.1 → cssrlib-1.3.2}/src/cssrlib/cssr_has.py

@@ -1,5 +1,4 @@
-"""
-Galileo HAS correction data decoder
+"""Galileo HAS correction data decoder.
 
 [1] Galileo High Accuracy Service Signal-in-Space
   Interface Control Document (HAS SIS ICD), Issue 1.0, May 2022
@@ -15,7 +14,10 @@ from binascii import unhexlify
 
 
 class cssr_has(cssr):
+    """Class to decode Galileo HAS Compact SSR messages."""
+
     def __init__(self, foutname=None):
+        """Initialize Galileo HAS SSR decoder."""
         super().__init__(foutname)
         self.MAXNET = 1
         self.cssrmode = sCSSRTYPE.GAL_HAS_SIS
@@ -29,13 +31,14 @@ class cssr_has(cssr):
         self.pb_scl = 0.01  # cycles
 
     def sval(self, u, n, scl):
-        """ calculate signed value based on n-bit int, lsb """
+        """Calculate signed value based on n-bit int, lsb."""
         invalid = -2**(n-1)
         dnu = 2**(n-1)-1
         y = np.nan if u == invalid or u == dnu else u*scl
         return y
 
     def decode_cssr_clk_sub(self, msg, i=0):
+        """Decode Galileo HAS clock subset correction message."""
         head, i = self.decode_head(msg, i)
         self.flg_net = False
         if self.iodssr != head['iodssr']:
@@ -58,12 +61,11 @@ class cssr_has(cssr):
         return i
 
     def decode_head(self, msg, i, st=-1):
-        """ decode header part of HAS messages """
-
+        """Decode header part of HAS messages."""
         if st == sCSSR.MASK:
-            ui = 0
+            udi = 0
         else:
-            ui = bs.unpack_from('u4', msg, i)[0]
+            udi = bs.unpack_from('u4', msg, i)[0]
             i += 4
 
         if self.tow0 >= 0:
@@ -71,13 +73,12 @@ class cssr_has(cssr):
             if self.week >= 0:
                 self.time = gpst2time(self.week, self.tow)
 
-        head = {'uint': ui, 'mi': 0, 'iodssr': self.iodssr}
+        head = {'udi': udi, 'mi': 0, 'iodssr': self.iodssr}
         return head, i
 
     def decode_cssr(self, msg, i=0):
-        """ decode HAS messages """
-
-        # Galileo HAS-SIS only defines MT=1, MT=4073 is CSSR message used in the
+        """Decode HAS messages."""
+        # Galileo HAS-SIS only defines MT=1; MT=4073 is CSSR used in the
         # Galileo HAS test dataset
         #
         if self.msgtype != 1 and self.msgtype != 4073:
@@ -131,9 +132,10 @@ class cssr_has(cssr):
 
 
 class cnav_msg():
-    """ class to handle Galileo E6 (CNAV) message """
+    """Class to handle Galileo E6 (CNAV) message."""
 
     def __init__(self):
+        """Initialize Galileo C/NAV HAS message decoder."""
         self.mid_decoded = []
         self.mid_ = -1
         self.ms_ = -1
@@ -148,10 +150,11 @@ class cnav_msg():
         self.monlevel = 0
 
     def load_gmat(self, file_gm):
+        """Load Reed-Solomon generator matrix for Galileo HAS decoding."""
         self.gMat = np.genfromtxt(file_gm, dtype="u1", delimiter=",")
 
     def decode_has_header(self, buff, i):
-        """ decode header of HAS pages (obsoleted) """
+        """Decode header of HAS pages (obsoleted)."""
         if bs.unpack_from('u24', buff, i)[0] == 0xaf3bc3:
             return 0, 0, 0, 0, 0
 
@@ -160,7 +163,7 @@ class cnav_msg():
         return hass, mt, mid, ms, pid
 
     def decode_has_page(self, idx, has_pages, gMat, ms):
-        """ HPVRS decoding for RS(255,32,224) """
+        """Decode HPVRS Reed-Solomon RS(255,32,224) page."""
         HASmsg = bytes()
         k = len(idx)
         if k >= ms:
@@ -172,8 +175,7 @@ class cnav_msg():
         return HASmsg
 
     def decode_cnav(self, tow, vi):
-        """ decode Galileo CNAV message """
-
+        """Decode Galileo CNAV message."""
         HASmsg = None
 
         for vn in vi:

{cssrlib-1.2.1 → cssrlib-1.3.2}/src/cssrlib/cssr_mdc.py

@@ -1,9 +1,8 @@
-"""
-QZSS MADOCA-PPP correction data decoder
+"""QZSS MADOCA-PPP correction data decoder.
 
 [1] Quasi-Zenith Satellite System Interface Specification Multi-GNSS
     Advanced Orbit and Clock Augmentation - Precise Point Positioning
-    (IS-QZSS-MDC-004), 2025
+    (IS-QZSS-MDC-005draft), October, 2026
 
 """
 
@@ -11,10 +10,14 @@ import numpy as np
 import bitstruct as bs
 from cssrlib.cssrlib import cssr, sCSSRTYPE, sCSSR, local_corr, sCType
 from cssrlib.gnss import gpst2time, time2str, uGNSS, prn2sat, rCST, sat2id
+from cssrlib.gnss import pos2ecef, gtime_t
 
 
 class areaInfo():
+    """Class to store MADOCA-PPP service area information."""
+
     def __init__(self, sid, latr, lonr, p1=0, p2=0):
+        """Initialize MADOCA-PPP area information."""
         self.sid = sid
         self.latr = latr
         self.lonr = lonr
@@ -26,7 +29,10 @@ class areaInfo():
 
 
 class ionoCorr():
+    """Class to store MADOCA-PPP regional ionospheric corrections."""
+
     def __init__(self):
+        """Initialize MADOCA-PPP ionospheric correction data."""
         self.t0 = None
         self.qi = []
         self.iodssr = 0
@@ -36,7 +42,10 @@ class ionoCorr():
 
 
 class cssr_mdc(cssr):
+    """Class to decode QZSS MADOCA-PPP L6D/L6E messages."""
+
     def __init__(self, foutname=None):
+        """Initialize MADOCA-PPP L6 message decoder."""
         super().__init__(foutname)
 
         self.cssrmode = sCSSRTYPE.QZS_MADOCA
@@ -44,6 +53,8 @@ class cssr_mdc(cssr):
 
         self.pnt = {}
         self.ci = {}
+        self.reg_alrt = {}
+        self.reg_iod = {}
 
         self.area = -1
         self.reg = -1
@@ -51,6 +62,10 @@ class cssr_mdc(cssr):
         self.narea_t = {1: 8, 2: 16, 3: 5, 4: 1, 5: 8}
         self.MAXNET = np.sum(list(self.narea_t.values()))
 
+        self.tmax = {sCType.CLOCK: 60.0, sCType.ORBIT: 60.0,
+                     sCType.CBIAS: 60.0, sCType.PBIAS: 60.0,
+                     sCType.TROP: 120.0, sCType.STEC: 300.0}
+
         for inet in range(self.MAXNET+1):
             self.lc.append(local_corr())
             self.lc[inet].inet = inet
@@ -60,35 +75,40 @@ class cssr_mdc(cssr):
             self.lc[inet].t0 = {}
 
     def find_grid_index(self, pos):
-        """ find region/area   """
-
+        """Find region/area."""
         latd, lond = np.rad2deg(pos[0:2])
+        rr = pos2ecef(pos)
         self.reg = self.area = -1
+        min_dist = 1e5
 
         for reg in self.pnt.keys():
+            if self.reg_alrt.get(reg, 0) != 0:
+                continue
             for area in self.pnt[reg].keys():
                 p = self.pnt[reg][area]
+                ref_ecef = pos2ecef(
+                    np.array([p.latr*rCST.D2R, p.lonr*rCST.D2R, 0.0]))
+                dist = np.linalg.norm(rr - ref_ecef)*1e-3  # [km]
+                if dist > min_dist:
+                    continue
                 if p.sid == 0:  # rectangle shape:
                     if p.latr-p.lats <= latd and latd <= p.latr+p.lats and \
                             p.lonr-p.lons <= lond and lond <= p.lonr+p.lons:
                         self.reg = reg
                         self.area = area
-                        break
+                        min_dist = dist
                 elif p.sid == 1:  # circle
-                    dn = (latd - p.latr)*rCST.D2R*rCST.RE_WGS84
-                    de = (lond - p.lonr)*rCST.D2R * \
-                        rCST.RE_WGS84*np.cos(p.lonr*rCST.D2R)
-                    r = np.sqrt(dn**2+de**2)
-                    if r <= p.rng:
+                    if dist <= p.rng:
                         self.reg = reg
                         self.area = area
-                        break
+                        min_dist = dist
 
         inet = self.get_inet(self.reg, self.area)
+        self.inet = inet
         return inet
 
     def get_dpos(self, pos):
-        """ calculate position offset from reference """
+        """Calculate position offset from reference."""
         if self.reg < 0 or self.area < 0:
             print("get_dpos: region, area not defined.")
             return 0, 0
@@ -101,7 +121,7 @@ class cssr_mdc(cssr):
         return dlat, dlon
 
     def get_stec(self, dlat=0.0, dlon=0.0):
-        """ calculate STEC correction by interporation """
+        """Calculate STEC correction by interporation."""
         if self.inet < 0:
             print("get_stec: region, area not defined.")
             return np.zeros(self.nsat_n)
@@ -118,14 +138,19 @@ class cssr_mdc(cssr):
         return stec
 
     def decode_mdc_stec_area(self, buff, i=0):
-        """ decoder for MT1 - STEC Coverage Message """
+        """Decode MT1 STEC coverage message."""
         tow, uid, mi, iodssr = bs.unpack_from('u20u4u1u4', buff, i)
         i += 29
+        self.tow = tow
         self.tow0 = tow//3600*3600
+        if self.week >= 0:
+            self.time = gpst2time(self.week, self.tow)
         reg, alrt, len_, narea = bs.unpack_from('u8u1u16u5', buff, i)
         i += 30
         self.reg = reg
         self.alrt = alrt
+        self.reg_alrt[reg] = alrt
+        self.reg_iod[reg] = iodssr
 
         if reg not in self.pnt:
             self.pnt[reg] = {}
@@ -152,8 +177,7 @@ class cssr_mdc(cssr):
         return i
 
     def get_inet(self, reg, area):
-        """ region, area to inet conversion """
-
+        """Convert region and area ID to network index."""
         if reg < 0 or area < 0:
             return -1
 
@@ -168,9 +192,17 @@ class cssr_mdc(cssr):
         return inet
 
     def decode_mdc_stec_corr(self, buff, i=0):
-        """ decoder for MT2 - STEC Correction Message """
+        """Decode MT2 STEC correction message."""
         dtow, uid, mi, iodssr = bs.unpack_from('u12u4u1u4', buff, i)
         i += 21
+        if self.tow0 >= 0:
+            tow = self.tow0 + dtow
+            if self.tow >= 0 and tow < self.tow - 1800:
+                tow += 3600
+            self.tow = tow
+        t0 = gpst2time(self.week, self.tow) if self.week >= 0 else gtime_t()
+        if self.week >= 0:
+            self.time = t0
         reg, area, stype_ = bs.unpack_from('u8u5u2', buff, i)
         i += 15
 
@@ -184,7 +216,6 @@ class cssr_mdc(cssr):
 
         inet = self.get_inet(reg, area)
 
-        self.lc[inet].nsat_n = np.sum(nsat)
         ci = {}
         qi = {}
         stype = {}
@@ -193,39 +224,29 @@ class cssr_mdc(cssr):
         for gnss in range(5):
             sys = sys_t[gnss]
             for k in range(nsat[gnss]):
-                c01 = c10 = c11 = c02 = c20 = 0.0
-
-                prn, qi_, c00 = bs.unpack_from('u6u6s14', buff, i)
-                i += 26
+                prn, cls_, val = bs.unpack_from('u6u3u3', buff, i)
+                i += 12
 
                 if sys == uGNSS.QZS:
                     prn += 192
                 sat = prn2sat(sys, prn)
-                qi[sat] = qi_
+                ci_sat, i = self.decode_cssr_stec_coeff(buff, stype_, i)
+                if np.isnan(ci_sat).any():
+                    continue
+                qi[sat] = self.quality_idx(cls_, val)
                 stype[sat] = stype_
-
-                if stype_ > 0:
-                    c01, c10 = bs.unpack_from('s12s12', buff, i)
-                    i += 24
-                if stype_ > 1:
-                    c11 = bs.unpack_from('s10', buff, i)[0]
-                    i += 10
-                if stype_ > 2:
-                    c02, c20 = bs.unpack_from('s8s8', buff, i)
-                    i += 16
-
                 sat_ += [sat]
-                ci[sat] = np.array(
-                    [c00*0.05, c01*0.02, c10*0.02, c11*0.02,
-                     c02*5e-3, c20*5e-3])
+                ci[sat] = ci_sat
+                self.set_t0(inet, sat, sCType.STEC, t0)
                 j += 1
 
+        self.lc[inet].nsat_n = len(sat_)
+        self.lc[inet].flg_stec = 2
         self.lc[inet].stype = stype
         self.lc[inet].sat_n = sat_
         self.lc[inet].ci = ci
         self.lc[inet].stec_quality = qi
 
-        t0 = gpst2time(self.week, self.tow0+dtow)
         self.set_t0(inet, 0, sCType.STEC, t0)
         self.lc[inet].cstat |= (1 << sCType.STEC)
 
@@ -236,7 +257,7 @@ class cssr_mdc(cssr):
     def out_log(self):
         # if self.msgtype not in (1, 2):
         #    return super(cssr_mdc, self).out_log()
-
+        """Write decoded MADOCA-PPP STEC messages to log file."""
         sz_t = [1, 3, 4, 6]
 
         if self.time == -1:
@@ -276,9 +297,11 @@ class cssr_mdc(cssr):
         self.fh.flush()
 
     def decode_cssr(self, msg, i=0):
-        """decode Compact SSR message with MADOCA-PPP extension """
+        """Decode Compact SSR message with MADOCA-PPP extension."""
         df = {'msgtype': 4073}
         while df['msgtype'] in [1, 2, 4073]:
+            if len(msg)*8 < i + 16:
+                return 0
             df = bs.unpack_from_dict('u12u4', ['msgtype', 'subtype'], msg, i)
             i += 16
             if df['msgtype'] not in [1, 2, 4073]:

{cssrlib-1.2.1 → cssrlib-1.3.2}/src/cssrlib/cssr_pvs.py

@@ -1,11 +1,10 @@
-"""
-PVS (PPP via SouthPAN) correction data decoder
+"""PVS (PPP via SouthPAN) correction data decoder.
 
 [1] Service Definition Document for Open Services, SBAS-STN-0001,
-    Revision 02, December 2022
+    Revision 04, April 2026
 
 [2] Service Definition Document for Data Access Services, SBAS-STN-0002,
-    Revision 01, November 2023
+    Revision 02, April 2026
 
 For use of SSR corrections see section A 3.9.2 in
 
@@ -24,7 +23,7 @@ from binascii import unhexlify
 
 
 def decode_sinca_line(line):
-    """ SINCA (SISNeT compression algorithm) decoder  """
+    """Decode SINCA (SISNeT compression algorithm) line."""
     if line[0:4] != '*MSG':
         return None
     v = line.split(',')
@@ -47,7 +46,10 @@ def decode_sinca_line(line):
 
 
 class cssr_pvs(cssr):
+    """Class to decode SouthPAN PVS PPP messages."""
+
     def __init__(self, foutname=None):
+        """Initialize PVS PPP decoder."""
         super().__init__(foutname)
         self.MAXNET = 1
         self.cssrmode = sCSSRTYPE.PVS_PPP
@@ -67,6 +69,7 @@ class cssr_pvs(cssr):
         self.nav_mode = {uGNSS.GPS: 0, uGNSS.GAL: 1}
 
     def slot2sat(self, slot):
+        """Convert satellite slot index to satellite number."""
         sat = 0
         if slot >= 1 and slot <= 32:
             prn = slot
@@ -86,14 +89,16 @@ class cssr_pvs(cssr):
         return sat
 
     def adjust_time_week(self, time, time0):
+        """Adjust time-of-week for week rollover."""
         dt = timediff(time, time0)
-        if dt > rCST.DAY_SEC/2:
-            time = timeadd(time, -rCST.DAY_SEC)
-        elif dt < -rCST.DAY_SEC/2:
-            time = timeadd(time,  rCST.DAY_SEC)
+        if dt > rCST.HALFWEEK_SEC:
+            time = timeadd(time, -rCST.WEEK_SEC)
+        elif dt < -rCST.HALFWEEK_SEC:
+            time = timeadd(time,  rCST.WEEK_SEC)
         return time
 
     def check_validity(self, time):
+        """Check validity time of PVS SSR corrections."""
         for sat in self.sat_n:
             if timediff(time, self.lc[0].t0[sat][sCType.CLOCK]) > \
                     self.tmax[sCType.CLOCK]:
@@ -104,7 +109,7 @@ class cssr_pvs(cssr):
                 self.lc[0].dorb[sat] = np.zeros(3)
 
     def decode_cssr_orb(self, msg, i, inet=0):
-        """ Types 32 clock-ephemeris correction and covariance matrix """
+        """Decode Type 32 clock-ephemeris correction and covariance matrix."""
         slot, iodn = bs.unpack_from('u8u10', msg, i)
         i += 18
         dx, dy, dz, db = bs.unpack_from('s11s11s11s12', msg, i)
@@ -158,7 +163,7 @@ class cssr_pvs(cssr):
         return i
 
     def decode_cssr(self, msg, i=0):
-
+        """Decode PVS PPP DFMC SBAS message."""
         _, mt = bs.unpack_from('u4u6', msg, i)
         i += 10