anchor.diff:      write (VTKFile,*) HED
       write (VTKFile,"(A5)") 'ASCII'
       write (VTKFile,"(A25)") 'DATASET UNSTRUCTURED_GRID'
  +    write (VTKFile,*) "POINTS",NUMNP,"double"
  +    do i = 1, NUMNP
  +        write (VTKFile,*) DA(NP(2)+i-1), DA(NP(3)+i-1), DA(NP(4)+i-1)           !X(i), Y(i), Z(i)
  +    end do
       
   case (2)                                                            !Called in elcal.f90, subroutine ELCAL
       select case (NPAR(1))
  
positive.code:     rewind (VTKTmpFile)
      do i = 1 , NLCASE
          read (VTKTmpFile) string(1:25), Dat(1:2)                    !Fetch Displacements of Load Cases
          write (VTKFile,*) string(1:25), Dat(1:2), "double"
          do j = 1, NUMNP
              read (VTKTmpFile) Dat1(1:Dat(1))
              write (VTKFile,*) Dat1(1:Dat(1))
  
negative.code:      end select
       
  case (3)                                                            !Called in stap.f90, STAP at solution phase IND=3
  !>>>>>>>>>>>>>>>>>>>>>>>>>>HOW TO RESOLVE CONNECTION POINTS BETWEEN TWO ELEMENT GROUPS??? THESE POINTS ARE REPEATED
      write (VTKFile,*) "CELLS ", NEL, NCONECT                        !Sum up all elements to generate a global picture
      rewind (VTKNodeTmp)
      do i = 1 , NEL
  