anchor.diff:  
     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)
  +  subroutine do_gf_oscillator_strength_Phoenix(itrans,ichunk,iso,nswap,nlines,energies,Jrot,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
  +     real(rk)  :: gf,tranfreq,energyf,energyi,gamma1,gamma2,n1,n2,ji,jf,abscoeff
  +     integer(ik),intent(in) :: gfunit
        integer(ik) :: i,ileveli,ilevelf,jp,jpp,wl2ind,il,ib,current_blk,irec_length_ph
        !
  -     real(rk) :: lambda,ratiolog,g_u
  +     real(rk) :: lambda,ratiolog
        !
        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)
  
positive.code:        oneline(:)%igw = -1
         !oneline(:)%ign = -1
         !
         do il = 1,blksize
           !
           if (ib+il==2) cycle
           i = ib-2+il
           !
           if (i>nswap) exit
           !
           ilevelf = ilevelf_ram(i)
           ileveli = ileveli_ram(i)
  
negative.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,1x,f16.6)") ielion,wl2ind,iener,iloggf,ilog_gamma1,ilog_n1,&
                                      ilog_gamma2,ilog_n2,tranfreq,abscoeff,energyi,lambda
           !
           !write(out,"(f16.6,1x,e16.6,1x,f16.6)") tranfreq,abscoeff,energyi
           !
         enddo loop_block
         !
         current_blk = current_blk+1
         !
  