anchor.diff:            halfwidth=dpwcoef*tranfreq
           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+500.0*halfwidth.or.tranfreq<freql-500.0*halfwidth) cycle
  +        if (tranfreq>freqr+ioffset*halfwidth.or.tranfreq<freql-ioffset*halfwidth) cycle
           !
           x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
           !
  
positive.code:             enddo
              !$omp end parallel do 
              !
          case ('MAX');
            !
            ipoint =  max(nint( ( tranfreq-freql)/dfreq )+1,1)
  
negative.code:    use  input
     !
     integer(ik) :: info,ipoint,nlevels,i,itemp,enunit,tunit,sunit,j,ilog,ib,ie,j0,ilevelf,ileveli,indexi,indexf,iline,maxitems,kitem,l
     real(rk)    :: beta,ln2,dtemp,dfreq,temp0,beta0,intband,dpwcoef,x0,tranfreq,tranfreq_i,abscoef,dfreq_,xp,xm,de,lor,b,lor2,dfreq_2,halfwidth0
     real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,elow,jf,ji,acoef,j0rk,alpha
     character(len=cl) :: ioname
     character(wl) :: string_tmp
     !
     real(rk),allocatable :: freq(:),intens(:),jrot(:),pf(:,:),energies(:),Asum(:)
     integer(ik),allocatable :: gtot(:),indices(:)
     character(len=20),allocatable :: quantum_numbers(:,:)
     !
     character(len=9) :: npoints_fmt  !text variable containing formats for reads/writes
  