Fixed typos

src/read_input_file.f90

@@ -250,11 +250,11 @@ do i=0,params%nmat
 
    mat(i)%plasticity_ss_type_coh='tot_strain'
    if (iproc==0) call scanfile (params%infile,'plasticity_ss_type_coh'//cm(il:3),mat(i)%plasticity_ss_type_coh,ires)
-   call mpi_bcast(mat(i)%plasticity_ss_type_coh,12,mpi_character,0,mpi_comm_world,ierr)
+   call mpi_bcast(mat(i)%plasticity_ss_type_coh,16,mpi_character,0,mpi_comm_world,ierr)
 
    mat(i)%plasticity_ss_type_phi='tot_strain'
    if (iproc==0) call scanfile (params%infile,'plasticity_ss_type_phi'//cm(il:3),mat(i)%plasticity_ss_type_phi,ires)
-   call mpi_bcast(mat(i)%plasticity_ss_type_phi,12,mpi_character,0,mpi_comm_world,ierr)
+   call mpi_bcast(mat(i)%plasticity_ss_type_phi,16,mpi_character,0,mpi_comm_world,ierr)
 
    mat(i)%plasticity_parameters=0.d0
    mat(i)%plasticity_parameters(3)=-1.d0

src/update_cloud_structure.f90

@@ -380,7 +380,7 @@ do i=1,np
       ww=h(k)/dxyz**3
       w(ic)=w(ic)+ww
       strain(ic)=strain(ic)+ww*cl%strain(i)
-      e2dp(ic)=e2dp(uc)+ww*cl%e2dp(i)
+      e2dp(ic)=e2dp(ic)+ww*cl%e2dp(i)
    enddo
 enddo
 

src/vrm.f90

@@ -281,7 +281,8 @@ select case (trim(plasticity_type))
     strain_soft_in_c = plasticity_parameters(6)
     if (strain_soft_in_phi.gt.0.d0 .or. strain_soft_in_c.gt.0.d0) then
       call strain_soften (c,phi,plasticity_ss_type_coh,plasticity_ss_type_phi, &
-                          plasticity_parameters,straintot,e2d,flag_vrm_pb)
+                          plasticity_parameters,straintot,e2d,e2dprev,         &
+                          flag_vrm_pb)
     endif
     friction_angle = phi*180.d0/pi
     sin_phi = sin(phi)
@@ -326,7 +327,7 @@ select case (trim(plasticity_type))
     if (strain_soft_in_c.gt.0.d0) then
       call strain_soften (yield,phi,plasticity_ss_type_coh,                    &
                           plasticity_ss_type_phi,plasticity_parameters,        &
-                          straintot,e2d,flag_vrm_pb)
+                          straintot,e2d,e2dprev,flag_vrm_pb)
     endif
     cos_theta = cos(theta)
     fail  = 2.d0*viscosity*e2d*cos_theta-yield 
@@ -350,7 +351,8 @@ end subroutine vrm
 !-------------------------------------------------------------------------------
 
 subroutine strain_soften (c,phi,plasticity_ss_type_coh,plasticity_ss_type_phi, &
-                          plasticity_parameters,straintot,e2d,flag_vrm_pb)
+                          plasticity_parameters,straintot,e2d,e2dprev,         &
+                          flag_vrm_pb)
 
 use constants
 
@@ -368,7 +370,7 @@ use constants
 implicit none
 character (len=16) plasticity_ss_type_coh,plasticity_ss_type_phi
 double precision plasticity_parameters(9)
-double precision :: c,phi,straintot,e2d
+double precision :: c,phi,straintot,e2d,e2dprev
 logical :: flag_vrm_pb
 
 !------------------------------------------------------------------------------|