From 991eed13ab54fdba1f8912d81648d40da54df4e8 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Mon, 21 Dec 2009 15:14:22 +0000
Subject: [PATCH] Added read/initialization of zisodisp array
---
src/define_cloud.f90 | 12 +++++++++++-
1 file changed, 11 insertions(+), 1 deletion(-)
diff --git a/src/define_cloud.f90 b/src/define_cloud.f90
index a35540fc..67081bf6 100644
--- a/src/define_cloud.f90
+++ b/src/define_cloud.f90
@@ -16,7 +16,7 @@
!------------------------------------------------------------------------------|
!------------------------------------------------------------------------------|
-subroutine define_cloud (cl,params)
+subroutine define_cloud (cl,params,zisodisp)
!------------------------------------------------------------------------------|
!(((((((((((((((( Purpose of the routine ))))))))))))))))))))))))))))))))))))))
@@ -38,6 +38,7 @@ implicit none
type (cloud) cl
type (parameters) params
+double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
!------------------------------------------------------------------------------|
!(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
@@ -69,6 +70,9 @@ if (params%irestart.eq.0) then
allocate (cl%temp(1), stat=err) ; if (err.ne.0) call stop_run ('Error alloc cl%temp in define_cloud$')
allocate (cl%press(1), stat=err) ; if (err.ne.0) call stop_run ('Error alloc cl%press in define_cloud$')
allocate (cl%tag(1), stat=err) ; if (err.ne.0) call stop_run ('Error alloc cl%tag in define_cloud$')
+ if (params%isobc) then
+ zisodisp=0.d0
+ endif
else
open (19,file=trim(params%restartfile),status='old',form='unformatted')
@@ -142,6 +146,12 @@ else
cl%tag(i), &
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
+
close (19)
endif
--
GitLab