anchor.diff:          pGst=p(m)
           uGst(:)=-u(:,m)
           TGst=temp(m)
  +        rhoGst=rho(m)
  +        rhouGst=-rhou(:,m)
  +        rhoEGst=rhoE(m)
           if(iBC(n)>0.and..false.)then
             select case(bc%t(iBC(n)))
             case(BC_WALL_TEMP,BC_WALL_TEMP_UDF) ! wall temperature boundary
  
positive.code:         p(n)=pGst
          u(:,n)=uGst(:)
          temp(n)=TGst
          rho(n)=rho(m)!pGst/Rgas/TGst
          rhou(:,n)=-rhou(:,m)!rho(n)*uGst(:)
          rhoE(n)=rhoE(m)!rho(n)*(1d0/(gamm-1d0)*Rgas*TGst+0.5d0*dot_product(uGst,uGst))
        end if
      end do
    end subroutine
  
negative.code:     call associateVector(newVector,y)
      
      rhoE(1:grid%nC)=x(1:grid%nC)
      forall(i=1:grid%nC)
        temp(i)=(rhoE(i)/rho(i)-0.5d0*dot_product(u(:,i),u(:,i)))/(1d0/(gamm-1d0))/Rgas
      end forall
      call setBC()
  