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:
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"
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_co21 = glob.glob(dir_data + "data/" + galname + "*image*")[0]
### step 2/10: Kelvin to Jansky conversion
print("### step 2/10: Kelvin to Jansky conversion")
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, restfreq)
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")
#pbimage = glob.glob(dir_data+"data/"+galname+"*pb.image")[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
fitsimages = glob.glob(dir_data + "data/" + galname + "*.fits")
for i in range(len(fitsimages)):
myim2.eazy_importfits(fitsimages[i])
# find imported images
image_co10 = glob.glob(dir_data + "data/" + galname + "*co10*image*")[0]
image_co21 = glob.glob(dir_data + "data/" + galname + "*co21*image*")[0]
### step 2/10: Kelvin to Jansky conversion
print("### step 2/10: Kelvin to Jansky conversion")
bunit = imhead(image_co10, "list")["bunit"]
synsbeam10 = imhead(image_co10, "list")["beammajor"]["value"]
if bunit == "K": # CO(1-0) conversion if bunit = K
myim2.easy_K2Jy(image_co10, synsbeam10, 115.27120)
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")
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")
### 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+"images/"+galname+"*32*smooth "\
+dir_data+galname+"/"+galname+"_co32_"+suffix+".cube")
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_co10 = glob.glob(dir_data + "data/" + galname + "*co10*image*")[0]
### step 2/10: Kelvin to Jansky conversion
print("### step 2/10: Kelvin to Jansky conversion")
bunit = imhead(image_co10, "list")["bunit"]
synsbeam10 = imhead(image_co10, "list")["beammajor"]["value"]
if bunit == "K": # CO(2-1) conversion if bunit = K
myim2.easy_K2Jy(image_co10, synsbeam10, 1900.53690000 / (1 + 4.69))
image_co10 = image_co10 + ".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_co10,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
