anchor.diff:          Cnew(i) = 0.5 * (-b + SQRT(b*b - 4*c))
         enddo
         !$omp end parallel do
  -    end if     
       else if (gSelect == 5) then
         !$omp parallel do shared(k,Cnew,aux,M,Csol,Mnew)
         do i = 1,n
  
positive.code:         Cnew(i) = 0.5 * (-b + SQRT(b*b - 4*c))
        enddo
        !$omp end parallel do
      end if 
      else if (gSelect == 4) then
        !$omp parallel do private(b,c,f) shared(k,Cnew,aux,M,Csol,Mnew,s,r)
        do i = 1,n
  
negative.code:         Cnew(i) = (k - aux*M(i))*Csol(i) / (aux * Mnew(i) + k)
        enddo
        !$omp end parallel do
      else if (gSelect == 3) then
        !$omp parallel do private(b,c) shared(k,Cnew,aux,M,Csol,Mnew)
        do i = 1,n
  