updated for output of elemental pressure, viscosity and is_plastic

src/define_cloud.f90

@@ -43,7 +43,7 @@ double precision zisodisp(2**params%levelmax_oct+1,2**params%levelmax_oct+1)
 !------------------------------------------------------------------------------|
 !(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
 !------------------------------------------------------------------------------|
-double precision x,y,z,u,v,w,p,r,s,t,e2d,xlsf,xn,yn,zn,don,con
+double precision x,y,z,u,v,w,p,r,s,t,e2d,xlsf,xn,yn,zn,don,con,wp,ev,ip,epr,espr
 double precision current_time,activation_time
 logical inf,wlif
 integer err,iproc,nproc,ierr,i,j,k,nl,nf,nn,nr,kx,ky,kz,kt,icon
@@ -100,6 +100,7 @@ else
                 u,               &
                 v,               &
                 w,               &
+                wp,              &
                 (xlsf,j=1,nlsf), &
                 t,               &
                 p,               &
@@ -109,7 +110,7 @@ else
                 kz,              &
                 kt,              &
                 i=1,nnode)
-   read (19) ((icon,k=1,8),con,e2d,wlif,i=1,nleaves)
+   read (19) ((icon,k=1,8),epr,espr,con,e2d,ev,ip,wlif,i=1,nleaves)
    read (19) (ioc,k=1,noctree)                    
    read (19) ((iface,k=1,9),i=1,nface)
    read (19) (nn,nl,nf,nr,inf,i=1,nnode)
@@ -145,10 +146,8 @@ else
                 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
+   if (params%isobc) read (19) ((zisodisp(i,j),j=1,2**params%levelmax_oct+1),&
+                               i=1,2**params%levelmax_oct+1)
 
    close (19)
 endif