From 1fbea68bb105d9908505e890c5b83dff60d1ed6c Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 14 Apr 2010 20:25:59 +0000
Subject: [PATCH] Edited post processor to use new isostasy setup and output
isostasy velocities as well
---
src/VTK/input_of_outputs.txt | 1 +
src/VTK/post.f90 | 17 ++++++++++++-----
2 files changed, 13 insertions(+), 5 deletions(-)
diff --git a/src/VTK/input_of_outputs.txt b/src/VTK/input_of_outputs.txt
index 11c01455..1bd1241f 100644
--- a/src/VTK/input_of_outputs.txt
+++ b/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
diff --git a/src/VTK/post.f90 b/src/VTK/post.f90
index d92ae834..4022fda1 100644
--- a/src/VTK/post.f90
+++ b/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
--
GitLab