diff --git a/Convert_Matlab/inst_eddies_v6.py b/Convert_Matlab/inst_eddies_v6.py index dd6ef7350bca52f09ad4a969291893abec021362..b4a125577d9ff8ede9bba569ccdd60fe00d9e6d9 100755 --- a/Convert_Matlab/inst_eddies_v6.py +++ b/Convert_Matlab/inst_eddies_v6.py @@ -133,20 +133,28 @@ for f in factories: while True: - os.symlink(adt_file, "adt.mat") - pty.spawn(["matlab", "-nojvm", "-r", "inst_eddies"]) - matlab_data = sio.loadmat("adt_v6.mat", squeeze_me = True, - variable_names = ("date_num", "X", "Y", - "Anticyclonic_Cell", - "Cyclonic_Cell", "ADT")) - # matlab_data is a dictionary. - - assert datetime.date.fromordinal(matlab_data["date_num"] - 366) == my_date - - for f in factories: - write(f["writers"], matlab_data[f["cell_name"]], f["cyclone"]) - - os.remove("adt.mat") + if os.access(adt_file, os.R_OK): + os.symlink(adt_file, "adt.mat") + pty.spawn(["matlab", "-nojvm", "-r", "inst_eddies"]) + matlab_data = sio.loadmat("adt_v6.mat", squeeze_me = True, + variable_names = ("date_num", "X", "Y", + "Anticyclonic_Cell", + "Cyclonic_Cell", "ADT")) + # matlab_data is a dictionary. + + assert datetime.date.fromordinal(matlab_data["date_num"] - 366) \ + == my_date + + for f in factories: + write(f["writers"], matlab_data[f["cell_name"]], f["cyclone"]) + + os.remove("adt.mat") + else: + print("Missing file:", adt_file) + + for f in factories: + write(f["writers"], [], f["cyclone"]) + my_date += datetime.timedelta(1) if my_date > final_date: break date_index += 1