anchor.diff:     ErrStat = ErrID_None
      ErrMsg  = ""
      
  -   
  -   
  +         
      IF(p%analysis_type == 2) THEN
          CALL BD_GA2( t, n, u, utimes, p, x, xd, z, OtherState, ErrStat, ErrMsg )
      ELSEIF(p%analysis_type == 1) THEN
  
positive.code:    p%IniDisp(:) = x%q(:)
     p%IniVelo(:) = x%dqdt(:)
  
  ! Acutator
     p%torq = .TRUE. 
     p%pitchK = 2.0D+07 
     p%pitchC = 5.0D+05
     p%pitchJ = 2.0D+02
     TmpDCM(:,:) = MATMUL(u%RootMotion%Orientation(:,:,1),TRANSPOSE(u%HubMotion%Orientation(:,:,1)))
     temp_CRV(:) = EulerExtract(TmpDCM)
     xd%thetaP = -temp_CRV(3)    
  
negative.code:    TYPE(BD_ContinuousStateType)                       :: x_tmp      ! Holds temporary modification to x
     TYPE(BD_OtherStateType     )                       :: OS_tmp     ! Holds temporary modification to x
     TYPE(BD_InputType)                                 :: u_interp   ! interpolated value of inputs
     INTEGER(IntKi)                                     :: ErrStat2   ! Temporary Error status
     CHARACTER(ErrMsgLen)                               :: ErrMsg2    ! Temporary Error message
     CHARACTER(*), PARAMETER                            :: RoutineName = 'BD_GA2'
  