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:   double precision,allocatable::condQ(:) !< heat conduction into cell [W]
    double precision,allocatable::carrier(:,:) !< {rho,rhou,rhoH} of each cell [*]
    double precision,allocatable::advRho(:) !< advection rate of rho [kg/s]
    double precision,allocatable::flowRho(:) !< flow rate of rho through pairs [kg/s]
    double precision,allocatable::fluxRhou(:,:,:) !< flux of rhou [kg/m/s^2]
    double precision,allocatable::advRhou(:,:) !< advection rate of rhou into cell [kg*m/s^2]
    double precision,allocatable::fluxRhoH(:,:) !< flux of rhoH [J/m^2/s]
  