1 files changed, 40 insertions, 0 deletions

diff --git a/src/mt/srfindtab.f b/src/mt/srfindtab.f
new file mode 100644
index 0000000..ad5fdd4
--- /dev/null
+++ b/src/mt/srfindtab.f
@@ -0,0 +1,40 @@
+ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+c
+c HTCd - Copyright (C) 1998-2006 Henrik Rydberg
+c
+c This program is free software; you can redistribute it and/or modify
+c it under the terms of the GNU General Public License as published by
+c the Free Software Foundation; either version 2 of the License, or
+c (at your option) any later version.
+c
+c This program is distributed in the hope that it will be useful,
+c but WITHOUT ANY WARRANTY; without even the implied warranty of
+c MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+c GNU General Public License for more details.
+c
+c You should have received a copy of the GNU General Public License
+c along with this program; if not, write to the Free Software
+c Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
+c
+
+      integer function srfindtab(s,sn,c,cn,tab)
+      byte s(0:1),c(0:1),tab(0:1)
+      integer i,j,m,sn,cn
+      m=mod(cn,4)
+      do 10 i=sn-1,0,-1
+         do 20 j=0,m-1
+            if(tab(s(i+j)).ne.tab(c(j)))go to 50
+ 20      continue
+         do 30 j=m,cn-1,4
+            if(tab(s(i+j)).ne.tab(c(j)))go to 50
+            if(tab(s(i+j+1)).ne.tab(c(j+1)))go to 50
+            if(tab(s(i+j+2)).ne.tab(c(j+2)))go to 50
+            if(tab(s(i+j+3)).ne.tab(c(j+3)))go to 50
+ 30      continue
+ 50      if(j.eq.cn)go to 100
+ 10   continue
+      srfindtab=-1
+      return
+ 100  srfindtab=i
+      return
+      end