Modified for output of dilatation rate

src/define_ov.f90

@@ -28,7 +28,6 @@ subroutine define_ov (ov,params,threadinfo)
 ! leveluniform_oct is the minimum/uniform octree level
 ! irestart is a flag to decide if this is a restart job or not
 ! restartfile is the name of the restart file if it is needed
- 
 
 !------------------------------------------------------------------------------|
 !((((((((((((((((  declaration of the subroutine arguments  ))))))))))))))))))))
@@ -49,7 +48,7 @@ type (thread) threadinfo
 
 integer iproc, nproc, err, ierr, levelv
 integer nface, nlsf, i, j, k, np, kfix
-double precision s,e2d,xlsf,crit,wiso,evisc
+double precision s,e2d,xlsf,crit,wiso,evisc,dilatr
 integer ioctree_number_of_elements
 external ioctree_number_of_elements
 double precision current_time,activation_time
@@ -67,7 +66,6 @@ call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 ov%noctree=params%noctreemax
 allocate (ov%octree(ov%noctree),stat=threadinfo%err) ;  call heap (threadinfo,'ov%octree','main',size(ov%octree),'int',+1)
 
-
 if (params%irestart.eq.0) then
    ! creates the octree
    call octree_init (ov%octree,ov%noctree)
@@ -169,6 +167,7 @@ else
               e2d,                      &
               evisc,                    & 
               is_plas,                  &
+              dilatr,                   &
               ov%whole_leaf_in_fluid(i),&
               i=1,ov%nleaves)
    read (9) (ov%octree(i),i=1,ov%octree(3))