iraf.net

Well thanks anyway. I did know about hte data types, and I was careful about
that. I was trying to use dtype=3 (short integers). It isn't exactly
filling every other pixl - it is every other column. For example if I
create an array:1 2 3 4 5
1 2 3 4 5
1 2 3 4 5
1 2 3 4 5
1 2 3 4 5The resultant image will be:0 1 0 2 0
0 1 0 2 0
0 1 0 2 0
0 1 0 2 0
0 1 0 2 0MOST Aggravating! In addition, any real values I stick into the header are
wrong, and it doesn't give a value for i_pixfile - so it can't find the pixels
and I have to edit afterwards (inconveninet, but liveable). ARGH! Who do
you suggest I talk to on Monday? ( I was hoping to get all the models run
this weekend, but not much point if I can't write the images! Too bad - Marc
Buie needs them for our HST proposal by Monday). I think I will go swim.Thanks.Karen.FYI - my test program: integer ixsz,iysz,axlen(7),naxis,dtype,ier,acmode,im,ibit,line(800),
1 l,ic,intensity(800,800)
logical bval
real bscale,bzero,bmin,bmax,platescale
character image*80,objectname*20,errmsg*80 bval = .true.
c
c - - - Initialize arrays
c
do 2 i=1,800
do 1 j=1,800
intensity(i,j) = i
1 continue
2 continue
print 5
5 format('Enter output image name: ',$)
read(*,14) image
15 format(a)
14 format(a80)
print 45
45 format('Image size in pixels x,y (max=800x800): ',$)
read*,ixsz,iysz
print 50
50 format('What is the desired plate scale in arcsec/pixel: ',$)
read*,platescale
print 65
65 format('Enter objectname (20 char max): ',$)
read(*,15) objectname
c
c - - - Setup parameters for writing out fits image
c
bzero = 1.234
bscale = 1000.1
bmax = 54.2
bmin = 1.2
ibit = 16
axlen(1) = ixsz
axlen(2) = iysz
naxis = 2
dtype = 3
acmode = 3
do 140 i=3,7
140 axlen(i) = 0
call imcrea(image,axlen,naxis,dtype,ier)
call imopen(image,acmode,im,ier)
c
c - - - Write header information
c
call imakwb(im,'SIMPLE',bval,' ',ier)
call imakwi(im,'BITPIX',ibit,'short integers',ier)
call imakwi(im,'NAXIS',naxis,'number of dimensions',ier)
call imakwi(im,'NAXIS1',ixsz,'number of pixels per row',ier)
call imakwi(im,'NAXIS2',iysz,'number of rows',ier)
call imakwr(im,'BSCALE',bscale,'scale factor',ier)
call imakwr(im,'BZERO',bmin,'zero level',ier)
call imakwr(im,'IRAF-MAX',bmax,'data max',ier)
call imakwr(im,'IRAF-MIN',bmin,'data min',ier)
call imakwc(im,'ORIGIN','cometimage','monte carlo program ver. 1',ier)
call imakwc(im,'OBJECT',objectname,comm,ier)
call imakwc(im,'IMAGETYP','OBJECT',comm,ier)
call imakwr(im,'SECPPIX',platescale,'arcsec / pixel',ier)
c
c - - - Write out CCD image
c
do 150 l=1,iysz
do 145 ic=1,ixsz
line(ic) = intensity(ic,l)
Print*,'line value = ',line(ic)
145 continue
call impl2s(im,line,l,ier)
if (ier .ne. 0) goto 500
150 continue
call imclos(im,ier)
stop
500 call imemsg(ier,errmsg)
write (*,fmt=501) errmsg
501 format('Error: ',a80)
stop
end

Bill,Again, sorry about the delay; this time I was on vacation. V2.10 will
be released in the next couple of months. If you want the fix and have the
necessary permissions or support to change the IRAF system here is what
you do.0. Login as IRAF and startup the CL. Load noao to define the logical
directories.1. Configure NOAO package for the appropriate architecture: no> cd noao
no> mkpkg <arch> Where <arch> would be sparc for a Sparc machine, ffpa for a Sun3
with float point accellerator, f68881 for a Sun3 with F68881
floating point processor. If these don't match your situation
let me know.2. Change onedspec$t_standard.x: no> cd onedspec
no> ed t_standard.xold:
w = wstart + x1 * wpc
call intrp (EXT_LOOKUP, Memr[EXT_WAVES(std)], Memr[EXT_EXTNS(std)],
EXT_NWAVES(std), w, extn, ierr)
flux = flux + (i1-x1) * spec[i1] * 10.0 ** (0.4 * extn * airmass) w = wstart + x2 * wpc
call intrp (EXT_LOOKUP, Memr[EXT_WAVES(std)], Memr[EXT_EXTNS(std)],
EXT_NWAVES(std), w, extn, ierr)
flux = flux + (1-(i2-x2)) * spec[i2] * 10.0 ** (0.4*extn*airmass)new:
if (i1 == i2) {
w = wstart + x1 * wpc
call intrp (EXT_LOOKUP, Memr[EXT_WAVES(std)], Memr[EXT_EXTNS(std)],
EXT_NWAVES(std), w, extn, ierr)
flux = flux + (x2-x1) * spec[i1] * 10.0 ** (0.4 * extn * airmass)
} else {
w = wstart + x1 * wpc
call intrp (EXT_LOOKUP, Memr[EXT_WAVES(std)], Memr[EXT_EXTNS(std)],
EXT_NWAVES(std), w, extn, ierr)
flux = flux + (i1-x1) * spec[i1] * 10.0 ** (0.4 * extn * airmass) w = wstart + x2 * wpc
call intrp (EXT_LOOKUP, Memr[EXT_WAVES(std)], Memr[EXT_EXTNS(std)],
EXT_NWAVES(std), w, extn, ierr)
flux = flux + (1-(i2-x2)) * spec[i2] * 10.0 ** (0.4*extn*airmass)
}3. Recompile and link: no> mkpkg -p noao updateSend mail to iraf@noao.edu if you have questions or problems with this.Cheers,
Frank