anchor.diff:       open(unit=gfunit,file=trim(output)//trim(t_fmt)//"_"//trim(c_fmt)//".bin",action='write',status='replace',&
                         form='unformatted',access='direct',recl=irec_length_ph)
        !
  -     write(gfunit,iostat=istat,rec=1) nlines_3,iblksize_3,nblock,ibsu
  +     if (verbose>=3) write(gfunit,iostat=istat,rec=1) nlines_3,iblksize_3,nblock,ibsu
  +     !
  +     if (verbose>=3) write(out,"('lines = ',i8,' chunks = ',i8,' block = ',i5)") nswap,nswap / iblksize_3,iblksize_3
        !
        do ib = 1,nswap,blksize
          !
  
positive.code: 
    end subroutine do_stick
  
    subroutine do_gf_oscillator_strength_Phoenix(itrans,ichunk,iso,nswap,nlines,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,gf_ram,abscoeff_ram,gfunit)
       !
       implicit none
       !
       integer(ik),intent(in) :: itrans,ichunk,iso,nswap,nlines,ilevelf_ram(nswap),ileveli_ram(nswap)
       real(rk),intent(in) :: gf_ram(nswap),jrot(nlines),energies(nlines)
       real(rk),intent(in) :: abscoeff_ram(nswap)
       real(rk)  :: gf,acoef,tranfreq,energyf,energyi,gamma1,gamma2,n1,n2,ji,jf,abscoeff
       integer(ik),intent(in) :: gtot(nlines),gfunit
       integer(ik) :: i,ileveli,ilevelf,jp,jpp,wl2ind,il,ib,current_blk,irec_length_ph
       !
       real(rk) :: lambda,ratiolog,g_u
       !
       integer(kind=2) :: iloggf,ilog_n1,ilog_n2,ilog_gamma1,ilog_gamma2,ielion,iener
       !type gausslines
       ! sequence
       ! integer(kind=4) :: iwl
       ! integer(kind=2) :: ielion,ielo,igflog,igr,igs,igw,ign
       !end type
       !
       integer, parameter :: blksize=65536
       type(gauss_linesT) :: oneline(blksize)
  
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)
  