anchor.diff:  !$omp parallel do  private(i)  
           do i=1,NL
              tau_s(i)=tau_s_l(i)
  -           j_s(i,3)=j_s_l(i,3)
  +           j_s(i,1:3)=j_s_l(i,1:3)
           end do        
           
        else if(Sym /= 1)then
  
positive.code: !$omp parallel do  private(i)  
          do i=1,NL
             tau_s(i)=tau_s_l(i)+tau_s_l(itable_sym(5,i))
             j_s(i,3)=j_s_l(i,3)+j_s_l(itable_sym(5,i),3)
          end do        
          
  
negative.code: !$omp parallel do  private(i)  
          do i=1,NL
             tau_s_l(i)=tau_s(i)+tau_s(itable_sym(2,i))
             j_s_l(i,1)=j_s(i,1)-j_s(itable_sym(2,i),2)
             j_s_l(i,2)=j_s(i,2)-j_s(itable_sym(2,i),1)
             j_s_l(i,3)=j_s(i,3)+j_s(itable_sym(2,i),3)
          end do        
  
  