anchor.diff:                  call do_lorentz(tranfreq,abscoef,dfreq,freq,halfwidth,offset,freql,intens_omp(:,iomp))
                   !
                 enddo
  -              !              
  +              !
               enddo
               !$omp end parallel do
               !
  
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:                 call do_rect(tranfreq,abscoef,dfreq,freq,offset,freql,intens_omp(:,iomp))
                  !
                enddo
                !
              enddo
              !$omp end parallel do
              !
  