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

Rename `coord_extr` to `extr_coord`

Rename `coord_extr` to `extr_coord` and `coord_extr_proj` to
`extr_coord_proj`, to have the same names as in procedure
`set_max_speed`.
parent a52041a8
No related branches found
No related tags found
No related merge requests found
......@@ -5,7 +5,7 @@ module get_1_outerm_m
contains
subroutine get_1_outerm(out_cont, cont_list, cont_list_proj, n_cont, &
cyclone, coord_extr, innermost_level_2, outside_points, ssh, corner, &
cyclone, extr_coord, innermost_level_2, outside_points, ssh, corner, &
step)
! This procedure gets one outermost good contour with sufficient
......@@ -80,7 +80,7 @@ contains
! number of good contours found and stored, 0 <= n_cont <= n_max_cont - 1
logical, intent(in):: cyclone
real, intent(in):: coord_extr(:) ! (2) longitude and latitude, in rad
real, intent(in):: extr_coord(:) ! (2) longitude and latitude, in rad
real, intent(in):: innermost_level_2
! ssh level of the innermost contour to consider around the target
......@@ -106,18 +106,18 @@ contains
logical mask(size(ssh, 1), size(ssh, 2))
integer n_max_cont ! >= 3
integer i
real corner_proj(2), coord_extr_proj(2)
real corner_proj(2), extr_coord_proj(2)
real outside_points_proj(2, size(outside_points, 2))
!-----------------------------------------------------------------
corner_proj = (corner - corner_whole) / step + 1.
coord_extr_proj = (coord_extr - corner_whole) / step + 1.
extr_coord_proj = (extr_coord - corner_whole) / step + 1.
outside_points_proj = convert_to_ind(outside_points, corner_whole, step)
n_max_cont = size(cont_list_proj)
n_cont = 0
cont_list_proj(1)%polyline = good_contour(corner_proj, ssh, &
innermost_level_2, coord_extr_proj, outside_points_proj)
innermost_level_2, extr_coord_proj, outside_points_proj)
test_n_points: if (.not. cont_list_proj(1)%closed) then
! cont_list_proj(1)%n_points == 0
......@@ -131,7 +131,7 @@ contains
call assert(merge(level_try > level_good, level_try < level_good, &
cyclone), "get_1_outerm level_try")
cont_list_proj(2)%polyline = good_contour(corner_proj, ssh, level_try, &
coord_extr_proj, outside_points_proj)
extr_coord_proj, outside_points_proj)
if (cont_list_proj(2)%closed) then
! cont_list_proj(2)%n_points /= 0
......@@ -147,7 +147,7 @@ contains
do while (abs(level_bad - level_good) > precision)
level_try = (level_good + level_bad) / 2.
cont_list_proj(n_cont + 1)%polyline = good_contour(corner_proj, &
ssh, level_try, coord_extr_proj, outside_points_proj)
ssh, level_try, extr_coord_proj, outside_points_proj)
if (cont_list_proj(n_cont + 1)%closed) then
! cont_list_proj(n_cont + 1)%n_points /= 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