From d31310ba737615b6f4671e6a20edea293d61e784 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 3 Feb 2010 03:32:23 +0000
Subject: [PATCH] modified to use new zi (zisodisp) structure

---
 src/isostasy.f90 | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/src/isostasy.f90 b/src/isostasy.f90
index 31c69c75..17989b3f 100644
--- a/src/isostasy.f90
+++ b/src/isostasy.f90
@@ -1,4 +1,4 @@
-subroutine isostasy (params,weightel,ov,surface,mat,zisodisp)
+subroutine isostasy (params,weightel,ov,surface,mat,zi)
 
 ! first go at imposing isostasy
 ! this routine assumes that the original surface are flat
@@ -15,7 +15,8 @@ type(octreev) ov
 double precision weightel(ov%nleaves)
 type(sheet) surface(params%ns)
 type(material) mat(0:params%nmat)
-double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
+type(ziso) zi
+!double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
 
 double precision xc,yc,zc,eps
 integer i,j,ii,jj,kk,ie,dlev,is,iproc,nproc,ierr,ic
@@ -261,7 +262,7 @@ enddo
   endif
 
 ! Store the total displacement of the basal node plane 
-if (params%isobc) zisodisp=zisodisp-work
+if (params%isobc) zi%zisodisp=zi%zisodisp-work
 
 ! transforms the displacements into velocities and add this contribution to the velocity
 ! solution of the Stokes equation before interpolation onto the surfaces and their advection
-- 
GitLab