diff --git a/src/phase_mechanical.f90 b/src/phase_mechanical.f90
index bfbaa1133..f4e400c38 100644
--- a/src/phase_mechanical.f90
+++ b/src/phase_mechanical.f90
@@ -481,7 +481,7 @@ function integrateStress(F,subFp0,subFi0,Delta_t,ph,en) result(broken)
           !crystallite_converged(ipc,ip,el) = .true.  !> See "phase_mechanical_constitutive" and "homogenization_mechanical_response"
           !crystallite_todo(ipc,ip,el) = .false. !> Can't find this
           ! _converged = .not. broken
-          subFp0 = matmul(deltaFp,subFp0) 
+          !subFp0 = matmul(deltaFp,subFp0) 
 ! subFp0 is input need to change "phase_mechanical_Fp(ph)%data(1:3,1:3,en) = Fp_new / math_det33(Fp_new)**(1.0_pReal/3.0_pReal)"
 
           !plasticState(ph)%state()
diff --git a/src/phase_mechanical_plastic_phenopowerlaw.f90 b/src/phase_mechanical_plastic_phenopowerlaw.f90
index 75390e6be..d58d0d812 100644
--- a/src/phase_mechanical_plastic_phenopowerlaw.f90
+++ b/src/phase_mechanical_plastic_phenopowerlaw.f90
@@ -281,7 +281,6 @@ do ph = 1, phases%length
   allocate(geom(ph)%IPareaNormal(3,nIPneighbors,Nmembers))
   allocate(geom(ph)%IParea(nIPneighbors,Nmembers))
   allocate(geom(ph)%IPcoordinates(3,Nmembers))
-  call storeGeometry(ph)                                                                     !Achal delete
 
 !--------------------------------------------------------------------------------------------------
 ! state aliases and initialization
@@ -702,7 +701,7 @@ associate(prm => param(ph), stt => state(ph))
   where(tau_tw > 0.0_pReal)
     dot_gamma_tw = (1.0_pReal-sum(stt%gamma_tw(:,en)/prm%gamma_char)) &                           ! only twin in untwinned volume fraction
                  * prm%dot_gamma_0_tw*(abs(tau_tw)/stt%xi_tw(:,en))**prm%n_tw
-    fdot_twin = (0.05_pReal*(abs(tau_tw)/state(ph)%xi_tw(:,en))**param(ph)%n_tw)/param(ph)%gamma_char
+    fdot_twin = (0.05_pReal*(abs(tau_tw)/stt%xi_tw(:,en))**prm%n_tw)/prm%gamma_char   !Achal 0.05 is constant
   else where
     dot_gamma_tw = 0.0_pReal
   end where
@@ -719,40 +718,4 @@ end associate
 
 end subroutine kinetics_tw
 
-!--------------------------------------------------------------------------------------------------
-!> Achal: Added from nonlocal
-!--------------------------------------------------------------------------------------------------
-
-subroutine storeGeometry(ph)
-
-  integer, intent(in) :: ph
-
-  integer :: ce, co, nCell
-  real(pReal), dimension(:), allocatable :: V
-  integer, dimension(:,:,:), allocatable :: neighborhood
-  real(pReal), dimension(:,:), allocatable :: area, coords
-  real(pReal), dimension(:,:,:), allocatable :: areaNormal
-
-  nCell = product(shape(IPvolume))
-
-  V = reshape(IPvolume,[nCell])
-  neighborhood = reshape(IPneighborhood,[3,nIPneighbors,nCell])
-  area = reshape(IParea,[nIPneighbors,nCell])
-  areaNormal = reshape(IPareaNormal,[3,nIPneighbors,nCell])
-  coords = reshape(discretization_IPcoords,[3,nCell])
-
-  do ce = 1, size(material_homogenizationEntry,1)
-    do co = 1, homogenization_maxNconstituents
-      if (material_phaseID(co,ce) == ph) then
-        geom(ph)%V_0(material_phaseEntry(co,ce)) = V(ce)
-        geom(ph)%IPneighborhood(:,:,material_phaseEntry(co,ce)) = neighborhood(:,:,ce)
-        geom(ph)%IParea(:,material_phaseEntry(co,ce)) = area(:,ce)
-        geom(ph)%IPareaNormal(:,:,material_phaseEntry(co,ce)) = areaNormal(:,:,ce)
-        geom(ph)%IPcoordinates(:,material_phaseEntry(co,ce)) = coords(:,ce)
-      end if
-    end do
-  end do
-
-end subroutine
-
 end submodule phenopowerlaw