diff --git a/src/source_damage_anisoBrittle.f90 b/src/source_damage_anisoBrittle.f90
index 7473524e7..3f8f4b13d 100644
--- a/src/source_damage_anisoBrittle.f90
+++ b/src/source_damage_anisoBrittle.f90
@@ -75,8 +75,10 @@ subroutine source_damage_anisoBrittle_init
   do p = 1, size(config_phase)
     source_damage_anisoBrittle_instance(p) = count(phase_source(:,1:p) == SOURCE_DAMAGE_ANISOBRITTLE_ID)
     do sourceOffset = 1, phase_Nsources(p)
-      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ANISOBRITTLE_ID) &
+      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ANISOBRITTLE_ID) then
         source_damage_anisoBrittle_offset(p) = sourceOffset
+        exit
+      endif
     enddo
 
     if (all(phase_source(:,p) /= SOURCE_DAMAGE_ANISOBRITTLE_ID)) cycle
diff --git a/src/source_damage_anisoDuctile.f90 b/src/source_damage_anisoDuctile.f90
index c4df0f686..b09b5404d 100644
--- a/src/source_damage_anisoDuctile.f90
+++ b/src/source_damage_anisoDuctile.f90
@@ -70,8 +70,10 @@ subroutine source_damage_anisoDuctile_init
   do p = 1, size(config_phase)
     source_damage_anisoDuctile_instance(p) = count(phase_source(:,1:p) == SOURCE_DAMAGE_ANISODUCTILE_ID)
     do sourceOffset = 1, phase_Nsources(p)
-      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ANISODUCTILE_ID) &
+      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ANISODUCTILE_ID) then
         source_damage_anisoDuctile_offset(p) = sourceOffset
+        exit
+      endif
     enddo
 
     if (all(phase_source(:,p) /= SOURCE_DAMAGE_ANISODUCTILE_ID)) cycle
diff --git a/src/source_damage_isoBrittle.f90 b/src/source_damage_isoBrittle.f90
index b2f83a0f3..6e0e4d3fd 100644
--- a/src/source_damage_isoBrittle.f90
+++ b/src/source_damage_isoBrittle.f90
@@ -65,8 +65,10 @@ subroutine source_damage_isoBrittle_init
   do p = 1, size(config_phase)
     source_damage_isoBrittle_instance(p) = count(phase_source(:,1:p) == SOURCE_DAMAGE_ISOBRITTLE_ID)
     do sourceOffset = 1, phase_Nsources(p)
-      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ISOBRITTLE_ID) &
+      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ISOBRITTLE_ID) then
         source_damage_isoBrittle_offset(p) = sourceOffset
+        exit
+      endif
     enddo
 
     if (all(phase_source(:,p) /= SOURCE_DAMAGE_ISOBRITTLE_ID)) cycle
diff --git a/src/source_damage_isoDuctile.f90 b/src/source_damage_isoDuctile.f90
index f891faa81..a92b2cf91 100644
--- a/src/source_damage_isoDuctile.f90
+++ b/src/source_damage_isoDuctile.f90
@@ -64,8 +64,10 @@ subroutine source_damage_isoDuctile_init
   do p = 1, size(config_phase)
     source_damage_isoDuctile_instance(p) = count(phase_source(:,1:p) == SOURCE_DAMAGE_ISODUCTILE_ID)
     do sourceOffset = 1, phase_Nsources(p)
-      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ISODUCTILE_ID) &
+      if (phase_source(sourceOffset,p) == SOURCE_DAMAGE_ISODUCTILE_ID) then
         source_damage_isoDuctile_offset(p) = sourceOffset
+        exit
+      endif
     enddo
 
     if (all(phase_source(:,p) /= SOURCE_DAMAGE_ISODUCTILE_ID)) cycle
diff --git a/src/source_thermal_dissipation.f90 b/src/source_thermal_dissipation.f90
index 3efd9a762..acaf30068 100644
--- a/src/source_thermal_dissipation.f90
+++ b/src/source_thermal_dissipation.f90
@@ -54,8 +54,10 @@ subroutine source_thermal_dissipation_init
   do p = 1, size(config_phase)
     source_thermal_dissipation_instance(p) = count(phase_source(:,1:p) == SOURCE_THERMAL_DISSIPATION_ID)
     do sourceOffset = 1, phase_Nsources(p)
-      if (phase_source(sourceOffset,p) == SOURCE_THERMAL_DISSIPATION_ID) &
+      if (phase_source(sourceOffset,p) == SOURCE_THERMAL_DISSIPATION_ID) then
         source_thermal_dissipation_offset(p) = sourceOffset
+        exit
+      endif
     enddo
 
     if (all(phase_source(:,p) /= SOURCE_THERMAL_DISSIPATION_ID)) cycle
diff --git a/src/source_thermal_externalheat.f90 b/src/source_thermal_externalheat.f90
index 3cd4f9b32..1cf69b8a5 100644
--- a/src/source_thermal_externalheat.f90
+++ b/src/source_thermal_externalheat.f90
@@ -58,8 +58,10 @@ subroutine source_thermal_externalheat_init
   do p = 1, size(config_phase)
     source_thermal_externalheat_instance(p) = count(phase_source(:,1:p) == SOURCE_thermal_externalheat_ID)
     do sourceOffset = 1, phase_Nsources(p)
-      if (phase_source(sourceOffset,p) == SOURCE_thermal_externalheat_ID) &
+      if (phase_source(sourceOffset,p) == SOURCE_thermal_externalheat_ID) then
         source_thermal_externalheat_offset(p) = sourceOffset
+        exit 
+      endif
     enddo
 
     if (all(phase_source(:,p) /= SOURCE_thermal_externalheat_ID)) cycle