anchor.diff:        !
       endif
       !
  -    if (use_resolving_power) then 
  +    if (use_resolving_power) then
         !
         npoints0 = nint(real((log(freqr)-log(freql))/resolving_f,rk))+1
         !
         if (npoints0>npoints) then
  -        write(out,"('For the resolving power of ',f15.1,' and range of',2f12.2,' npoints(max) must be > ',i18,'; please increase npoints = ',i15)") resolving_power,freql,freqr,npoints0
  +        write(out,"('For the resolving power of ',f15.1,' and range of',2f12.2)") resolving_power,freql,freqr
  +        write(out,"('npoints(max) must be > ',i18)") npoints0
           write(out,"('Consider increasing npoints > ',i15)") npoints
  -        write(out,"('Npoints(max) = ln(nu2/nu1)/ln(1+1/R)+1')") 
  +        write(out,"('Npoints(max) = ln(nu2/nu1)/ln(1+1/R)+1')")
           stop "Too small number of points for resolving_power and range given!"
         endif
         !
         npoints = npoints0
         !
  -      if (freql<small_) then 
  +      if (freql<small_) then
           write(out,"('use_resolving_power cannot be used for range starting at zero',2f11.3)") freql,freqr
           stop "Illegal use_resolving_power with zero nu1"
  -      endif 
  +      endif
         !
       endif
       !
       !   half width for Doppler profiling
       if (proftype(1:3)=='VOI'.or.proftype(1:3)=='PSE'.or.proftype(1:3)=='LOR'.and.Nspecies>0) then
         !
  -      halfwidth = 0 
  +      halfwidth = 0
         !
         lineprofile_do = .true.
         !
  
positive.code:       call ArrayStart('gtot',info,size(gtot),kind(gtot))
        call ArrayStart('indices',info,size(indices),kind(indices))
        !
        if (trim(specttype)=='LIFETIME') THEN 
          !
          ! allocate the matrix for the sum of the A-coeffs 
          !
          allocate(Asum(nlines),stat=info)
          call ArrayStart('Asum',info,size(Asum),kind(Asum))
  
negative.code:        !        0.000000000248206_rk, 0.000000000000126_rk/)
         !   return
         !  !
         !endif
         !
         ! do only half
         ! m=(n+1)/2
         !
         m=n
  