Проблема в том, что ваши значения в z больше 25. Если вы хотите ограничить значения шкалы 25 или меньше, вам нужно либо обрезать значения в z, которые больше 25, либо удалить их, прежде чем сделатьplot.
Вот пример усечения значений от z до 25 или меньше:
import numpy
import numpy.random
import matplotlib.pyplot as plt
znew = z
znew[z > 25] = 25
t = numpy.linspace(0, 1, 20)
x, y = numpy.meshgrid(t, t)
# this is supposed to be the maximum value of the plot
max_value = 25 # or 100
fig = plt.figure()
axis = fig.add_subplot(1, 1, 1)
plot = axis.contourf(x, y, znew, 100, cmap=None, vmin=0, vmax=max_value)
cbar = plt.colorbar(plot)
plt.show(block=False)

В качестве альтернативы выможно просто удалить эти значения, установив их на NaN
:
znew = z
znew[z > 25] = numpy.nan
t = numpy.linspace(0, 1, 20)
x, y = numpy.meshgrid(t, t)
# this is supposed to be the maximum value of the plot
max_value = 25 # or 100
fig = plt.figure()
axis = fig.add_subplot(1, 1, 1)
plot = axis.contourf(x, y, znew, 100, cmap=None, vmin=0, vmax=max_value)
cbar = plt.colorbar(plot)
plt.show(block=False)
25 removed">