anchor.diff:      ! 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.
  
positive.code:    REAL(BDKi)                  :: b11(3,3)
     REAL(BDKi)                  :: b12(3,3)
     REAL(BDKi)                  :: alpha(3,3)
     INTEGER(IntKi)              :: i, j
  
     INTEGER(IntKi)              :: idx_qp      !< index of current quadrature point
     
  
negative.code: 
     TYPE(BD_InputType)                            :: u_interp                     ! temporary copy of inputs, transferred to BD local system
     REAL(BDKi)                                    :: ScaleFactor                  ! Factor for scaling applied loads at each step
     INTEGER(IntKi)                                :: j                            ! Generic counters
     INTEGER(IntKi)                                :: piter
     REAL(BDKi)                                    :: gravity_temp(3)
  