Added output of zisodisp surface

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