<< CSPlotter.m

Contact structure tight at infinity

rmin = 0 + 0.01;
rmax = 3  Pi + 0.05;
dr = (rmax - rmin)/16;
rmax2 = rmax + 5 dr;
scale = dr;
delta = scale 1.5;
conormalright[r_, θ_, z_, 
          nn_] = {0, 0, 1} Cos[Pi/rmax nn r] + 
          Sin[Pi/rmax nn r]/r {-r Sin[θ], r Cos[θ], 
              0};
cotight[r_, θ_, z_] = {0, 0, 1} + 
          r {-r Sin[θ], r Cos[θ], 
              0} 0.80;

outs = 20;

n = 1;
twisted = 
        Table[contactelement[{r Cos[θ], r Sin[θ], 0}, 
            conormalright[r, θ, 0, n], scale], {r, rmin, rmax, 
            dr}, {θ, 0, Pi 2, 2 Pi/outs}];
tight = 
        Table[contactelement[{r Cos[θ], r Sin[θ], 0}, 
            cotight[(r - rmax - delta/2), θ, 0], scale], {r, 
            rmax + delta, rmax2, dr}, {θ, 0, Pi 2, 
            2  Pi/outs}];
out = 
      Show[tight, twisted, PlotRange -> All, 
        Boxed -> False]

Export["TInf.eps", out];
Export["TInf.png", out];
Export["TInf.pdf",  out];

n = 2;
twisted = 
        Table[contactelement[{r Cos[θ], r Sin[θ], 0}, 
            conormalright[r, θ, 0, n], scale], {r, rmin, rmax, 
            dr}, {θ, 0, Pi 2, 2 Pi/outs}];
tight = 
        Table[contactelement[{r Cos[θ], r Sin[θ], 0}, 
            cotight[(r - rmax - delta/2), θ, 0], scale], {r, 
            rmax + delta, rmax2, dr}, {θ, 0, Pi 2, 
            2  Pi/outs}];
out = 
      Show[tight, twisted, PlotRange -> All, 
        Boxed -> False]

Export["T2Inf.eps", out];
Export["T2Inf.png", out];
Export["T2Inf.pdf",  out];