From cd20f2e91716e3b632de592796e08b72944ab762 Mon Sep 17 00:00:00 2001
From: Jan Polcher <jan.polcher@lmd.jussieu.fr>
Date: Sat, 25 Jan 2020 23:47:15 +0100
Subject: [PATCH] A bug in the delta longitude and latitude computed to get the
boundary of the Hydro grid.
---
HydroGrid.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/HydroGrid.py b/HydroGrid.py
index 3327490..8a108c2 100644
--- a/HydroGrid.py
+++ b/HydroGrid.py
@@ -33,8 +33,8 @@ def corners(lon, lat) :
cornerspoly = []
cornersll = []
index = []
- hdlon=(np.max(lon[0,:])-np.min(lon[0,:]))/iim/2.0
- hdlat=(np.max(lat[:,0])-np.min(lat[:,0]))/jjm/2.0
+ hdlon=np.mean(np.abs(np.diff(lon[0,:])))/2.0
+ hdlat=np.mean(np.abs(np.diff(lat[:,0])))/2.0
#
for i in range(iim) :
for j in range(jjm) :
--
GitLab