Repatch (from linux) of OSX IRAF

1 files changed, 39 insertions, 0 deletions

diff --git a/pkg/tbtables/fitsio/ftg3db.f b/pkg/tbtables/fitsio/ftg3db.f
new file mode 100644
index 00000000..be3b48a6
--- /dev/null
+++ b/pkg/tbtables/fitsio/ftg3db.f
@@ -0,0 +1,39 @@
+C--------------------------------------------------------------------------
+        subroutine ftg3db(ounit,group,nulval,dim1,dim2,nx,ny,nz,
+     &                    array,anyflg,status)
+
+C       Read a 3-d cube of byte values from the primary array.
+C       Data conversion and scaling will be performed if necessary
+C       (e.g, if the datatype of the FITS array is not the same
+C       as the array being read).
+
+C       ounit   i  Fortran output unit number
+C       group   i  number of the data group, if any
+C       nulval  c*1  undefined pixels will be set to this value (unless = 0)
+C       dim1    i  actual first dimension of ARRAY
+C       dim2    i  actual second dimension of ARRAY
+C       nx      i  size of the cube in the x direction
+C       ny      i  size of the cube in the y direction
+C       nz      i  size of the cube in the z direction
+C       array   c*1  the array of values to be read
+C       anyflg  l  set to true if any of the image pixels were undefined
+C       status  i  returned error stataus
+
+C       written by Wm Pence, HEASARC/GSFC, June 1991
+
+        integer ounit,group,dim1,dim2,nx,ny,nz,status
+        character*1 array(dim1,dim2,*),nulval
+        logical anyflg,ltemp
+        integer fpixel,row,band
+
+        anyflg=.false.
+        fpixel=1
+        do 20 band=1,nz
+        do 10 row = 1,ny
+                call ftgpvb(ounit,group,fpixel,nx,nulval,
+     &              array(1,row,band),ltemp,status)
+                if (ltemp)anyflg=.true.
+                fpixel=fpixel+nx
+10      continue
+20      continue
+        end