dot state in loop

src/crystallite.f90

@@ -1505,15 +1505,26 @@ subroutine integrateStateRKCK45
            nonlocalBroken = .true.
          if(.not. crystallite_todo(g,i,e)) cycle
 
+         call constitutive_collectDotState(crystallite_S(1:3,1:3,g,i,e), &
+                                           crystallite_partionedF0, &
+                                           crystallite_Fi(1:3,1:3,g,i,e), &
+                                           crystallite_partionedFp0, &
+                                           crystallite_subdt(g,i,e)*CC(stage), g,i,e)
+         crystallite_todo(g,i,e) = all(.not. IEEE_is_NaN(plasticState(p)%dotState(:,c)))
+         do s = 1, phase_Nsources(p)
+           crystallite_todo(g,i,e) = crystallite_todo(g,i,e) .and. all(.not. IEEE_is_NaN(sourceState(p)%p(s)%dotState(:,c)))
+         enddo
+         if(.not. (crystallite_todo(g,i,e) .or. crystallite_localPlasticity(g,i,e))) &
+           nonlocalBroken = .true.
+         if(.not. crystallite_todo(g,i,e)) cycle
+
        endif
    enddo; enddo; enddo
    !$OMP END PARALLEL DO
 
-  if(nonlocalBroken) where(.not. crystallite_localPlasticity) crystallite_todo = .false.
-    call update_dotState(CC(stage))
-
  enddo
 
+  if(nonlocalBroken) where(.not. crystallite_localPlasticity) crystallite_todo = .false.
 
 !--------------------------------------------------------------------------------------------------
 ! --- STATE UPDATE WITH ERROR ESTIMATE FOR STATE ---