Initial commit

1 files changed, 26 insertions, 0 deletions

diff --git a/pkg/tbtables/fitsio/ftpi1b.f b/pkg/tbtables/fitsio/ftpi1b.f
new file mode 100644
index 00000000..dc97b8b4
--- /dev/null
+++ b/pkg/tbtables/fitsio/ftpi1b.f
@@ -0,0 +1,26 @@
+C----------------------------------------------------------------------
+        subroutine ftpi1b(ounit,nvals,incre,chbuff,status)
+
+C       Write an array of Integer*1 bytes to the output FITS file.
+
+        integer nvals,incre,ounit,status,i,offset
+        character*1 chbuff(nvals)
+
+C       ounit   i  fortran unit number
+C       nvals   i  number of pixels in the i2vals array
+C       incre   i  byte increment between values
+C       chbuff  c*1 array of input byte values
+C       status  i  output error status
+
+        if (incre .le. 1)then
+                call ftpcbf(ounit,0,nvals,chbuff,status)
+        else
+C               offset is the number of bytes to move between each value
+                offset=incre-1
+                call ftpcbf(ounit,0,1,chbuff,status)
+                do 10 i=2,nvals
+                        call ftmoff(ounit,offset,.true.,status)
+                        call ftpcbf(ounit,0,1,chbuff(i),status)
+10              continue
+        end if
+        end