anchor.diff:     call ArrayStart('swap:ileveli_RAM',info,size(ileveli_RAM),kind(ileveli_RAM))
      call ArrayStart('swap:ilevelf_RAM',info,size(ilevelf_RAM),kind(ilevelf_RAM))
      !
  +   if (lineprofile_do) then
  +     allocate(gamma_RAM(N_to_RAM),stat=info)
  +     call ArrayStart('swap:gamma_RAM',info,size(gamma_RAM),kind(gamma_RAM))
  +   endif
  +   !
      nswap = 0
      Noffset = max(nint(offset/dfreq,ik),1)
      !
  
positive.code:       !
        halfwidth = 0 
        !
        do i=1,Nspecies
          !
          halfwidth =  halfwidth + species(i)%ratio*species(i)%gamma*(species(i)%T0/Temp)**species(i)%N*pressure/species(i)%P0
  
negative.code:             !
              intens_omp = 0
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
                !
                do iswap = iomp,nswap,N_omp_procs
                  !
                  abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  halfwidth = gamma_ram(iswap)
                  !
                  call do_pseud(tranfreq,abscoef,dfreq,freq,halfwidth,offset,freql,dpwcoef,intens_omp(:,iomp))
                  !
  