### Define Parameters
#####################
#####################
### Main Procedure
#####################
os.system("rm -rf "+dir_data+galname+"/*"+suffix+"*")
### step 1/10: importfits
print("### step 1/10: importfits")
fitsimages = glob.glob(dir_data+"data/"+galname+"*.fits")
for i in range(len(fitsimages)):
myim2.eazy_importfits(fitsimages[i],defaultaxes=False)
# find imported images
image_co32 = glob.glob(dir_data+"data/"+galname+"*co32*image*")[0]
### step 2/10: Kelvin to Jansky conversion
print("### step 2/10: Kelvin to Jansky conversion")
bunit = imhead(image_co32,"list")["bunit"]
synsbeam21 = imhead(image_co32,"list")["beammajor"]["value"]
if bunit == "K": # CO(2-1) conversion if bunit = K
myim2.easy_K2Jy(image_co32,synsbeam21,230.53800)
image_co32 = image_co32 + ".jy"
else:
thres_masking = 8.5 # threshold s/n ratio for masking
thres_mom = 1.0 # threshold s/n ratio for immoments
chans = "4~70" # channel selection for immoments
pbcut = 0.50 # pirmary beam cut
#####################
### Main Procedure
#####################
os.system("rm -rf " + dir_data + galname + "/*" + suffix + "*")
### step 1/10: importfits
print("### step 1/10: importfits")
fitsimages = glob.glob(dir_data + "images/" + galname + "*.fits")
for i in range(len(fitsimages)):
myim2.eazy_importfits(fitsimages[i])
# find imported images
image_co32 = glob.glob(dir_data + "images/" + galname + "*.image")[0]
### step 2/10: Kelvin to Jansky conversion
print("### step 2/10: Kelvin to Jansky conversion")
bunit = imhead(image_co32, "list")["bunit"]
synsbeam10 = imhead(image_co32, "list")["beammajor"]["value"]
if bunit == "K": # CO(1-0) conversion if bunit = K
myim2.easy_K2Jy(image_co32, synsbeam10, 115.27120)
image_co32 = image_co32 + ".jy"
else:
print("# skip easy_K2Jy for the CO(3-2) data")