File:Gaps in gas poor protoplanetary disk 1 r 1.png

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

1. include "stars.inc"

global_settings { assumed_gamma 1.0 }

camera {

 location  <0,100,100>*1.5

 look_at   <0.0, 0.0, 0.0> 
 angle 40

}

/*

// stars sphere { 0,1

    texture{ Starfield2 scale 0.25
              } // 

scale 10000

}

*/  
  
  

1. declare sun1= sphere {

   <0,0,0>,2
   texture {
   
   pigment {
       color White 
       transmit 0
   }

   finish {ambient 3}
   }   

}

light_source {

   <0, 0, 0>

   color Gray75 *50
   fade_distance 5
   fade_power 2  
   
 looks_like {sun1}

}

1. declare fn_Pigm=function {

   pigment {
     agate
     color_map {
       [0 color rgb 0]
       [1 color rgb 1]
     }
   }
 }

1. declare motikka= isosurface {

//function { f_sphere(x, y, z, 0.5) - f_noise3d(x, y, z) } // function { f_sphere(x, y, z, 0.5) - f_noise3d(x * 5, y * 5, z * 5) * 0.5 }

 function { f_sphere(x, y, z, 0.5)-fn_Pigm(x/2, y/2, z/2).gray*0.5 }
 threshold 0
 accuracy 0.005

texture { pigment {color Gray }

       finish{
       diffuse 0.7
       ambient 1.0
       specular 0.0
       reflection 0.0
     }

} }

//#declare numasteroids=0;

//#declare rmax=0; //#declare rmin=0;

1. macro asteroid_zone(rmin, rmax, numasteroids)

1. declare rdelta=rmax-rmin;
2. declare rmid=(rmax+rmin)/2;

1. declare fiidelta=10;

#declare nn=10;

1. local nn=0;

1. declare theta=0.0;
2. declare fii=0.0;

1. while (nn<numasteroids)

   object 
   {
       motikka    
       

       #declare seed1=seed(nn*1);
       #declare seed2=seed(nn*15143);
       #declare seed3=seed(nn*76359);
       #declare seed4=seed(nn*8987547);
       #declare seed5=seed(nn*1543); 
       
       #declare seed6=seed(nn*76923439);
       #declare seed7=seed(nn*87639);
       #declare seed8=seed(nn*76312); 
       
       #declare theta=radians( rand(seed1)*360 );
 
       #declare a1=Rand_Normal( 1.0, 0.33, seed2 );

       #declare rc1=a1*rdelta+rmid;   
        
       #declare a2=Rand_Normal( 1.0, 0.5, seed3 );
        
       
       #declare zc1= rand(seed4)*rc1/10-rc1/20; 
       
       #declare zc2=zc1*a2;
                                  
                                    
       #declare x1=cos(theta)*rc1;
       #declare y1=sin(theta)*rc1;
       #declare z1=zc2;

       #declare de1=rand(seed5);
       
       // fictional size distribution
       #declare skale1=pow(de1,-0.111)*0.2;  
         
       #declare rot1=rand(seed6)*360;
       #declare rot2=rand(seed7)*360;     
       #declare rot3=rand(seed8)*360;
         
     
       #scale skale1
        scale 1
       
       rotate <rot1,rot2, rot3>  
 
       translate <x1,y1,z1>
    }
     
 #local nn=nn+1;

1. end

1. end

1. macro planet1(distance1, radius1, rotation1, color1, color2)

object { sphere {0,1}

scale 0.2*radius1

texture { pigment { gradient y frequency 3 sine_wave color_map { [0 rgb color1] [1 rgb color2] } } finish {ambient 0.1 diffuse 0.5} }

translate x*distance1 rotate z*rotation1

}

1. end

#declare fog1=object {

 sphere // transparent sphere containing media
{ 0,1 pigment { rgbt 1 } hollow
  interior
  { media
    {

samples 40,80

    // emission 1

     scattering {
     2,<1,1/2,1/4>*1/1000
     extinction 0.0001
     }

      density
      {

       //   function { -exp((x*x+y*y+z*z)/1)/10 }
       function {
       f_wrinkles(x,y,z)*-exp((x*x+y*y+z*z)/1)/10
       }

       /*
       //    spherical
       density_map

        { [0 rgb 0]
          [0.5 rgb <1,0,0>]
          [0.8 rgb <1,0.5,0>]
          [1 rgb <1,1,0.5>]
        }
        */

      }
    }
 }

}

scale 200

  scale z/3

}

// scene ...

object {fog1}

1. declare numasteroids=1000;

1. declare rmin=5;
2. declare rmax=10;

asteroid_zone(rmin, rmax, numasteroids)

1. declare numasteroids=1000;
2. declare rmin=20;
3. declare rmax=25;

asteroid_zone(rmin, rmax, numasteroids)

1. declare numasteroids=1000;
2. declare rmin=35;
3. declare rmax=40;

asteroid_zone(rmin, rmax, numasteroids)

planet1(20,15, 30, <1,0.5,0.1>, <1,1,1> )

planet1(35,7, -120, <0,1,1>, <0,0.5,1> )