code: def stepceil_with_units(param, step, unit):
      
      counter = 0 * unit
      while counter < param.to(unit):
          counter += step * unit
      return counter
docs: This function returns the smallest multiple of 'step' greater than or
      equal to 'param' and outputs the result in Pint units.
      This function is unit-aware and functions without requiring translation
      so long as 'param' and 'unit' are of the same dimensionality.