From 134bfc6c4b363d4b1ccc21aec60d53e1731b43bf Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Fri, 16 Apr 2010 03:13:34 +0000
Subject: [PATCH] Fixed zisodisp tracking to work only when isobc is true
---
src/DOUAR.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/DOUAR.f90 b/src/DOUAR.f90
index 31878ebd..ea426e21 100644
--- a/src/DOUAR.f90
+++ b/src/DOUAR.f90
@@ -1233,7 +1233,7 @@ do while (istep.le.params%nstep)
call show_time (total,step,inc,1,'Compute isostasy and adjust vertical velocity$')
!allocate(ov%wpreiso(ov%nnode),stat=err) ; if (err.ne.0) call stop_run ('Error alloc ov%wpreiso in main$')
call isostasy (params,weightel,ov,surface,mat,0,zi)
- zi%zisodisp=zi%zisodisp-zi%zisoinc
+ if (params%isobc) zi%zisodisp=zi%zisodisp-zi%zisoinc
do is=1,params%ns
if (current_time+tiny(current_time).ge. surface(is)%activation_time) then
--
GitLab