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

Remove interpolat field from shapefiles

parent a24e8dfd
No related branches found
No related tags found
No related merge requests found
......@@ -26,7 +26,6 @@ def define_fields(writers):
writers["extr"].field("ssh", "N", 13, 6)
writers["extr"].field("days_1950", "N", 5)
writers["extr"].field("eddy_index", "N", 5)
writers["extr"].field("interpolat", "N", 1)
writers["extr"].field("valid", "N", 1)
writers["extr"].field("speed", "N", 13, 6)
......@@ -67,9 +66,9 @@ def write(writers, cell, cyclone):
i = np.argwhere(matlab_data["X"] == eddy[0]).item()
j = np.argwhere(matlab_data["Y"] == eddy[1]).item()
speed = 1e4 if np.isnan(eddy[12]) else eddy[12]
writers["extr"].record(ssh = matlab_data["ADT"][i, j], interpolat = 0,
valid = 1, days_1950 = d,
eddy_index = eddy_index, speed = speed)
writers["extr"].record(ssh = matlab_data["ADT"][i, j], valid = 1,
days_1950 = d, eddy_index = eddy_index,
speed = speed)
if cyclone:
ssh = matlab_data["ADT"][i, j] + eddy[7]
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment