anchor.diff:        omega_N_O_N = matmul(u%Mesh%Orientation(:,:,1),u%Mesh%RotationVel(:,1))
         alpha_N_O_N = matmul(u%Mesh%Orientation(:,:,1),u%Mesh%RotationAcc(:,1))
   
  -   ! Start of the proposed change : Change the TMD DOF Mode to the specific number {1: TMD_X_DOF; 2: TMD_Y_DOF; 3: TMD_XY_DOF} by sm 2015-0922 
         IF (p%TMD_DOF_MODE == ControlMode_None) THEN
        
         ! tmd external forces of dependent degrees:
  
positive.code:    p%K_Y = InputFileData%TMD_Y_K
     p%C_Y = InputFileData%TMD_Y_C
  
     ! Start of the proposed change : Add TMD Mass_XY by sm 2015-0922 
     p%M_XY = InputFileData%TMD_XY_M
     ! End of the proposed change : Add TMD Mass_XY by sm 2015-0922 
     
       ! vector parameters
     ! stop positions
  
negative.code:       B_Y = - rddot_N_N(2) + a_G_N(2) + 1 / p%M_XY * ( OtherState%F_ext(2) + OtherState%F_stop(2) - OtherState%F_table(2) )
        END IF
        
  
         IF (p%TMD_DOF_MODE == ControlMode_None) THEN
              dxdt%tmd_x(1) = 0
              dxdt%tmd_x(2) = 0
  