Similarly, mths = np.arange(12) pts = np.arange(numpts) out_array[mths, pts] = array[mths, 0, r, c] Does not work either... -- View this message in context: http://old.nabble.com/removing-for-loops...-tp28633477p28819859.html Sent from the Scipy-User mailing list archive at Nabble.com.