anchor.diff:      allocate(presF(DIMS,grid%nC))
       allocate(presF1(DIMS,grid%nC))
       allocate(condQ(grid%nC))
  +    allocate(flowRho(grid%nP))
       allocate(advRho(grid%nC))
       allocate(advRho1(grid%nC))
       allocate(advRhou(DIMS,grid%nC))
  
positive.code:     use ieee_arithmetic
      use modNumerics
      use modGradient
      use modDiffusion
      use modRhieChow
      type(C_PTR),value::oldVector !< old N_Vector
  
negative.code:     end if
      tmpP(:)=p(1:grid%nC)
      call setBC()
      call findMassFlow(grid,rhou,flowRho)
      call pressureEq%solve(tmpP,info=info)
      needRetry=info<0
      call pressureEq%getNIt(n)
  