1 files changed, 71 insertions, 0 deletions
diff --git a/pkg/tbtables/fitsio/ftgtbs.f b/pkg/tbtables/fitsio/ftgtbs.f
new file mode 100644
index 00000000..2a659b66
--- /dev/null
+++ b/pkg/tbtables/fitsio/ftgtbs.f
@@ -0,0 +1,71 @@
+C----------------------------------------------------------------------
+        subroutine ftgtbs(iunit,frow,fchar,nchars,svalue,status)
+
+C       read a consecutive string of characters from an ascii or binary
+C       table. This will span multiple rows of the table if NCHARS+FCHAR is
+C       greater than the length of a row.
+
+C       iunit   i  fortran unit number
+C       frow    i  starting row number (1st row = 1)
+C       fchar   i  starting character/byte in the row to read (1st character=1)
+C       nchars  i  number of characters/bytes to read (can span multiple rows)
+C       svalue  c  returned string of characters
+C       status  i  output error status
+C
+C       written by Wm Pence, HEASARC/GSFC, June 1991
+
+        integer iunit,frow,fchar,nchars,status
+        character*(*) svalue
+
+C       COMMON BLOCK DEFINITIONS:--------------------------------------------
+        integer nf,nb,ne
+        parameter (nb = 20)
+        parameter (nf = 3000)
+        parameter (ne = 200)
+        integer bufnum,chdu,hdutyp,maxhdu,hdstrt,hdend,nxthdr,dtstrt
+        integer nxtfld
+        logical wrmode
+        common/ft0001/bufnum(199),chdu(nb),hdutyp(nb),maxhdu(nb),
+     &  wrmode(nb),hdstrt(nb,ne),hdend(nb),nxthdr(nb),dtstrt(nb),nxtfld
+        integer tfield,tstart,tbcol,rowlen,tdtype,trept,tnull,scount
+        integer theap,nxheap
+        double precision tscale,tzero
+        common/ft0002/tfield(nb),tstart(nb),tbcol(nf),rowlen(nb),
+     &  tdtype(nf),trept(nf),tscale(nf),tzero(nf),tnull(nf),scount(nb)
+     &  ,theap(nb),nxheap(nb)
+C       END OF COMMON BLOCK DEFINITIONS-----------------------------------
+
+        integer ibuff,bstart,nget
+        
+        if (status .gt. 0)return
+
+        ibuff=bufnum(iunit)
+
+C       check for errors
+        if (nchars .le. 0)then
+C               zero or negative number of character requested
+                return
+        else if (frow .lt. 1)then
+C               error: illegal first row number
+                status=307
+                return
+        else if (fchar .lt. 1)then
+C               error: illegal starting character
+                status=308
+                return
+        end if
+        
+C       move the i/o pointer to the start of the sequence of characters
+        bstart=dtstrt(ibuff)+(frow-1)*rowlen(ibuff)+fchar-1
+        call ftmbyt(iunit,bstart,.false.,status)
+
+C       get the string of characters, (up to the length of the input string)
+        if (len(svalue) .ne. 1)then
+            svalue=' '
+            nget=min(nchars,len(svalue))
+        else
+C           assume svalue was dimensioned as: character*1 svalue(nchars)
+            nget=nchars
+        end if
+        call ftgcbf(iunit,1,nget,svalue,status)
+        end