Support numpy>=1.24 (#318)

* Fix value error with numpy>=1.24

* Fix value error with numpy>=1.24

Author: dnerini

pysteps/noise/fftgenerators.py

@@ -507,8 +507,8 @@ def initialize_nonparam_2d_ssft_filter(field, **kwargs):
     # SSFT algorithm
 
     # prepare indices
-    idxi = np.zeros((2, 1), dtype=int)
-    idxj = np.zeros((2, 1), dtype=int)
+    idxi = np.zeros(2, dtype=int)
+    idxj = np.zeros(2, dtype=int)
 
     # number of windows
     num_windows_y = np.ceil(float(dim_y) / win_size[0]).astype(int)
@@ -793,8 +793,8 @@ def generate_noise_2d_ssft_filter(F, randstate=None, seed=None, **kwargs):
     cN = np.zeros(dim)
     sM = np.zeros(dim)
 
-    idxi = np.zeros((2, 1), dtype=int)
-    idxj = np.zeros((2, 1), dtype=int)
+    idxi = np.zeros(2, dtype=int)
+    idxj = np.zeros(2, dtype=int)
 
     # get the window size
     win_size = (float(dim_y) / F.shape[0], float(dim_x) / F.shape[1])

pysteps/nowcasts/sseps.py

@@ -466,8 +466,8 @@ def estimator(precip, parsglob=None, idxm=None, idxn=None):
     M, N = precip.shape[1:]
     n_windows_M = np.ceil(1.0 * M / win_size[0]).astype(int)
     n_windows_N = np.ceil(1.0 * N / win_size[1]).astype(int)
-    idxm = np.zeros((2, 1), dtype=int)
-    idxn = np.zeros((2, 1), dtype=int)
+    idxm = np.zeros(2, dtype=int)
+    idxn = np.zeros(2, dtype=int)
 
     if measure_time:
         starttime = time.time()
@@ -677,8 +677,8 @@ def worker(j):
 
             # then the local steps
             if n_windows_M > 1 or n_windows_N > 1:
-                idxm = np.zeros((2, 1), dtype=int)
-                idxn = np.zeros((2, 1), dtype=int)
+                idxm = np.zeros(2, dtype=int)
+                idxn = np.zeros(2, dtype=int)
                 precip_l = np.zeros((M, N), dtype=float)
                 M_s = np.zeros((M, N), dtype=float)
                 for m in range(n_windows_M):