anchor.diff:                 !
                case ('GF')
                  !
  -               ! oscillator strength gf
  +               ! oscillator strength gf; abscoef is to check agains the intensity thresholds 
                  !
  -               abscoef=gfcoef*acoef*gtot(ilevelf)/(vellgt*tranfreq)**2
  +               abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-beta*energyi)*(1.0_rk-exp(-beta*tranfreq))/(tranfreq**2*partfunc)
  +               acoef=gfcoef*acoef*gtot(ilevelf)/(vellgt*tranfreq)**2
                  !
                case ('LIFETIME')
                  !
  
positive.code:          endif
           ilog_gamma2 = int(log(gamma2)*(1.0_rk/(log(10.0_rk)*0.001_rk))) + 2**14
           !
           write(gfunit,*) &
                       ielion,wl2ind,iener,iloggf,ilog_gamma1,ilog_gamma2,ilog_n
       enddo
       !
       !ielion, wl, gf-value, xi, igr, igs, igw
  
negative.code:            !
             if (nchars_>15) cycle
             !
             write(sunit,b_fmt,advance="no") trim(quantum_numbers(kitem,ileveli))
             !
           enddo
           !
           if ( mod(nint(2*Jf),2)==0 ) then
             !
             write(sunit,'(i7,1x,a1,1x,a1)',advance="no") nint(Jf),trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
             write(sunit,'(i7,1x,a1,1x,a1)',advance="no") nint(Ji),trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
             !
           else
             !
             write(sunit,'(f7.1,1x,a1,1x,a1)',advance="no") Jf,trim(quantum_numbers(1,ilevelf)),trim(quantum_numbers(2,ilevelf))
             write(sunit,'(f7.1,1x,a1,1x,a1)',advance="no") Ji,trim(quantum_numbers(1,ileveli)),trim(quantum_numbers(2,ileveli))
             !
           endif
           !
  