from astropy.table import Table
import matplotlib.pyplot as plt
import numpy as np
import matplotlib.cm as cm

#fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_90pc_props.fits.bz2'
#fp = '/data/kant/0/sun.1608/PHANGS/ALMA/v3p4-CPROPS/STv1p5/90pc_homogenized/ngc3621_12m+7m+tp_co21_90pc_props.fits.bz2'
# 120pc_homogenized  150pc_matched     90pc_homogenized  README
# 120pc_matched      60pc_homogenized  90pc_matched
# 150pc_homogenized  60pc_matched      native

res = [60,90,120,150]
for i in res:
    fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_'+str(res[i])+'pc_props.fits.bz2'



tab = Table.read(fp)
x = np.array(tab['XCTR_DEG'])
y = np.array(tab['YCTR_DEG'])
ys = y*np.cos(x)
vel = tab['VCTR_KMS']
plt.scatter(x, y, c=vel, cmap='RdBu')
plt.xlabel('R.A.')
plt.ylabel('DEC')
cbar = plt.colorbar(cmap='RdBu')
cbar.set_label('Velocity km/s')

#plt.colorbar()
plt.show()

92711c037dffc201fdd35cea1d27807b30dbda99
<mask token>

<mask token>
for i in res:
    fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_' + str(res[i]
        ) + 'pc_props.fits.bz2'
<mask token>
plt.scatter(x, y, c=vel, cmap='RdBu')
plt.xlabel('R.A.')
plt.ylabel('DEC')
<mask token>
cbar.set_label('Velocity km/s')
plt.show()

<mask token>
res = [60, 90, 120, 150]
for i in res:
    fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_' + str(res[i]
        ) + 'pc_props.fits.bz2'
tab = Table.read(fp)
x = np.array(tab['XCTR_DEG'])
y = np.array(tab['YCTR_DEG'])
ys = y * np.cos(x)
vel = tab['VCTR_KMS']
plt.scatter(x, y, c=vel, cmap='RdBu')
plt.xlabel('R.A.')
plt.ylabel('DEC')
cbar = plt.colorbar(cmap='RdBu')
cbar.set_label('Velocity km/s')
plt.show()

from astropy.table import Table
import matplotlib.pyplot as plt
import numpy as np
import matplotlib.cm as cm
res = [60, 90, 120, 150]
for i in res:
    fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_' + str(res[i]
        ) + 'pc_props.fits.bz2'
tab = Table.read(fp)
x = np.array(tab['XCTR_DEG'])
y = np.array(tab['YCTR_DEG'])
ys = y * np.cos(x)
vel = tab['VCTR_KMS']
plt.scatter(x, y, c=vel, cmap='RdBu')
plt.xlabel('R.A.')
plt.ylabel('DEC')
cbar = plt.colorbar(cmap='RdBu')
cbar.set_label('Velocity km/s')
plt.show()

from astropy.table import Table
import matplotlib.pyplot as plt
import numpy as np
import matplotlib.cm as cm

#fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_90pc_props.fits.bz2'
#fp = '/data/kant/0/sun.1608/PHANGS/ALMA/v3p4-CPROPS/STv1p5/90pc_homogenized/ngc3621_12m+7m+tp_co21_90pc_props.fits.bz2'
# 120pc_homogenized  150pc_matched     90pc_homogenized  README
# 120pc_matched      60pc_homogenized  90pc_matched
# 150pc_homogenized  60pc_matched      native

res = [60,90,120,150]
for i in res:
    fp = '/Users/josh/projects/intro/ngc3621_12m+7m+tp_co21_'+str(res[i])+'pc_props.fits.bz2'



tab = Table.read(fp)
x = np.array(tab['XCTR_DEG'])
y = np.array(tab['YCTR_DEG'])
ys = y*np.cos(x)
vel = tab['VCTR_KMS']
plt.scatter(x, y, c=vel, cmap='RdBu')
plt.xlabel('R.A.')
plt.ylabel('DEC')
cbar = plt.colorbar(cmap='RdBu')
cbar.set_label('Velocity km/s')

#plt.colorbar()
plt.show()
