Problem with contour plots in PythonPlot.jl

With reference to my earlier post here, I have stumbled upon a bad issue.

As in the earlier post, my code is as follows:

using PythonPlot, LaTeXStrings, Measures, PlotThemes

f(X, Y) = (100/(3*s*k^2))*(((1-(X/Y)^2)^3))*(Y^4)

mₕ, v, mₜ, Nₜ = 125, 246, 172.76, 3
k= pi*mₕ*v
s = 1 - (Nₜ/(3*k^2))mₜ^4
X, Y = range(0, 200, 1000), range(0, 500, 2000)
Z = @. f(X', Y)
zc = [4.2, 5, 10, 30, 70, 100]

contour(X, Y, Z, levels=zc, lw=1, color = :turbo, clims=extrema(zc), cbar = false, clabels = true, framestyle = :box, fontfamily = "cmr10",
        ylims = (150, 480), xlabel = "M [GeV]", ylabel = L"$m_{\Phi}$ [GeV]", title = L"\Delta \lambda/ \lambda^{SM}_{hhh}\ (\%)")

The output I was getting earlier is:
Screenshot 2024-08-21 122157

However, I was running the code again after adding an additional point for z=4.2. The code is the same - the output this time is:

I don’t understand this at all - the code seems to completely ignore the labels, titles, axes, everything.

Moreover, if I import PlotThemes.jl and run theme(:bright), the statement shoots up an error:

UndefVarError: `theme` not defined

Can anyone help me with this? It really doesn’t matter to me if the solution is not based in PythonPlot at all - I just want to create the figure.

EDIT:
Sometimes I get the warning:

sys:1: UserWarning: The following kwargs were not used by contour: 'lw', 'color', 'clims', 'cbar', 'clabels', 'framestyle', 'fontfamily', 'ylims', 'xlabel', 'ylabel', 'title'

You seem to be mixing Plots.jl with PythonPlot.jl code.
Using Plots.jl’s pythonplot() backend, your example runs just fine.

Plots.jl pythonplot() code
using LaTeXStrings, Measures, Plots
pythonplot(dpi=100)

f(X, Y) = (100/(3*s*k^2))*(((1-(X/Y)^2)^3))*(Y^4)

mₕ, v, mₜ, Nₜ = 125, 246, 172.76, 3
k= pi*mₕ*v
s = 1 - (Nₜ/(3*k^2))mₜ^4
X, Y = range(0, 200, 1000), range(0, 500, 2000)
Z = @. f(X', Y)
zc = [4.2, 5, 10, 30, 70, 100]

theme(:bright)
contour(X, Y, Z, levels=zc, lw=1, color=:turbo, clims=extrema(zc), cbar=false, clabels=true, framestyle=:box, fontfamily="cmr10",
        ylims=(150, 480), xlabel="M [GeV]", ylabel=L"$m_{\Phi}$ [GeV]", title=L"\Delta \lambda/ \lambda^{SM}_{hhh}\ (\%)")

Yes, well, that is the problem - for some reason, parts of this code don’t seem to run. The theme part, for example.

I stopped importing PythonPlot.jl and ran the code entirely with Plots.jl with the default backend. Turns out that works almost perfectly except for the fact that the contour labels (5, 10, 30, 70, 100) are written twice on their corresponding contour lines. In the end, I switched off the label and printed the contour plot and added the labels by hand.