diff --git a/src/VTK/post.f90 b/src/VTK/post.f90
index fe731590c46f411937722a7b4c11a0470cf01176..d92ae834ccb31603fc60b370e1c86e090fc29562 100644
--- a/src/VTK/post.f90
+++ b/src/VTK/post.f90
@@ -78,6 +78,7 @@ integer :: output_ps
 integer :: output_cubes
 integer :: output_rivers
 integer :: output_regular
+integer :: output_zisodisp
 integer iproc,nproc,ierr,nnn,nnne,nnnep,levelmax,ix,iy,iz,nlinks,iordermin
 integer iconmin,iconmax
 integer,dimension(:),allocatable::nstrain,donor,order
@@ -188,6 +189,7 @@ read(77,*) output_ps
 read(77,*) output_cubes
 read(77,*) output_rivers
 read(77,*) output_regular
+read(77,*) output_zisodisp
 
 write(*,*) 'output u      field       ->',(output_u_field==1)
 write(*,*) 'output v      field       ->',(output_v_field==1)
@@ -213,6 +215,7 @@ write(*,*) 'output principal stresses ->',(output_ps==1)
 write(*,*) 'output elemental fields   ->',(output_cubes==1)
 write(*,*) 'output rivers             ->',(output_rivers==1)
 write(*,*) 'output regular            ->',(output_regular==1)
+write(*,*) 'output zisodisp           ->',(output_zisodisp==1)
 write(*,*) '--------------------------------------------------------------------------'
 
 !==============================================================================!
@@ -349,6 +352,8 @@ end if
 
 end do
 
+if (output_zisodisp==1) then
+
    read (7) (dumpdp,      & 
                 dumpdp,    &
                 dumpdp,    &
@@ -367,6 +372,7 @@ end do
    dxy=1.d0/real(nb)
    allocate(zisodisp(nb+1,nb+1),zisoslx(nb+1,nb+1))
    if (1==1) read (7) ((zisodisp(i,j),j=1,nb+1),i=1,nb+1)
+endif
 
 close(7)
 
@@ -1239,7 +1245,7 @@ endif
 !======[produce zisodisp.vtk file]=============================================
 !==============================================================================
 
-if (1==1) then
+if (output_zisodisp==1) then
   open(unit=30,file='zisodisp.vtk')
   write(30,'(a)')'# vtk DataFile Version 3.0'
   write(30,'(a)')'zisodisp'