anchor.diff:     ! local variables
      INTEGER(IntKi)          :: i                ! do-loop counter
      INTEGER(IntKi)          :: j                ! do-loop counter
  -   INTEGER(IntKi)          :: idx_qp           !< index of current quadrature point in loop
      INTEGER(IntKi)          :: member_first_kp
      INTEGER(IntKi)          :: member_last_kp
  -   INTEGER(IntKi)          :: temp_id2
      REAL(BDKi)              :: eta
      REAL(BDKi)              :: temp_POS(3)
      REAL(BDKi)              :: temp_CRV(3)
  
positive.code:     ! local variables
      INTEGER(IntKi)                                 :: i
      INTEGER(IntKi)                                 :: idx_dof
      REAL(BDKi), allocatable                        :: RHS_m(:,:), RHS_p(:,:)
      CHARACTER(*), PARAMETER                        :: RoutineName = 'BD_FD_Stat'
  
      ! zero out the local matrices.
  
negative.code:    REAL(BDKi)                  :: b11(3,3)
     REAL(BDKi)                  :: b12(3,3)
     REAL(BDKi)                  :: alpha(3,3)
  
     INTEGER(IntKi)              :: idx_qp      !< index of current quadrature point
     
  