// genesis
// copied from targenesis_resonator/Gran_layer_2D_cluster

   include defaults
   include Stell_layer_const.g


int   i
float xpos, ypos, zpos, dx, dy, dz
float randfactor = 3


   echo randomizing stellate cell positions
   for (i = {number_stellate_cells}; i > 0; i = i - 1)

       ce /molecular_layer/Stellate[{i-1}]

       xpos = {getfield . x}
       ypos = {getfield . y}
       zpos = {getfield . z}

       dx = {rand {- stc_separation / {randfactor}} {stc_separation / {randfactor}} } 
       dz = {rand {- stc_separation / 2} {stc_separation / 2} }
       dy = {rand {- stc_separation * 0.5 * {sqrt 3} / {randfactor}} {stc_separation * 0.5 * {sqrt 3} / {randfactor}} }

       position . {xpos + dx} {ypos + dy} {zpos + dz}

       echo {xpos} {ypos} {zpos} {dx} {dy} {dz} {xpos + dx} {ypos + dy} {zpos + dz}

   end