anchor.diff:      ! Find number of cells
       geo%n_cells = size(x)
       geo%n_masked = geo%n_cells
  -
  -    print *,geo%n_cells
   
       ! Allocate cells
       allocate(geo%cells(geo%n_cells))
  
positive.code:           geo%mask_map(iv) = ic
         end if
      end do
      
      
      where(geo%volume < 0._dp) geo%volume = 0._dp
  
      geo%n_dim = 3
  
negative.code:     call mp_table_read_column_auto(group, 'cells', 'neighbours', neighbors)
      call mp_table_read_column_auto(group, 'cells', 'bmin', bmin)
      call mp_table_read_column_auto(group, 'cells', 'bmax', bmax)
  
      allocate(x(size(points,2)))
      allocate(y(size(points,2)))
      allocate(z(size(points,2)))
  
      x = points(1,:)
      y = points(2,:)
      z = points(3,:)
  