diff --git a/miscmaths.cc b/miscmaths.cc
index e9a092b76debaa23302ea478e6a7b4ccabfbf140..aba20f3940d0b258daec63b6b651e30072825397 100644
--- a/miscmaths.cc
+++ b/miscmaths.cc
@@ -1455,6 +1455,24 @@ ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2)
   res.Release();
   return res;
 }
+ReturnMatrix geqt(const Matrix& mat1,const float a) 
+{
+  int ncol = mat.Ncols();
+  int nrow = mat.Nrows();
+  Matrix res(nrow,ncol);
+  res=0.0;
+
+  for (int ctr1 = 1; ctr1 <= nrows; ctr1++) {
+    for (int ctr2 =1; ctr2 <= ncols; ctr2++) {
+      if( mat1(ctr1,ctr2) >= a){
+	res(ctr1,ctr2) = 1.0;
+      }
+    }
+  }
+
+  res.Release();
+  return res;
+}
 
 
 ReturnMatrix leqt(const Matrix& mat1,const Matrix& mat2) 
diff --git a/miscmaths.h b/miscmaths.h
index ddd9c48c57a549b6ec9e920d978d927b7d5c8e8a..84164f5c7c7e0899cffaeb68aef9a4f4ffbcb042 100644
--- a/miscmaths.h
+++ b/miscmaths.h
@@ -193,7 +193,8 @@ namespace MISCMATHS {
   ReturnMatrix min(const Matrix& mat);
   ReturnMatrix gt(const Matrix& mat1,const Matrix& mat2); 
   ReturnMatrix lt(const Matrix& mat1,const Matrix& mat2); 
-  ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2); 
+  ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2);  
+  ReturnMatrix geqt(const Matrix& mat1,const float a); 
   ReturnMatrix leqt(const Matrix& mat1,const Matrix& mat2); 
   ReturnMatrix eq(const Matrix& mat1,const Matrix& mat2); 
   ReturnMatrix neq(const Matrix& mat1,const Matrix& mat2);