find_refant includes case of missing/flagged antenna.

e-merlin · Sep 25, 2017 · e262761 · e262761
1 parent 489230d
commit e262761
Showing 1 changed file with 7 additions and 4 deletions.
diff --git a/functions/eMERLIN_CASA_functions.py b/functions/eMERLIN_CASA_functions.py
@@ -621,10 +621,13 @@ def find_refant(msfile, field, antennas='', spws='', scan=''):
                 snrj = []
                 for basel in antennas:
                     if basel != anten:
-                        d2 = visstat(msfile, antenna='{}&{}'.format(anten, basel),
-                                     field=field, axis='amplitude', scan=scan,
-                                     spw=str(spw)+channels, correlation=corr)
-                        snrj.append(d2['DATA']['median']/d2['DATA']['stddev'])
+                        try:
+                            d2 = visstat(msfile, antenna='{}&{}'.format(anten, basel),
+                                         field=field, axis='amplitude', scan=scan,
+                                         spw=str(spw)+channels, correlation=corr)
+                            snrj.append(d2['DATA']['median']/d2['DATA']['stddev'])
+                        except:
+                            snrj.append(0.0)
                 snr[i, j, k] = np.mean(snrj)
 
     # Sorting antennas by average S/N, averaged in polarization and spw.