Added dummy read of number of elements along edge of octree

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