diff --git a/Projections.py b/Projections.py index 4cb5feac326501f911a4485d96251b3e95fbead7..eb0287cfbbcc7782f4a9f5e13809c9fa931c87f9 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