anchor.diff:  
     ! Now read in the geometry information
     read(iu, iostat=ierr) lnspecies
  -  if ( lnspecies < ispecie ) stop 'Error in reading data, not allocated, nspecies<ispecie'
  +  if ( ispecie < 1 ) call iostat_update(-6)
  +  if ( lnspecies < ispecie ) call iostat_update(-6)
     call iostat_update(ierr)
     allocate(label(lnspecies))
     allocate(zval(lnspecies))
  
positive.code: 
    ! Read overall data
    read(iu, iostat=ierr) lna_u, lno_u, lnspin, nspecies, nsc
    if ( nspecies < ispecie ) call iostat_update(-6)
    call iostat_update(ierr)
  
    read(iu, iostat=ierr) !ucell, Ef, qtot, temp
  
negative.code: 
    do is = 1, nspecies
      do io = 1 , no(is)
        read(iu, iostat=ierr) !n(is,io), l(is,io), zeta(is,io)
        call iostat_update(ierr)
      end do
    end do
  