Python pcolormesh() alternative in Julia

Please explain in detail. I also have to also add contours of psi[1,:,:] and
alfsurf[1,:,:] on plot. rho[t]*v_p[t] is also 270x100 matrix. How to put it in color?