Merge pull request #352 from volodia99/rel-v0.16

REL: prep release 0.16

Author: volodia99

CITATION.cff

@@ -8,6 +8,6 @@ authors:
   given-names: "Clément"
   orcid: "https://orcid.org/0000-0001-8629-7068"
 title: "nonos"
-version: 0.15.1
+version: 0.16
 date-released: 2024-04-08
 url: "https://github.com/volodia99/nonos"

nonos/__version__.py

@@ -1 +1 @@
-__version__ = "0.15.1"
+__version__ = "0.16"

nonos/api/analysis.py

@@ -787,8 +787,8 @@ def latitudinal_projection(self, theta=None, *, operation_name=None) -> "GasFiel
             km = find_nearest(self.coords.theta, self.coords.theta.min())
             kp = find_nearest(self.coords.theta, self.coords.theta.max())
             if theta is not None:
-                km = find_nearest(self.coords.theta, np.pi / 2 + theta)
-                kp = find_nearest(self.coords.theta, np.pi / 2 - theta)
+                kp = find_nearest(self.coords.theta, np.pi / 2 + theta)
+                km = find_nearest(self.coords.theta, np.pi / 2 - theta)
             ret_data = (
                 np.sum(
                     (

tests/test_load_data.py

@@ -119,6 +119,14 @@ def test_npy_radial_at_r(test_data_dir, tmp_path):
     assert list(dsnpy.keys()) == ["RHO"]
 
 
+def test_latitudinal_projection(test_data_dir):
+    os.chdir(test_data_dir / "idefix_spherical_planet3d")
+    theta = np.arctan(5 * 0.1)
+    gflpt = GasDataSet(500)["RHO"].latitudinal_projection(theta)
+    gflp = GasDataSet(500)["RHO"].latitudinal_projection()
+    np.testing.assert_equal(gflpt.data, gflp.data)
+
+
 def test_radial_average_interval_vmin_vmax(test_data_dir):
     os.chdir(test_data_dir / "idefix_spherical_planet3d")
     with pytest.raises(