anchor.diff:     use  input
      use VoigtKampff_module
      !
  -   integer(ik) :: info,ipoint,nlevels,i,itemp,enunit,tunit,sunit,bunit,j,j0,ilevelf,ileveli,indexi,indexf,iline,maxitems,k,gfunit
  +   integer(ik) :: info,ipoint,nlevels,i,itemp,enunit,tunit,sunit,bunit,j,j0,ilevelf,ileveli,indexi,indexf,iline,maxitems,gfunit
      integer(ik) :: indexf_,indexi_,kitem,nlines,ifilter
      real(rk)    :: beta,ln2,ln22,dtemp,dfreq,temp0,beta0,intband,dpwcoef,tranfreq,abscoef,halfwidth0,tranfreq0
      real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,jf,ji,acoef,j0rk,gfcoef
      real(rk)    :: acoef_,cutoff
  -   integer(ik) :: Jmax,Jp,Jpp,Kp,Kpp,J_,Noffset,Nspecies_,Nvib_states,ivib1,ivib2,ivib,JmaxAll
  +   integer(ik) :: Jmax,Jp,Jpp,J_,Noffset,Nspecies_,Nvib_states,ivib1,ivib2,ivib,JmaxAll
      real(rk)    :: gamma_,n_,gamma_s,ener_vib,ener_rot
      character(len=cl) :: ioname
      !
  
positive.code:          oneline(il)%igw    = ilog_gamma2
           !oneline(il)%ign = ilog_n2
           !
           if (verbose>=4) write(out,"(i7,i12,1x,6(1x,i7),f16.6,1x,e16.6,1x,f16.6)") ielion,wl2ind,iener,iloggf,ilog_gamma1,ilog_n1,&
                                      ilog_gamma2,ilog_n2,tranfreq,abscoeff,energyi
           !
           !write(out,"(f16.6,1x,e16.6,1x,f16.6)") tranfreq,abscoeff,energyi
           !
         enddo
         !
         current_blk = current_blk+1
         !
  
negative.code:        oneline(:)%igw = -1
         !oneline(:)%ign = -1
         !
         if (verbose>=4) write(out,"(' block = ',i8)") current_blk
         !
         loop_block :  do il = 1,blksize
           !
           if (ib+il==2) cycle
           !
           ! the wavenumbers are sorted in increasing order, therefore we start from the last 
           ! to get the libe list sorted with wavelength increasing 
           !
           i = nswap - (ib-1+il)
           !
           if (i<1) exit
           !
           ilevelf = ilevelf_ram(i)
           ileveli = ileveli_ram(i)
  