From 12b3b1c5b5e99a77e77744ec3e44fb454e8d6d32 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 10 Mar 2010 03:16:11 +0000
Subject: [PATCH] Added dummy read of number of elements along edge of octree

---
 src/define_cloud.f90 | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

diff --git a/src/define_cloud.f90 b/src/define_cloud.f90
index a541ab84..7d8dadd7 100644
--- a/src/define_cloud.f90
+++ b/src/define_cloud.f90
@@ -48,7 +48,7 @@ double precision x,y,z,u,v,w,p,r,s,t,e2d,xlsf,xn,yn,zn,don,con,wp,ev,epr,espr
 double precision current_time,activation_time
 logical inf,wlif,ip
 integer err,iproc,nproc,ierr,i,j,k,nl,nf,nn,nr,kx,ky,kz,kt,icon
-integer ioc,iface,ilsf,noctree,nnode,nface,nlsf,nleaves,ns,nt
+integer ioc,iface,ilsf,noctree,nnode,nface,nlsf,nleaves,ns,nt,dumpi
 
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
@@ -151,8 +151,11 @@ else
                 i=1,cl%np)
 
    ! read isostasy basal displacement array - dwhipp 11/09
-   if (params%isobc) read (19) ((zi%zisodisp(i,j),j=1,2**params%levelmax_oct+1),&
-                               i=1,2**params%levelmax_oct+1)
+   if (params%isobc) then
+     read (19) dumpi
+     read (19) ((zi%zisodisp(i,j),j=1,2**params%levelmax_oct+1),&
+                                  i=1,2**params%levelmax_oct+1)
+   endif
 
 ! correct for vertical exaggeration
 
-- 
GitLab