Polar Stream Plot in Mathematica

Example of PolarStreamPlot in Mathematica

Needs["VectorAnalysis`"]
Clear[field, r, θ, ϕ];
m = Transpose[
     Transpose[JacobianMatrix[#, Spherical @@ #]]/
      ScaleFactors[Spherical @@ #]] &@{r, θ, ϕ};
field[r_, θ_, ϕ_] = 
  Simplify[m.{3/(4 r^4) (3 Cos[θ]^2 - 1), 
     3/(4 r^4) Sin[2 θ], 0}];

StreamPlot[

 Delete[

 field @@ CoordinatesFromCartesian[{x, 0, z}, Spherical], 2], 

 {x, -3, 3}, {z, -3, 3}]