anchor.diff:    use modCondition
     use modUDF
     use modSUNDIALS
  +  use modSparse
     use iso_c_binding
     
     public
  
positive.code:   
    double precision,allocatable::massFlowSensU(:,:,:) !< sensitivity of mass flow rate on velocity
    double precision,allocatable::massFlowSensP(:,:) !< sensitivity of mass flow rate on pressure
    double precision,allocatable::pbcDiagPrec(:) !< diagonal preconditioner
    
    ! scales and scaling factors
    double precision,allocatable::rhoScale !< density scale [kg/m^3]
  
negative.code:     integer(C_INT)::pbcPSet
      double precision,pointer::x(:),xScale(:),f(:),fScale(:)
      type(C_PTR)::foo
      double precision,parameter::PBC_PREC_DROPTOL=1e-6
      
      call associateVector(c_x,x)
      call associateVector(c_xScale,xScale)
  