light_source{<-0,1000,-100> color rgb 1.*<0.8,0.8,1>} 
light_source{<-200,300,500> color red 0.2}  
light_source{<200,300,-50> color green 0.2} 
#include "functions.inc"          
#include "glass.inc"
#include "consts.inc"
#include "metals.inc"

background{rgb 1}  

//global_settings { max_trace_level 20 }

//fog {distance 200 color rgb <0.2,0.3,1>}

plane {y,-22 pigment{checker rgb 0 rgb 1 scale 50}}

camera {
        location < 0, 115,-215>
        look_at  < 1, 0.5,  0>
        }

//cylinder{0, 1000*y, 20 }

#declare RRR=seed(160823);

//#declare n1=0.25;n*1.5/25-0.25,<(rand(RRR)-0.5)*10,(rand(RRR)-0.5)*10,(rand(RRR)-0.5)*10>
//#declare splinepoints=50;

#declare Norm=<0,0,0>;
#declare numarraypoints=1000;

#declare POINTS=array[numarraypoints+3];
#declare RADIUS=array[numarraypoints+3];

#declare UUU=3;

#declare n=0;
#while(n<numarraypoints+3)

        //#declare POINTS[n]=30*(<cos(2*pi*n/numarraypoints*UUU),0,sin(2*pi*n/numarraypoints*UUU)>)+(n/numarraypoints)*1000*y; 
        #declare POINTS[n]=vaxis_rotate(80*(1-n/numarraypoints)*x, y, UUU*360*n/numarraypoints)+(n/numarraypoints)*100*y;
            
        #declare RADIUS[n]= 1+8*(numarraypoints-n)/numarraypoints+0.1;
        
#declare n=n+1;
#end



#declare MyTrack =
spline {

   cubic_spline
   
   #declare n=-1;
   #while(n<numarraypoints+2)
   
        n/numarraypoints, POINTS[n+1]
           
   #declare n=n+1;
   #end
   
 }

#declare MyWaist =
spline {

   cubic_spline
   
   #declare n=-1;
   #while(n<numarraypoints+2)
   
        n/numarraypoints, RADIUS[n+1]
   
   #declare n=n+1;
   #end
   
 }

//display points       
#declare Q_Display_spheres=no;
#if(Q_Display_spheres)
        #declare n=-1;
        #while(n<2*numarraypoints+2)
   
   //sphere{ POINTS[n+1],.28 pigment{red 1} 

        sphere{      MyTrack(n/numarraypoints/2),
                        MyWaist(n/numarraypoints/2).x*2.1 
   
                        pigment{red 1} 
                        }
                
        #declare n=n+1;
        #end             
#end

#declare npoints=125;  //number of polygons or segments on rail

#declare numrotes=25;

#declare MESHY=
         
mesh2 {
  vertex_vectors {
       6*npoints*numrotes,
       
       #declare n=0;
       #while(n<npoints) 
               #declare nclock=n/npoints+0.0;
               #declare p0=MyTrack(nclock);  
               #declare arr0=vlength(MyWaist(nclock).x);
               
               #declare p1=MyTrack(nclock+1/npoints);
               #declare arr1=vlength(MyWaist(nclock+1/npoints).x);
               #declare p2=MyTrack(nclock+2/npoints);
               

               #declare dr1=p1-p0;  
               //#debug "\n"               #debug str(vlength(p0),6,4)               #debug str(nclock,6,4)               #debug "\n"
               
               
               
               #declare drx1=arr0*vnormalize(vcross(vrotate(dr1,40),dr1));
               #declare dry1=arr0*vnormalize(vcross(drx1-0.000,dr1));                 
                                                                
               #declare dr2=p2-p1;
               #declare drx2=arr1*vnormalize(vcross(vrotate(dr1,40),dr2));
               #declare dry2=arr1*vnormalize(vcross(drx2-0.000,dr2));
               
               
               #declare m=0;
               #while(m<numrotes)
               
               
                        p0+cos(m/numrotes*2*pi)*drx1    +sin(m/numrotes*2*pi)*dry1,
                        p0+cos((m+1)/numrotes*2*pi)*drx1+sin((m+1)/numrotes*2*pi)*dry1,
                        p1+cos((m+1)/numrotes*2*pi)*drx2+sin((m+1)/numrotes*2*pi)*dry2,
                                                                              
                                                                              
                        p0+cos(m/numrotes*2*pi)*drx1    +sin(m/numrotes*2*pi)*dry1,
                        p1+cos((m+0)/numrotes*2*pi)*drx2+sin((m+0)/numrotes*2*pi)*dry2,
                        p1+cos((m+1)/numrotes*2*pi)*drx2+sin((m+1)/numrotes*2*pi)*dry2
                        
                #declare m=m+1;
                #end                        
                              
                                                     
                                                     

               #if(n<npoints-3)
                       ,
               #end
     
                             
       #declare n=n+1;
       #end
  }
  
  normal_vectors 
      {
     6*npoints*numrotes,
     #declare n=0;
     #while(n<npoints) 
              #declare nclock=n/npoints+0.0;
              #declare p0=MyTrack(nclock);  
              #declare arr0=MyWaist(nclock);
              
              #declare p1=MyTrack(nclock+1/npoints);
              #declare arr1=MyWaist(nclock+1/npoints);
              #declare p2=MyTrack(nclock+2/npoints);   
              #declare arr2=MyWaist(nclock+2/npoints);
               

               #declare dr1=p1-p0;  
               //#debug "\n"               #debug str(vlength(p0),6,4)               #debug str(nclock,6,4)               #debug "\n"
               
               
               
               #declare drx1=arr1*vnormalize(vcross(vrotate(dr1,40),dr1));
               #declare dry1=arr1*vnormalize(vcross(drx1-0.0001,dr1));                 
                                                                
               #declare dr2=p2-p1;
               #declare drx2=arr2*vnormalize(vcross(vrotate(dr1,40),dr2));
               #declare dry2=arr2*vnormalize(vcross(drx2-0.0001,dr2));
               
               
               #declare m=0;
               #while(m<numrotes)
               
               
                       /* ( p0+cos(m/numrotes*2*pi)*drx1    +sin(m/numrotes*2*pi)*dry1      )/3+
                        ( p0+cos((m+1)/numrotes*2*pi)*drx1+sin((m+1)/numrotes*2*pi)*dry1  )/3+ 
                        ( p1+cos((m+1)/numrotes*2*pi)*drx2+sin((m+1)/numrotes*2*pi)*dry2  )/3  ,
                         */                                                                         
                                                                                                  
                        p0+cos(m/numrotes*2*pi)*drx1    +sin(m/numrotes*2*pi)*dry1,
                        p0+cos((m+1)/numrotes*2*pi)*drx1+sin((m+1)/numrotes*2*pi)*dry1,
                        p1+cos((m+1)/numrotes*2*pi)*drx2+sin((m+1)/numrotes*2*pi)*dry2,
                                                                              
                                                                              
                        p0+cos(m/numrotes*2*pi)*drx1    +sin(m/numrotes*2*pi)*dry1,
                        p1+cos((m+0)/numrotes*2*pi)*drx2+sin((m+0)/numrotes*2*pi)*dry2,
                        p1+cos((m+1)/numrotes*2*pi)*drx2+sin((m+1)/numrotes*2*pi)*dry2
                        
                #declare m=m+1;
                #end                        
                              
                                                     
                                                     

               #if(n<npoints-3)
                       
               #end
     
                             
       #declare n=n+1;
       #end
      }
  
  face_indices {
     2*npoints*numrotes,
      

               #declare n=0;
               #while (n<6*npoints*numrotes -1)      
     
                       <n,n+1,n+2>,
                       
               #declare n=n+3;
               #end

              
        //<2*npoints*numrotes-3,2*npoints*numrotes-2,2*npoints*numrotes-1>
             
               
               }
  /*
  normal_indices 
      {
        6*npoints*numrotes,
      

               #declare n=0;
               #while (n<2*npoints*numrotes -0)      
     
                       <n,n+1,n+2>
               
                        #if(n<6*npoints*numrotes-3)
                         
                        #end
               
                       
                       
               #declare n=n+3;
               
               
               #end

              
        //<2*npoints*numrotes-3,2*npoints*numrotes-2,2*npoints*numrotes-1>
             
              
      }
    */  
  
  pigment {rgb 0.5}       finish {ambient 0.2 reflection 0.5 specular 0.5}//phong 0.2}   
  //rotate 210*y
  
  //texture { T_Glass1 }//  interior { ior Flint_Glass_Ior }
  //texture{T_Copper_1B}
  //rotate -180*x // *clock
  //translate 3*z 
}



#object{MESHY}