anchor.diff:       !
     end subroutine do_pse_R
         !
  -  subroutine  do_pse_L(tranfreq,abscoef0,dfreq,freq,halfwidth,offset,freql,dpwcoef,intens)
  +  subroutine  do_pse_L(tranfreq,abscoef,dfreq,freq,halfwidth,offset,freql,dpwcoef,intens)
        !
        implicit none
        !
  -     real(rk),intent(in) :: tranfreq,abscoef0,dfreq,halfwidth,offset,freql,dpwcoef
  -     real(rk),intent(in) :: freq(:)
  -     real(rk),intent(out) :: intens(:)
  -     real(rk) :: alpha,abscoef,dfreq_,de,xp,xm,ln2,b,dnu_half,bnorm,f,eta1,eta2,intens2,halfwidth0
  +     real(rk),intent(in) :: tranfreq,abscoef,dfreq,halfwidth,offset,freql,dpwcoef
  +     real(rk),intent(in) :: freq(npoints)
  +     real(rk),intent(inout) :: intens(npoints)
  +     real(rk) :: alpha,dfreq_,de,xp,xm,ln2,b,dnu_half,bnorm,f,eta1,eta2,intens2,halfwidth0
        real(rk) :: d,x0,intens1
        integer(ik) :: ib,ie,ipoint
        !
  
positive.code:           !
            if (mod(npoints,2)==0) npoints = npoints + 1
            !
          case ("ABSORPTION","EMISSION","LIFETIME","PHOENIX")
            !
            specttype = trim(w)
            !
         case ("MEM","MEMORY")
            !
           call readf(memory_limit)
            !
          case ("ISO","ISOTOPE")
            !
  
negative.code:                !
                 abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-beta*energyi)*(1.0_rk-exp(-beta*tranfreq))/(tranfreq**2*partfunc)
                 !
               case ('emission','EMISSION','EMISS','emiss')
                 !
                 ! emission coefficient [Ergs/mol/Sr]
                 !
                 abscoef=emcoef*acoef*gtot(ilevelf)/real(2*ji+1,rk)*real(2*jf+1,rk)*exp(-beta*energyf)*tranfreq/(partfunc)
                 !
               case ('LIFETIME')
                 !
                 if (Asum(ilevelf)<0) Asum(ilevelf) = 0 
                 !
  