code: def _geoid_radius(latitude: float) -> float:
      
      lat = deg2rad(latitude)
      return sqrt(1/(cos(lat) ** 2 / Rmax_WGS84 ** 2 + sin(lat) ** 2 / Rmin_WGS84 ** 2))
docs: Calculates the GEOID radius at a given latitude
  
      Parameters
      ----------
      latitude : float
          Latitude (degrees)
  
      Returns
      -------
      R : float
          GEOID Radius (meters)