anchor.diff:            !
             call readi(verbose)
             !
  -        case ("PARTFUNC","PARTITION-FUNCTION","COOLING")
  +        case ("PARTFUNC","PARTITION-FUNCTION")
             !
             specttype = trim(w)
             !
  
positive.code:                !
                 !print*,ilevelf,indexf,indexi,acoef
                 !
                 cycle 
                 !
               end select
               !
               if (abscoef<abscoef_thresh) cycle
               !
               cutoff = thresh
               !
  
negative.code:      !
    end subroutine do_pse_R
        !
    subroutine  do_pse_L(tranfreq,abscoef,dfreq,freq,halfwidth,offset,freql,dpwcoef,intens)
       !
       implicit none
       !
       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
       !
  