I tried slightly different functions inside the ColorFunction. Basically I compare #1 to zero using <, >, <=, >=, and == (#1 < 0 &, #1 <= 0 &, etc.). Three of the five look wrong to me. I tried increasing PlotPoints, changed all the “0” to “0.0”. Nothing helped. The use of ColorFunction is copied straight from the documentation of ContourPlot (last example under Options > ColorFunction):
Make everything red above a contour at z==1/2:
ContourPlot[Cos[x] + Cos[y], {x, 0, 4 Pi}, {y, 0, 4 Pi}, Contours -> {1/2}, ColorFunction -> (If[#1 > 1/2, Red, White] &)]
I am using Mathematica 11.0.0.0, Student Edition. I hope I am wrong.
Thanks. it is good to know I was missing something simple. I had run into that scaling thing before, but did not anticipate its effect could be so drastic. It is in that borderline where you could almost call it a bug (at my work we would), but I can live with it. Thanks.