diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp
index c3b17d61..59c305db 100644
--- a/models/ColorModel.cpp
+++ b/models/ColorModel.cpp
@@ -606,6 +606,10 @@ void ScaLBL_ColorModel::Run(){
 						fprintf(kr_log_file,"%i %.5g %.5g %.5g %.5g %.5g %.5g ",timestep-analysis_interval+20,muA,muB,5.796*alpha,Fx,Fy,Fz);
 						fprintf(kr_log_file,"%.5g %.5g %.5g %.5g %.5g %.5g %.5g %.5g\n",volA,volB,vA_x,vA_y,vA_z,vB_x,vB_y,vB_z);
 						fclose(kr_log_file);
+						// Add exit criteria based on relative permeability ratio
+						if (volA*sqrt(vA_x*vA_x + vA_y*vA_y + vA_z*vA_z)/(volB*sqrt(vB_x*vB_x + vB_y*vB_y + vB_z*vB_z)) < 0.05){
+							timestep = timestepMax;
+						}
 
 						printf("  Measured capillary number %f \n ",Ca);
 					}