ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
cccc                     08/24/04                          cccc
cccc                    mkhead.for                         cccc
cccc   input - list of images and important header info    cccc
cccc   output - ascii file to be executed in iraf to       cccc 
cccc   update headers of images - superfast                cccc
cccc   need to update #images, and any fields that may     cccc
cccc   not be included in the observation file             cccc
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc  

	implicit none

	character*50 input, output
	character*10 image(87)
	character*10 field(87),filter(87),UT(87),LST(87),am(87)
	character*10 ra(87),dec(87)
	integer k
	
 200	format(a8,2x,a5,1x,a5,1x,a5,2x,a5,2x,a5)	
 300	format(a50)

	write(*,*), 'Enter name of observation log: '
	read(*,300), input
	
	write(*,*), 'Enter name of output ascii file: '
	read(*,300), output
	
c	initialize arrays to zero
	do k=1,87
		image(k)='a'
		field(k)='a'
		filter(k)='a'
		UT(k)='a'
		LST(k)='a'
		am(k)='a'
	enddo
	write(*,*), 'Initialization Complete'
	
c	read in observation log	
	open(unit=2,name=input,type='old',access='sequential')
	do k=1,87
		read(2,*),image(k),field(k),filter(k),UT(k),LST(k),am(k),ra(k),dec(k)
	enddo
	close(unit=2)
	write(*,*), 'Input Complete'
	
c	write out iraf header creation commands into ascii file
	open(unit=3,name=output,type='new',access='sequential')
	do k=1,87
		if (field(k).eq.'flat') then
			write(3,*), 'hedit ',image(k),' OBJECT ',field(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' FILTERS ',filter(k),' add+ verify-'
			write(3,*), 'hedit ',image(k),' CCDTYPE flat add+ verify-'
			write(3,*), 'hedit ',image(k),' IMAGETYP flat add+ verify-'
		elseif (field(k).eq.'bias') then
			write(3,*), 'hedit ',image(k),' OBJECT ',field(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' CCDTYPE bias add+ verify-'
			write(3,*), 'hedit ',image(k),' IMAGETYP zero add+ verify-'
		elseif (field(k).eq.'dark') then
                	write(3,*), 'hedit ',image(k),' OBJECT ',field(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' CCDTYPE dark add+ verify-'
			write(3,*), 'hedit ',image(k),' IMAGETYP dark add+ verify-'
                else
			write(3,*), 'hedit ',image(k),' AIRMASS ',am(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' FILTERS ',filter(k),' add+ verify-'
			write(3,*), 'hedit ',image(k),' OBJECT ',field(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' TELRA ',ra(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' TELDEC ',dec(k),' update+ verify-'
			write(3,*), 'hedit ',image(k),' CCDTYPE object add+ verify-'
			write(3,*), 'hedit ',image(k),' IMAGETYP object add+ verify-'
		endif
	enddo
	write(*,*), 'Output Complete'
	
	stop
	end