anchor.diff:             fact = .FALSE.
          ENDIF
   
  -       CALL BD_GenerateDynamicElementGA2(uuN0,rrN0,uuNf,vvNf,aaNf,            &
  +       CALL BD_GenerateDynamicElementGA2(uuNf,vvNf,aaNf,            &
                                            Stif0,Mass0,gravity,u,damp_flag,beta,&
                                            elem_total,node_elem,dof_node,ngp,   &
                                            quadrature,gw,hhx,hpx,Jaco,uu0,E10,fact,&
  
positive.code:                    y%BldMotion%RotationAcc(3,temp_id2) = temp_cc(1)
             ENDDO
         ENDDO
         CALL BD_DynamicSolutionForce(p%uuN0,p%rrN0,x_tmp%q,x_tmp%dqdt,OS_tmp%Acc,                        &
                                      p%Stif0_GL,p%Mass0_GL,p%gravity,u_tmp,                       &
                                      p%damp_flag,p%beta,                                          &
                                      p%node_elem,p%dof_node,p%elem_total,p%dof_total,p%node_total,&
                                      p%ngp,p%quadrature,p%GLw,p%Shp,p%Der,p%Jacobian,p%uu0,p%E10,&
                                      temp_Force,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     ELSEIF(p%analysis_type .EQ. 1) THEN
         CALL BD_StaticSolutionForce( p%uuN0,p%rrN0,x%q,x%dqdt,p%Stif0_GL,p%Mass0_GL,&
                 p%gravity,u_tmp,p%node_elem,p%dof_node,p%elem_total,p%dof_total,p%node_total,&
                 p%ngp,p%quadrature,p%GLw,p%Shp,p%Der,p%Jacobian,p%uu0,p%E10,&
                 temp_Force,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code:          call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
  
     ! find x, acc, and xcc at t+dt
     CALL BD_DynamicSolutionGA2( x%q,x%dqdt,OtherState%acc,OtherState%xcc,&
                                 p%Stif0_GL,p%Mass0_GL,p%gravity,u_interp,              &
                                 p%damp_flag,p%beta,                                    &
                                 p%node_elem,p%dof_node,p%elem_total,p%dof_total,       &
  