diff --git a/Inst_eddies/Analysis/eddy_dump.py b/Inst_eddies/Analysis/eddy_dump.py index bca355615243bc675ee8045bd861ce8ce2a9c55e..5f7720c1700e739470883d8dd732127948edb2b5 100755 --- a/Inst_eddies/Analysis/eddy_dump.py +++ b/Inst_eddies/Analysis/eddy_dump.py @@ -26,6 +26,7 @@ SHPC = util_eddies.SHPC_class(args.directory, def_orient=args.orientation) if args.ishape: ishape = args.ishape + i_slice = 0 else: if args.node: with open("e_overestim.txt") as f: @@ -37,7 +38,8 @@ else: eddy_index = int(reply[1]) assert eddy_index >= 1 - ishape = SHPC.comp_ishape(date, eddy_index, 0, args.orientation) + i_slice = SHPC.get_slice(date) + ishape = SHPC.comp_ishape(date, eddy_index, i_slice, args.orientation) print("ishape =", ishape) @@ -45,9 +47,9 @@ for layer in ["extremum", "outermost_contour", "max_speed_contour", "centroid"]: print() try: - shapeRec = SHPC.get_reader(0, args.orientation, layer).shapeRecord( - ishape - ) + shapeRec = SHPC.get_reader( + i_slice, args.orientation, layer + ).shapeRecord(ishape) except AttributeError: pass else: