CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA.

We present new HI spectral line images of the nearby low-mass galaxy NGC 5238, acquired with the Karl G. Jansky Very Large Array (VLA). Located at a distance of 4.51+/-0.04 Mpc, NGC 5238 is an actively star-forming galaxy with widespread H-alpha and UV continuum emission. The source is included in many ongoing and recent nearby galaxy surveys, but until this work the spatially resolved qualities of its neutral interstellar medium have remained unstudied. Our HI images resolve the disk on physical scales of ~400 pc, allowing us to undertake a detailed comparative study of the gaseous and stellar components. The HI disk is asymmetric in the outer regions, and the areas of high HI mass surface density display a crescent-shaped morphology that is slightly offset from the center of the stellar populations. The HI column density exceeds 10^21 cm^-2 in much of the disk. We quantify the degree of co-spatiality of dense HI gas and sites of ongoing star formation as traced by far-UV and H-alpha emission. The neutral gas kinematics are complex; using a spatially-resolved position-velocity analysis, we infer a rotational velocity of 31+/-5 km/s. We place NGC 5238 on the baryonic Tully-Fisher relation and contextualize the system amongst other low-mass galaxies.

We analyze the relationships between atomic, neutral hydrogen (HI) and star formation (SF) in the 12 low-mass SHIELD galaxies. We compare high spectral (~0.82 km/s/channel) and spatial resolution (physical resolutions of 170 pc - 700 pc) HI imaging from the VLA with H\alpha and far-ultraviolet imaging. We quantify the degree of co-spatiality between star forming regions and regions of high HI column densities. We calculate the global star formation efficiencies (SFE, $\Sigma_{\rm SFR}$ / $\Sigma_{\rm HI}$), and examine the relationships among the SFE and HI mass, HI column density, and star formation rate (SFR). The systems are consuming their cold neutral gas on timescales of order a few Gyr. While we derive an index for the Kennicutt-Schmidt relation of N ~ 0.68 $\pm$ 0.04 for the SHIELD sample as a whole, the values of N vary considerably from system to system. By supplementing SHIELD results with those from other surveys, we find that HI mass and UV-based SFR are strongly correlated over five orders of magnitude. Identification of patterns within the SHIELD sample allows us to bin the galaxies into three general categories: 1) mainly co-spatial HI and SF regions, found in systems with highest peak HI column densities and highest total HI masses, 2) moderately correlated HI and SF regions, found in systems with moderate HI column densities, and 3) obvious offsets between HI and SF peaks, found in systems with the lowest total HI masses. SF in these galaxies is dominated by stochasticity and random fluctuations in their ISM.

We present kinematic analyses of the 12 galaxies in the "Survey of HI in Extremely Low-mass Dwarfs" (SHIELD). We use multi-configuration interferometric observations of the HI 21cm emission line from the Karl G. Jansky Very Large Array (VLA) to produce image cubes at a variety of spatial and spectral resolutions. Both two- and three-dimensional fitting techniques are employed in an attempt to derive inclination-corrected rotation curves for each galaxy. In most cases, the comparable magnitudes of velocity dispersion and projected rotation result in degeneracies that prohibit unambiguous circular velocity solutions. We thus make spatially resolved position-velocity cuts, corrected for inclination using the stellar components, to estimate the circular rotation velocities. We find circular velocities <30 km/s for the entire survey population. Baryonic masses are calculated using single-dish HI fluxes from Arecibo and stellar masses derived from HST and Spitzer imaging. Comparison is made with total dynamical masses estimated from the position-velocity analysis. The SHIELD galaxies are then placed on the baryonic Tully-Fisher relation. There exists an empirical threshold rotational velocity <15 km/s, below which current observations cannot differentiate coherent rotation from pressure support. The SHIELD galaxies are representative of an important population of galaxies whose properties cannot be described by current models of rotationally-dominated galaxy dynamics.