#!/common/casa/casa-release-5.4.1-31.el6/bin/casa
from suncasa.eovsa import eovsa_diskmodel as ed
from glob import glob
eofiles = glob('eovsa_????????.spw??-??.tb.fits')
eofiles = sorted(eofiles)
eofiles_new = []
diskxmlfile = '/data1/eovsa/fits/UDBms_slfcaled/201910/UDB20191001.ms.SOLDISK.xml'
for idx, eofile in enumerate(eofiles):
    eomap_disk, tb_disk, eofile_new = ed.image_adddisk(eofile, diskxmlfile)
    eofiles_new.append(eofile_new)
ed.plt_eovsa_image(eofiles_new)