function s = h3d_create_struct(angular_cell_size_deg, rho_cell_size, rho_min, rho_max)
	s.angular_cell_size_deg = angular_cell_size_deg;
	s.rho_min = rho_min;
	s.rho_max = rho_max;
	s.rho_ncells = ceil( (s.rho_max-s.rho_min)/rho_cell_size );
	s.rho_cell_size = rho_cell_size;
	s.points = [];
	
	% Cube data structure
	
	% one side of the cube is about pi/2 rad
	
	s.cube_ncells =  ceil( (pi/2) / deg2rad(angular_cell_size_deg));
	
	for face=1:6
	for i=1:s.cube_ncells
	for j=1:s.cube_ncells
	
		[face, u, v] = h3d_coords_cell_to_avg_cube(face, i, j, s.cube_ncells);
		normal = h3d_coords_cube_to_s2(face, u, v);
		
		s.cube_normals{face, i, j} = normal;
	end
	end
	end
	
	
	s.cube_ht = zeros(6, s.cube_ncells, s.cube_ncells, s.rho_ncells);
	s.cube_hs = zeros(6, s.cube_ncells, s.cube_ncells);