We develop a method that allows numerical ocean models to include drag from icebergs, even for icebergs smaller than the model grid scale. This builds upon previous models that have either neglected iceberg drag or required higher resolution to model individual icebergs. We test our model against higher-resolution models, as well as models without iceberg drag, and show that including drag from icebergs is important for capturing realistic ocean circulation, temperature, and ice melt rates.