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

Polish

parent efdfee04
No related branches found
No related tags found
No related merge requests found
......@@ -232,12 +232,9 @@ print("Iterating on edges...")
for edge in g.edges():
source_node = edge.source()
target_node = edge.target()
latitude = (pos_last[source_node][1] + pos_first[target_node][1]) / 2
# Because of the wrapping issue (360° wrapping incorrectly to 0°),
# we check for that here:
lon_diff = abs(pos_last[source_node][0] - pos_first[target_node][0])
if lon_diff > math.radians(300):
lon_diff = 2 * math.pi - lon_diff
......@@ -246,7 +243,7 @@ for edge in g.edges():
+ (pos_last[source_node][1] - pos_first[target_node][1]) ** 2
)
# Rossbies:
# Rossby numbers:
if first_av_ros[target_node] and last_av_ros[source_node]:
Delta_Ro = last_av_ros[source_node] - first_av_ros[target_node]
else:
......
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