diff --git a/chkhist.f b/chkhist.f
new file mode 100644
index 000000000..22522dbb7
--- /dev/null
+++ b/chkhist.f
@@ -0,0 +1,23 @@
+      subroutine chkhist(mrsym,nmax,ipk)
+
+      integer mrsym(63)
+      integer hist(0:63)
+
+      do i=0,63
+         hist(i)=0
+      enddo
+      do j=1,63
+         i=mrsym(j)
+         hist(i)=hist(i)+1
+      enddo
+
+      nmax=0
+      do i=0,63
+         if(hist(i).gt.nmax) then
+            nmax=hist(i)
+            ipk=i+1
+         endif
+      enddo
+
+      return
+      end