anchor.diff:  
     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:   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)
        call iostat_update(ierr)
      end do
    end do
  
negative.code:       read(iu, iostat=ierr) (n_specie(io), l_specie(io), zeta_specie(io), io=1,no_specie)
        call iostat_update(ierr)
      else
        read(iu, iostat=ierr) !n(is,io), l(is,io), zeta(is,io)
        call iostat_update(ierr)
      end if
    end do
  