include "fitsio.h"

procedure fsg3di(ounit,group,nulval,dim1,dim2,nx,ny,nz,
                    array,anyflg,status)

# Read a 3-d cube of i*2 values from the primary array.
# Data conversion and scaling will be performed if necessary
# (e.g, if the datatype of the FITS array is not the same
# as the array being read).

int     ounit           # i output file pointer
int     group           # i group number
short   nulval          # i value for undefined pixels
int     dim1            # i size of 1st dimension
int     dim2            # i size of 2nd dimension
int     nx              # i size of x axis
int     ny              # i size of y axis
int     nz              # i size of z axis
short   array[ARB]      # o array of values
bool    anyflg          # o any null values?
int     status          # o error status

begin

call ftg3di(ounit,group,nulval,dim1,dim2,nx,ny,nz,
                    array,anyflg,status)
end