Updated for output of e2dp on cloud

src/VTK/post.f90

@@ -23,7 +23,8 @@ MODULE definitions
    end type sheet
 
    type cloud
-      double precision,dimension(:),pointer::x,y,z,x0,y0,z0,strain,lsf0,temp,press
+      double precision,dimension(:),pointer :: x,y,z,x0,y0,z0,strain,lsf0,temp
+      double precision,dimension(:),pointer :: press,e2dp
       integer,dimension(:),pointer::tag
    end type cloud
 
@@ -343,7 +344,7 @@ allocate(ov%icon(8,ov%nleaves))
 allocate(ov%octree(ov%noctree))
 allocate(cl%x(npcl),cl%y(npcl),cl%z(npcl),cl%x0(npcl),cl%y0(npcl),cl%z0(npcl))
 allocate(cl%strain(npcl),cl%lsf0(npcl),cl%temp(npcl),cl%press(npcl))
-allocate(cl%tag(npcl))
+allocate(cl%e2dp(npcl),cl%tag(npcl))
 
 !========================
 !=====[nodal values]=====
@@ -487,6 +488,7 @@ read (7) (cl%x(i),      &
           cl%lsf0(i),   &
           cl%temp(i),   &
           cl%press(i),  &
+          cl%e2dp(i),   &
           cl%tag(i),    &
           i=1,npcl)
 
@@ -1669,6 +1671,11 @@ if (output_cloud==1) then
   do i=1,npcl
     if (cl%lsf0(i).le.0.d0) write(iunit,'(f16.11)') cl%press(i)
   enddo
+  write(30,'(a)')'SCALARS e2dp float 1'
+  write(30,'(a)')'LOOKUP_TABLE default'
+  do i=1,npcl
+    if (cl%lsf0(i).le.0.d0) write(iunit,'(f16.11)') cl%e2dp(i)
+  enddo
   write(30,'(a)')'SCALARS tag integer 1'
   write(30,'(a)')'LOOKUP_TABLE default'
   do i=1,npcl