A9.13  Appendix: Computer Programs

*		        PROGRAM I.16				    *

*******************************************************************************
*                                FACTOR                                       *
* THIS IS THE PROGRAM THAT READS THE INPUT ANGLES OF THE ESCAPING             *
* TRAJECTORIES AND COMPUTES G FOR A GIVEN ENERGY.  IT HAS TO LOOP THROUGH     *
* ALL DETECTOR AREA ELEMENTS, AS WELL AS EACH ESCAPING TRAJECTORY.  NOTE      *
* THAT FOR MY SITUATION, I USED DELPHI = 2 DEGREES, AND DELTHETA = 2 DEGREES. *
* NOTE ALSO THAT I RENAMED THE *.E*IN FILES 1.DAT, 2.DAT,...ETC.              *
*******************************************************************************

	program factor

	real*8 G,theta,omega
	real*8 area(13)
	integer int,n,npas,j
	character*2 fname
	parameter (delth = .0349066,delphi = .0349066)

	G = 0.
	do n=1,13
	  int = n
	  write (fname,'(i2)') int
	  write (6,*) fname
	  open (unit=1,status='old',access='sequential',file=fname)
 	  read (1,*) area(n)
	  write (6,*) area(n)
          read (1,*) npas
	  write (6,*) npas
	  omega = 0.
	  do j=1,npas
	    read (1,*) theta
c	    write (6,*) theta
	    theta = theta*.0174533
	    omega = omega + delth*delphi*sin(theta + delth/2)
c	    write (6,*) 'omega',omega
	  end do
	  if (n.lt.6) then
	    G = G + area(n)*omega
	  else
	    G = G + 2*area(n)*omega
	  end if
	  close (1)
	end do
	write (6,*) 'G=',G
	end