C   This program was NOT used to generate the system file - the source of
C    that file was lost. However, this program does manage to repeat it
C    fairly accurately. Note that this throughput applies to measurements
C    made through a 0.5 arcsec aperture!

	real otax(10000), otay(10000)
	real ccdx(10000), ccdy(10000)
	real wfpcx(10000), wfpcy(10000)

        open(1,file='ota',status='old')

	i=1
1	read(1,*,end=99) otax(i),otay(i)
        i=i+1
	goto 1
99	nota=i-1
	close(1)

        open(1,file='ccd.avg',status='old')
	i=1
101	read(1,*,end=199) ccdx(i),ccdy(i)
        i=i+1
	goto 101
199	nccd=i-1
	close(1)

        open(1,file='/home/avalon2/holtz/wfpc2/response/wfpc2.mod2',
     &       status='old')
	i=1
201	read(1,*,end=299) wfpcx(i),wfpcy(i)
        i=i+1
	goto 201
299	nwfpc=i-1
	close(1)

C   take out OTA obscurations
	pi=3.14159
        area=pi*1200.**2
        obssec=pi*396.**2
        obsspi=4.*(1200.-396.)*12.827*2.
        obspad=3.*pi*75.184**2-75.184*4.*12.827
        clearap=area-obssec-obsspi-obspad
        otathru=clearap/area

C   put in WFPC2 obscurations
        area=pi*1200.**2
        obssec=pi*(0.42*1200.)**2
        obsspi=1.*(1200.-0.42*1200.)*12.827*2.+
     &             3.*(1200.-0.42*1200.)*12.827*2.*2.
        obspad=3.*pi*75.184**2-75.184*4.*12.827
        clearap=area-obssec-obsspi-obspad
        wfpcthru=clearap/area

	do i=2000,12000,10
	  call linint(otax,otay,nota,float(i),ota)
	  call linint(wfpcx,wfpcy,nwfpc,float(i),wfpc)
	  call linint(ccdx,ccdy,nccd,float(i),ccd)
	  print *, float(i), ota*wfpc*ccd*wfpcthru/otathru
	end do


	stop
	end