From cc1fe809237cd880b191da1748b87a04f0cb58b1 Mon Sep 17 00:00:00 2001
From: POLCHER Jan <jan.polcher@lmd.jussieu.fr>
Date: Wed, 15 Jul 2020 10:36:08 +0200
Subject: [PATCH] Correction of a small error on the resolution calculation.
---
HydroGrid.py | 10 +++++++---
1 file changed, 7 insertions(+), 3 deletions(-)
diff --git a/HydroGrid.py b/HydroGrid.py
index 07db97b..2f1c2f4 100644
--- a/HydroGrid.py
+++ b/HydroGrid.py
@@ -187,10 +187,14 @@ class HydroGrid :
return indices, polylist
def resolution(self,i,j) :
- if i+1 >= self.iim or j+1 >= self.jjm :
- inx=i-1; jnx=j-1
+ if i+1 >= iim :
+ inx=i-1
else :
- inx=i+1; jnx=j+1
+ inx=i+1
+ if j+1 >= jjm :
+ jnx=j-1
+ else :
+ jnx=j+1
dlon = np.radians(np.abs(self.lon[j,i]-self.lon[j,inx]))
alon = np.cos(np.radians(self.lat[j,i]))**2*(np.sin(dlon/2))**2
clon = 2*np.arctan2(np.sqrt(alon), np.sqrt(1-alon))
--
GitLab