Updated output for new values (eviscosity, is_plastic, etc.)

src/VTK/post.f90

@@ -8,11 +8,11 @@ MODULE definitions
       double precision,dimension(:),pointer::x,y,z
       double precision,dimension(:),pointer::u,v,w,wpreiso,temp,pressure,strain
       double precision,dimension(:),pointer::nodal_pressure,spressure
-      double precision,dimension(:),pointer::e2d
+      double precision,dimension(:),pointer::e2d,eviscosity
       double precision,dimension(:),pointer::crit
       integer,dimension(:,:),pointer::icon,iconr
       double precision,dimension(:,:),pointer:: lsf
-      integer,dimension(:),pointer::on(:)
+      integer,dimension(:),pointer::on(:),is_plastic
       logical, dimension(:), pointer :: whole_leaf_in_fluid
    end type structure
 
@@ -63,9 +63,12 @@ integer :: output_disp_field
 integer :: output_press_fieldn 
 integer :: output_press_fielde 
 integer :: output_smooth_press_fielde
+integer :: output_countnode_field
 integer :: output_temp_field
 integer :: output_e2d_fielde 
-integer :: output_e2d_fieldn 
+integer :: output_e2d_fieldn
+integer :: output_eviscosity_fielde 
+integer :: output_is_plastic_fielde
 integer :: output_crit_field
 integer :: output_strain_field
 integer :: output_lsf_field
@@ -170,9 +173,12 @@ read(77,*) output_disp_field
 read(77,*) output_press_fieldn 
 read(77,*) output_press_fielde 
 read(77,*) output_smooth_press_fielde
+read(77,*) output_countnode_field
 read(77,*) output_temp_field 
 read(77,*) output_e2d_fielde 
 read(77,*) output_e2d_fieldn
+read(77,*) output_eviscosity_fielde 
+read(77,*) output_is_plastic_fielde
 read(77,*) output_crit_field 
 read(77,*) output_strain_field 
 read(77,*) output_lsf_field 
@@ -192,9 +198,12 @@ write(*,*) 'output disp   field       ->',(output_disp_field==1)
 write(*,*) 'output node press field   ->',(output_press_fieldn==1) 
 write(*,*) 'output raw press field    ->',(output_press_fielde==1) 
 write(*,*) 'output smooth press field ->',(output_smooth_press_fielde==1)
+write(*,*) 'output countnode field    ->',(output_countnode_field==1)
 write(*,*) 'output temp   field       ->',(output_temp_field==1)
 write(*,*) 'output elemental e2d field->',(output_e2d_fielde==1) 
 write(*,*) 'output nodal e2d field    ->',(output_e2d_fieldn==1)
+write(*,*) 'output elem eff viscosity ->',(output_eviscosity_fielde==1) 
+write(*,*) 'output element is plastic ->',(output_is_plastic_fielde==1)
 write(*,*) 'output crit   field       ->',(output_crit_field==1)
 write(*,*) 'output strain field       ->',(output_strain_field==1)
 write(*,*) 'output lsf    field       ->',(output_lsf_field==1)
@@ -215,7 +224,7 @@ allocate(ov%on(nn))
 allocate(ov%x(nn),ov%y(nn),ov%z(nn))
 allocate(ov%u(nn),ov%v(nn),ov%w(nn),ov%wpreiso(nn))
 allocate(ov%pressure(ov%nleaves),ov%strain(nn))
-allocate(ov%e2d(ov%nleaves))
+allocate(ov%e2d(ov%nleaves),ov%eviscosity(ov%nleaves),ov%is_plastic(ov%nleaves))
 allocate(ov%whole_leaf_in_fluid(ov%nleaves))
 allocate(ov%crit(ov%nleaves))
 allocate(ov%temp(nn))
@@ -252,6 +261,8 @@ read(7) (ov%icon(1:8,ie),            &
          ov%spressure(ie),           &
          ov%crit(ie),                &
          ov%e2d(ie),                 &
+         ov%eviscosity(ie),          & 
+         ov%is_plastic(ie),          & 
          ov%whole_leaf_in_fluid(ie), & 
          ie=1,ov%nleaves)
 
@@ -371,8 +382,6 @@ do i=1,ov%nnode
    end if
 end do
 
-deallocate(countnode)
-
 !==============================================================================
 !==============================================================================
 
@@ -785,6 +794,14 @@ if(output_press_fieldn==1) then
    enddo
 end if
 
+if(output_countnode_field==1) then 
+   write(iunit,'(a)')'SCALARS countnode float 1' 
+   write(iunit,'(a)')'LOOKUP_TABLE default' 
+   do i=1,nnn 
+      write(iunit,'(e11.4)') countnode(rtf(i)) 
+   enddo 
+end if 
+ 
 if(output_e2d_fieldn==1) then 
    write(iunit,'(a)')'SCALARS e2dn float 1'
    write(iunit,'(a)')'LOOKUP_TABLE default'
@@ -899,11 +916,32 @@ if(output_smooth_press_fielde==1) then
    enddo 
 end if
 
+if(output_eviscosity_fielde==1) then  
+   write(iunit,'(a)')'SCALARS eff_visc_e float 1' 
+   write(iunit,'(a)')'LOOKUP_TABLE default' 
+   do ie=1,ov%nleaves 
+      if (elvoid(ie).eq.0) then 
+         write(iunit,'(e11.4)') ov%eviscosity(ie) 
+      endif 
+   enddo 
+end if 
+	 
+if(output_is_plastic_fielde==1) then  
+   write(iunit,'(a)')'SCALARS is_plastic_e integer 1' 
+   write(iunit,'(a)')'LOOKUP_TABLE default' 
+   do ie=1,ov%nleaves 
+      if (elvoid(ie).eq.0) then 
+         write(iunit,'(I10)') ov%is_plastic(ie) 
+      endif 
+   enddo 
+end if
+
 write(iunit,'(a)')'done visu.vtk'
 close(iunit)
 write(*,*) '--------------------------------------------------------------------------'
 write(*,*)'opla I am done producing visu.vtk '
 
+deallocate(countnode)
 
 !==============================================================================
 !======[produce norm<i>.vtk file]==============================================
@@ -1179,4 +1217,4 @@ endif
 
 write(*,*) '**************************************************************************'
 
-end
+end
\ No newline at end of file