anchor.diff:      call associateVector(newVector,y)
       
       rhoE(1:grid%nC)=x(1:grid%nC)
  -    forall(i=1:grid%nE)
  +    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()
  
positive.code:   
    !> set the boundary conditions
    subroutine setBC()
      double precision::pGst,TGst,uGst(DIMS),Mach,pP(DIMS),Tt,pt
      
      do i=1,grid%nP
        m=grid%iEP(1,i)
  
negative.code:         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
  