anchor.diff:    double precision,allocatable::tempScale !< temperature scale [K]
     double precision,allocatable::pbcXFact(:) !< pbc solution scaling factor vector
     double precision,allocatable::pbcRFact(:) !< pbc residual scaling factor vector
  +  double precision,allocatable::fullXFact(:) !< full system solution scaling factor vector
  +  double precision,allocatable::fullRFact(:) !< full system residual scaling factor vector
     
     ! data for algebraic solver
     type(NewtonKrylov)::pbcEq !< momentum and pressure equations (isothermal) solved by Newton-GMRES
     type(fixPt)::energyEq !< energy equation as a fix point problem
  -  integer::nItPBC,nItEnergy,nItOuter !< number of iterations
  +  type(NewtonKrylov)::fullEq !< full NS equations solved by Newton-GMRES
  +  integer::nItPBC,nItEnergy,nItOuter,nItFull !< number of iterations
     
   contains
     
  
positive.code:     end if
    end function
    
  end module
  
negative.code:     allocate(condQ(grid%nC))
      allocate(pbcXFact((DIMS+1)*grid%nC))
      allocate(pbcRFact((DIMS+1)*grid%nC))
      allocate(fullXFact((DIMS+2)*grid%nC))
      allocate(fullRFact((DIMS+2)*grid%nC))
      do i=1,grid%nE
        if(udfIc/=0)then
          pE(:)=grid%p(:,i)
  