anchor.diff:    double precision,allocatable::rho0(:),rhou0(:,:),rhoE0(:),p0(:),u0(:,:),temp0(:),massFrac0(:,:)
     
     ! state at the beginning of an outer iteration
  -  double precision,allocatable::rho1(:)
  -  double precision,allocatable::p1(:)
  +  double precision,allocatable::rho1(:),p1(:)
     
     double precision,allocatable::visc(:) !< viscosity [Pa*s]
     double precision,allocatable::cond(:) !< thermal conductivity [W/m/K]
  
positive.code:     call associateVector(c_fScale,fScale)
      
      ! construct preconditioning matrix
      ! factor preconditioning matrix
      pbcPSet=0
      foo=dat
    end function
  
negative.code:   
    !> momentum and pressure preconditioning matrix setup/factor
    function pbcPSet(c_x,c_xScale,c_f,c_fScale,dat)
      use iso_c_binding
      use modSUNDIALS
      use modAdvection
      use modPressure
      type(C_PTR),value::c_x,c_xScale,c_f,c_fScale,dat
      integer(C_INT)::pbcPSet
      double precision,pointer::x(:),xScale(:),f(:),fScale(:)
  