def azimuth_direction(azstart, azend):
  """
  Direction for shortest distance between two azimuths

  Returns the direction and the distance
  """
  delta_az = azend-azstart
  if delta_az > 180:
    return 'ccw', 360-delta_az
  elif delta_az < -180:
    return 'cw', 360+delta_az
  elif delta_az < 0:
    return 'ccw', -delta_az
  else:
    return 'cw', delta_az