cube = fits.open(cubename+"_sub.fits")
# kernel = ((2.5*60)**2 - 50.**2)**0.5 / sqrt(8*log(2)) = 60 arcsec
# 60 arcsec / 15 "/pixel = 4
cubesm2 = FITS_tools.cube_regrid.gsmooth_cube(cube[0].data, [5,4,4], use_fft=True, psf_pad=False, fft_pad=False)
cubesm = smooth_cube(cube[0].data, kernelwidth=4, interpolate_nan=True)
cube[0].data = cubesm
cube.writeto(cubename+"_sub_smoothtoCband.fits",clobber=True)
cube[0].data = cubesm2
cube.writeto(cubename+"_sub_smoothtoCband_vsmooth.fits",clobber=True)
#makecube.make_taucube(cubename,cubename+"_continuum.fits",etamb=0.886)
#makecube.make_taucube(cubename,cubename+"_continuum.fits",etamb=0.886, suffix="_sub_smoothtoCband.fits")
# -0.4 is the most negative point in the continuum map...
makecube.make_taucube(cubename,
cubename+"_continuum.fits",
etamb=0.886,
suffix="_sub_smoothtoCband_vsmooth.fits",
outsuffix="_smoothtoCband_vsmooth.fits",
TCMB=2.7315+0.4)
makecube.make_taucube(cubename,
cubename+"_continuum.fits",
etamb=0.886,
suffix="_sub_smoothtoCband.fits",
outsuffix="_smoothtoCband.fits",
TCMB=2.7315+0.4)
makecube.make_taucube(cubename,
cubename+"_continuum.fits",
etamb=0.886,
suffix="_sub.fits",
outsuffix=".fits",
velocityrange=[-400,400],excludefitrange=[-225,250],
diagnostic_plot_name=fn.replace('.fits','_data_scrubbed.png'),
smoothto=2)
#makecube.add_file_to_cube(fn,
# cubename_discrete+'.fits',
# nhits=cubename_discrete+'_nhits.fits',wcstype='V',
# add_with_kernel=False,
# velocityrange=[-400,400],
# excludefitrange=[-225,250],
# smoothto=2)
import os
os.system('./LimaBean_H2CO11_cube_starlink.sh')
makecube.make_flats(cubename,vrange=[-20,60],noisevrange=[250,300])
makecube.make_taucube(cubename,cubename+"_continuum.fits",etamb=0.98)
# SAMPLERS:
# [('A10', 5500000000.0),
# ('A13', 4829659400.0),
# ('A14', 5500000000.0),
# ('A9', 4829659400.0),
# ('B17', 5200000000.0),
# ('B18', 5008992000.0),
# ('B21', 5200000000.0),
# ('B22', 5008992000.0),
# ('C25', 4593089000.0),
# ('C26', 4630000000.0),
# ('C29', 4593089000.0),
# ('C30', 4630000000.0),
