Modified to read dilatation rate

src/define_surface.f90

@@ -47,9 +47,10 @@ double precision total, step, inc, current_time
      
 integer err,ierr,iproc,nproc,iface,ioc,i,j,k,is,noctree,nnode,nleaves,nface,nt
 integer icon,inn,inl,nlsf,ncl,nf,inr,kx,ky,kz,kt
-logical inf,wlif
+logical inf,wlif,is_plas
 character*7 ic
-double precision s,p,t,x,y,z,u,v,w,wiso,xlsf,e2d,don,con,activation_time
+double precision s,p,t,x,y,z,u,v,w,wiso,xlsf,e2d,don,con,activation_time,dilatr
+double precision epr,espr,evisc
 
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
@@ -153,7 +154,7 @@ else
                kz,             &
                kt,             &
                i=1,nnode)
-   read (9) ((icon,k=1,8),con,e2d,wlif,i=1,nleaves)
+   read (9) ((icon,k=1,8),epr,espr,con,e2d,evisc,is_plas,dilatr,wlif,i=1,nleaves)
    read (9) (ioc,k=1,noctree)                    
    read (9) ((iface,k=1,9),i=1,nface)
    read (9) (inn, &