anchor.diff:         dset % translations     = translations
          dset % wyckoffs         = wyckoffs
          dset % equivalent_atoms = equivalent_atoms
  +       dset % crystallographic_orbits = crystallographic_orbits
          dset % std_types        = std_types
          dset % std_positions    = std_positions
   
  
positive.code:        allocate( dset % translations    (3,    n_operations))
         allocate( dset % wyckoffs        (n_atoms))
         allocate( dset % equivalent_atoms(n_atoms))
         allocate( dset % std_types       (n_std_atoms))
         allocate( dset % std_positions   (3, n_std_atoms))
  
  
negative.code:      integer(c_int), allocatable :: wyckoffs(:)
       character(len=7), allocatable :: site_symmetry_symbols(:)
       integer(c_int), allocatable :: equivalent_atoms(:) !Beware mapping refers to positions starting at 0
       integer(c_int), allocatable :: crystallographic_orbits(:) !Beware mapping refers to positions starting at 0
       real(c_double) :: primitive_lattice(3,3)
       integer(c_int), allocatable :: mapping_to_primitive(:) !Beware mapping refers to positions starting at 0
       integer(c_int) :: n_std_atoms
       real(c_double) :: std_lattice(3,3)
       integer(c_int), allocatable :: std_types(:)
       real(c_double), allocatable :: std_positions(:,:)
       real(c_double)  :: std_rotation_matrix(3,3)
       integer(c_int), allocatable :: std_mapping_to_primitive(:) !Beware mapping refers to positions starting at 0
       character(len=6) :: pointgroup_symbol
       integer(kind(SPGLIB_SUCCESS)) :: spglib_error
  
  