We consider a passive scalar field under the action of pumping, diffusion and advection by a smooth flow with a Lagrangian chaos. We present theoretical arguments showing that scalar statistics is not conformal invariant and formulate new effective semi-analytic algorithm to model the scalar turbulence. We then carry massive numerics of passive scalar turbulence with the focus on the statistics of nodal lines. The distribution of contours over sizes and perimeters is shown to depend neither on the flow realization nor on the resolution (diffusion) scale $r_d$ for scales exceeding $r_d$. The scalar isolines are found fractal/smooth at the scales larger/smaller than the pumping scale $L$. We characterize the statistics of bending of a long isoline by the driving function of the L\"owner map, show that it behaves like diffusion with the diffusivity independent of resolution yet, most surprisingly, dependent on the velocity realization and the time of scalar evolution.