File:Kepler 186 planet system 1 1 1 1.png

1. include "functions.inc"
2. include "rand.inc"

camera { location <0,0,10>*5 look_at 0 up -3/10

   right -10/10

}

light_source { <-10000,0,10000>*10000 color rgb <1,1,1>*1.5 }

1. declare samples1=2;

1. macro star2(xloc1, yloc1, zloc1,size1,color1, color2)

union {

 sphere {
 0, 1
 pigment { color rgb color2 }
 finish {ambient color2*2}
    scale 0.25
   }

sphere {
 0, 1
 pigment { color rgbf 1 }
 hollow
 interior {
   media {
     emission 10
     method 3
     intervals 2
     samples samples1
    // emission_type 2
    // emission_extinction 2
     density {

// spherical

//function {exp(-sqrt(x*x+y*y+z*z))}

           function { 1/(3000)*pow( sqrt(x*x+y*y+z*z),-4 ) *f_granite(atan(x/y),atan(y/x),atan(z/x))
            *pow( sqrt(x*x+y*y+z*z),-4 )
           }    
               
 //  function {pow(1*f_spherical(x,y,z),4)*0.33 }

// scallop_wave

//    function {pow( (2*(sqrt(x*x+y*y+z*z))) ,-1) } 
     //    function {pow(f_spherical(x,y,z),2)*f_wrinkles(x*1000,y*1000,z*1000) }  
         //      function {f_wrinkles(x*10,y*10,z*10) }  

       color_map {
         [ 0.0 rgb 0.0 ] 
      //      [ 0.2 rgb color1*0.05   ]
       //  [ 0.6 rgb  color1*0.5 ]
         [ 1.0 rgb  color2*2 ]
       }
     }
   }

 // absorbing stuff#declare samples1=30;

 }

} scale size1

	translate x*xloc1

translate y*yloc1 translate z*zloc1 } // juunion

1. end

1. macro print_string(str1, x1, y1)

 text {
   ttf "timrom.ttf" str1 0.1, 0
   pigment { color rgb <1,1,1> } 
   finish {ambient 1}
   scale 3
   translate x*x1
   translate y*y1
   translate z*0.2 
   
 }

1. end

1. macro print_caption(str1, x1, y1, skale1)

 text {
   ttf "timrom.ttf" str1 0.1, 0
   pigment { color rgb <1,1,1> } 
       finish {ambient 1}
   scale 3*skale1
   translate x*x1
   translate y*y1
   translate z*0.3 
   
 }

1. end

1. macro star01 (xloc1, yloc1, zloc1, size1, name1, color1, color2)

object { sphere {0,1} texture { pigment { wrinkles scale 1/10 turbulence 0.5 color_map { [0 color rgb color1] [1 color rgb color2] } }

finish { diffuse 0.6 ambient 1} }

scale size1

translate x*xloc1 translate y*yloc1 translate z*zloc1

}

1. end

1. macro planet01 (xloc1, yloc1,zloc1, size1, name1, color1, color2)

object { sphere {0,1} texture { pigment { wrinkles scale y/3 color_map { [0 color rgb color1] [1 color rgb color2] } }

finish { diffuse 0.6 ambient 1}

}

scale size1 translate x*xloc1 translate y*yloc1 translate z*zloc1 }

1. end

1. macro planet_base_1 (xloc1, yloc1,zloc1, size1, type1)

#declare randi1= SRand(0) ;
#declare randi2= SRand(1) ;
#declare randi3= SRand(2) ;
#declare randi4= SRand(3) ;

object { sphere {0,1} #if(type1="base")

texture {

pigment {

wrinkles scale 1/3

turbulence randi2/5 rotate z*randi2*90

color_map { [0 color rgb <0.741176, 0.486275, 0.423529>] [1 color rgb <0.305882, 0.121569, 0.0862745>]

} }

finish { diffuse 0.6 ambient 0}

}

#end

#if(type1="rock")

texture {

pigment {

wrinkles warp { turbulence 0.2} scale 1/3 octaves 60

turbulence randi2/5 rotate z*randi2*90

color_map { [0 color rgb <0.521569, 0.431373, 0.337255>+randi3/10] [1 color rgb <0.807843, 0.588235, 0.415686>*1.0-randi2/10]

} }

finish { diffuse 0.6 ambient 0 phong 0.05 specular 0.1 }

}

#end

#if(type1="water")

texture {

pigment {

wrinkles scale y/3 turbulence randi2/5 rotate z*randi2*90

color_map { [0 color rgb <0.239216, 0.380392, 0.568627>] [1 color rgb <0.960784, 0.992157, 0.992157>] } }

finish { diffuse 0.6 ambient 0}

}

#end

#if(type1="gas")

texture {

pigment {

wrinkles

turbulence randi2/5 rotate z*randi2*90

scale y/5

color_map { [0 color rgb <0.764706, 0.615686, 0.466667>] [1 color rgb <0.960784, 0.960784, 0.992157>] } }

finish { diffuse 0.6 ambient 0}

}

#end

// objekt modifiers

scale size1 translate x*xloc1 translate y*yloc1 translate z*zloc1 }

1. end

///////// parameters of planet system

//## kepler-186

1. declare numplanets1=5;

1. declare systemname1="Kepler-186";
2. declare names1=array[numplanets1]{"b", "c", "d", "e", "f" }
3. declare distances1=array[numplanets1]{0.0378 , 0.0574 ,0.0861 , 0.1216 , 0.432 }
4. declare radiuses1=array[numplanets1]{1.08,1.25, 1.39, 1.33, 1.17}
5. declare masses1=array[numplanets1]{1.24, 2.1, 2.54, 2.15, 1.71}
6. declare types1=array[numplanets1]{"rock", "rock" ,"rock" ,"rock" ,"rock" }
7. declare radiuscoeff1=2;

//////////////////////////////////////7

1. declare starx1=-35;
2. declare dmax1=45 ; //abs(starx1)*2-30;

1. declare starsize1=5;
2. declare planetoffset1=10.0;

//star01(starx1,0.0,1.0,starsize1,"star1",<1,1,0>, <1,1,0.5>)

//#declare starcolor1 = <1,1,0>; //#declare starcolor2 = <1,1,0.5>;

1. declare starcolor1 = <1.0,0.637,0.34>;
2. declare starcolor2 = <1.0,0.637,0.34>*0.95;

1. declare starsize1 =5;

star2(starx1, 0,0, starsize1,starcolor1, starcolor2)

1. declare beginx1=distances1[0];
2. declare endx1=distances1[numplanets1-1];

1. declare dfx1=log(endx1)-log(beginx1);
2. declare dkx1=dmax1/dfx1;

1. declare color1=<0.368627, 0.317647, 0.258824>;
2. declare color2=<0.956863, 0.701961, 0.470588>;

print_caption(systemname1, -10, 8,1.5)

1. for (nn,0,(numplanets1-1) )

#declare planetx1=distances1[nn]; #declare planetr1=radiuses1[nn]*radiuscoeff1;

#declare dpx1=log(planetx1)-log(beginx1); #declare sitex1=starx1+starsize1+planetoffset1+dpx1*dkx1;

//planet01(sitex1,0,0, planetr1,"b",color1, color2)

#declare type1=types1[nn];

planet_base_1 (sitex1, 0,0, planetr1, type1)

print_string(names1[nn], sitex1-1, -8)

1. end