From 68af944734c7d3b64e15cf5a6518838dc3795522 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 3 Feb 2010 03:31:47 +0000
Subject: [PATCH] Added output of zisodisp surface
---
src/VTK/post.f90 | 20 +++++++++++---------
1 file changed, 11 insertions(+), 9 deletions(-)
diff --git a/src/VTK/post.f90 b/src/VTK/post.f90
index 38e3a31c..d84fa123 100644
--- a/src/VTK/post.f90
+++ b/src/VTK/post.f90
@@ -363,8 +363,10 @@ end do
i=1,npcl)
! read isostasy basal displacement array - dwhipp 11/09
- dxy=abs(ov%x(1)-ov%x(2))
- nb=1.d0/dxy
+ !dxy=abs(ov%x(1)-ov%x(2))
+ dxy=1.d0/128.d0
+ !nb=int(1.d0/dxy)
+ nb=128
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)
@@ -1245,28 +1247,28 @@ if (1==1) then
write(30,'(a)')'zisodisp'
write(30,'(a)')'ASCII'
write(30,'(a)')'DATASET UNSTRUCTURED_GRID'
- write(30,'(a7,i10,a6)')'POINTS ',(nb+1)**2.d0,' float'
+ write(30,'(a7,i10,a6)')'POINTS ',int((nb+1)**2.d0),' float'
do j=1,nb+1
do i=1,nb+1
- write(30,'(3f16.11)') ((i-1)*dxy)/nb,((j-1)*dxy)/nb,zisodisp(i,j)
+ write(30,'(3f16.11)') (i-1)*dxy,(j-1)*dxy,zisodisp(i,j)
enddo
enddo
- write(30,'(a6, 2I10)') 'CELLS ',nb**2.d0,(4+1)*nb**2.d0
+ write(30,'(a6, 2I10)') 'CELLS ',int(nb**2.d0),int((4+1)*nb**2.d0)
ptcnt=0
do j=1,nb
do i=1,nb
+ write(30,'(9I10)')4,ptcnt,ptcnt+1,ptcnt+1+(nb+1),ptcnt+(nb+1)
ptcnt=ptcnt+1
- write(30,'(9I10)')4,ptcnt,ptcnt+1,ptcnt+1+nb,ptcnt+nb
enddo
ptcnt=ptcnt+1
enddo
- write(30,'(A11, I10)') 'CELL_TYPES ',nb**2.d0
+ write(30,'(A11, I10)') 'CELL_TYPES ',int(nb**2.d0)
do j=1,nb
do i=1,nb
write(30,'(I2)')9 ! VTK quad (4 nodes)
enddo
enddo
- write(30,'(a11,i10)') 'POINT_DATA ',(nb+1)**2.d0
+ write(30,'(a11,i10)') 'POINT_DATA ',int((nb+1)**2.d0)
write(30,'(a)')'SCALARS zisodisp float 1'
write(30,'(a)')'LOOKUP_TABLE default'
do j=1,nb+1
@@ -1274,7 +1276,7 @@ if (1==1) then
write(30,'(e11.4)') zisodisp(i,j)
enddo
enddo
- write(30,'(a11,i10)') 'POINT_DATA ',(nb+1)**2.d0
+ !write(30,'(a11,i10)') 'POINT_DATA ',int((nb+1)**2.d0)
write(30,'(a)')'SCALARS zisoslx float 1'
write(30,'(a)')'LOOKUP_TABLE default'
do j=1,nb+1
--
GitLab