Merge branch 'develop' into enh/net-thrust

Author: kevin-alcaniz

CHANGELOG.md

@@ -32,18 +32,36 @@ Attention: The newest changes should be on top -->
 
 ### Added
 
+
+### Changed
+
+
+### Fixed
+
+
+## v1.9.0 - 2025-03-24
+
+### Added
+
 - ENH: Parallel mode for monte-carlo simulations 2 [#768](https://github.com/RocketPy-Team/RocketPy/pull/768)
 - DOC: ASTRA Flight Example [#770](https://github.com/RocketPy-Team/RocketPy/pull/770)
+- ENH: Add Eccentricity to Stochastic Simulations [#792](https://github.com/RocketPy-Team/RocketPy/pull/792)
+- ENH: Introduce the StochasticAirBrakes class [#785](https://github.com/RocketPy-Team/RocketPy/pull/785)
 
 ### Changed
 
+- DEP: Remove Pending Deprecations and Add Warnings Where Needed [#794](https://github.com/RocketPy-Team/RocketPy/pull/794)
 - DOCS: reshape docs (closes #659) [#781](https://github.com/RocketPy-Team/RocketPy/pull/781)
 - MNT: EmptyMotor class inherits from Motor(ABC) [#779](https://github.com/RocketPy-Team/RocketPy/pull/779)
 
 ### Fixed
 
 - BUG: do not allow drawing rockets with no aerodynamic surface [#774](https://github.com/RocketPy-Team/RocketPy/pull/774)
 - BUG: update flight simulation logic to include burn start time [#778](https://github.com/RocketPy-Team/RocketPy/pull/778)
+- BUG: fixes get_instance_attributes for Flight objects containing a Rocket object without rail buttons [#786](https://github.com/RocketPy-Team/RocketPy/pull/786)
+- BUG: fixed AGL altitude print for parachutes with lag [#788](https://github.com/RocketPy-Team/RocketPy/pull/788)
+- BUG: fix the wind velocity factors usage and better visualization of uniform distributions in Stochastic Classes [#783](https://github.com/RocketPy-Team/RocketPy/pull/783)
+
 
 ## [v1.8.0] - 2025-01-20
 

Makefile

@@ -21,6 +21,7 @@ install:
 	$(PYTHON) -m pip install --upgrade pip
 	pip install -r requirements.txt
 	pip install -r requirements-optional.txt
+	pip install -r requirements-tests.txt
 	pip install -e .
 
 format:

docs/conf.py

@@ -27,7 +27,7 @@
 author = "RocketPy Team"
 
 # The full version, including alpha/beta/rc tags
-release = "1.8.0"
+release = "1.9.0"
 
 
 # -- General configuration ---------------------------------------------------

docs/user/installation.rst

@@ -19,7 +19,7 @@ If you want to choose a specific version to guarantee compatibility, you may ins
 
 .. code-block:: shell
 
-    pip install rocketpy==1.8.0
+    pip install rocketpy==1.9.0
 
 
 Optional Installation Method: ``conda``

pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "rocketpy"
-version = "1.8.0"
+version = "1.9.0"
 description="Advanced 6-DOF trajectory simulation for High-Power Rocketry."
 dynamic = ["dependencies"]
 readme = "README.md"

rocketpy/__init__.py

@@ -44,6 +44,7 @@
 from .sensors import Accelerometer, Barometer, GnssReceiver, Gyroscope
 from .simulation import Flight, MonteCarlo
 from .stochastic import (
+    StochasticAirBrakes,
     StochasticEllipticalFins,
     StochasticEnvironment,
     StochasticFlight,

rocketpy/environment/environment.py

@@ -1653,38 +1653,6 @@ def process_wyoming_sounding(self, file):  # pylint: disable=too-many-statements
         # Save maximum expected height
         self._max_expected_height = data_array[-1, 1]
 
-    def process_noaaruc_sounding(self, file):  # pylint: disable=too-many-statements
-        """Import and process the upper air sounding data from `NOAA
-        Ruc Soundings` database (https://rucsoundings.noaa.gov/) given as
-        ASCII GSD format pages passed by its url to the file parameter. Sets
-        pressure, temperature, wind-u, wind-v profiles and surface elevation.
-
-        Parameters
-        ----------
-        file : string
-            URL of an upper air sounding data output from `NOAA Ruc Soundings`
-            in ASCII GSD format.
-
-            Example:
-
-            https://rucsoundings.noaa.gov/get_raobs.cgi?data_source=RAOB&latest=latest&start_year=2019&start_month_name=Feb&start_mday=5&start_hour=12&start_min=0&n_hrs=1.0&fcst_len=shortest&airport=83779&text=Ascii%20text%20%28GSD%20format%29&hydrometeors=false&start=latest
-
-
-        See also
-        --------
-        This method is deprecated and will be fully deleted in version 1.8.0.
-
-        Returns
-        -------
-        None
-        """
-        warnings.warn(
-            "NOAA RUC models are no longer available. "
-            "This method is deprecated and will be fully deleted in version 1.8.0.",
-            DeprecationWarning,
-        )
-        return file
-
     def process_forecast_reanalysis(self, file, dictionary):  # pylint: disable=too-many-locals,too-many-statements
         """Import and process atmospheric data from weather forecasts
         and reanalysis given as ``netCDF`` or ``OPeNDAP`` files.
@@ -2259,26 +2227,6 @@ def select_ensemble_member(self, member=0):
         self.calculate_speed_of_sound_profile()
         self.calculate_dynamic_viscosity()
 
-    def load_international_standard_atmosphere(self):  # pragma: no cover
-        """Defines the pressure and temperature profile functions set
-        by `ISO 2533` for the International Standard atmosphere and saves
-        them as ``Environment.pressure_ISA`` and ``Environment.temperature_ISA``.
-
-        Notes
-        -----
-        This method is **deprecated** and will be removed in version 1.6.0. You
-        can access :meth:`rocketpy.Environment.pressure_ISA` and
-        :meth:`rocketpy.Environment.temperature_ISA` directly without the need
-        to call this method.
-        """
-        warnings.warn(
-            "load_international_standard_atmosphere() is deprecated in version "
-            "1.5.0 and will be removed in version 1.7.0. This method is no longer "
-            "needed as the International Standard Atmosphere is already calculated "
-            "when the Environment object is created.",
-            DeprecationWarning,
-        )
-
     @funcify_method("Height Above Sea Level (m)", "Pressure (Pa)", "spline", "natural")
     def pressure_ISA(self):
         """Pressure, in Pa, as a function of height above sea level as defined
@@ -2619,6 +2567,11 @@ def geodesic_to_utm(
         EW : string
             Returns "W" for western hemisphere and "E" for eastern hemisphere
         """
+        warnings.warn(
+            "This function is deprecated and will be removed in v1.10.0. "
+            "Please use the new method `tools.geodesic_to_utm` instead.",
+            DeprecationWarning,
+        )
         return geodesic_to_utm_tools(lat, lon, semi_major_axis, flattening)
 
     @staticmethod
@@ -2656,6 +2609,11 @@ def utm_to_geodesic(
         lon : float
             latitude of the analyzed point
         """
+        warnings.warn(
+            "This function is deprecated and will be removed in v1.10.0. "
+            "Please use the new method `tools.utm_to_geodesic` instead.",
+            DeprecationWarning,
+        )
         return utm_to_geodesic_tools(x, y, utm_zone, hemis, semi_major_axis, flattening)
 
     @staticmethod

rocketpy/environment/fetchers.py

@@ -5,7 +5,6 @@
 
 import re
 import time
-import warnings
 from datetime import datetime, timedelta, timezone
 
 import netCDF4
@@ -328,33 +327,6 @@ def fetch_wyoming_sounding(file):
     return response
 
 
-@exponential_backoff(max_attempts=5, base_delay=2, max_delay=60)
-def fetch_noaaruc_sounding(file):  # pragma: no cover
-    """Fetches sounding data from a specified file using the NOAA RUC soundings.
-
-    Parameters
-    ----------
-    file : str
-        The URL of the file to fetch.
-
-    Returns
-    -------
-    str
-        The content of the fetched file.
-
-    Raises
-    ------
-    ImportError
-        If unable to load the specified file or the file content is too short.
-    """
-    warnings.warn(
-        "The NOAA RUC soundings are deprecated since September 30th, 2024. "
-        "This method will be removed in version 1.8.0.",
-        DeprecationWarning,
-    )
-    return file
-
-
 @exponential_backoff(max_attempts=5, base_delay=2, max_delay=60)
 def fetch_gefs_ensemble():
     """Fetches the latest GEFS (Global Ensemble Forecast System) dataset from

rocketpy/environment/tools.py

@@ -443,9 +443,44 @@ def get_interval_date_from_time_array(time_array, units=None):
 
 
 def geodesic_to_utm(lat, lon, semi_major_axis=6378137.0, flattening=1 / 298.257223563):  # pylint: disable=too-many-locals,too-many-statements
-    # NOTE: already documented in the Environment class.
-    # TODO: deprecated the static method from the environment class, use only this one.
+    """Function which converts geodetic coordinates, i.e. lat/lon, to UTM
+    projection coordinates. Can be used only for latitudes between -80.00°
+    and 84.00°
 
+    Parameters
+    ----------
+    lat : float
+        The latitude coordinates of the point of analysis, must be contained
+        between -80.00° and 84.00°
+    lon : float
+        The longitude coordinates of the point of analysis, must be
+        contained between -180.00° and 180.00°
+    semi_major_axis : float
+        The semi-major axis of the ellipsoid used to represent the Earth,
+        must be given in meters (default is 6,378,137.0 m, which corresponds
+        to the WGS84 ellipsoid)
+    flattening : float
+        The flattening of the ellipsoid used to represent the Earth, usually
+        between 1/250 and 1/150 (default is 1/298.257223563, which
+        corresponds to the WGS84 ellipsoid)
+
+    Returns
+    -------
+    x : float
+        East coordinate, always positive
+    y : float
+        North coordinate, always positive
+    utm_zone : int
+        The number of the UTM zone of the point of analysis, can vary
+        between 1 and 60
+    utm_letter : string
+        The letter of the UTM zone of the point of analysis, can vary
+        between C and X, omitting the letters "I" and "O"
+    hemis : string
+        Returns "S" for southern hemisphere and "N" for Northern hemisphere
+    EW : string
+        Returns "W" for western hemisphere and "E" for eastern hemisphere
+    """
     # Calculate the central meridian of UTM zone
     if lon != 0:
         signal = lon / abs(lon)
@@ -529,9 +564,37 @@ def geodesic_to_utm(lat, lon, semi_major_axis=6378137.0, flattening=1 / 298.2572
 def utm_to_geodesic(  # pylint: disable=too-many-locals,too-many-statements
     x, y, utm_zone, hemis, semi_major_axis=6378137.0, flattening=1 / 298.257223563
 ):
-    # NOTE: already documented in the Environment class.
-    # TODO: deprecate the static method from the environment class, use only this one.
+    """Function to convert UTM coordinates to geodesic coordinates
+    (i.e. latitude and longitude).
 
+    Parameters
+    ----------
+    x : float
+        East UTM coordinate in meters
+    y : float
+        North UTM coordinate in meters
+    utm_zone : int
+        The number of the UTM zone of the point of analysis, can vary
+        between 1 and 60
+    hemis : string
+        Equals to "S" for southern hemisphere and "N" for Northern
+        hemisphere
+    semi_major_axis : float
+        The semi-major axis of the ellipsoid used to represent the Earth,
+        must be given in meters (default is 6,378,137.0 m, which corresponds
+        to the WGS84 ellipsoid)
+    flattening : float
+        The flattening of the ellipsoid used to represent the Earth, usually
+        between 1/250 and 1/150 (default is 1/298.257223563, which
+        corresponds to the WGS84 ellipsoid)
+
+    Returns
+    -------
+    lat : float
+        latitude of the analyzed point
+    lon : float
+        latitude of the analyzed point
+    """
     if hemis == "N":
         y = y + 10000000