Skip to content
Snippets Groups Projects
Commit 83f21057 authored by Lionel GUEZ's avatar Lionel GUEZ
Browse files

Do not re-create `adt_file` for first iteration

We already have it: it is `sys.argv[1]`.
parent a5074514
No related branches found
No related tags found
No related merge requests found
......@@ -89,7 +89,8 @@ def write(writers, cell, cyclone):
n_shapes = len(writers["extr"])
writers["ishape_last"].write(str(n_shapes - 1) + "\n")
dirname, basename = path.split(sys.argv[1])
adt_file = sys.argv[1]
dirname, basename = path.split(adt_file)
my_date = datetime.datetime.strptime(basename, "adt_%Y-%m-%d_v6.mat").date()
if len(sys.argv) == 3:
......@@ -121,9 +122,7 @@ for f in factories:
file = path.join(f["SHP_coll"], "ishape_last.txt")
f["writers"]["ishape_last"] = open(file, "w")
while my_date <= final_date:
adt_file = my_date.strftime("adt_%Y-%m-%d_v6.mat")
adt_file = path.join(dirname, adt_file)
while True:
matlab_data = sio.loadmat(adt_file, squeeze_me = True,
variable_names = ("date_num", "X", "Y",
"Anticyclonic_Cell",
......@@ -138,6 +137,9 @@ while my_date <= final_date:
write(f["writers"], matlab_data[f["cell_name"]], f["cyclone"])
my_date += datetime.timedelta(1)
if my_date > final_date: break
adt_file = my_date.strftime("adt_%Y-%m-%d_v6.mat")
adt_file = path.join(dirname, adt_file)
nml = f90nml.Namelist()
nml["grid_nml"] = dict(corner_deg = [matlab_data["X"][0], matlab_data["Y"][0]],
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment