Urban land surface processes need to be represented to inform future urban climate and building energy projections. Here, the single layer urban canopy model Town Energy Balance (TEB) is coupled to the Weather Research and Forecasting (WRF) model to create WRF-TEB. The coupling method is described generically, implemented into software, and the code and data are released with a Singularity image to address issues of scientific reproducibility. The coupling is implemented modularly and verified by an integration test. Results show no detectable errors in the coupling. Separately, a meteorological evaluation is undertaken using observations from Toulouse, France. The latter evaluation, during an urban canopy layer heat island episode, shows reasonable ability to estimate turbulent heat flux densities and other meteorological quantities. We conclude that new model couplings should make use of integration tests as meteorological evaluations by themselves are insufficient, given that errors are difficult to attribute because of the interplay between observational errors and multiple parameterization schemes (e.g., radiation, microphysics, and boundary layer).