Bug fix: do not try to distinguish orientation

An SHPC now contains a single orientation (since after v0.12).

Inst_eddies/Analysis/stat.py

@@ -9,9 +9,7 @@ parser = argparse.ArgumentParser()
 parser.add_argument("directory", help = "containing the shapefiles")
 args = parser.parse_args()
 
-n_valid_cycl = 0
 n_valid = 0
-n_minima = 0
 n_radius4 = np.zeros(3, dtype = int)
 n_outer = 0
 n_points_valid = 0
@@ -26,7 +24,6 @@ with shapefile.Reader(path.join(args.directory, "extremum")) as extremum, \
     
     for rec_extr, shape_rec_outer in zip(extremum.iterRecords(),
                                          outermost_contour):
-        if rec_extr.cyclone == 1: n_minima += 1
         l = len(shape_rec_outer.shape.points)
         n_points += l
 
@@ -34,7 +31,6 @@ with shapefile.Reader(path.join(args.directory, "extremum")) as extremum, \
             n_valid += 1
             n_points_valid += l
             
-        if rec_extr.cyclone == 1 and rec_extr.valid == 1: n_valid_cycl += 1
         if shape_rec_outer.record.r_eq_area >= 0: n_outer += 1
         
         if shape_rec_outer.record.radius4 == 0:
@@ -55,14 +51,10 @@ with shapefile.Reader(path.join(args.directory, "extremum")) as extremum, \
             n_valid_speed += 1
 
 print("number of extrema of SSH:", n_extr)
-print("number of minima of SSH:", n_minima)
-print("number of maxima of SSH:", n_extr - n_minima)
 print()
 
 print("number of valid eddies (outermost contour found with sufficient area):",
       n_valid)
-print("number of valid cyclones:", n_valid_cycl)
-print("number of valid anticyclones:", n_valid - n_valid_cycl)
 print()
 
 print("number of extrema with radius4 == 0 (no outermost contour or "