anchor.diff:                    angle = TwoPi*(k-1.0_ReKi)/NumSegments1
                     xyz(1) = r*COS(angle)
                     xyz(2) = r*SIN(angle)
  -                  WRITE(Un,VTK_AryFmt) M%Position(:,i) + M%TranslationDisp(:,i) + MATMUL(M%Orientation(:,:,i),xyz)
  +                  WRITE(Un,VTK_AryFmt) M%Position(:,i) + M%TranslationDisp(:,i) + MATMUL(xyz,M%Orientation(:,:,i))
                     !WRITE(Un,VTK_AryFmt) M%Position(:,i) + M%TranslationDisp(:,i) + xyz
                  END DO            
               END DO
  
positive.code:       WRITE(Un,'(A)')         '      <Points>'
        WRITE(Un,'(A)')         '        <DataArray type="Float32" NumberOfComponents="3" format="ascii">'
        
        if (present(verts)) then
           DO i=1,M%Nnodes
              DO j=1,NumSegments1
                 delX = verts(1,j,i)
                 delY = verts(2,j,i)
                 WRITE(Un,VTK_AryFmt) M%Position(:,i) + M%TranslationDisp(:,i) + delX*M%Orientation(1,:,i) + delY*M%Orientation(2,:,i)
              END DO
           END DO         
        else               
           DO i=1,M%Nnodes
              DO j=1,NumSegments1
                 angle = TwoPi*(j-1.0_ReKi)/NumSegments1
                 delX = radius(i)*COS(angle)
                 delY = radius(i)*SIN(angle)
                 WRITE(Un,VTK_AryFmt) M%Position(:,i) + M%TranslationDisp(:,i) + delX*M%Orientation(1,:,i) + delY*M%Orientation(2,:,i)
              END DO
           END DO
        end if
  
negative.code:    INTEGER(IntKi)                        :: i,j           ! loop counters
     CHARACTER(1024)                       :: FileName
     REAL(SiKi)                            :: angle
     REAL(SiKi)                            :: xyz(3)
  
     INTEGER(IntKi)                        :: firstPntEnd, firstPntStart, secondPntStart, secondPntEnd  ! node indices for forming rectangle 
     INTEGER(IntKi)                        :: NumSegments1
  