anchor.diff:               
         if (allocated(map_mod1_mod2%dm%mi))   call wrmatrix(map_mod1_mod2%dm%mi,  99,Fmt,'Mi')
         if (allocated(map_mod1_mod2%dm%fx_p)) call wrmatrix(map_mod1_mod2%dm%fx_p,99,Fmt,'fx_p')
  -      if (allocated(map_mod1_mod2%dm%tv))   call wrmatrix(map_mod1_mod2%dm%tv,  99,Fmt,'tv')
  -      if (allocated(map_mod1_mod2%dm%ta1))  call wrmatrix(map_mod1_mod2%dm%ta1, 99,Fmt,'ta1')
  -      if (allocated(map_mod1_mod2%dm%ta2))  call wrmatrix(map_mod1_mod2%dm%ta2, 99,Fmt,'ta2')
  -            
  -      if (allocated(map_mod2_mod1%dm%li ))  call wrmatrix(map_mod2_mod1%dm%li,  99,Fmt,'li')
  -      if (allocated(map_mod2_mod1%dm%M_uD)) call wrmatrix(map_mod2_mod1%dm%M_uD,99,Fmt,'M_uD')
  -      if (allocated(map_mod2_mod1%dm%M_uS)) call wrmatrix(map_mod2_mod1%dm%M_uS,99,Fmt,'m_uS')
  -      if (allocated(map_mod2_mod1%dm%M_f))  call wrmatrix(map_mod2_mod1%dm%M_f, 99,Fmt,'m_f')      
  -      
  -      
  -      write(*,ErrTxtFmt) 'Field','#','Absolute Err','Relative Err','Max perturb:D','Max OP:S','Max delta:S','Max delta:S','Max delta:S','Max delta:S'
  -      write(*,ErrTxtFmt) '-----','-','------------','------------','-------------','--------','-----------','-----------','-----------','-----------'
  +      !if (allocated(map_mod1_mod2%dm%tv_uD))   call wrmatrix(map_mod1_mod2%dm%tv_uD,  99,Fmt,'tv_uD')
  +      !if (allocated(map_mod1_mod2%dm%tv_uS))   call wrmatrix(map_mod1_mod2%dm%tv_uS,  99,Fmt,'tv_uD')
  +      !if (allocated(map_mod1_mod2%dm%ta_uD))  call wrmatrix(map_mod1_mod2%dm%ta_uD, 99,Fmt,'ta_uD')
  +      !if (allocated(map_mod1_mod2%dm%ta_uS))  call wrmatrix(map_mod1_mod2%dm%ta_uS, 99,Fmt,'ta_uS')
  +      !if (allocated(map_mod1_mod2%dm%ta_rv))  call wrmatrix(map_mod1_mod2%dm%ta_rv, 99,Fmt,'ta_rv')
  +      !      
  +      !if (allocated(map_mod2_mod1%dm%li ))  call wrmatrix(map_mod2_mod1%dm%li,  99,Fmt,'li')
  +      !if (allocated(map_mod2_mod1%dm%M_uD)) call wrmatrix(map_mod2_mod1%dm%M_uD,99,Fmt,'M_uD')
  +      !if (allocated(map_mod2_mod1%dm%M_uS)) call wrmatrix(map_mod2_mod1%dm%M_uS,99,Fmt,'m_uS')
  +      !if (allocated(map_mod2_mod1%dm%M_f))  call wrmatrix(map_mod2_mod1%dm%M_f, 99,Fmt,'m_f')      
  +      
  +      
  +      write(*,ErrTxtFmt) 'Field','#','Absolute Err','Percent Err','Max perturb:D','Max OP:S','Max delta:S','Max delta:S','Max delta:S','Max delta:S'
  +      write(*,ErrTxtFmt) '-----','-','------------','-----------','-------------','--------','-----------','-----------','-----------','-----------'
         ! ..................
         ! Rotational displacement:
         ! ..................
         
  -      !call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
  +      call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
  +      call meshcopy( Mesh2_O_op, Mesh2_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
         do n1=1,15
            
            ! perturb x^S:
            call getRotationPerturb(LinVec_1)
                     
            LinVec_2_a = matmul( map_mod1_mod2%dm%mi, LinVec_1 ) ! approximate delta theta^D         
  -             
  -         !call wrmatrix(LinVec_1,  99,Fmt,'LinVec_1')
  -         !call wrmatrix(LinVec_2_a,  99,Fmt,'LinVec_2_a')
  -         
  +           
            do i=1,Mesh1_O%nnodes
               j = (i-1)*3 + 1
   
                  ! delta theta^S
  -            !Orientation = EulerConstruct(LinVec_1(j:j+2))
               call SmllRotTrans( 'orientation', LinVec_1(j  ) &
                                               , LinVec_1(j+1) &
                                               , LinVec_1(j+2) & 
  
positive.code:       if (allocated(LinVec_1_a)) deallocate(LinVec_1_a)
        if (allocated(LinVec_1_b)) deallocate(LinVec_1_b)
        if (allocated(LinVec_1_c)) deallocate(LinVec_1_c)      
        if (allocated(LinVec_2)) deallocate(LinVec_2)      
        if (allocated(LinVec_2_a)) deallocate(LinVec_2_a)      
        if (allocated(LinVec_2_b)) deallocate(LinVec_2_b)      
  
negative.code:             j = (i-1)*3 + 1
              LinVec_2(j:j+2) = Mesh2_I%RotationAcc(:,i) - Mesh2_I_op%RotationAcc(:,i)          
           end do
               
           call WrErrorLine('Rotational Acceleration', LinVec_2, LinVec_2_a, Mesh1_O_op%RotationAcc, LinVec_1)
  
        end do            
        write(*,*) ! blank line
        
        ! ..................
        ! Translational acceleration:
        ! ..................
  
        call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
        call meshcopy( Mesh2_O_op, Mesh2_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
  