anchor.diff:                                      betaC,Gd,Xd,Yd,ErrStat2,ErrMsg2)
                 CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
          ENDIF
  -       CALL BD_GyroForce(mEta,rho,uuu,vvv,Fb,ErrStat2,ErrMsg2)
  +       CALL BD_GyroForce(mEta,rho,vvv,Fb,ErrStat2,ErrMsg2)
             CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
   
          DO i=1,node_elem
  
positive.code:    ErrMsg  = ""
  
  
     CALL BD_GenerateDynamicElementForce(uuN0,rrN0,uuN,vvN,aaN,     &
                                         Stif0,Mass0,gravity,u,&
                                         damp_flag,beta,&
                                         elem_total,node_elem,dof_node,ngp,&
                                         quadrature,gw,hhx,hpx,Jaco,uu0,E10,&
  
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,       &
  