anchor.diff:          !
           !   half width for Doppler profiling
           if (proftype(1:5)=='DOPPL') then
  +          !
  +          if (tranfreq<small_) cycle
  +          !
             halfwidth=dpwcoef*tranfreq
  +        endif
  +        !
  +        if (offset<0) offset = ioffset*halfwidth
  +        !
  +        x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
  +        !
  +        !   half width for Doppler profiling
  +        if (proftype(1:5)=='VOIGT'.or.proftype(1:5)=='PSEUD') then
  +          if (tranfreq<small_) cycle
  +          halfwidth0=dpwcoef*tranfreq
  +          x0 = sqrt(ln2)/halfwidth0*dfreq*0.5_rk
           endif 
           !
  -        !   half width for Doppler profiling
  -        if (proftype(1:5)=='VOIGT') then
  -          halfwidth0=dpwcoef*tranfreq
  -        endif 
  -        !
           !   if transition frequency is out of selected range
  -        if (tranfreq>freqr+ioffset*halfwidth.or.tranfreq<freql-ioffset*halfwidth) cycle
  -        !
  -        x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
  +        if (tranfreq>freqr+offset.or.tranfreq<freql-offset) cycle
           !
           intband = intband + abscoef
           !
  
positive.code:         !
        endif 
        !
        indices(itemp) = i
        !
        j = 2*int(jrot(i))+1
        energy = energies(i)
  
negative.code:                intens(ipoint)=intens(ipoint)+voigt_humlicek(tranfreq_i,tranfreq,halfwidth,halfwidth0)*abscoef
                 !
              enddo
              !$omp end parallel do      
              !
          case ('MAX');
            !
            ipoint =  max(nint( ( tranfreq-freql)/dfreq )+1,1)
  