add awareness of polarization basis and some initial capabilities

Smithsonian · Apr 23, 2024 · a0ead2f · a0ead2f
1 parent 0de11e2
commit a0ead2f
Show file tree

Hide file tree

Showing 3 changed files with 190 additions and 124 deletions.
diff --git a/ngehtsim/const_def.py b/ngehtsim/const_def.py
@@ -76,6 +76,13 @@
 # default feed angle, if not otherwise known
 feed_angle = 0.0
 
+# default polarization basis, if not otherwise known
+pol_basis = 'circular'
+
+# polarization basis conversion matrices
+circ_to_lin = np.array([[1.0,1.0],[-1.0j,1.0j]])/np.sqrt(2.0)
+lin_to_circ = np.array([[1.0,1.0j],[1.0,-1.0j]])/np.sqrt(2.0)
+
 ###################################################
 # physical constants
 
@@ -354,12 +361,12 @@
 ###################################################
 # pull antenna properties from table
 
-known_stations, tlcs, diam_arr, surf_arr, mnts_arr, fa_arr, altnames = np.loadtxt(path_to_tsm,
-                                                                                  delimiter=',',
-                                                                                  skiprows=1,
-                                                                                  usecols=(0, 1, 7, 8, 9, 10, 11),
-                                                                                  dtype=str,
-                                                                                  unpack=True)
+known_stations, tlcs, diam_arr, surf_arr, polbasis_arr, mnts_arr, fa_arr, altnames = np.loadtxt(path_to_tsm,
+                                                                                     delimiter=',',
+                                                                                     skiprows=1,
+                                                                                     usecols=(0, 1, 7, 8, 9, 10, 11, 12),
+                                                                                     dtype=str,
+                                                                                     unpack=True)
 lat_arr, lon_arr, elev_arr = np.loadtxt(path_to_tsm,
                                         delimiter=',',
                                         skiprows=1,
@@ -391,6 +398,11 @@
     if (elev_arr[i] != ''):
         known_elevations[known_stations[i]] = elev_arr[i]
 
+known_polbases = {}
+for i in range(len(known_stations)):
+    if (polbasis_arr[i] != ''):
+        known_polbases[known_stations[i]] = polbasis_arr[i]
+
 known_mount_types = {}
 for i in range(len(known_stations)):
     if (mnts_arr[i] != ''):