anchor.diff:  
      ! specify time increment; currently, all modules will be time integrated with this increment size
   !   dt_global = 1.0D-03
  -   dt_global = 5.0D-03
  +   dt_global = 5.0D-04
   
      n_t_final = ((t_final - t_initial) / dt_global )
   
  
positive.code:    Integer(IntKi)                     :: temp_count
     INTEGER(IntKi),PARAMETER:: QiHUnit = 30
     INTEGER(IntKi),PARAMETER:: LoadUnit = 40
     REAL(ReKi):: InputLoad(20001,2)
  
  
  
  
negative.code: BD_Input(3)%DistrLoad%Force(:,:)  = 0.0D0
  BD_Input(3)%DistrLoad%Moment(:,:) = 0.0D0
  
  BD_Input(1)%PointLoad%Force(2,BD_Parameter%node_total) = InputLoad(n_t_global+1,2)
  BD_Input(2)%PointLoad%Force(2,BD_Parameter%node_total) = InputLoad(n_t_global+2,2)
  BD_Input(3)%PointLoad%Force(2,BD_Parameter%node_total) = InputLoad(n_t_global+3,2)
  
       CALL BeamDyn_CalcOutput( t_global, BD_Input(1), BD_Parameter, BD_ContinuousState, BD_DiscreteState, &
                               BD_ConstraintState, &
  