Edited post processor to use new isostasy setup and output isostasy velocities as well

src/VTK/input_of_outputs.txt

@@ -3,6 +3,7 @@
 0  ! output w scalar field 
 1  ! output velocity vectors 
 1  ! output pre-isostasy velocity vectors
+1  ! output isostasy (only) velocity vectors
 0  ! output disp=sqrt(u**2+v**2+w**2) field
 1  ! output nodal pressure field 
 1  ! output raw elemental pressure field 

src/VTK/post.f90

@@ -6,7 +6,7 @@ MODULE definitions
       integer noctree,nnode,nleaves
       integer nelemr,nlsf 
       double precision,dimension(:),pointer::x,y,z
-      double precision,dimension(:),pointer::u,v,w,wpreiso,temp,pressure,strain
+      double precision,dimension(:),pointer::u,v,w,wiso,temp,pressure,strain
       double precision,dimension(:),pointer::nodal_pressure,spressure
       double precision,dimension(:),pointer::e2d,eviscosity
       double precision,dimension(:),pointer::crit
@@ -225,7 +225,7 @@ read (7) ov%noctree,ov%nnode,ov%nleaves,ner,ov%nlsf,npcl,current_time
 nn=ov%nnode
 allocate(ov%on(nn))
 allocate(ov%x(nn),ov%y(nn),ov%z(nn))
-allocate(ov%u(nn),ov%v(nn),ov%w(nn),ov%wpreiso(nn))
+allocate(ov%u(nn),ov%v(nn),ov%w(nn),ov%wiso(nn))
 allocate(ov%pressure(ov%nleaves),ov%strain(nn))
 allocate(ov%e2d(ov%nleaves),ov%eviscosity(ov%nleaves),ov%is_plastic(ov%nleaves))
 allocate(ov%whole_leaf_in_fluid(ov%nleaves))
@@ -246,7 +246,7 @@ allocate(ov%octree(ov%noctree))
 
 read(7)(ov%x(i),ov%y(i),ov%z(i),   &
         ov%u(i),ov%v(i),ov%w(i),   &
-        ov%wpreiso(i),             &
+        ov%wiso(i),             &
         ov%lsf(i,1:ov%nlsf),       &
         ov%temp(i),                &
         ov%nodal_pressure(i),      &
@@ -422,7 +422,7 @@ write(*,'(a,2f30.20)') 'z                 range : ',minval(ov%z(1:nn)),      max
 write(*,'(a,2f30.20)') 'u                 range : ',minval(ov%u(1:nn)),      maxval(ov%u(1:nn))
 write(*,'(a,2f30.20)') 'v                 range : ',minval(ov%v(1:nn)),      maxval(ov%v(1:nn))
 write(*,'(a,2f30.20)') 'w                 range : ',minval(ov%w(1:nn)),      maxval(ov%w(1:nn))
-write(*,'(a,2f30.20)') 'w (preiso)        range : ',minval(ov%wpreiso(1:nn)),maxval(ov%wpreiso(1:nn))
+write(*,'(a,2f30.20)') 'w (preiso)        range : ',minval(ov%wiso(1:nn)),maxval(ov%wiso(1:nn))
 ! Line below added by dwhipp - 12/09 
 write(*,'(a,2f30.20)') 'pressure          range : ',minval(ov%pressure),     maxval(ov%pressure) 
 ! Line below added by dwhipp - 12/09 
@@ -881,7 +881,14 @@ endif
 if(output_preiso_velo_vect==1) then 
    write(iunit,'(a)')'VECTORS preiso-velo float' 
    do i=1,nnn 
-      write(iunit,'(3e15.4)')ov%u(rtf(i)),ov%v(rtf(i)),ov%wpreiso(rtf(i)) 
+      write(iunit,'(3e15.4)')ov%u(rtf(i)),ov%v(rtf(i)),ov%w(rtf(i))-ov%wiso(rtf(i)) 
+   enddo 
+endif
+
+if(output_iso_velo_vect==1) then 
+   write(iunit,'(a)')'VECTORS iso-velo float' 
+   do i=1,nnn 
+      write(iunit,'(3e15.4)')0.d0*ov%u(rtf(i)),0.d0*ov%v(rtf(i)),ov%wiso(rtf(i)) 
    enddo 
 endif