Updated alpha.py

punch-mission · Jun 8, 2023 · aa0d313 · aa0d313
1 parent 6e407c7
commit aa0d313
Showing 1 changed file with 5 additions and 7 deletions.
diff --git a/solpolpy/alpha.py b/solpolpy/alpha.py
@@ -17,7 +17,7 @@ def radial90(shape):
     return np.fliplr(np.rot90(np.arctan2(yy, xx), k=1))*u.radian
 
 
-def radial_north(shape):
+def radial_west(shape):
     '''
     assumes solar north is up
     creates radial polarization map
@@ -26,11 +26,9 @@ def radial_north(shape):
     x = np.arange(-x_size // 2, x_size // 2)
     y = np.arange(-y_size // 2, y_size // 2)
     xx, yy = np.meshgrid(x, y)
-    # return np.fliplr(np.arctan2(yy, xx))*u.radian
-    return np.rot90(np.fliplr(np.arctan2(yy, xx) + np.pi), k=1) * u.radian
-
+    return np.flipud(np.fliplr(np.arctan2(yy, xx) + np.pi))*u.radian
 
-def radial_west(shape):
+def radial_north(shape):
     '''
     assumes solar north is up
     creates radial polarization map
@@ -39,8 +37,8 @@ def radial_west(shape):
     x = np.arange(-x_size // 2, x_size // 2)
     y = np.arange(-y_size // 2, y_size // 2)
     xx, yy = np.meshgrid(x, y)
-    return np.fliplr(np.arctan2(yy, xx) + np.pi)*u.radian
-    # return np.rot90(np.fliplr(np.arctan2(yy, xx) + np.pi), k=1) * u.radian
+    # return np.fliplr(np.arctan2(yy, xx))*u.radian
+    return np.flipud(np.rot90(np.fliplr(np.arctan2(yy, xx) + np.pi), k=1))* u.radian
 #
 
 def zeros(shape):