dir_casa_region = dir_data1 + "casa_region/"
done = glob.glob(dir_casa_region)
if not done:
os.mkdir(dir_casa_region)
### Nyquist sampling with varying aperture size
for i in range(len(apertures)):
print("##### sampling with aperture = " + str(apertures[i]))
# define sampling grid at the given aperture size
if galname == "ngc3627":
stp_ra, stp_dec, rng_ra, rng_dec = Nyq.def_step2(apertures[i], fov)
else:
stp_ra, stp_dec, rng_ra, rng_dec = Nyq.def_step(apertures[i], fov)
# rms calculation at the given aperture size
Sa_co, Sb_co = Nyq.def_area(m0_co21, apertures[i], beam_co)
rms_apt_co10 = rms_co10 * velres * sqrt(nchan) / sqrt(Sa_co) #average
rms_apt_co21 = rms_co21 * velres * sqrt(nchan) / sqrt(Sa_co) #average
Sa_wise, Sb_wise = Nyq.def_area(image_w3, apertures[i], beam_wise)
rms_apt_w1 = rms_w1 / sqrt(Sa_wise) #average
rms_apt_w2 = rms_w2 / sqrt(Sa_wise) #average
rms_apt_w3 = rms_w3 / sqrt(Sa_wise) #average
# setup for imval txt output
if type(apertures[i]) == int:
name_size = "{0:02d}".format(apertures[i])
elif type(apertures[i]) == float:
name_size = str(apertures[i]).replace(".", "p").replace("0", "")
product_file = dir_data1+galname+"_flux_"+suf_wise+"_"\
+str(name_size)+".txt"
if not done:
os.mkdir(dir_casa_region)
### Nyquist sampling with varying aperture size
for i in range(len(apertures)):
print("##### sampling with aperture = " + str(apertures[i]))
# define sampling grid at the given aperture size
if galname == "ngc3627":
stp_ra,stp_dec,rng_ra,rng_dec=Nyq.def_step2(apertures[i],fov)
else:
stp_ra,stp_dec,rng_ra,rng_dec=Nyq.def_step(apertures[i],fov)
# rms calculation at the given aperture size
Sa_co,Sb_co = Nyq.def_area(image_co21,apertures[i],beam_size)
rms_apt_co10 = rms_co10*velres*sqrt(nchan)/sqrt(Sa_co) #average
rms_apt_co21 = rms_co21*velres*sqrt(nchan)/sqrt(Sa_co) #average
# setup for imval txt output
if type(apertures[i]) == int:
name_size = "{0:02d}".format(apertures[i])
elif type(apertures[i]) == float:
name_size = str(apertures[0]).replace(".","p")
product_file = dir_data+galname+"_flux_"+suffix+"_"\
+str(name_size)+".txt"
os.system("rm -rf "+product_file)
f = open(product_file,"a")
f.write(header_output+"\n")
f.close()
done = glob.glob(dir_casa_region)
if not done:
os.mkdir(dir_casa_region)
# Nyquist sampling with varying aperture size
for i in range(len(apertures)):
print("##### sampling with aperture = " + str(apertures[i]))
# define sampling grid at the given aperture size
stp_ra,stp_dec,rng_ra,rng_dec=Nyq.def_step(apertures[i],fov)
#ra_dgr = ra_dgr + stp_ra
#dec_dgr = dec_dgr - stp_dec
# rms calculation at the given aperture size
Sa_co,Sb_co = Nyq.def_area(m0_co21,apertures[i],beam)
rms_apt_co10 = rms_co10*velres*sqrt(nchan)/sqrt(Sa_co) #average
rms_apt_co21 = rms_co21*velres*sqrt(nchan)/sqrt(Sa_co) #average
# setup for imval txt output
if type(apertures[i]) == int:
name_size = "{0:02d}".format(apertures[i])
elif type(apertures[i]) == float:
name_size = str(apertures[0]).replace(".","p")
product_file = dir_data1+galname+"_fluxin_"+suffix+"_"\
+str(name_size)+".txt"
os.system("rm -rf "+product_file)
f = open(product_file,"a")
f.write(header_output+"\n")
f.close()