From fa4f187f7d02a1aff28522840aebaa0281fc27e6 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 17 Feb 2010 03:01:19 +0000
Subject: [PATCH] Added flag for zisodisp output
---
src/VTK/post.f90 | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/src/VTK/post.f90 b/src/VTK/post.f90
index fe731590..d92ae834 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'
--
GitLab