anchor.diff:              !
           case ('PHOEN')
               !
  -            call do_gf_oscillator_strength_Phoenix(iso,nswap,nlines,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,acoef_ram,abscoef_ram,gfunit)
  +            call do_gf_oscillator_strength_Phoenix(i,ichunk,iso,nswap,nlines,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,acoef_ram,abscoef_ram,gfunit)
               !
               call TimerStop('Calc')
               cycle loop_tran
  
positive.code:    endif
     !
     if ( trim(proftype)=='PHOENIX' ) then
  
       write(ioname, '(a)') 'Phoenix'
       call IOstart(trim(ioname),gfunit)
       open(unit=gfunit,file=trim(output)//".gf",action='write',status='replace',BLOCKSIZE=65536)
       write(gfunit,"('    iel          wv   ener     gf gamma1 gamma2      n')") 
       !
       if (trim(specttype)/="GF") then 
         write(out,"('Change the spectral type to GF (oscillator strength) for Phoenix'/)")
  
negative.code:              energyi = energies(ileveli)
               !
               if (energyf-energyi<-1e1) then
                 write(out,"('Error Ei>Ef: i,f,indi,indf,Aif,Ef,Ei = ',4i12,2x,3es16.8)") ilevelf,ileveli,indexf,indexi,acoef,energyf,energyi
                 stop 'wrong order of indices'
                 cycle loop_swap
               elseif (energyf-energyi<small_) then
  