diff --git a/Convert_Matlab/inst_eddies_v6.py b/Convert_Matlab/inst_eddies_v6.py index 11bd4cfd026b20126219817ce9b675a80083bd68..f8111189f7d39d1f7a300a34d2e238c4d3dda85b 100755 --- a/Convert_Matlab/inst_eddies_v6.py +++ b/Convert_Matlab/inst_eddies_v6.py @@ -42,7 +42,7 @@ def define_fields(writers): writers["max_speed"].field("days_1950", "N", 5) writers["max_speed"].field("eddy_index", "N", 5) -def write(writers, cell, X, Y, ADT, cyclone): +def write_to_shp(writers, cell, X, Y, ADT, cyclone): """Write to shapefiles eddies with a given orientation, at a given date. @@ -173,16 +173,16 @@ while True: == my_date for factory in factories: - write(factory["writers"], matlab_data[factory["cell_name"]], - matlab_data["X"], matlab_data["Y"], matlab_data["ADT"], - factory["cyclone"]) + write_to_shp(factory["writers"], matlab_data[factory["cell_name"]], + matlab_data["X"], matlab_data["Y"], matlab_data["ADT"], + factory["cyclone"]) os.remove("adt.mat") else: print("Missing file:", adt_file) for factory in factories: - write(factory["writers"], [], factory["cyclone"]) + write_to_shp(factory["writers"], [], factory["cyclone"]) my_date += datetime.timedelta(1) if my_date > final_date: break