diff --git a/src/define_cloud.f90 b/src/define_cloud.f90
index a541ab845e3d8d271a172c59b22c87a4c129b145..7d8dadd705a2dc189565dc68537112822bb36839 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