anchor.diff:      end if
       tmpP(:)=p(1:grid%nC)
       call setBC()
  -    call findMassFlow(grid,rhou,flowRho)
  +    call findMassFlow(grid,rhou,flowRho1)
       call pressureEq%solve(tmpP,info=info)
       needRetry=info<0
       call pressureEq%getNIt(n)
  
positive.code:     integer(C_INT)::pressureRHS !< error code
      double precision,pointer::x(:) !< fortran pointer associated with oldVector
      double precision,pointer::y(:) !< fortran pointer associated with newVector
      double precision,allocatable,save::tmpFlowRho(:),tmpAdvRho(:)
      
      call associateVector(oldVector,x)
      call associateVector(newVector,y)
  
negative.code:   double precision,allocatable::rho0(:),rhou0(:,:),rhoE0(:),p0(:),u0(:,:),temp0(:),Y0(:,:)
    
    ! auxiliary variables for the PISO iteration
    double precision,allocatable::rho1(:),p1(:),laP1(:),presF1(:,:),temp1(:),rhou1(:,:),flowRho1(:)
    
    ! scales
    double precision,allocatable::rhoScale !< density scale [kg/m^3]
  