anchor.diff:     x_tmp%dqdt(4:6,1) = m%u%Rootmotion%RotationVel(:,1)
   
         ! Root velocities/angular velocities and accelerations/angular accelerations
  -      
  -   
  -      
  -      ! Calculate Quadrature point values needed for BldForce results 
  +
  +
  +
  +      ! Calculate Quadrature point values needed for BldForce results
      CALL BD_QuadraturePointData( p,x_tmp,m )   ! Calculate QP values uuu, uup, RR0, kappa, E1
  -   
  - 
  +
  +
   
      IF(p%analysis_type .EQ. BD_DYNAMIC_ANALYSIS) THEN
   
  
positive.code:       ENDIF
  
           ! Apply accelerations using F=ma ?  Is that what this routine does?
           ! Calculate Quadrature point values needed 
        CALL BD_QuadraturePointData( p,x,m )         ! Calculate QP values uuu, uup, RR0, kappa, E1
        CALL BD_GenerateDynamicElementGA2( x, OtherState, u, p, m,fact)
  
  
negative.code: 
                    ! possible type conversions here:
                DCM = BDrot_to_FASTdcm(p%uuN0(4:6,j,i),p)
  
                    ! set the reference position and orientation for each node.
                temp_id = (i-1)*p%nodes_per_elem+j
                CALL MeshPositionNode ( Mesh    = y%BldMotion   &
  