diff --git a/Inst_eddies/Analysis/stat.py b/Inst_eddies/Analysis/stat.py index 55423673c5b7e9bebac1122479dd56ad8816d7d7..d2d95ebf1799b315b6c117093981c96bf7c3c133 100755 --- a/Inst_eddies/Analysis/stat.py +++ b/Inst_eddies/Analysis/stat.py @@ -32,15 +32,16 @@ with shapefile.Reader(path.join(args.directory, "extremum")) as extremum, \ n_points_valid += l if shape_rec_outer.record.r_eq_area >= 0: n_outer += 1 - - if shape_rec_outer.record.radius4 == 0: - n_radius4[0] += 1 - elif shape_rec_outer.record.radius4 == 1: - n_radius4[1] += 1 - elif shape_rec_outer.record.radius4 >= 2: - n_radius4[2] += 1 - - if rec_extr.speed == 1e4: + + if "radius4" in shape_rec_outer.record: + if shape_rec_outer.record.radius4 == 0: + n_radius4[0] += 1 + elif shape_rec_outer.record.radius4 == 1: + n_radius4[1] += 1 + elif shape_rec_outer.record.radius4 >= 2: + n_radius4[2] += 1 + + if rec_extr.speed == 1e4 and "radius4" in shape_rec_outer.record: if shape_rec_outer.record.radius4 == 0: n_missing_speed[0] += 1 elif shape_rec_outer.record.radius4 == 1: