def CompareUV(uvFile1, uvFile2, err, verbose=True): """ Compare UV FITS files *uvFile1* and *uvFile2*. * uvFile1, uvFile2 = UV FITS files to be compared * err = OErr object """ uv1 = UV.newPFUV('uv1', uvFile1, 0, True, err) uv2 = UV.newPFUV('uv2', uvFile2, 0, True, err) rms = UV.PUtilVisCompare(uv1, uv2, err) if verbose: print("UV data set 1: " + uvFile1) print("UV data set 2: " + uvFile2) print("RMS of real, imag differences / amplitude of data set 2 =\n%6i" % (rms))
Input['Type'] = 0 # CC #replace data with model? #Input['REPLACE'] = True # Subtract #print "Subtract" #print "Replace" #SkyModel.PSubUV(err, Input) #OErr.printErrMsg(err, "Error subtracting") # Divide print "Divide" SkyModel.PDivUV(err, Input) OErr.printErrMsg(err, "Error dividing") # Compare with master lie [rms diff] masterData = UV.newPFUV("Master UVData", masterFile, masterDisk, True, err) diff = UV.PUtilVisCompare(outData, masterData, err) print "Comparison with master lie, fractional RMS R,I difference", diff # Say something #print "Subtracted",inModel,"From", inFile,"to",outFile print "Divided FT of", inModel, "into", inFile, "to", outFile #def PUtilUVUtilVisCompare (in1UV, in2UV, err): # returns RMS # Shutdown Obit OErr.printErr(err) OSystem.Shutdown(ObitSys)