add script for tuning demixing

Closes RAP-90, script: tune_demixing_parameters.py

python tune_demixing_parameters.py --MS '/some/path/data/L99*MS' --sky_model /path/to/A-team-model.txt

Will examine all patches in given sky model, and output witch patches to actually use in demixing.

Further tuning with number of iterations, frequency resolution is being developed.

Requirements: LINC, lsmtool, DP3

scripts/tune_demixing_parameters.py

@@ -346,8 +346,12 @@ class DemixingTuner:
           self.AIC[index]=rss*self.N*self.N+np.sum(chosen_dirs)*self.N
 
     def report(self):
-        index=np.argmin(self.AIC[1:])+1
-        chosen_dirs=self.scalar_to_kvec(index)
+        indx=np.argsort(self.AIC[1:])+1
+        chosen_dirs=np.zeros(self.K)
+        for index in indx[:3]:
+            chosen_dirs+=self.scalar_to_kvec(index)
+        # remove -ve elevation
+        chosen_dirs[self.elevation<=5]=0
         chosen_patches=itertools.compress(self.ateam_names,chosen_dirs)
         print('***************** Result:')
         for patch in chosen_patches: