anchor.diff:              !
           case ('PHOEN')
               !
  -            call do_gf_oscillator_strength_Phoenix(iso,nswap,nlines,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,abscoef_ram,gfunit)
  +            call do_gf_oscillator_strength_Phoenix(iso,nswap,nlines,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,acoef_ram,abscoef_ram,gfunit)
               !
               call TimerStop('Calc')
               cycle loop_tran
  
positive.code:    !
     cmcoef = cmcoef*abundance
     emcoef = emcoef*abundance
     gfcoef = gfcoef*abundance
     !
     !   half width for Doppler profiling
     !
  
negative.code:            !
             if (nchars_>15) cycle
             !
             write(sunit,b_fmt,advance="no") trim(quantum_numbers(kitem,ilevelf))
             !
           enddo
           !
  