anchor.diff:        p%torqM(1,2) =  p%pitchJ * p%dt
         p%torqM(2,2) =  p%pitchJ
         denom        =  p%pitchJ + p%pitchC*p%dt + p%pitchK*p%dt**2
  -      if (EqualRealNos(denom,0.0_ReKi)) then
  +      if (EqualRealNos(denom,0.0_BDKi)) then
            call SetErrStat(ErrID_Fatal,"Cannot invert matrix for pitch actuator: J+c*dt+k*dt^2 is zero.",ErrStat,ErrMsg,RoutineName)
         else         
            p%torqM(:,:) =  p%torqM / denom
  
positive.code:    RETURN
  
  END SUBROUTINE BD_StaticElementMatrixForce
  
  SUBROUTINE BD_GA2(t,n,u,utimes,p,x,xd,z,OtherState,ErrStat,ErrMsg)
  !-------------------------------------------------------
  ! This subroutine performs time marching from t_i to t_f
  
negative.code:            temp_CRV(3) = temp_POS(1)
             
             CALL BD_CrvMatrixR(temp_CRV,TmpDCM,ErrStat2,ErrMsg2)
                CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )  
                
                ! possible type conversions here:
             DCM = TRANSPOSE(TmpDCM)
             
             temp_id = (i-1)*p%node_elem+j
             CALL MeshPositionNode ( Mesh    = y%BldMotion   &
                                    ,INode   = temp_id       &
                                    ,Pos     = Pos           &
                                    ,ErrStat = ErrStat2      &
                                    ,ErrMess = ErrMsg2       &
                                    ,Orient  = DCM           )
  