anchor.diff:              j = (i-1)*3 + 1
               LinVec_2(j:j+2) = Mesh2_I%TranslationAcc(:,i) - Mesh2_I_op%TranslationAcc(:,i)          
            end do
  -                  
  -         e=TwoNorm( LinVec_2 - LinVec_2_a )               
  -         !if (.not. equalrealNos(e,0.0_ReKi)) &
  -            write(*,ErrFmt) 'Translational Acceleration: ', n1, e, e/max(mmin, TwoNorm(LinVec_2)), maxval(abs(LinVec_2)), &
  -               mm, maxval(abs(LinVec_1)), maxval(abs(LinVec_1_a)), maxval(abs(LinVec_1_b)), maxval(abs(LinVec_1_c)) 
  -                              
  +
  +         
  +         ! delta u^D
  +         do i=1,Mesh2_I%nnodes
  +            j = (i-1)*3 + 1
  +            LinVec_2_b(j:j+2) = Mesh2_I%TranslationDisp(:,i) - Mesh2_I_op%TranslationDisp(:,i)          
  +         end do
  +         
  +         LinVec_2_a = matmul( map_mod1_mod2%dm%ta_uD,  LinVec_2_b ) + matmul( map_mod1_mod2%dm%ta_uS,  LinVec_1_a ) &
  +                    + matmul( map_mod1_mod2%dm%ta_rv,  LinVec_1_b ) &
  +                    + matmul( map_mod1_mod2%dm%mi,     LinVec_1_c ) + matmul( map_mod1_mod2%dm%fx_p, LinVec_1_d )! approximate delta theta^D         
  +                               
  +         call WrErrorLine('Translational Acceleration', LinVec_2, LinVec_2_a, Mesh1_O_op%TranslationAcc, LinVec_1, LinVec_1_a, LinVec_1_b, LinVec_1_c)
  +                                       
         end do      
         write(*,*) ! blank line
         
  
positive.code:             j = (i-1)*3 + 1
              LinVec_2(j:j+2) = Mesh2_I%RotationVel(:,i) - Mesh2_I_op%RotationVel(:,i)          
           end do
                    
           e=TwoNorm( LinVec_2 - LinVec_2_a )               
           !if (.not. equalrealNos(e,0.0_ReKi)) &
              write(*,ErrFmt) 'Rotational Velocity: ', n1, e, e/max(mmin, TwoNorm(LinVec_2)), maxval(abs(LinVec_2)), &
                 mm , maxval(abs(LinVec_1)) 
        end do            
        
        write(*,*) ! blank line
  
negative.code:             j = (i-1)*3 + 1
              LinVec_1(j:j+2) = Mesh1_I%Force(:,i) - Mesh1_I_op%Force(:,i)   
           end do
            
           call WrErrorLine('Forces', LinVec_1, LinVec_1_a, Mesh2_O_op%Force, LinVec_2)
           
        end do            
        write(*,*) ! blank line
        
  