Currently, base_models.physical_base.deposition_velocity uses lots of matlib operations (e.g. repmat). We should remove these dependencies and make better use of numpy broadcasting. Moreover, the iterative solver for the settling velocity could likely be improved in terms of readability and computational efficiency.
Currently,
base_models.physical_base.deposition_velocityuses lots ofmatliboperations (e.g.repmat). We should remove these dependencies and make better use of numpy broadcasting. Moreover, the iterative solver for the settling velocity could likely be improved in terms of readability and computational efficiency.