From 78cc96fa36006fad87d4c9a80c1b271c2bee7b97 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Mon, 21 Dec 2009 16:17:38 +0000
Subject: [PATCH] updated for output of elemental pressure, viscosity and
is_plastic
---
src/define_cloud.f90 | 11 +++++------
1 file changed, 5 insertions(+), 6 deletions(-)
diff --git a/src/define_cloud.f90 b/src/define_cloud.f90
index 8d550ad0..ebd06389 100644
--- a/src/define_cloud.f90
+++ b/src/define_cloud.f90
@@ -43,7 +43,7 @@ double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
!------------------------------------------------------------------------------|
!(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
!------------------------------------------------------------------------------|
-double precision x,y,z,u,v,w,p,r,s,t,e2d,xlsf,xn,yn,zn,don,con
+double precision x,y,z,u,v,w,p,r,s,t,e2d,xlsf,xn,yn,zn,don,con,wp,ev,ip,epr,espr
double precision current_time,activation_time
logical inf,wlif
integer err,iproc,nproc,ierr,i,j,k,nl,nf,nn,nr,kx,ky,kz,kt,icon
@@ -100,6 +100,7 @@ else
u, &
v, &
w, &
+ wp, &
(xlsf,j=1,nlsf), &
t, &
p, &
@@ -109,7 +110,7 @@ else
kz, &
kt, &
i=1,nnode)
- read (19) ((icon,k=1,8),con,e2d,wlif,i=1,nleaves)
+ read (19) ((icon,k=1,8),epr,espr,con,e2d,ev,ip,wlif,i=1,nleaves)
read (19) (ioc,k=1,noctree)
read (19) ((iface,k=1,9),i=1,nface)
read (19) (nn,nl,nf,nr,inf,i=1,nnode)
@@ -145,10 +146,8 @@ else
i=1,cl%np)
! read isostasy basal displacement array - dwhipp 11/09
- if (params%isobc) then
- read (19) ((zisodisp(i,j),j=1,2**params%levelmax_oct+1),i=1,&
- 2**params%levelmax_oct+1)
- endif
+ if (params%isobc) read (19) ((zisodisp(i,j),j=1,2**params%levelmax_oct+1),&
+ i=1,2**params%levelmax_oct+1)
close (19)
endif
--
GitLab