From 9f76e0025c6f182ff585bb25bfaf76505946b26e Mon Sep 17 00:00:00 2001
From: POLCHER Jan <jan.polcher@lmd.jussieu.fr>
Date: Sun, 2 Feb 2020 01:09:06 +0100
Subject: [PATCH] Small corrections on the llij function for the Lambert
projection.
---
Projections.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/Projections.py b/Projections.py
index 4cb5fea..eb0287c 100644
--- a/Projections.py
+++ b/Projections.py
@@ -69,14 +69,14 @@ def LC_llij(polyll, hemi, truelat1, truelat2, polei, polej, rebydx, stdlon, cone
for lon,lat in polyll :
- deltalon = lon - self.stdlon
+ deltalon = lon - stdlon
if (deltalon > +180.) : deltalon = deltalon - 360.
if (deltalon < -180.) : deltalon = deltalon + 360.
tl1r = truelat1 * rad_per_deg
ctl1r = np.cos(tl1r)
- rm = self.rebydx * ctl1r/cone * (np.tan((90.*hemi-lat)*rad_per_deg/2.) / np.tan((90.*hemi-truelat1)*rad_per_deg/2.))**cone
+ rm = rebydx * ctl1r/cone * (np.tan((90.*hemi-lat)*rad_per_deg/2.) / np.tan((90.*hemi-truelat1)*rad_per_deg/2.))**cone
arg = cone*(deltalon*rad_per_deg)
i = polei + hemi * rm * np.sin(arg)
@@ -127,7 +127,7 @@ class LambertC:
def llij(self, polyll) :
- polyij = LC_llij(polyij, self.hemi, self.truelat1, self.truelat2, self.polei, self.polej, self.rebydx, self.stdlon, self.cone)
+ polyij = LC_llij(polyll, self.hemi, self.truelat1, self.truelat2, self.polei, self.polej, self.rebydx, self.stdlon, self.cone)
return polyij
--
GitLab