From 5773a7a4720505fc30469a8d496c3a02e2828d50 Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ipsl.fr>
Date: Thu, 4 Jul 2024 18:13:04 +0200
Subject: [PATCH] Simplifiy computation of `corner_proj`
---
Inst_eddies/Tests/test_set_max_speed.f90 | 9 +--------
1 file changed, 1 insertion(+), 8 deletions(-)
diff --git a/Inst_eddies/Tests/test_set_max_speed.f90 b/Inst_eddies/Tests/test_set_max_speed.f90
index 52912cc8..7873ab0c 100644
--- a/Inst_eddies/Tests/test_set_max_speed.f90
+++ b/Inst_eddies/Tests/test_set_max_speed.f90
@@ -36,13 +36,9 @@ program test_set_max_speed
character(len = :), allocatable:: path
! Window around the extremum:
-
integer llc(2) ! indices in global grid of lower left corner
integer urc(2) ! indices in global grid of upper right corner
- real corner_window(2) ! longitude and latitude of the window around
- ! the extremum, in rad
-
integer, parameter:: n_max_cont = 31
type(polyline) cont_list(n_max_cont)
@@ -132,14 +128,11 @@ program test_set_max_speed
urc(2) = min(urc(2), nlat)
if (.not. periodic) urc(1) = min(urc(1), nlon)
- corner_window = corner + (llc - 1) * step
-
call create_cont_list
call set_max_speed(e%speed_cont, e%max_speed, cont_list, cont_list_proj, &
n_cont, e%extr%coord, real(e%extr%coord_proj), &
ssh(llc(1):urc(1), llc(2):urc(2)), u(llc(1):urc(1), llc(2):urc(2)), &
- v(llc(1):urc(1), llc(2):urc(2)), &
- (corner_window - corner) / step + 1., step)
+ v(llc(1):urc(1), llc(2):urc(2)), real(llc), step)
call close_cont_list
call shpc_create(hshpc, shpc_dir = "SHPC", cyclone = cyclone, slice = 0, &
grid_lon_lat = .true.)
--
GitLab