#version unofficial megapov 1.21;

#include "stdinc.inc"
#include "metals.inc"
#include "textures.inc"
#include "glass.inc"
#include "math.inc"
#include "rad_def.inc"
#declare Interior_Glass =
  interior {
    ior 1.5
}
//#include "triangulation.inc"

global_settings {

   max_trace_level 10
}

/*
#declare dome =
	pigment {
		image_map {
                  jpeg "night04.jpg"
			map_type 1 
			interpolate 2
		}
	}
*/	
	
sky_sphere {
	//pigment {image_map {jpeg "night04.jpg" map_type 1 interpolate 2 }} 
	pigment {rgb 0}
	rotate <180,120,0>
	translate <0,0,0>
}



#declare stream=seed(225111);
#declare stream2=seed(435111);
#declare stream3=seed(4145111);
#declare stream4 = seed(1231);


union {
  sphere { <0,0,0>, 40 }
  sphere { <0,0,0>, 7 }
  cylinder{
     <0,-40,0>,<0,-6.57,0>,4 
     open
  }
  texture {
      pigment { color rgbf <0.82,0.82,0.82,0.9> }
      finish { F_Glass3 }
   }
  interior { Interior_Glass caustics 1 }
}

  
   #macro pulse(col,pos)
          glow { 
	   	   location pos
		   color rgb col
		   size 0.75/20 
		   type 0
		   fade_power 1.5 
	        } 
	   
   #end 
   
   #macro pulseline(col,spaceing,pos1,pos2,amp,freq)
        #local ii = 1;
        #local pos = pos1;
        #while (ii<vlength(pos2-pos1)/spaceing)
          #local pos = pos1 + ii*(pos2-pos1)/(vlength(pos2-pos1)/spaceing) ;
          #local perp1 = VPerp_To_Vector(pos);
          #local perp2 = VPerp_To_Plane(pos,perp1);
          #local perturb = SRand(stream4)/2;
          #local locamp = (perturb + amp)*vlength(pos)/vlength(pos2-pos1);
          #local pos = pos + locamp*vnormalize(perp1)*sin(2*pi*freq*vlength(pos)) + locamp*vnormalize(perp2)*cos(2*pi*freq*vlength(pos));
          glow { 
	   	   location pos
		   color rgb col
		   size 0.75/20 
		   type 0
		   fade_power 1.5 
	        } 
	      //sphere {pos,0.3 pigment {rgb <0,0,1>} finish {ambient 1 diffuse 0}}  
	  #local ii = ii + 1;
	#end 	   
   #end 
   
   #declare bluep=(0.8)*<0.4,0.4,1>;

#declare i = 0;
#while (i<=10)
 #declare p = vnormalize(<SRand(stream),SRand(stream2),SRand(stream3)>);
 pulseline(bluep,0.1,p*6.5,p*40,1,0.1)
 #declare i = i+1;
#end


camera {
        perspective
        up <0,1,0>
        right -x*image_width/image_height
        location <0,0,600>
        look_at <0,-30,0>
        angle 18 // horizontal FOV angle
        //rotate <0,0,0> //roll
        rotate <-4,0,0> //pitch
        //rotate <0,45,0> //yaw
        //translate <-0.3125,0,0>
} 
light_source {
  <0,10,42> //light position
   color rgb <1,1,1>*.5
    parallel
    point_at <0.0,0>
    //rotate <0,0,0> //roll
    rotate <-25,0,0> //elevation
    rotate <0,45,0> //rotation
}


cylinder{<0,-350,0>,<0,-40,0>,130  pigment {rgb <0.3,0.3,0.3>} finish{reflection {0.1}} normal {agate  agate_turb 0.5 scale <0.0025,0.00005,0.0025>} 
}

cylinder{<0,-350,0>,<0,-40.01,0>,130.1  pigment {rgb 0}
}