anchor.diff:     subroutine CreateOutputMeshes_Test5Orient
      
         real(reki) :: z
  -      real(reki) :: DCM(3,3)
  +      real(Dbki) :: DCM(3,3)
         
         
         Mesh1Type = ELEMENT_LINE2
  
positive.code:       CALL MeshMapCreate( Mesh1_O, Mesh2_I,     Map_Mod1_Mod2, ErrStat, ErrMsg );       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))   
        CALL MeshMapCreate( Mesh2_O, Mesh1_I,     Map_Mod2_Mod1, ErrStat, ErrMsg );       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        
        CALL MeshMapCreate( Mesh1_I, Mesh1_I_1PT, Map_Mod1_I_1PT,  ErrStat, ErrMsg );       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        CALL MeshMapCreate( Mesh2_O, Mesh2_O_1PT, Map_Mod2_O_1PT,  ErrStat, ErrMsg );       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
              
                    
        ! ..............................................................................................................................   
        ! Map the outputs to inputs:
        ! ..............................................................................................................................   
        IF (Mesh1Type == ELEMENT_POINT ) THEN
           IF ( Mesh2Type == ELEMENT_POINT ) THEN
              CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
              CALL Transfer_Point_to_Point( Mesh2_O, Mesh1_I, Map_Mod2_Mod1, ErrStat, ErrMsg, Mesh2_I, Mesh1_O );   IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))   
           ELSEIF ( Mesh2Type == ELEMENT_LINE2) THEN                                                                
              CALL Transfer_Point_to_Line2( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
              CALL Transfer_Line2_to_Point( Mesh2_O, Mesh1_I, Map_Mod2_Mod1, ErrStat, ErrMsg, Mesh2_I, Mesh1_O );   IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))         
  
negative.code:       ! ..............................................................................................................................   
        ! Write results to file(s)
        ! ..............................................................................................................................   
        call WriteMappingTransferToFile(Mesh1_I, Mesh1_O, Mesh2_I, Mesh2_O, Map_Mod1_Mod2, Map_Mod2_Mod1, BinOutputName)                       
        
        if (TestNumber == 13) then
  
           !write( 17, '(A)' )
           DO i=1,mesh2_o%nnodes
              n1 = mesh1_o%ElemTable(ELEMENT_LINE2)%Elements(Map_Mod1_Mod2%MapMotions(i)%OtherMesh_Element)%ElemNodes(1)
              n2 = mesh1_o%ElemTable(ELEMENT_LINE2)%Elements(Map_Mod1_Mod2%MapMotions(i)%OtherMesh_Element)%ElemNodes(2)
        
              WRITE (17, '(6(F15.5))' ) atan2( mesh2_i%refOrientation(1,2,i ) , mesh2_i%refOrientation(1,1,i ) ),&
                                        atan2( mesh2_i%Orientation(   1,2,i ) , mesh2_i%Orientation(   1,1,i ) ),&
                                        atan2( mesh1_o%refOrientation(1,2,n1) , mesh1_o%refOrientation(1,1,n1) ),&
                                        atan2( mesh1_o%Orientation(   1,2,n1) , mesh1_o%Orientation(   1,1,n1) ),&      
                                        atan2( mesh1_o%refOrientation(1,2,n2) , mesh1_o%refOrientation(1,1,n2) ),&
                                        atan2( mesh1_o%Orientation(   1,2,n2) , mesh1_o%Orientation(   1,1,n2) )
           end do                              
     
  
  
        end if
       
  if ( TestLinearization ) then      
        ! ..............................................................................................................................   
        ! Linearize about the operating points:
        ! ..............................................................................................................................   
        
        call AllocAry( LinVec_1, Mesh1_O%nnodes*3, 'LinVec_1', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_2, Mesh2_O%nnodes*3, 'LinVec_2', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_theta, Mesh1_O%nnodes*3, 'LinVec_1_theta', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_2_theta, Mesh2_O%nnodes*3, 'LinVec_2_theta', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_a, Mesh1_O%nnodes*3, 'LinVec_1_a', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_b, Mesh1_O%nnodes*3, 'LinVec_1_b', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        call AllocAry( LinVec_1_c, Mesh2_O%nnodes*3, 'LinVec_1_c', ErrStat, ErrMsg);       IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        
        
        CALL MeshCopy (        SrcMesh          = mesh1_O          &
                             , DestMesh         = mesh1_O_op       &
                             , CtrlCode         = MESH_NEWCOPY     &
                             , ErrStat          = ErrStat          &
                             , ErrMess          = ErrMsg           )
              IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        
        CALL MeshCopy (        SrcMesh          = mesh2_O          &
                             , DestMesh         = mesh2_O_op       &
                             , CtrlCode         = MESH_NEWCOPY     &
                             , ErrStat          = ErrStat          &
                             , ErrMess          = ErrMsg           )
              IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        
              
        CALL MeshCopy (        SrcMesh          = mesh1_I          &
                             , DestMesh         = mesh1_I_op       &
                             , CtrlCode         = MESH_NEWCOPY     &
                             , ErrStat          = ErrStat          &
                             , ErrMess          = ErrMsg           )
              IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
        
        CALL MeshCopy (        SrcMesh          = mesh2_I          &
                             , DestMesh         = mesh2_I_op       &
                             , CtrlCode         = MESH_NEWCOPY     &
                             , ErrStat          = ErrStat          &
                             , ErrMess          = ErrMsg           )
              IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
              
        
  write(99, *) '---------------------------------------------------------------'
  write(99, *) '   Test ', TestNumber
  write(99, *) '---------------------------------------------------------------'
        
        IF (Mesh1Type == ELEMENT_POINT ) THEN
           IF ( Mesh2Type == ELEMENT_POINT ) THEN            
        
              CALL Linearize_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
              CALL Linearize_Point_to_Point( Mesh2_O, Mesh1_I, Map_Mod2_Mod1, ErrStat, ErrMsg, Mesh2_I, Mesh1_O );   IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))   
                       
              call wrmatrix(map_mod1_mod2%dm%mi,  99,Fmt,'Mi')
              call wrmatrix(map_mod1_mod2%dm%fx_p,99,Fmt,'fx_p')
              call wrmatrix(map_mod1_mod2%dm%tv,  99,Fmt,'tv')
              call wrmatrix(map_mod1_mod2%dm%ta1, 99,Fmt,'ta1')
              call wrmatrix(map_mod1_mod2%dm%ta2, 99,Fmt,'ta2')
              
  
              call wrmatrix(map_mod2_mod1%dm%li,99,Fmt,'li')
              call wrmatrix(map_mod2_mod1%dm%M1,99,Fmt,'m1')
              call wrmatrix(map_mod2_mod1%dm%M2,99,Fmt,'m2')      
              
              
           ELSEIF ( Mesh2Type == ELEMENT_LINE2) THEN                                                                
              !CALL Linearize_Point_to_Line2( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
              !CALL Linearize_Line2_to_Point( Mesh2_O, Mesh1_I, Map_Mod2_Mod1, ErrStat, ErrMsg, Mesh2_I, Mesh1_O );   IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))         
           END IF                                                                                                   
        ELSEIF ( Mesh1Type == ELEMENT_LINE2 ) THEN                                                                  
           IF ( Mesh2Type == ELEMENT_LINE2 ) THEN                                                                   
              !CALL Linearize_Line2_to_Line2( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))   
              !if (TestNumber == 5 .or. TestNumber == 12 ) call InitTest5Loads()
              !if (TestNumber == 13 ) call InitTest13Loads()
              !CALL Linearize_Line2_to_Line2( Mesh2_O, Mesh1_I, Map_Mod2_Mod1, ErrStat, ErrMsg, Mesh2_I, Mesh1_O );   IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))        
           ELSEIF ( Mesh2Type == ELEMENT_POINT ) THEN        
              !CALL Linearize_Line2_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))                  
              !CALL Linearize_Point_to_Line2( Mesh2_O, Mesh1_I, Map_Mod2_Mod1, ErrStat, ErrMsg, Mesh2_I, Mesh1_O );   IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
           END IF            
        END IF
        
        
        ! src=Mesh1_O, dest=Mesh2_I, Map_Mod1_Mod2
     if (Mesh1Type == ELEMENT_POINT  .and.  Mesh2Type == ELEMENT_POINT ) THEN            
        
        
        !........................ 
        ! get Mesh1_O_op orientation parameters, i.e., LinVec_1_theta = theta^S
        !........................
        j=1
        do i=1,Mesh1_O%nnodes
           LinVec_1_theta(j:j+2) = GetSmllRotAngs(Mesh1_O_op%Orientation(:,:,i),ErrStat,ErrMsg)
           j = j+3
           !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
        end do
        
        !call wrmatrix(LinVec_1_theta,  99,Fmt,'LinVec_1_theta')
           
        ! ..................
        ! Rotational displacement:
        ! ..................
        
        do n1=1,15
           
           ! perturb x^S:
           call random_number( LinVec_1 )
           LinVec_1 = LinVec_1/real(n1**3)               
           
           LinVec_2_theta = matmul( map_mod1_mod2%dm%mi, LinVec_1 ) ! approximate delta theta^D         
               
           !call wrmatrix(LinVec_1,  99,Fmt,'LinVec_1')
           !call wrmatrix(LinVec_2_theta,  99,Fmt,'LinVec_2_theta')
           
           j=1
           do i=1,Mesh1_O%nnodes
  
                 ! Mesh1_O%Orientation = theta^S|_op + delta theta^S
              call SmllRotTrans( 'orientation', LinVec_1_theta(j  )+LinVec_1(j  ) &
                                              , LinVec_1_theta(j+1)+LinVec_1(j+1) &
                                              , LinVec_1_theta(j+2)+LinVec_1(j+2) & 
                                              , Mesh1_O%Orientation(:,:,i), ErrStat=ErrStat, ErrMsg=ErrMsg)            
                                      
              j = j+3            
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
           end do
  
           !M( theta^S|_op + delta theta^S )
           CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
  
           
           j=1
           do i=1,Mesh2_O%nnodes
              LinVec_2(j:j+2) = GetSmllRotAngs(Mesh2_I%Orientation(:,:,i),ErrStat,ErrMsg) 
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
              
              LinVec_2(j:j+2) = LinVec_2(j:j+2) - GetSmllRotAngs(Mesh2_I_op%Orientation(:,:,i),ErrStat,ErrMsg)
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
              
              j = j+3
           end do
           !call wrmatrix(LinVec_2-LinVec_2_theta,  99,Fmt,'LinVec_2')
           
           
           print *, 'Rotational Displacement: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )
              
        end do
        
        ! ..................
        ! Translational displacement:
        ! ..................
        
        call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
        do n1=1,15
           
           ! perturb x^S:
           call random_number( LinVec_1 )
           LinVec_1 = LinVec_1/real(n1**3)               
           call random_number( LinVec_1_a )
           LinVec_1_a = LinVec_1_a/real(n1**3)               
           
           LinVec_2_theta = matmul( map_mod1_mod2%dm%mi, LinVec_1 ) + matmul( map_mod1_mod2%dm%fx_p, LinVec_1_a ) ! approximate delta theta^D         
                        
           j=1
           do i=1,Mesh1_O%nnodes
  
                 ! Mesh1_O%TranslationDisp = theta^S|_op + delta theta^S
              Mesh1_O%TranslationDisp(:,i) = Mesh1_O_op%TranslationDisp(:,i) + LinVec_1(j:j+2)
              
                 ! Mesh1_O%Orientation = theta^S|_op + delta theta^S
              call SmllRotTrans( 'orientation', LinVec_1_theta(j  )+LinVec_1_a(j  ) &
                                              , LinVec_1_theta(j+1)+LinVec_1_a(j+1) &
                                              , LinVec_1_theta(j+2)+LinVec_1_a(j+2) & 
                                              , Mesh1_O%Orientation(:,:,i), ErrStat=ErrStat, ErrMsg=ErrMsg)            
                                      
              j = j+3            
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
                          
           end do
           
  
           !M( theta^S|_op + delta theta^S )
           CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
  
           ! delta theta^D = M( theta^S|_op + delta theta^S ) - x^D|_op
           j=1
           do i=1,Mesh2_O%nnodes
              LinVec_2(j:j+2) = Mesh2_I%TranslationDisp(:,i) - Mesh2_I_op%TranslationDisp(:,i)          
              j = j+3
           end do
                    
           print *, 'Translational Displacement: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )
                                
        end do      
        
        
        ! ..................
        ! Rotational velocity:
        ! ..................
        
        call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
        do n1=1,15
           
           ! perturb x^S:
           call random_number( LinVec_1 )
           LinVec_1 = LinVec_1/real(n1**3)               
           
           LinVec_2_theta = matmul( map_mod1_mod2%dm%mi, LinVec_1 )  ! approximate delta theta^D         
                        
           j=1
           do i=1,Mesh1_O%nnodes
  
                 ! Mesh1_O%RotationVel = theta^S|_op + delta theta^S
              Mesh1_O%RotationVel(:,i) = Mesh1_O_op%RotationVel(:,i) + LinVec_1(j:j+2)
                                                  
              j = j+3            
                          
           end do
           
  
           !M( theta^S|_op + delta theta^S )
           CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
  
           ! delta theta^D = M( theta^S|_op + delta theta^S ) - x^D|_op
           j=1
           do i=1,Mesh2_O%nnodes
              LinVec_2(j:j+2) = Mesh2_I%RotationVel(:,i) - Mesh2_I_op%RotationVel(:,i)          
              j = j+3
           end do
                    
           print *, 'Rotational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )                              
        end do            
        
        
        ! ..................
        ! Translational velocity:
        ! ..................
        
        call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
        do n1=1,15
           
           ! perturb x^S:
           call random_number( LinVec_1 )
           LinVec_1 = LinVec_1/real(n1**3)               
           call random_number( LinVec_1_a )
           LinVec_1_a = LinVec_1_a/real(n1**3)               
           call random_number( LinVec_1_b )
           LinVec_1_b = LinVec_1_b/real(n1**3)               
           
           LinVec_2_theta = matmul( map_mod1_mod2%dm%tv, LinVec_1 ) + matmul( map_mod1_mod2%dm%mi, LinVec_1_a ) &
                          + matmul( map_mod1_mod2%dm%fx_p, LinVec_1_b )! approximate delta theta^D         
                        
           j=1
           do i=1,Mesh1_O%nnodes
              
                 ! Mesh1_O%Orientation = theta^S|_op + delta theta^S
              call SmllRotTrans( 'orientation', LinVec_1_theta(j  )+LinVec_1(j  ) &
                                              , LinVec_1_theta(j+1)+LinVec_1(j+1) &
                                              , LinVec_1_theta(j+2)+LinVec_1(j+2) & 
                                              , Mesh1_O%Orientation(:,:,i), ErrStat=ErrStat, ErrMsg=ErrMsg)            
  
                          
                 ! Mesh1_O%TranslationVel = theta^S|_op + delta theta^S
              Mesh1_O%TranslationVel(:,i) = Mesh1_O_op%TranslationVel(:,i) + LinVec_1_a(j:j+2)
              
                 ! Mesh1_O%RotationVel = theta^S|_op + delta theta^S
              Mesh1_O%RotationVel(:,i) = Mesh1_O_op%RotationVel(:,i) + LinVec_1_b(j:j+2)
              
              j = j+3            
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
                          
           end do
           
  
           !M( theta^S|_op + delta theta^S )
           CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
  
           ! delta theta^D = M( theta^S|_op + delta theta^S ) - x^D|_op
           j=1
           do i=1,Mesh2_O%nnodes
              LinVec_2(j:j+2) = Mesh2_I%TranslationVel(:,i) - Mesh2_I_op%TranslationVel(:,i)          
              j = j+3
           end do
                    
           print *, 'Translational Velocity: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )
                                
        end do      
        
        ! ..................
        ! Rotational acceleration:
        ! ..................
        
        call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
        do n1=1,15
           
           ! perturb x^S:
           call random_number( LinVec_1 )
           LinVec_1 = LinVec_1/real(n1**3)               
           
           LinVec_2_theta = matmul( map_mod1_mod2%dm%mi, LinVec_1 )  ! approximate delta theta^D         
                        
           j=1
           do i=1,Mesh1_O%nnodes
  
                 ! Mesh1_O%RotationAcc = theta^S|_op + delta theta^S
              Mesh1_O%RotationAcc(:,i) = Mesh1_O_op%RotationAcc(:,i) + LinVec_1(j:j+2)
                                                  
              j = j+3            
                          
           end do
           
  
           !M( theta^S|_op + delta theta^S )
           CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
  
           ! delta theta^D = M( theta^S|_op + delta theta^S ) - x^D|_op
           j=1
           do i=1,Mesh2_O%nnodes
              LinVec_2(j:j+2) = Mesh2_I%RotationAcc(:,i) - Mesh2_I_op%RotationAcc(:,i)          
              j = j+3
           end do
                    
           print *, 'Rotational RotationAcc: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )                              
        end do            
                    
        ! ..................
        ! Translational acceleration:
        ! ..................
        
        call meshcopy( Mesh1_O_op, Mesh1_O, MESH_UPDATECOPY, ErrStat, ErrMsg) 
        do n1=1,15
           
           ! perturb x^S:
           call random_number( LinVec_1 )
           LinVec_1 = LinVec_1/real(n1**3)               
           call random_number( LinVec_1_a )
           LinVec_1_a = LinVec_1_a/real(n1**3)               
           call random_number( LinVec_1_b )
           LinVec_1_b = LinVec_1_b/real(n1**3)               
           call random_number( LinVec_1_c )
           LinVec_1_c = LinVec_1_c/real(n1**3)               
           
           LinVec_2_theta = matmul( map_mod1_mod2%dm%ta1, LinVec_1   ) + matmul( map_mod1_mod2%dm%ta2,  LinVec_1_a ) &
                          + matmul( map_mod1_mod2%dm%mi,  LinVec_1_b ) + matmul( map_mod1_mod2%dm%fx_p, LinVec_1_c )! approximate delta theta^D         
                        
           j=1
           do i=1,Mesh1_O%nnodes
              
                 ! Mesh1_O%Orientation = theta^S|_op + delta theta^S
              call SmllRotTrans( 'orientation', LinVec_1_theta(j  )+LinVec_1(j  ) &
                                              , LinVec_1_theta(j+1)+LinVec_1(j+1) &
                                              , LinVec_1_theta(j+2)+LinVec_1(j+2) & 
                                              , Mesh1_O%Orientation(:,:,i), ErrStat=ErrStat, ErrMsg=ErrMsg)            
  
                          
                 ! Mesh1_O%RotationVel = theta^S|_op + delta theta^S
              Mesh1_O%RotationVel(:,i) = Mesh1_O_op%RotationVel(:,i) + LinVec_1_a(j:j+2)
              
                 ! Mesh1_O%TranslationAcc = theta^S|_op + delta theta^S
              Mesh1_O%TranslationAcc(:,i) = Mesh1_O_op%TranslationAcc(:,i) + LinVec_1_b(j:j+2)
              
                 ! Mesh1_O%RotationAcc = theta^S|_op + delta theta^S
              Mesh1_O%RotationAcc(:,i) = Mesh1_O_op%RotationAcc(:,i) + LinVec_1_c(j:j+2)
              
              j = j+3            
              !IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))
                          
           end do
           
  
           !M( theta^S|_op + delta theta^S )
           CALL Transfer_Point_to_Point( Mesh1_O, Mesh2_I, Map_Mod1_Mod2, ErrStat, ErrMsg );                     IF (ErrStat /= ErrID_None) CALL WrScr("*******"//TRIM(ErrMsg))     
  
           ! delta theta^D = M( theta^S|_op + delta theta^S ) - x^D|_op
           j=1
           do i=1,Mesh2_O%nnodes
              LinVec_2(j:j+2) = Mesh2_I%TranslationAcc(:,i) - Mesh2_I_op%TranslationAcc(:,i)          
              j = j+3
           end do
                    
           print *, 'Translational Acceleration: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )
                                
        end do      
        
        
  end if
  
  
  else      
        
  end if ! linearization or mapping test
  
        ! ..............................................................................................................................   
        ! Destroy data structures:
        ! ..............................................................................................................................   
  