anchor.diff:     if (ascii_in) then
         open(lu, file=file, iostat=iostat)
         call check_open
  -      write(lu,'(a,i0.1)') '# sphere_tesselate version ', st_version
  +      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
  
positive.code: integer, parameter, private :: lu = 10
  ! File format version descriptor
  integer, parameter, private :: st_version = 1
  
  ! Tolerance for duplicate point searching; this must be less than the closest
  ! expected point spacing.  It is the 3D distance between the points
  
negative.code:    else
        open(lu, file=file, form='unformatted', iostat=iostat)
        call check_open
        write(string,'(a,i0.1)') '# sphere_tesselate version ', st_neighbours_version
        write(lu) string
        write(lu) st_neighbours_version
        write(lu) lnbrs%level
        write(lu) lnbrs%np
        write(lu) lnbrs%nt
  