diff --git a/src/core/system_parameters/gnss_almanac.cc b/src/core/system_parameters/gnss_almanac.cc
index 2c609a53d..923416832 100644
--- a/src/core/system_parameters/gnss_almanac.cc
+++ b/src/core/system_parameters/gnss_almanac.cc
@@ -63,7 +63,7 @@ double Gnss_Almanac::predicted_doppler(double rx_time_s,
     // Position in EFEF
     const std::vector<double> pos_rx = {(v + h) * cosp * cosl, (v + h) * cosp * sinl, (v * (1.0 - e2) + h) * sinp};
 
-    // Velovity in EFEF
+    // Velocity in EFEF
     const double t = cosp * vu - sinp * vn;
     const std::vector<double> vel_rx = {cosl * t - sinl * ve, sinl * t + cosl * ve, sinp * vu + cosp * vn};
 
diff --git a/src/core/system_parameters/gnss_almanac.h b/src/core/system_parameters/gnss_almanac.h
index a8e28aeeb..a8c265a80 100644
--- a/src/core/system_parameters/gnss_almanac.h
+++ b/src/core/system_parameters/gnss_almanac.h
@@ -38,7 +38,8 @@ public:
      */
     Gnss_Almanac() = default;
 
-    /* \brief Computes prediction of the Doppler shift for a given time and receiver's position and velocity.
+    /*!
+     * \brief Computes prediction of the Doppler shift for a given time and receiver's position and velocity.
      * \f[
      * f_{d} = - \mathbf{v} \frac{\mathbf{x}^{T}}{\left| \mathbf{x} \right| } \frac{f_{L}}{c}
      * \f]
@@ -72,6 +73,9 @@ public:
         double vu,
         int band) const;
 
+    /*!
+     * \brief Computes satellite Position and Velocity, in ECEF, for a given time (expressed in seconds of week)
+     */
     void satellitePosVelComputation(double transmitTime, std::array<double, 7>& pos_vel_dtr) const;
 
     uint32_t PRN{};     //!< SV PRN NUMBER
diff --git a/src/core/system_parameters/gnss_ephemeris.cc b/src/core/system_parameters/gnss_ephemeris.cc
index 9505242ca..d95d14a28 100644
--- a/src/core/system_parameters/gnss_ephemeris.cc
+++ b/src/core/system_parameters/gnss_ephemeris.cc
@@ -59,7 +59,7 @@ double Gnss_Ephemeris::predicted_doppler(double rx_time_s,
     // Position in EFEF
     const std::vector<double> pos_rx = {(v + h) * cosp * cosl, (v + h) * cosp * sinl, (v * (1.0 - e2) + h) * sinp};
 
-    // Velovity in EFEF
+    // Velocity in EFEF
     const double t = cosp * vu - sinp * vn;
     const std::vector<double> vel_rx = {cosl * t - sinl * ve, sinl * t + cosl * ve, sinp * vu + cosp * vn};