anchor.diff:  
         call MeshMapDestroy(Map_Mod1_Mod2, ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
         call MeshMapDestroy(Map_Mod2_Mod1, ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      
  -
  -      
  -      CALL MeshDestroy( mesh_Motion_1PT, ErrStat, ErrMsg ); IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      CALL MeshDestroy( mesh1_I_1PT, ErrStat, ErrMsg );     IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      CALL MeshDestroy( mesh2_O_1PT, ErrStat, ErrMsg );     IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -   
  -      call MeshMapDestroy(Map_Mod1_I_1PT, ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      call MeshMapDestroy(Map_Mod2_O_1PT, ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -   
  -      
  -#ifdef MESH_DEBUG  
  -      CALL MeshDestroy( mesh2_O_1PT_augmented,      ErrStat, ErrMsg );   IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      CALL MeshDestroy( mesh2_O_1PT_lumped,         ErrStat, ErrMsg );   IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      CALL MeshDestroy( mesh1_I_1PT_lumped,         ErrStat, ErrMsg );   IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -
  -      call MeshMapDestroy(Map_Mod2_O_1PT_augmented, ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      call MeshMapDestroy(Map_Mod2_O_1PT_lumped,    ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -      call MeshMapDestroy(Map_Mod1_I_1PT_lumped,    ErrStat, ErrMsg);IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -#endif      
  -      
  -      ! ..............................................................................................................................   
  -      ! open files:
  -      ! ..............................................................................................................................   
  -      close( un_out )
  -
  -   
  +                    
  +      if (allocated(LinVec_1)) deallocate(LinVec_1)
  +      if (allocated(LinVec_2)) deallocate(LinVec_2)      
  +      if (allocated(LinVec_1_theta)) deallocate(LinVec_1_theta)
  +      if (allocated(LinVec_2_theta)) deallocate(LinVec_2_theta)      
  +      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)      
  +      
      end do
      
      
      
      
  -contains
  -   ! ..............................................  
  -   subroutine CreateTotalLoadsPointMeshes( )
  -   
  -      CALL MeshCreate( BlankMesh       = mesh1_I_1PT        &
  -                     , IOS              = COMPONENT_INPUT   &
  -                     , NNodes           = 1                 &
  -                     , Force            = .TRUE.            &
  -                     , Moment           = .TRUE.            &
  -                     , ErrStat          = ErrStat           &
  -                     , ErrMess          = ErrMsg            )
  -      
  -      IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -         
  -         
  -      position = (/0.0_ReKi, 0.0_ReKi, 0.0_ReKi/)
  -      CALL MeshPositionNode ( mesh1_I_1PT, 1, position, ErrStat, ErrMsg ) ; IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))             
  -
  -      CALL MeshConstructElement ( mesh1_I_1PT, ELEMENT_POINT, ErrStat, ErrMsg, 1 );                 IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))                     
  -                  
  -      CALL MeshCommit ( mesh1_I_1PT, ErrStat, ErrMsg )   ;      IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg)) 
  -      
  -      
  -      CALL MeshCopy( mesh1_I_1PT, mesh_Motion_1PT, MESH_SIBLING, ErrStat, ErrMsg &
  -                     , IOS              = COMPONENT_OUTPUT  &
  -                     , TranslationDisp  = .TRUE.            ) ;      IF (ErrStat /= ErrID_None) CALL WrScr(TRIM(ErrMsg))
  -                     
  -      !.....                 
  -      CALL MeshCopy( mesh1_I_1PT, mesh2_O_1PT, MESH_NEWCOPY, ErrStat, ErrMsg )  ! This thinks it's for input, but really it's for output. I don't think it matters...
  -
  -      
  -      
  -   end subroutine CreateTotalLoadsPointMeshes
  -   
  +contains   
      ! ..............................................   
      subroutine CreateOutputMeshes_Test1()   
         ! this is a point-to-point mapping, with one point going to many.
  
positive.code:    subroutine CreateOutputMeshes_Test5Orient
     
        real(reki) :: z
        real(reki) :: DCM(3,3)
        
        
        Mesh1Type = ELEMENT_LINE2
  
negative.code:    
        WRITE(BinOutputName,'(A,A,A)') 'Test', TRIM(Num2LStr(TestNumber)),'Meshes.bin'
        
     
        ! ..............................................................................................................................   
        ! Create sibling input meshes: 
  