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)
