anchor.diff:        write(Label,'(A,A)') 'NW.Bld', i2ABC(iW)
         Filename = TRIM(FileRootName)//'.'//trim(Label)//'.'//Tstr//'.vtk'
         if (m%FirstCall) then ! Small Hack - At t=0, NW not set, but first NW panel is the LL panel
  -         call WrVTK_Lattice(FileName, m%r_LL(1:3,:,1:2,iW), m%Gamma_LL(:,iW:iW))
  +         allocate(dxdt_0(3, size(m%dxdt_NW,2) , m%nNW+1)); dxdt_0=0.0_ReKi
  +         call WrVTK_Lattice(FileName, m%r_LL(1:3,:,1:2,iW), m%Gamma_LL(:,iW:iW),dxdt_0)
  +         deallocate(dxdt_0)
         else
  -         call WrVTK_Lattice(FileName, x%r_NW(1:3,:,1:m%nNW+1,iW), x%Gamma_NW(:,1:m%nNW,iW))
  +         call WrVTK_Lattice(FileName, x%r_NW(1:3,:,1:m%nNW+1,iW), x%Gamma_NW(:,1:m%nNW,iW), m%dxdt_NW(:,:,1:m%nNW+1,iW))
         endif
      enddo
      ! --------------------------------------------------------------------------------}
  
positive.code:       call vtk_cell_data_scalar(SegEpsilon,'Epsilon')
  !       call vtk_cell_data_scalar(real(SegConnct(3,:), ReKi),'Age')
        !call vtk_cell_data_scalar(real(SegConnct(4,:), ReKi),'Span')
        !call vtk_point_data_init()
        !call vtk_point_data_vector(Sgmt%UconvP(1:3,1:Sgmt%nP_Storage),'Uconv')
        call vtk_close_file()
     endif
  end subroutine
  
  subroutine WrVTK_Lattice(filename, LatticePoints, LatticeGamma)
     use VTK ! for all the vtk_* functions
     character(len=*), intent(in)                         :: filename
     real(Reki), dimension(:,:,:), intent(in  )           :: LatticePoints !< Array of points 3 x nSpan x nDepth
     real(Reki), dimension(:,:), intent(in  )             :: LatticeGamma  !< Array of            nSpan x nDepth
     !real(Reki), dimension(:,:,:), intent(in  ), optional :: LatticeData3d !< Array of n x nSpan x nDepth KEEP ME
     !
     integer(IntKi), dimension(:,:), allocatable :: Connectivity
     real(ReKi), dimension(:,:), allocatable     :: Points
  
negative.code:    logical              :: bMirror
     integer(IntKi)       :: ErrStat2
     character(ErrMsgLen) :: ErrMsg2
     real(Reki), dimension(:,:,:), allocatable :: dxdt_0 !<
  
     if (DEV_VERSION) then
        print*,'------------------------------------------------------------------------------'
  