#declare Earth_Radius=1;

#declare Media_Intensity=50;

//#declare Media_Emission=0.065*10;

1. declare Media_Eccentricity=0.56;

//#declare Media_Eccentricity=0.3;

1. declare Atmosphere_Top=0.04;

//#declare Cloud_Brightness=0.42; //#declare Light_Intensity=4.1;

1. declare atm_samples= 20; // must be 20+

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

// pov-ray 3.8 source code

// giant planet 

// 27.9.2022 v 0000.0002

default { finish { ambient 0.000000 diffuse 0.7 } }    
 

camera {

 location  <0,0,-1>*5
 angle 35
 look_at   0

}

light_source {

<-3/3,-0,-1>*1000000

color rgb <1,1,1> *1.5 }

1. declare dark1= texture {

     pigment {
    wrinkles  scale 0.05
    // granite
 //  gradient y   
  // frequency 2 
   scale y*0.2 
  
  // sine_wave    
   
   
   turbulence 0.2
 //  omega 0.1
  // lambda 0.1

/*

   0.97, 1.00, 1.00  // white
  0.72, 0.50, 0.27 // ora
 0.55, 0.38, 0.30  // red
 0.44, 0.42, 0.32 // gray

• /

   pigment_map {
   
        [0.0 color rgb <0.55, 0.38, 0.30 >*1.5  ]   
      [1.0 color rgb <0.91, 0.65, 0.45 >*1.05  ]  

// [0.0 color rgb <.85,.8,.725>*0.95 ]

 //      [1.0 color rgb <.85,.8,.725>*1.05  ]  

  //     [0.0 agate color_map{ [ 0 color rgb <.85,.8,.725>*0.9 ] [ 1 color rgb <.85,.8,.725>*1.0 ] } ]
   //    [1.0 granite  color_map{ [ 0 color rgb <.85,.8,.725>*0.9 ] [ 1 color rgb <.85,.8,.725>*1.1 ] } ]
    }         
    
    }

   normal {
  granite 
  //     gradient y   
//   frequency 2 
  // sine_wave  
    scale y*0.1
    turbulence 0.1
   //granite scale y*0.03 bump_size 0.05
      //    turbulence 0.5
//   omega 0.1
 //  lambda 0.1
   bump_size 0.01
   
    }    
   finish {ambient .00 diffuse .55 phong .075 phong_size 1.5 specular .025 roughness .01}

}

1. declare light1 = texture

{

 // pigment { color rgb <0.9,1,1> }     
        pigment {

// granite

 wrinkles   

scale 0.2

 // agate
  // gradient y   
   turbulence 0.2
//   frequency 2   
   scale y*0.15
 //  sine_wave

/*

   0.97, 1.00, 1.00  // white
  0.72, 0.50, 0.27 // ora
 0.55, 0.38, 0.30  // red
 0.44, 0.42, 0.32 // gray

• /

   pigment_map {
   
//   [0.0 color rgb  <1.00, 0.99, 0.94> ]
  //     [1.0  color rgb <1.00, 0.92, 0.82>*0.9  ]
   
     //     [0.0 agate color_map{ [ 0 color rgb <0.9,1,1>*1 ] [ 1 color rgb <0.9,1,1>*1.1 ] } ]
     //      [1.0  color rgb <0.9,1,1>*0.95  ]
      [1.0 granite  color_map{ [ 0 color rgb <0.97,1,1>*0.8 ] [ 1 color rgb <1.00, 0.99, 0.94>*1.2 ] } ]
          [1.0  color rgb <0.97,1,1>*1.05  ]
    }         
    
    }
  
  normal {
      granite
  // gradient y   
   turbulence 0.2
   frequency 2  
   scale y*0.05
  //granite scale y*0.03 bump_size 0.05
  //    gradient y   
 //  frequency 10 
 //  sine_wave
      bump_size -0.003
  }  
  finish {ambient .00 diffuse .65 phong .075 phong_size 1.5 specular .025 roughness .01}

}

#macro tex1_storms1()
                  
           

1. declare seed1 = seed(2531);

scale y/2

   #declare radi1=0.99;

    #declare anx1=radians(180+20);
   #declare any1=radians(20);
 
   #declare theta1=anx1;
   #declare fii1=any1;

   #declare x1=radi1*sin(theta1)*cos(fii1);
   #declare y1=radi1*sin(theta1)*sin(fii1);
   #declare z1=radi1*cos(theta1);

        warp {black_hole <x1,y1,z1>,0.35 falloff 1 strength 2 turbulence <1,1,1>*10 inverse } 
     #declare radi1=0.99;
 

1. for (n,1,10)

    #declare j1=rand(seed1)*360+20;
         #declare j2=rand(seed1)*90;
 #declare anx1=radians(180-j1);
   #declare any1=radians(j2);
 
   #declare theta1=anx1;
   #declare fii1=any1;

   #declare x1=radi1*sin(theta1)*cos(fii1);
   #declare y1=radi1*sin(theta1)*sin(fii1);
   #declare z1=radi1*cos(theta1);

        warp {black_hole <x1,y1,z1>,0.15*rand(1) falloff 2 strength 2.5 turbulence 0.4 inverse  } 
   #end

  //    scale y/2

1. end

#declare tex1= texture {


    function { ((sin(abs(y*2)*pi)/2)+.5)*f_wrinkles(0,y*2,0)*f_wrinkles(0,y*1,0)*1 }

scallop_wave

// sine_wave

  //
      //     frequency 2  

//

 //      turbulence .03    

          scale 3
       warp {turbulence <0.07, 0.3, 0.07>*1}
           scale 1/3
               scale 10
       warp {turbulence <0.1, 0.2, 0.1>*2}
           scale 1/10   
  
                   scale 30
       warp {turbulence <0.1, 0.2, 0.1>*1}
           scale 1/30                       

                       scale 100
       warp {turbulence <0.1, 0.2, 0.1>*1}
           scale 1/100                       

       texture_map {
       [0 dark1 ]
         [1 light1  ]   
      //    [1 dark1  ]   

       }
       

 tex1_storms1()     

    scale y/2

}

#declare tex2= texture {


       gradient y 

//   scale pi/2*0.01

sine_wave

// scallop_wave                   
   scale 2

translate y*1/2

// turbulence 0.4

       texture_map {  
             [0 pigment {rgbt <1,1,1,1> }      finish {ambient .00 diffuse .55 phong .075 phong_size 1.5 specular .025 roughness .01} ] 
             [1 pigment {rgbt <0.95, 0.85, 0.67,0.5 > }      finish {ambient .00 diffuse .55 phong .075 phong_size 1.5 specular .025 roughness .01} 

// normal {granite scale 0.1 bump_size 0.05 turbulence 0.4 } ]

       }
       

}

1. declare Density1=density { function

    {

    exp(-6.7*(sqrt(x*x+y*y+z*z) -1.00001)/0.05 )/1000 
    
    }   
    }
 
 
 

1. declare Mat_Atm =

 material {
   texture {
     pigment {
       color rgbt <1.0, 1.0, 1.0, 1.0>
     }
   }
   interior {
     media {
       method 3  
                      //   emission <0.73, 0.65*0.8, 0.46>*30    
//   emission  <0.73, 0.65*0.8, 0.46>
               scattering { 4 color rgb  <0.4, 0.6, 1>*Media_Intensity  }

          //    scattering { 5 color rgb  <0.73, 0.65*0.8, 0.46>*Media_Intensity eccentricity Media_Eccentricity }

      //    scattering { 5 color rgb  <0.2, 0.4, 1>*Media_Intensity eccentricity Media_Eccentricity }

         samples atm_samples  // increase = more precise
 //intervals 1
       
       density { 

        Density1

       }
     }
   }
 }
          
          
          

1. declare atmos1=difference {

 sphere {
   <0,0,0>, 1
 }
 sphere { 
   <0,0,0>, 1
   scale (Earth_Radius+0.001)/(Earth_Radius+Atmosphere_Top)
 }
 material { Mat_Atm }
 scale Earth_Radius+Atmosphere_Top
 hollow on

}

//background{rgb .1}

// giant planet

1. declare giantplanet1=sphere {0,1

      texture {
    
       tex1
      }

}

1. declare giantplanet2=sphere {0,1

      texture {
    
       tex2
      }

}

1. declare gigaplanet1= union {

   object { giantplanet1}

   object { giantplanet2 scale 1.002}

 object {atmos1}  
  //   rotate y*-35
 rotate x*40

   }

object {gigaplanet1}