anchor.diff:    call iostat_update(ierr)
     do is = 1, nspecies
       do io = 1 , no(is)
  -      read(iu, iostat=ierr) !nfio(is,io), lfio(is,io), zetafio(is,io)
  +      read(iu, iostat=ierr) !n(is,io), l(is,io), zeta(is,io)
         call iostat_update(ierr)
       end do
     end do
  
positive.code: 
    do is = 1, nspecies
      do io = 1 , no(is)
        read(iu, iostat=ierr) !nfio(is,io), lfio(is,io), zetafio(is,io)
        call iostat_update(ierr)
      end do
    end do
  
negative.code: 
    read(iu, iostat=ierr) isc, xa, isa!, lasto
    call iostat_update(ierr)
    deallocate(isc, xa)
  
    read(iu, iostat=ierr) (label(1:20,is), zval(is),no(is),is=1,nspecies)
    call iostat_update(ierr)
  