Hi Rafael. Thanks for your response. Strangely, this solution does not work when the Z dimension only has a single unique value. Here is an example:
using Plots
pyplot()
x = .1:0.05:2
y = .1:0.05:1
X = repeat(reshape(x, 1, :), length(y), 1)
Y = repeat(y, 1, length(x))
Z = fill(.655848484848, size(X))
p1 = contour(x, y, Z, fill = true, colorbar_ticks=[.65])
It is certainly an improvement because it does not have a long number, but is there a way to still show it somehow?