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:
print("# skip easy_K2Jy for the CO(2-1) data")
### step 3/10: imregrid
print("### step 3/10: imregrid")
pbimage = glob.glob(dir_data+"data/"+galname+"*pb.image")[0]
myim2.easy_imregrid(pbimage,image_co32,False) # pbimage
pbimage = pbimage + ".regrid"
### step 4/10: imsmooth
print("### step 4/10: imsmooth")
beam_mask = beam_size * increment_mask # beam size for the masking
myim2.easy_imsmooth(image_co32,beam_mask,False) # co32
### mv the cubes to the working directory
os.system("mkdir "+dir_data+galname)
os.system("mv "+dir_data+"data/"+galname+"*32*smooth "\
image_co10 = image_co10 + ".jy"
else:
print("# skip easy_K2Jy for the CO(1-0) data")
bunit = imhead(image_co21, "list")["bunit"]
synsbeam21 = imhead(image_co21, "list")["beammajor"]["value"]
if bunit == "K": # CO(2-1) conversion if bunit = K
myim2.easy_K2Jy(image_co21, synsbeam21, 230.53800)
image_co21 = image_co21 + ".jy"
else:
print("# skip easy_K2Jy for the CO(2-1) data")
### step 3/10: imregrid
print("### step 3/10: imregrid")
myim2.easy_imregrid(image_co10, image_co21) # co10
image_co10 = image_co10 + ".regrid"
pbimage = glob.glob(dir_data + "data/" + galname + "*.pb")[0]
myim2.easy_imregrid(pbimage, image_co21, False) # pbimage
pbimage = pbimage + ".regrid"
### step 4/10: imsmooth
print("### step 4/10: imsmooth")
beam_mask = beam_size * increment_mask # beam size for the masking
myim2.easy_imsmooth(image_co10, beam_mask, False) # co10
myim2.easy_imsmooth(image_co21, beam_mask, False) # co21
### mv the cubes to the working directory
os.system("mkdir " + dir_data + galname)