diff --git a/ChrMod_main.m b/ChrMod_main.m
index e0f06cb..6bd9ab1 100644
--- a/ChrMod_main.m
+++ b/ChrMod_main.m
@@ -99,7 +99,7 @@ function ChrMod_main(ChrNo,n)
     %Grid search for omega 
     lnL_om=GS_omega(ChrIF,weight,Ensemble,EstLen,Dim);
     omega=find(lnL_om==max(lnL_om));
-    omega=0.25+0.25*omega;
+    omega=0.25+0.25*mean(omega);
     OmegaRec(k)=omega;
     
 end