gnss-js 1.0.0 → 1.1.0

package/dist/index.cjs

@@ -206,10 +206,13 @@ __export(src_exports, {
   phiBPSK: () => phiBPSK,
   rad2deg: () => rad2deg,
   readString: () => readString,
+  reportDecodeError: () => reportDecodeError,
   resetGloFreqCache: () => resetGloFreqCache,
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   selectEphemeris: () => selectEphemeris,
+  setGloFreqNumber: () => setGloFreqNumber,
+  setRtcm3DebugHandler: () => setRtcm3DebugHandler,
   systemCmp: () => systemCmp,
   systemName: () => systemName,
   updateStationMeta: () => updateStationMeta,
@@ -298,9 +301,22 @@ var SYS_SHORT = {
 var C_LIGHT = 299792458;
 var FREQ = {
   G: { "1": 157542e4, "2": 12276e5, "5": 117645e4 },
-  R: { "3": 1202025e3 },
-  E: { "1": 157542e4, "5": 117645e4, "6": 127875e4, "7": 120714e4, "8": 1191795e3 },
-  C: { "1": 157542e4, "2": 1561098e3, "5": 117645e4, "6": 126852e4, "7": 120714e4 },
+  // R bands 4/6 are the CDMA L1OC/L2OC center frequencies
+  R: { "3": 1202025e3, "4": 1600995e3, "6": 124806e4 },
+  E: {
+    "1": 157542e4,
+    "5": 117645e4,
+    "6": 127875e4,
+    "7": 120714e4,
+    "8": 1191795e3
+  },
+  C: {
+    "1": 157542e4,
+    "2": 1561098e3,
+    "5": 117645e4,
+    "6": 126852e4,
+    "7": 120714e4
+  },
   J: { "1": 157542e4, "2": 12276e5, "5": 117645e4, "6": 127875e4 },
   I: { "5": 117645e4, "9": 2492028e3 },
   S: { "1": 157542e4, "5": 117645e4 }
@@ -315,11 +331,28 @@ var BAND_LABELS = {
   S: { "1": "L1", "5": "L5" }
 };
 var DUAL_FREQ_PAIRS = {
-  G: [["1", "2"], ["1", "5"]],
-  R: [["1", "2"], ["1", "3"]],
-  E: [["1", "5"], ["1", "7"], ["1", "6"]],
-  C: [["2", "7"], ["2", "6"], ["1", "5"]],
-  J: [["1", "2"], ["1", "5"]],
+  G: [
+    ["1", "2"],
+    ["1", "5"]
+  ],
+  R: [
+    ["1", "2"],
+    ["1", "3"]
+  ],
+  E: [
+    ["1", "5"],
+    ["1", "7"],
+    ["1", "6"]
+  ],
+  C: [
+    ["2", "7"],
+    ["2", "6"],
+    ["1", "5"]
+  ],
+  J: [
+    ["1", "2"],
+    ["1", "5"]
+  ],
   I: [["5", "9"]],
   S: [["1", "5"]]
 };
@@ -363,11 +396,11 @@ function buildGloChannelMap(slots) {
   return map;
 }
 function gloFreq(gloChannels, prn, band) {
+  if (band === "3" || band === "4" || band === "6") return FREQ.R[band];
   const k = gloChannels[prn];
   if (k === void 0) return void 0;
   if (band === "1") return GLO_F1_BASE + k * GLO_F1_STEP;
   if (band === "2") return GLO_F2_BASE + k * GLO_F2_STEP;
-  if (band === "3") return GLO_F3;
   return void 0;
 }
 function getFreq(gloChannels, prn, band) {
@@ -964,7 +997,10 @@ function vincenty(lat1, lon1, lat2, lon2) {
     const cosLat2 = Math.cos(lat2), sinLat2 = Math.sin(lat2);
     const dLon = lon2 - lon1;
     sigma = Math.acos(
-      Math.min(1, Math.max(-1, sinLat1 * sinLat2 + cosLat1 * cosLat2 * Math.cos(dLon)))
+      Math.min(
+        1,
+        Math.max(-1, sinLat1 * sinLat2 + cosLat1 * cosLat2 * Math.cos(dLon))
+      )
     );
     const distance2 = WGS84_SEMI_MAJOR_AXIS * sigma;
     const ib = Math.atan2(
@@ -1162,10 +1198,20 @@ function parseCrxDataLine(line, ntype) {
     if (token.length >= 3 && token[1] === "&") {
       const arcOrder = parseInt(token[0]);
       const value = parseInt(token.substring(2));
-      fields.push({ empty: false, init: true, arcOrder, value: isNaN(value) ? 0 : value });
+      fields.push({
+        empty: false,
+        init: true,
+        arcOrder,
+        value: isNaN(value) ? 0 : value
+      });
     } else {
       const value = parseInt(token);
-      fields.push({ empty: false, init: false, arcOrder: -1, value: isNaN(value) ? 0 : value });
+      fields.push({
+        empty: false,
+        init: false,
+        arcOrder: -1,
+        value: isNaN(value) ? 0 : value
+      });
     }
   }
   const flagStr = pos < line.length ? line.substring(pos) : "";
@@ -1175,7 +1221,10 @@ function crxDecompress(prev, field) {
   if (field.init) {
     const a = new Array(Math.max(field.arcOrder + 1, 1)).fill(0);
     a[0] = field.value;
-    return { state: { accum: a, order: 0, arcOrder: field.arcOrder }, result: field.value };
+    return {
+      state: { accum: a, order: 0, arcOrder: field.arcOrder },
+      result: field.value
+    };
   }
   if (!prev) {
     return { state: { accum: [0], order: 0, arcOrder: 0 }, result: 0 };
@@ -1188,7 +1237,10 @@ function crxDecompress(prev, field) {
   for (let k = order; k >= 1; k--) {
     accum[k - 1] += accum[k];
   }
-  return { state: { accum, order, arcOrder: prev.arcOrder }, result: accum[0] };
+  return {
+    state: { accum, order, arcOrder: prev.arcOrder },
+    result: accum[0]
+  };
 }
 
 // src/rinex/format.ts
@@ -1210,7 +1262,9 @@ var CHUNK_SIZE = 2 * 1024 * 1024;
 var yieldToMain = () => new Promise((r) => setTimeout(r, 0));
 var noYield = () => Promise.resolve();
 var SYSTEM_ORDER = ["G", "R", "E", "C", "J", "I", "S"];
-var SYSTEM_RANK = Object.fromEntries(SYSTEM_ORDER.map((s, i) => [s, i]));
+var SYSTEM_RANK = Object.fromEntries(
+  SYSTEM_ORDER.map((s, i) => [s, i])
+);
 function systemCmp(a, b) {
   const ra = SYSTEM_RANK[a.charAt(0)] ?? 99;
   const rb = SYSTEM_RANK[b.charAt(0)] ?? 99;
@@ -1396,7 +1450,11 @@ function parseEpochLine3(line) {
   const wholeSec = Math.floor(sc);
   const ms = Math.round((sc - wholeSec) * 1e3);
   const time = Date.UTC(yr, mo - 1, dy, hr, mn, wholeSec, ms);
-  return { time, flag: isNaN(flag) ? 0 : flag, numSats: isNaN(numSats) ? 0 : numSats };
+  return {
+    time,
+    flag: isNaN(flag) ? 0 : flag,
+    numSats: isNaN(numSats) ? 0 : numSats
+  };
 }
 function parseEpochLine2(line) {
   const yr2 = parseInt(line.substring(1, 3));
@@ -1418,7 +1476,12 @@ function parseEpochLine2(line) {
     const id = satPart.substring(i, i + 3).trim();
     if (id) satIds.push(id);
   }
-  return { time, flag: isNaN(flag) ? 0 : flag, numSats: isNaN(numSats) ? 0 : numSats, satIds };
+  return {
+    time,
+    flag: isNaN(flag) ? 0 : flag,
+    numSats: isNaN(numSats) ? 0 : numSats,
+    satIds
+  };
 }
 async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode) {
   const skipEpochs = !!workerMode;
@@ -1471,16 +1534,19 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
     if (!skipEpochs) {
       const snrBySys = {};
       for (const [sys, vals] of Object.entries(snrPerSystemAccum)) {
-        if (vals.length > 0) snrBySys[sys] = vals.reduce((a, b) => a + b, 0) / vals.length;
+        if (vals.length > 0)
+          snrBySys[sys] = vals.reduce((a, b) => a + b, 0) / vals.length;
       }
       const meanSnr = snrValues.length > 0 ? snrValues.reduce((a, b) => a + b, 0) / snrValues.length : null;
       const snrPerSat = {};
       for (const [prn, vals] of Object.entries(snrPerSatAccum)) {
-        if (vals.length > 0) snrPerSat[prn] = vals.reduce((a, b) => a + b, 0) / vals.length;
+        if (vals.length > 0)
+          snrPerSat[prn] = vals.reduce((a, b) => a + b, 0) / vals.length;
       }
       const snrPerSatBand = {};
       for (const [key, vals] of Object.entries(snrPerSatBandAccum)) {
-        if (vals.length > 0) snrPerSatBand[key] = vals.reduce((a, b) => a + b, 0) / vals.length;
+        if (vals.length > 0)
+          snrPerSatBand[key] = vals.reduce((a, b) => a + b, 0) / vals.length;
       }
       epochs.push({
         time: epochInfo.time,
@@ -1525,7 +1591,8 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
     const codes = header.obsTypes[sys] ?? [];
     if (onSatObs && epochInfo) {
       const values = new Array(codes.length);
-      for (let i = 0; i < codes.length; i++) values[i] = readObsValue(obsLine, i);
+      for (let i = 0; i < codes.length; i++)
+        values[i] = readObsValue(obsLine, i);
       if (!skipEpochs) {
         for (const { idx, band } of getSnrIndices(sys)) {
           const val = values[idx] ?? null;
@@ -1762,7 +1829,10 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
           v2SnrBandInfo = snrIndicesWithBandV2(header.obsTypes);
           v2SatIndex = 0;
           v2CurrentSatLine = 0;
-          const continuationLines = Math.max(0, Math.ceil((info.numSats - 12) / 12));
+          const continuationLines = Math.max(
+            0,
+            Math.ceil((info.numSats - 12) / 12)
+          );
           if (continuationLines > 0) {
             v2ContinuationSatsRemaining = continuationLines;
           } else {
@@ -1802,7 +1872,9 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
   const magic = new Uint8Array(await file.slice(0, 2).arrayBuffer());
   const isGz = magic[0] === 31 && magic[1] === 139;
   if (magic[0] === 31 && magic[1] === 157) {
-    throw new Error("Unix compress (.Z) files are not supported. Please decompress first (e.g. uncompress or gzip -d).");
+    throw new Error(
+      "Unix compress (.Z) files are not supported. Please decompress first (e.g. uncompress or gzip -d)."
+    );
   }
   if (isGz) {
     let bytesRead = 0;
@@ -1818,7 +1890,9 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
       if (done) break;
       bytesRead += value.byteLength;
       processChunkText(decoder.decode(value, { stream: true }), false);
-      onProgress?.(Math.min(99, Math.round(bytesRead / (file.size * 4) * 100)));
+      onProgress?.(
+        Math.min(99, Math.round(bytesRead / (file.size * 4) * 100))
+      );
       await yield_();
     }
     processChunkText(decoder.decode(), true);
@@ -1828,7 +1902,10 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
       const end = Math.min(offset + CHUNK_SIZE, file.size);
       const slice = file.slice(offset, end);
       const arrayBuf = await slice.arrayBuffer();
-      processChunkText(decoder.decode(arrayBuf, { stream: end < file.size }), false);
+      processChunkText(
+        decoder.decode(arrayBuf, { stream: end < file.size }),
+        false
+      );
       onProgress?.(Math.min(99, Math.round(end / file.size * 100)));
       await yield_();
     }
@@ -1836,19 +1913,27 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
       processChunkText("", true);
     }
   }
-  if (epochInfo && satLinesRemaining === 0 && crxPhase === "epoch") finishEpoch();
-  if (!header) throw new Error("No valid RINEX header found (missing END OF HEADER).");
+  if (epochInfo && satLinesRemaining === 0 && crxPhase === "epoch")
+    finishEpoch();
+  if (!header)
+    throw new Error("No valid RINEX header found (missing END OF HEADER).");
   onProgress?.(100);
-  return { header, epochs, stats: computeStats(header, epochs, satellitesSeen) };
+  return {
+    header,
+    epochs,
+    stats: computeStats(header, epochs, satellitesSeen)
+  };
 }
 function computeStats(header, epochs, satellitesSeen) {
   const n = epochs.length;
   const startTime = n > 0 ? new Date(epochs[0].time) : header.timeOfFirstObs;
   const endTime = n > 0 ? new Date(epochs[n - 1].time) : header.timeOfLastObs;
   let duration = null;
-  if (startTime && endTime) duration = (endTime.getTime() - startTime.getTime()) / 1e3;
+  if (startTime && endTime)
+    duration = (endTime.getTime() - startTime.getTime()) / 1e3;
   let interval = header.interval;
-  if (interval === null && n >= 2) interval = (epochs[1].time - epochs[0].time) / 1e3;
+  if (interval === null && n >= 2)
+    interval = (epochs[1].time - epochs[0].time) / 1e3;
   const uniqueSatsPerSystem = {};
   let totalUnique = 0;
   const systems = [];
@@ -1903,10 +1988,18 @@ var WarningAccumulator = class {
       );
     }
     if (!header.antType) {
-      this.add("MISSING_ANT_TYPE", "info", "Missing antenna type (ANT # / TYPE)");
+      this.add(
+        "MISSING_ANT_TYPE",
+        "info",
+        "Missing antenna type (ANT # / TYPE)"
+      );
     }
     if (!header.receiverType) {
-      this.add("MISSING_REC_TYPE", "info", "Missing receiver type (REC # / TYPE / VERS)");
+      this.add(
+        "MISSING_REC_TYPE",
+        "info",
+        "Missing receiver type (REC # / TYPE / VERS)"
+      );
     }
     if (!header.markerName) {
       this.add("MISSING_MARKER", "info", "Missing marker name (MARKER NAME)");
@@ -1915,7 +2008,11 @@ var WarningAccumulator = class {
       this.add("MISSING_INTERVAL", "info", "No observation interval in header");
     }
     if (Object.keys(header.obsTypes).length === 0) {
-      this.add("NO_OBS_TYPES", "error", "No observation types defined in header");
+      this.add(
+        "NO_OBS_TYPES",
+        "error",
+        "No observation types defined in header"
+      );
     }
   }
   setInterval(intervalMs) {
@@ -1982,7 +2079,12 @@ var WarningAccumulator = class {
     };
     const range = ranges[sys];
     if (!range) {
-      this.add("UNKNOWN_SYSTEM", "info", "Unknown satellite system letter", prn);
+      this.add(
+        "UNKNOWN_SYSTEM",
+        "info",
+        "Unknown satellite system letter",
+        prn
+      );
     } else if (num < range[0] || num > range[1]) {
       this.add(
         "PRN_OUT_OF_RANGE",
@@ -2072,7 +2174,15 @@ var R4_MSG_LINES = {
   EOP: 3,
   ION: 3
 };
-var SUPPORTED_EPH_MSGS = /* @__PURE__ */ new Set(["LNAV", "INAV", "FNAV", "D1", "D2", "SBAS", "FDMA"]);
+var SUPPORTED_EPH_MSGS = /* @__PURE__ */ new Set([
+  "LNAV",
+  "INAV",
+  "FNAV",
+  "D1",
+  "D2",
+  "SBAS",
+  "FDMA"
+]);
 function parseFloat19(s) {
   return parseFloat(s.trim().replace(/[dD]/g, "E"));
 }
@@ -2279,9 +2389,25 @@ function parseNavFile(text) {
       }
       if (dataLines2.length === numDataLines2) {
         if (sys2 === "R" || sys2 === "S") {
-          ephemerides.push(buildStateVectorEphemeris(sys2, prn2, parsed.date, parsed.values, dataLines2));
+          ephemerides.push(
+            buildStateVectorEphemeris(
+              sys2,
+              prn2,
+              parsed.date,
+              parsed.values,
+              dataLines2
+            )
+          );
         } else if (sys2 === "G" || sys2 === "E" || sys2 === "C" || sys2 === "J" || sys2 === "I") {
-          ephemerides.push(buildKeplerEphemeris(sys2, prn2, parsed.date, parsed.values, dataLines2));
+          ephemerides.push(
+            buildKeplerEphemeris(
+              sys2,
+              prn2,
+              parsed.date,
+              parsed.values,
+              dataLines2
+            )
+          );
         }
       }
       i++;
@@ -2317,9 +2443,13 @@ function parseNavFile(text) {
     }
     if (dataLines.length === numDataLines) {
       if (sys === "R" || sys === "S") {
-        ephemerides.push(buildStateVectorEphemeris(sys, prn, date, epochVals, dataLines));
+        ephemerides.push(
+          buildStateVectorEphemeris(sys, prn, date, epochVals, dataLines)
+        );
       } else if (sys === "G" || sys === "E" || sys === "C" || sys === "J" || sys === "I") {
-        ephemerides.push(buildKeplerEphemeris(sys, prn, date, epochVals, dataLines));
+        ephemerides.push(
+          buildKeplerEphemeris(sys, prn, date, epochVals, dataLines)
+        );
       }
     }
     i++;
@@ -2430,27 +2560,31 @@ function writeRinexNav(nav) {
 }
 
 // src/rtcm3/decoder.ts
+var debugHandler = null;
+function setRtcm3DebugHandler(fn) {
+  debugHandler = fn;
+}
+function reportDecodeError(context, error) {
+  if (debugHandler) debugHandler(context, error);
+}
 var BitReader = class {
   data;
   bitPos = 0;
   constructor(data) {
     this.data = data;
   }
-  /** Read `n` bits as unsigned integer (max 32). */
+  /** Read `n` bits as unsigned integer (max 53). */
   readU(n) {
     if (n === 0) return 0;
     let val = 0;
     for (let i = 0; i < n; i++) {
       const byteIdx = this.bitPos + i >> 3;
       const bitIdx = 7 - (this.bitPos + i & 7);
-      if (byteIdx < this.data.length) {
-        val = val << 1 | this.data[byteIdx] >> bitIdx & 1;
-      } else {
-        val <<= 1;
-      }
+      const bit = byteIdx < this.data.length ? this.data[byteIdx] >> bitIdx & 1 : 0;
+      val = val * 2 + bit;
     }
     this.bitPos += n;
-    return val >>> 0;
+    return val;
   }
   /** Read `n` bits as signed (two's complement) integer. */
   readS(n) {
@@ -2461,9 +2595,9 @@ var BitReader = class {
   /** Read `n` bits as sign-magnitude integer (MSB = sign, rest = magnitude). */
   readSM(n) {
     const val = this.readU(n);
-    const sign = val >> n - 1;
-    const mag = val & (1 << n - 1) - 1;
-    return sign ? -mag : mag;
+    const half = 2 ** (n - 1);
+    const mag = val % half;
+    return val >= half ? -mag : mag;
   }
   /** Skip `n` bits. */
   skip(n) {
@@ -2967,31 +3101,32 @@ function decodeEphemeris(frame) {
       default:
         return null;
     }
-  } catch {
+  } catch (err) {
+    reportDecodeError("ephemeris", err);
     return null;
   }
 }
 
 // src/rtcm3/msm.ts
-var GPS_L1 = 157542e4;
-var GPS_L2 = 12276e5;
-var GPS_L5 = 117645e4;
-var GLO_L1a = 1600995e3;
-var GLO_L2a = 124806e4;
-var GLO_L3 = 1202025e3;
-var GAL_E1 = 157542e4;
-var GAL_E5a = 117645e4;
-var GAL_E5b = 120714e4;
-var GAL_E5 = 1191795e3;
-var GAL_E6 = 127875e4;
-var BDS_B1 = 1561098e3;
-var BDS_B3 = 126852e4;
-var BDS_B2 = 120714e4;
-var BDS_B1C = 157542e4;
-var BDS_B2a = 117645e4;
-var BDS_B2b = 120714e4;
-var QZSS_L6 = 127875e4;
-var NAVIC_S = 2492028e3;
+var GPS_L1 = FREQ.G["1"];
+var GPS_L2 = FREQ.G["2"];
+var GPS_L5 = FREQ.G["5"];
+var GLO_L1a = FREQ.R["4"];
+var GLO_L2a = FREQ.R["6"];
+var GLO_L3 = FREQ.R["3"];
+var GAL_E1 = FREQ.E["1"];
+var GAL_E5a = FREQ.E["5"];
+var GAL_E5b = FREQ.E["7"];
+var GAL_E5 = FREQ.E["8"];
+var GAL_E6 = FREQ.E["6"];
+var BDS_B1 = FREQ.C["2"];
+var BDS_B3 = FREQ.C["6"];
+var BDS_B2 = FREQ.C["7"];
+var BDS_B1C = FREQ.C["1"];
+var BDS_B2a = FREQ.C["5"];
+var BDS_B2b = FREQ.C["7"];
+var QZSS_L6 = FREQ.J["6"];
+var NAVIC_S = FREQ.I["9"];
 var EMPTY = { code: "", freq: 0 };
 function sd(code, freq) {
   return { code, freq };
@@ -3243,43 +3378,6 @@ function signalTable(sys) {
       return GPS_SIGNALS;
   }
 }
-var BitReader2 = class {
-  data;
-  pos;
-  constructor(data, startBit = 0) {
-    this.data = data;
-    this.pos = startBit;
-  }
-  /** Read unsigned value of numBits bits */
-  u(numBits) {
-    let val = 0;
-    for (let i = 0; i < numBits; i++) {
-      const byteIdx = this.pos + i >>> 3;
-      const bitIdx = 7 - (this.pos + i & 7);
-      if (byteIdx < this.data.length) {
-        val = val << 1 | this.data[byteIdx] >>> bitIdx & 1;
-      } else {
-        val <<= 1;
-      }
-    }
-    this.pos += numBits;
-    return val;
-  }
-  /** Read signed value (two's complement) */
-  s(numBits) {
-    const raw = this.u(numBits);
-    const sign = 1 << numBits - 1;
-    return (raw & sign - 1) - (raw & sign);
-  }
-  /** Skip bits */
-  skip(n) {
-    this.pos += n;
-  }
-  /** Current bit position */
-  get position() {
-    return this.pos;
-  }
-};
 var gloFreqNum = new Int8Array(64).fill(-128);
 function gloWavelength(satIdx, freqMarker) {
   const k = gloFreqNum[satIdx];
@@ -3314,28 +3412,9 @@ function lockTimeSec(msmType, lti) {
   if (variant <= 5) {
     return LTI_TABLE_4BIT[lti] ?? 0;
   }
-  if (lti < 64) return lti;
-  if (lti < 96) return (lti - 64) * 2 + 64;
-  if (lti < 128) return (lti - 96) * 4 + 128;
-  if (lti < 160) return (lti - 128) * 8 + 256;
-  if (lti < 192) return (lti - 160) * 16 + 512;
-  if (lti < 224) return (lti - 192) * 32 + 1024;
-  if (lti < 256) return (lti - 224) * 64 + 2048;
-  if (lti < 288) return (lti - 256) * 128 + 4096;
-  if (lti < 320) return (lti - 288) * 256 + 8192;
-  if (lti < 352) return (lti - 320) * 512 + 16384;
-  if (lti < 384) return (lti - 352) * 1024 + 32768;
-  if (lti < 416) return (lti - 384) * 2048 + 65536;
-  if (lti < 448) return (lti - 416) * 4096 + 131072;
-  if (lti < 480) return (lti - 448) * 8192 + 262144;
-  if (lti < 512) return (lti - 480) * 16384 + 524288;
-  if (lti < 544) return (lti - 512) * 32768 + 1048576;
-  if (lti < 576) return (lti - 544) * 65536 + 2097152;
-  if (lti < 608) return (lti - 576) * 131072 + 4194304;
-  if (lti < 640) return (lti - 608) * 262144 + 8388608;
-  if (lti < 672) return (lti - 640) * 524288 + 16777216;
-  if (lti < 704) return (lti - 672) * 1048576 + 33554432;
-  return (lti - 704) * 2097152 + 67108864;
+  if (lti < 64) return lti / 1e3;
+  const n = Math.min(20, Math.floor((lti - 64) / 32));
+  return ((lti - 64 - 32 * n) * 2 ** (n + 1) + 2 ** (n + 6)) / 1e3;
 }
 function msmSystem(type) {
   if (type >= 1071 && type <= 1077) return "G";
@@ -3361,10 +3440,10 @@ function decodeMsmFull(frame) {
   const payload = frame.payload;
   if (!payload || payload.length < 10) return null;
   try {
-    const bits = new BitReader2(payload);
+    const bits = new BitReader(payload);
     bits.skip(12);
     bits.skip(12);
-    const epochMs = bits.u(30);
+    const epochMs = bits.readU(30);
     bits.skip(1);
     bits.skip(3);
     bits.skip(7);
@@ -3372,8 +3451,8 @@ function decodeMsmFull(frame) {
     bits.skip(2);
     bits.skip(1);
     bits.skip(3);
-    const satMaskHi = bits.u(32) >>> 0;
-    const satMaskLo = bits.u(32) >>> 0;
+    const satMaskHi = bits.readU(32) >>> 0;
+    const satMaskLo = bits.readU(32) >>> 0;
     const satIndices = [];
     for (let i = 0; i < 32; i++) {
       if (satMaskHi & 1 << 31 - i) satIndices.push(i + 1);
@@ -3383,7 +3462,7 @@ function decodeMsmFull(frame) {
     }
     const numSat = satIndices.length;
     if (numSat === 0) return null;
-    const sigMask = bits.u(32) >>> 0;
+    const sigMask = bits.readU(32) >>> 0;
     const sigIndices = [];
     for (let i = 0; i < 32; i++) {
       if (sigMask & 1 << 31 - i) sigIndices.push(i);
@@ -3393,7 +3472,7 @@ function decodeMsmFull(frame) {
     const numCells = numSat * numSig;
     const cellMask = [];
     for (let i = 0; i < numCells; i++) {
-      cellMask.push(bits.u(1) === 1);
+      cellMask.push(bits.readU(1) === 1);
     }
     let activeCells = 0;
     for (const c of cellMask) if (c) activeCells++;
@@ -3403,14 +3482,13 @@ function decodeMsmFull(frame) {
     const extsat = new Int8Array(numSat);
     const rdop = new Float64Array(numSat);
     if (variant === 4 || variant === 6) {
-      for (let j = 0; j < numSat; j++) rrint[j] = bits.u(8);
-      for (let j = 0; j < numSat; j++) extsat[j] = bits.u(4);
-      for (let j = 0; j < numSat; j++) rrmod[j] = bits.u(10) / 1024;
+      for (let j = 0; j < numSat; j++) rrint[j] = bits.readU(8);
+      for (let j = 0; j < numSat; j++) rrmod[j] = bits.readU(10) / 1024;
     } else {
-      for (let j = 0; j < numSat; j++) rrint[j] = bits.u(8);
-      for (let j = 0; j < numSat; j++) extsat[j] = bits.u(4);
-      for (let j = 0; j < numSat; j++) rrmod[j] = bits.u(10) / 1024;
-      for (let j = 0; j < numSat; j++) rdop[j] = bits.s(14) * 1;
+      for (let j = 0; j < numSat; j++) rrint[j] = bits.readU(8);
+      for (let j = 0; j < numSat; j++) extsat[j] = bits.readU(4);
+      for (let j = 0; j < numSat; j++) rrmod[j] = bits.readU(10) / 1024;
+      for (let j = 0; j < numSat; j++) rdop[j] = bits.readS(14) * 1;
     }
     if (sys === "R") {
       for (let j = 0; j < numSat; j++) {
@@ -3430,31 +3508,31 @@ function decodeMsmFull(frame) {
     cp.fill(-1e30);
     dop.fill(-1e30);
     if (variant === 4) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(15) / (1 << 24);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(22) / (1 << 29);
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(4);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(6);
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(6);
     } else if (variant === 5) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(15) / (1 << 24);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(22) / (1 << 29);
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(4);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(6);
-      for (let i = 0; i < activeCells; i++) dop[i] = bits.s(15) * 1e-4;
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(6);
+      for (let i = 0; i < activeCells; i++) dop[i] = bits.readS(15) * 1e-4;
     } else if (variant === 6) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(20) / (1 << 29);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(24) / (1 << 31);
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(10);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(10) / 16;
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(20) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(24) / 2 ** 31;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(10);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(10) / 16;
     } else {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(20) / (1 << 29);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(24) / (1 << 31);
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(10);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(10) / 16;
-      for (let i = 0; i < activeCells; i++) dop[i] = bits.s(15) * 1e-4;
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(20) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(24) / 2 ** 31;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(10);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(10) / 16;
+      for (let i = 0; i < activeCells; i++) dop[i] = bits.readS(15) * 1e-4;
     }
     const sigTable = signalTable(sys);
     const observations = [];
@@ -3488,11 +3566,11 @@ function decodeMsmFull(frame) {
           wavelength
         };
         const psrVal = psr[cellIdx];
-        if (psrVal > -1 / (1 << 10)) {
+        if (psrVal > -(2 ** -10)) {
           signal.pseudorange = psrVal * C_LIGHT / 1e3 + roughRange_m;
         }
         const cpVal = cp[cellIdx];
-        if (cpVal > -1 / (1 << 8) && wavelength > 0) {
+        if (cpVal > -(2 ** -8) && wavelength > 0) {
           signal.phase = cpVal * C_LIGHT / 1e3 / wavelength + roughRange_m / wavelength;
         }
         if (variant === 5 || variant === 7) {
@@ -3520,12 +3598,13 @@ function decodeMsmFull(frame) {
       system: sys,
       observations
     };
-  } catch {
+  } catch (err) {
+    reportDecodeError("msm", err);
     return null;
   }
 }
-var GPS_EPOCH = Date.UTC(1980, 0, 6);
-var MS_PER_WEEK = 6048e5;
+var GPS_EPOCH = START_GPS_TIME.getTime();
+var MS_PER_WEEK = MILLISECONDS_IN_WEEK;
 function msmEpochToDate(sys, epochMs, refTime = /* @__PURE__ */ new Date()) {
   const refMs = refTime.getTime();
   if (sys === "R") {
@@ -3540,11 +3619,12 @@ function msmEpochToDate(sys, epochMs, refTime = /* @__PURE__ */ new Date()) {
     return new Date(t2);
   }
   if (sys === "C") {
-    const BDS_EPOCH = Date.UTC(2006, 0, 1) - 14e3;
+    const BDS_EPOCH = START_BDS_TIME.getTime();
     const weeksSinceEpoch2 = Math.floor((refMs - BDS_EPOCH) / MS_PER_WEEK);
     let t2 = BDS_EPOCH + weeksSinceEpoch2 * MS_PER_WEEK + epochMs;
     if (t2 - refMs > MS_PER_WEEK / 2) t2 -= MS_PER_WEEK;
     else if (refMs - t2 > MS_PER_WEEK / 2) t2 += MS_PER_WEEK;
+    t2 -= 18e3;
     return new Date(t2);
   }
   const weeksSinceEpoch = Math.floor((refMs - GPS_EPOCH) / MS_PER_WEEK);
@@ -3554,6 +3634,11 @@ function msmEpochToDate(sys, epochMs, refTime = /* @__PURE__ */ new Date()) {
   t -= 18e3;
   return new Date(t);
 }
+function setGloFreqNumber(slot, k) {
+  if (slot >= 1 && slot <= 64 && k >= -7 && k <= 13) {
+    gloFreqNum[slot - 1] = k;
+  }
+}
 function resetGloFreqCache() {
   gloFreqNum.fill(-128);
 }
@@ -3687,7 +3772,8 @@ function updateStationMeta(meta, frame) {
       default:
         return false;
     }
-  } catch {
+  } catch (err) {
+    reportDecodeError("station", err);
     return false;
   }
 }
@@ -3875,6 +3961,9 @@ function updateStreamStats(stats, frames, rawBytes) {
     const eph = decodeEphemeris(frame);
     if (eph) {
       stats.ephemerides.set(eph.prn, eph);
+      if (eph.freqChannel !== void 0 && eph.prn.startsWith("R")) {
+        setGloFreqNumber(parseInt(eph.prn.slice(1), 10), eph.freqChannel);
+      }
     }
     updateStationMeta(stats.stationMeta, frame);
   }
@@ -3894,6 +3983,8 @@ var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
 var OMEGA_E = 72921151467e-15;
+var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
+var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
 var J2_GLO = 108263e-8;
 var TWO_PI = 2 * Math.PI;
@@ -3919,10 +4010,7 @@ function keplerPosition(eph, t) {
   }
   const sinE = Math.sin(Ek);
   const cosE = Math.cos(Ek);
-  const vk = Math.atan2(
-    Math.sqrt(1 - eph.e * eph.e) * sinE,
-    cosE - eph.e
-  );
+  const vk = Math.atan2(Math.sqrt(1 - eph.e * eph.e) * sinE, cosE - eph.e);
   const phik = vk + eph.omega;
   const sin2phi = Math.sin(2 * phik);
   const cos2phi = Math.cos(2 * phik);
@@ -3934,7 +4022,7 @@ function keplerPosition(eph, t) {
   const ik = eph.i0 + dik + eph.idot * tk;
   const xp = rk * Math.cos(uk);
   const yp = rk * Math.sin(uk);
-  const isBdsGeo = eph.system === "C" && Math.abs(eph.i0) < 0.1 && a > 4e7;
+  const isBdsGeo = eph.system === "C" && Math.abs(eph.i0) < BDS_GEO_MAX_INCLINATION_RAD && a > BDS_GEO_MIN_SEMIMAJOR_AXIS_M;
   if (isBdsGeo) {
     const omegak2 = eph.omega0 + eph.omegaDot * tk - OMEGA_E * eph.toe;
     const cosO2 = Math.cos(omegak2);
@@ -4085,7 +4173,10 @@ function invert4x4(m) {
   let det = a[0] * inv[0] + a[1] * inv[4] + a[2] * inv[8] + a[3] * inv[12];
   if (Math.abs(det) < 1e-20) return null;
   det = 1 / det;
-  const result = Array.from({ length: 4 }, () => new Array(4).fill(0));
+  const result = Array.from(
+    { length: 4 },
+    () => new Array(4).fill(0)
+  );
   for (let i = 0; i < 4; i++) {
     for (let j = 0; j < 4; j++) {
       result[i][j] = inv[i * 4 + j] * det;
@@ -4112,9 +4203,10 @@ function computeSatPosition(eph, timeMs) {
   if (eph.system === "R" || eph.system === "S") {
     return glonassPosition(eph, timeMs / 1e3);
   }
-  const GPS_EPOCH2 = Date.UTC(1980, 0, 6);
+  const GPS_EPOCH2 = START_GPS_TIME.getTime();
   const gpsSeconds = (timeMs - GPS_EPOCH2) / 1e3;
-  const tow = gpsSeconds % (7 * 86400);
+  let tow = gpsSeconds % (7 * 86400);
+  if (eph.system === "C") tow = (tow - 14 + 604800) % 604800;
   return keplerPosition(eph, tow);
 }
 function navTimesFromEph(ephs) {
@@ -4196,16 +4288,22 @@ function selectBest(ephs, timeMs) {
   }
   return best;
 }
-var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
-var BDS_EPOCH_MS = Date.UTC(2006, 0, 1);
+var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var BDS_EPOCH_MS = START_BDS_TIME.getTime();
 var GAL_EPOCH_MS = GPS_EPOCH_MS;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
-    if (info.x === void 0 || info.y === void 0 || info.z === void 0) return null;
-    if (info.vx === void 0 || info.vy === void 0 || info.vz === void 0) return null;
+    if (info.x === void 0 || info.y === void 0 || info.z === void 0)
+      return null;
+    if (info.vx === void 0 || info.vy === void 0 || info.vz === void 0)
+      return null;
     const now = new Date(info.lastReceived);
-    const utcMidnight = Date.UTC(now.getUTCFullYear(), now.getUTCMonth(), now.getUTCDate());
+    const utcMidnight = Date.UTC(
+      now.getUTCFullYear(),
+      now.getUTCMonth(),
+      now.getUTCDate()
+    );
     const tbSec = (info.tb ?? 0) * 900;
     const moscowOffset = 3 * 3600;
     const tocMs = utcMidnight + (tbSec - moscowOffset) * 1e3;
@@ -4230,7 +4328,11 @@ function ephInfoToEphemeris(info) {
       freqNum: info.freqChannel ?? 0
     };
   }
-  if (info.sqrtA === void 0 || info.eccentricity === void 0 || info.inclination === void 0 || info.omega0 === void 0 || info.argPerigee === void 0 || info.meanAnomaly === void 0 || info.toe === void 0 || info.week === void 0) return null;
+  if (info.sqrtA === void 0 || info.eccentricity === void 0 || info.inclination === void 0 || info.omega0 === void 0 || info.argPerigee === void 0 || info.meanAnomaly === void 0 || info.toe === void 0 || info.week === void 0)
+    return null;
+  if (info.eccentricity < 0 || info.eccentricity >= 1 || info.sqrtA < 3e3 || info.sqrtA > 9e3 || info.toe < 0 || info.toe > 604800) {
+    return null;
+  }
   let epochMs;
   const tocSec = info.toc ?? info.toe;
   if (sys === "C") {
@@ -4406,14 +4508,18 @@ async function ntripFetch(proxyUrl, path, headers, signal) {
   } catch (err) {
     if (err instanceof Error && err.name === "AbortError") throw err;
     const message = err instanceof Error ? err.message : "The service may be temporarily unavailable.";
-    throw new Error(`Could not reach the NTRIP proxy: ${message}`, { cause: err });
+    throw new Error(`Could not reach the NTRIP proxy: ${message}`, {
+      cause: err
+    });
   }
 }
 async function fetchSourcetable(proxyUrl, info, signal) {
   const headers = ntripHeaders(info);
   const res = await ntripFetch(proxyUrl, "/", headers, signal);
   if (res.status === 401) {
-    throw new Error("Authentication required. Please provide valid credentials.");
+    throw new Error(
+      "Authentication required. Please provide valid credentials."
+    );
   }
   if (!res.ok) {
     throw new Error(`Caster returned ${res.status} ${res.statusText}`);
@@ -4425,7 +4531,12 @@ async function connectToMountpoint(proxyUrl, info, signal) {
   const headers = ntripHeaders(info);
   const controller = new AbortController();
   const combinedSignal = signal ? AbortSignal.any([signal, controller.signal]) : controller.signal;
-  const res = await ntripFetch(proxyUrl, `/${info.mountpoint}`, headers, combinedSignal);
+  const res = await ntripFetch(
+    proxyUrl,
+    `/${info.mountpoint}`,
+    headers,
+    combinedSignal
+  );
   if (res.status === 401) {
     throw new Error("Authentication required for this mountpoint.");
   }
@@ -4436,7 +4547,9 @@ async function connectToMountpoint(proxyUrl, info, signal) {
     throw new Error(`Caster returned ${res.status} ${res.statusText}`);
   }
   if (!res.body) {
-    throw new Error("No response body \u2014 streaming not supported by this environment.");
+    throw new Error(
+      "No response body \u2014 streaming not supported by this environment."
+    );
   }
   return {
     reader: res.body.getReader(),
@@ -4530,8 +4643,18 @@ var MultipathAccumulator = class {
         time: t,
         mp: arc.rawMp[i] - mean
       }));
-      const rms = Math.sqrt(points.reduce((s, p) => s + p.mp * p.mp, 0) / points.length);
-      this.results.push({ prn, system: sys, band, refBand, label: label2, points, rms });
+      const rms = Math.sqrt(
+        points.reduce((s, p) => s + p.mp * p.mp, 0) / points.length
+      );
+      this.results.push({
+        prn,
+        system: sys,
+        band,
+        refBand,
+        label: label2,
+        points,
+        rms
+      });
     }
     state.arc = { times: [], rawMp: [] };
   }
@@ -4685,10 +4808,20 @@ var CycleSlipAccumulator = class {
           const sigKey = `${sys}:${signal}`;
           this.countEpoch(sigKey);
           const stddev = mwSmoothStddev(prevSmooth);
-          const threshold = Math.max(MW_MIN_THRESHOLD, MW_SIGMA_FACTOR * stddev);
+          const threshold = Math.max(
+            MW_MIN_THRESHOLD,
+            MW_SIGMA_FACTOR * stddev
+          );
           const deviation = Math.abs(mw - prevSmooth.mean);
           if (deviation > threshold) {
-            this.recordSlip(sigKey, time, prn, signal, deviation, /* @__PURE__ */ new Set([b1, b2]));
+            this.recordSlip(
+              sigKey,
+              time,
+              prn,
+              signal,
+              deviation,
+              /* @__PURE__ */ new Set([b1, b2])
+            );
             mwSmooth.set(pairKey, mwSmoothInit(mw));
           } else {
             mwSmooth.set(pairKey, mwSmoothUpdate(prevSmooth, mw));
@@ -4743,7 +4876,13 @@ var CycleSlipAccumulator = class {
         }
       }
     }
-    this.prev.set(prn, { time, mwSmooth, gf: gfMap, phase: phaseM, code: codeM });
+    this.prev.set(prn, {
+      time,
+      mwSmooth,
+      gf: gfMap,
+      phase: phaseM,
+      code: codeM
+    });
   };
   finalize() {
     const satSlipCounts = {};
@@ -4751,7 +4890,10 @@ var CycleSlipAccumulator = class {
       satSlipCounts[ev.prn] = (satSlipCounts[ev.prn] ?? 0) + 1;
     }
     const signalStats = [];
-    const allKeys = /* @__PURE__ */ new Set([...this.signalEpochs.keys(), ...this.signalSlips.keys()]);
+    const allKeys = /* @__PURE__ */ new Set([
+      ...this.signalEpochs.keys(),
+      ...this.signalSlips.keys()
+    ]);
     for (const sigKey of allKeys) {
       const [sys, ...rest] = sigKey.split(":");
       const signal = rest.join(":");
@@ -5321,32 +5463,33 @@ function parseNmeaFile(content) {
 }
 
 // src/signals/definitions.ts
-var FREQ_GPS_L1 = 1575.42;
-var FREQ_GPS_L2 = 1227.6;
-var FREQ_GPS_L5 = 1176.45;
-var FREQ_GLO_L1 = 1602;
-var FREQ_GLO_L2 = 1246;
-var DELTA_GLO_L1 = 0.5625;
-var DELTA_GLO_L2 = 0.4375;
-var FREQ_GLO_L1OC = 1600.995;
-var FREQ_GLO_L2OC = 1248.06;
-var FREQ_GLO_L3OC = 1202.025;
-var FREQ_GAL_E1 = 1575.42;
-var FREQ_GAL_E5 = 1191.795;
-var FREQ_GAL_E5a = 1176.45;
-var FREQ_GAL_E5b = 1207.14;
-var FREQ_GAL_E6 = 1278.75;
-var FREQ_BDS_B1I = 1561.098;
-var FREQ_BDS_B2I = 1207.14;
-var FREQ_BDS_B3I = 1268.52;
-var FREQ_BDS_B1A = 1575.42;
-var FREQ_BDS_B1C = 1575.42;
-var FREQ_BDS_B3A = 1268.52;
-var FREQ_QZS_L1 = 1575.42;
-var FREQ_QZS_L2 = 1227.6;
-var FREQ_QZS_L5 = 1176.45;
-var FREQ_QZS_L6 = 1278.75;
-var FREQ_NAVIC_L5 = 1176.45;
+var MHz = 1e6;
+var FREQ_GPS_L1 = FREQ.G["1"] / MHz;
+var FREQ_GPS_L2 = FREQ.G["2"] / MHz;
+var FREQ_GPS_L5 = FREQ.G["5"] / MHz;
+var FREQ_GLO_L1 = GLO_F1_BASE / MHz;
+var FREQ_GLO_L2 = GLO_F2_BASE / MHz;
+var DELTA_GLO_L1 = GLO_F1_STEP / MHz;
+var DELTA_GLO_L2 = GLO_F2_STEP / MHz;
+var FREQ_GLO_L1OC = FREQ.R["4"] / MHz;
+var FREQ_GLO_L2OC = FREQ.R["6"] / MHz;
+var FREQ_GLO_L3OC = FREQ.R["3"] / MHz;
+var FREQ_GAL_E1 = FREQ.E["1"] / MHz;
+var FREQ_GAL_E5 = FREQ.E["8"] / MHz;
+var FREQ_GAL_E5a = FREQ.E["5"] / MHz;
+var FREQ_GAL_E5b = FREQ.E["7"] / MHz;
+var FREQ_GAL_E6 = FREQ.E["6"] / MHz;
+var FREQ_BDS_B1I = FREQ.C["2"] / MHz;
+var FREQ_BDS_B2I = FREQ.C["7"] / MHz;
+var FREQ_BDS_B3I = FREQ.C["6"] / MHz;
+var FREQ_BDS_B1A = FREQ.C["1"] / MHz;
+var FREQ_BDS_B1C = FREQ.C["1"] / MHz;
+var FREQ_BDS_B3A = FREQ.C["6"] / MHz;
+var FREQ_QZS_L1 = FREQ.J["1"] / MHz;
+var FREQ_QZS_L2 = FREQ.J["2"] / MHz;
+var FREQ_QZS_L5 = FREQ.J["5"] / MHz;
+var FREQ_QZS_L6 = FREQ.J["6"] / MHz;
+var FREQ_NAVIC_L5 = FREQ.I["5"] / MHz;
 var GPS_COLOR = "#4ade80";
 var GPS_M_COLOR = "#ef4444";
 var GPS_P_COLOR = "#dbc51f";
@@ -6272,10 +6415,13 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   phiBPSK,
   rad2deg,
   readString,
+  reportDecodeError,
   resetGloFreqCache,
   rhumbLine,
   rtcm3Constellation,
   selectEphemeris,
+  setGloFreqNumber,
+  setRtcm3DebugHandler,
   systemCmp,
   systemName,
   updateStationMeta,