Add subtract of rotationmeasure from rotation

revoltek · Nov 6, 2024 · b0e1a17 · b0e1a17
1 parent 0a67206
commit b0e1a17
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/losoto/operations/residuals.py b/losoto/operations/residuals.py
@@ -3,7 +3,9 @@
 
 # Residual operation for LoSoTo
 
-# This operation subtracts/divides two tables or a clock/tec/tec3rd/rm from a phase
+# This operation subtracts/divides two tables, or subtracts a clock/tec/tec3rd/rm
+# from a phase, or subtracts a rotationmeasure from a rotation.
+
 # Operation is flag-only capable
 
 from losoto.lib_operations import *
@@ -21,7 +23,9 @@ def _run_parser(soltab, parser, step):
 
 def run( soltab, soltabsToSub, ratio=False ):
     """
-    Subtract/divide two tables or a clock/tec/tec3rd/rm from a phase.
+    Subtract/divide two tables of the same type, or subtract a clock/tec/tec3rd/rm
+    table from a phase table, or subtract a rotationmeasure table from a rotation
+    table.
 
     Parameters
     ----------
@@ -123,6 +127,10 @@ def run( soltab, soltabsToSub, ratio=False ):
 
                 vals = np.swapaxes(vals, 0, idxPol)
 
+            elif soltabsub.getType() == 'rotationmeasure' and soltab.getType() == 'rotation':
+                wav = 2.99792458e8 / freq
+                vals -= wav * wav * valsSub
+
             else:
                 if ratio:
                     vals = (vals - valsSub) / valsSub