anchor.diff:     else
         open(lu, file=file, form='unformatted', iostat=iostat)
         call check_open
  -      write(string,'(a,i0.1)') '# sphere_tesselate version ', st_version
  +      write(string,'(a,i0.1)') '# sphere_tesselate version ', st_neighbours_version
         write(lu) string
  -      write(lu) st_version
  +      write(lu) st_neighbours_version
         write(lu) lnbrs%level
         write(lu) lnbrs%np
         write(lu) lnbrs%nt
  
positive.code:       end subroutine check_open
  
        subroutine check_version
           if (version /= st_version) then
              write(0,'(2(a,i0.1),a)') 'st_load_tesselation: Error: Input file "' // &
                 trim(file) // '" has a different version (', version, &
                 ') to that expected (', st_version ,')'
              error stop
           endif
        end subroutine check_version
  
negative.code:    if (ascii_in) then
        open(lu, file=file, iostat=iostat)
        call check_open
        write(lu,'(a,i0.1)') '# sphere_tesselate version ', st_neighbours_version
        write(lu,'(i0.1," ",i0.1," ",i0.1)') lnbrs%level, lnbrs%np,lnbrs%nt ! Check on level and num points
        do i = 1, lnbrs%np
           write(lu,*) lnbrs%list(i)%idxs
  