# retrouver les catégories 
# première colonne de station = altitude
# intervalle = largeur d'un intervalle
catAlt = np.floor((stations[:,1]-stations[:,1].min()) / intervalle)

# le début du code suivant peut-être
pDispoSAlt = np.zeros(np.unique(catAlt).shape[0])
for a in np.unique(catAlt): # fonction pratique pour isoler les valeurs uniques
    [...]