Modified to output the elemental dilatation rate

src/write_global_output.f90

@@ -95,6 +95,7 @@ if (iproc.eq.0) then
                osolve%nlsf,      &
                cl%np,            &
                current_time
+
    ! write info on octree solve nodes (x,y,z,u,v,w,lsf,temp)
    write (9) (osolve%x(i),                         &
                 osolve%y(i),                         &  
@@ -112,18 +113,21 @@ if (iproc.eq.0) then
                 osolve%kfix((i-1)*3+3),              &
                 osolve%kfixt(i),                     &
                 i=1,osolve%nnode)
+
    ! write icon array
    write (9) ((osolve%icon(k,i),k=1,8),   &
-! Line below uncommented by dwhipp - 12/09 
+                 ! Line below uncommented by dwhipp - 12/09 
                  osolve%pressure(i),        & 
-! Line below added by dwhipp - 12/09 
+                 ! Line below added by dwhipp - 12/09 
                  osolve%spressure(i),       &
                  osolve%crit(i),            &
                  osolve%e2d(i),             &
                  osolve%eviscosity(i),      & 
-                 osolve%is_plastic(i),      & 
+                 osolve%is_plastic(i),      &
+                 osolve%dilatr(i),          &
                  ov%whole_leaf_in_fluid(i), &
                  i=1,osolve%nleaves)
+
    ! write octree information
    write (9) (osolve%octree(i),i=1,osolve%octree(3))
    ! write bad face information