1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
C--------------------------------------------------------------------------
subroutine ftp3di(ounit,group,dim1,dim2,nx,ny,nz,array,status)
C Write a 3-d cube of i*2 values into 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 written).
C ounit i Fortran output unit number
C group i number of the data group, if any
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 i*2 the array of values to be written
C status i returned error stataus
C written by Wm Pence, HEASARC/GSFC, June 1991
integer ounit,group,dim1,dim2,nx,ny,nz,status
integer*2 array(dim1,dim2,*)
integer fpixel,row,band
fpixel=1
do 20 band=1,nz
do 10 row = 1,ny
call ftppri(ounit,group,fpixel,nx,array(1,row,band),status)
fpixel=fpixel+nx
10 continue
20 continue
end
|
