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

Rename `R_Vmax` to radius

parent aab21362
No related branches found
No related tags found
No related merge requests found
......@@ -50,14 +50,14 @@ def calculate_radii_rossby(list_eddies, e_overestim, handlers, array_d_init):
lat_in_deg = shapeRec.shape.points[0][1] # in degrees
f = 2 * Omega * math.sin(math.radians(lat_in_deg)) # in s-1
V_max = shapeRec.record[4] # in m/s
R_Vmax = handlers[i_SHPC]["readers"]["max_speed_contour"]\
radius = handlers[i_SHPC]["readers"]["max_speed_contour"]\
.record(ishape)['r_eq_area'] * 1000 # in m
if abs(V_max) < 100:
avg_Rossby += V_max / (f * R_Vmax)
avg_Rossby += V_max / (f * radius)
n_valid_Rossby += 1
avg_rad += R_Vmax # in m
avg_rad += radius # in m
avg_rad /= len(list_eddies)
......
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