Accurate simulations of atmospheric boundary layer flow are vital for predicting dispersion of contaminant releases, particularly in densely populated urban regions where first responders must react within minutes and the consequences of forecast errors are potentially disastrous. Current mesoscale models do not account for urban effects, and conversely urban scale models do not account for mesoscale weather features or atmospheric physics. The ultimate goal of this research is to develop and implement an immersed boundary method (IBM) along with a surface roughness parameterization into the mesoscale Weather Research and Forecasting (WRF) model. IBM will be used in WRF to represent the complex boundary conditions imposed by urban landscapes, while still including forcing from regional weather patterns and atmospheric physics. This document details preliminary results of this research, including the details of three distinct implementations of the immersed boundary method. Results for the three methods are presented for the case of a rotation influenced neutral atmospheric boundary layer over flat terrain.