anchor.diff:        pbcXFact((DIMS+1)*(i-1)+(DIMS+1))=1d0/pScale
         pbcRFact((DIMS+1)*(i-1)+1:(DIMS+1)*(i-1)+DIMS)=1d0/rhouScale
         pbcRFact((DIMS+1)*(i-1)+(DIMS+1))=1d0/rhoScale
  +      fullXFact((DIMS+2)*(i-1)+1:(DIMS+2)*(i-1)+DIMS)=1d0/uScale
  +      fullXFact((DIMS+2)*(i-1)+(DIMS+1))=1d0/pScale
  +      fullXFact((DIMS+2)*(i-1)+(DIMS+2))=1d0/tempScale
  +      fullRFact((DIMS+2)*(i-1)+1:(DIMS+2)*(i-1)+DIMS)=1d0/rhouScale
  +      fullRFact((DIMS+2)*(i-1)+(DIMS+1))=1d0/rhoScale
  +      fullRFact((DIMS+2)*(i-1)+(DIMS+2))=1d0/rhoEScale
       end forall
       call pbcEq%setScale(pbcXFact,pbcRFact)
  +    call fullEq%setScale(fullXFact,fullRfact)
       
       nItPBC=0
       nItEnergy=0
  +    nItFull=0
       
       write(*,'(a,g12.6,a,g12.6)')'[i] starting step, t: ',t,' dt: ',dt
     end subroutine
  
positive.code:     end if
    end function
    
  end module
  
negative.code:   integer,parameter::MAXIT_PBC=50 !< max number of momentum and pressure equation iterations
    integer,parameter::MAXIT_ENERGY=50 !< max number of energy equation iterations
    integer,parameter::MAXIT_OUTER=20 !< max number of outer iterations
    integer,parameter::MAXIT_FULL=50 !< max number of full NS equation iterations
    
    integer,parameter::BC_WALL_TEMP=0 !< wall boundary with prescribed temperature
    integer,parameter::BC_WALL_TEMP_UDF=5 !< wall boundary with prescribed temperature by UDF
  