diff --git a/miscmaths.cc b/miscmaths.cc
index 92d58046f28808eb25fe19b6463800a56ab9c878..7b8cb6324f51a703dde663306ddec0ed249a0851 100644
--- a/miscmaths.cc
+++ b/miscmaths.cc
@@ -1497,6 +1497,24 @@ ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2)
   res.Release();
   return res;
 }
+ReturnMatrix geqt(const Matrix& mat,const float a) 
+{
+  int ncols = mat.Ncols();
+  int nrows = mat.Nrows();
+  Matrix res(nrows,ncols);
+  res=0.0;
+
+  for (int ctr1 = 1; ctr1 <= nrows; ctr1++) {
+    for (int ctr2 =1; ctr2 <= ncols; ctr2++) {
+      if( mat(ctr1,ctr2) >= a){
+	res(ctr1,ctr2) = 1.0;
+      }
+    }
+  }
+
+  res.Release();
+  return res;
+}
 
 
 ReturnMatrix leqt(const Matrix& mat1,const Matrix& mat2)