C  PROGRZMMED BY BOB SNYDER  4/15/67  FOR  I.B.M. 1620

C PROGRAM TO GENERATE A TABLE OF POWERS OF FILM FACTORS MULTIPLIED BY RE
C INTENSITY UNITS TO BE USED IN SCALING INTENSITIES TO THE FIRST FILM (F
C MULTIPLE FILM TECHNIQUES INVOLVING UP TO 10 FILMS)

C  INPUT DATA

C***********************************************************************
C  CARD 1 TITLE CARD

C     COL1-80  ANY ALPHAMERIC INFORMATION

C***********************************************************************
C  CARD 2  (I5,F5.0,I1)

C     COL 1-5 THE HIGHEST VALUE(IN RELATIVE INTENSITY UNITS) ON YOUR FIL
C     COL5-10 THE INTERVAL DESIRED BETWEEN EACH VALUE OF INTENSITY(USUAL
C     COL11 =1 MEANS READ ANEW TITLE AND NEW DATA CARDS

C           =2 MEANS END OF JOB

C***********************************************************************
C  CARD 3  (5F10.4)   FILM FACTORS

C     COL1-10 THE FILM FACTOR BETWEEN FILM 1 AND FILM 2

C     COL11-20 FF BETWEEN FILM 2 AND FILM 3 ETC.

C***********************************************************************
C  CARD4  (4F10.4)    FILM FACTORS

C     COL1-10 THE FILM FACTOR BETWEEN FILM 6 AND FILM 7

C     COL11-20  FF BETWEEN FILM 7 AND FILM 8  ETC.

C***********************************************************************
C***********************************************************************
C                 SAMPLE  DATA

C  FILM FACTOR TABLE FOR ZERO LAYER DATA  PYHBRCL2

C 150 1.0

C  2.75     2.60     2.60      2.60     2.60

C 2.60     0.0        0.0      0.0

C***********************************************************************
C***********************************************************************
C  PROGRAM

	DIMENSION TITLE(80)

60	READ61,TITLE

	PRINT70

	PRINT12

	PRINT61,TITLE

	PRINT12

	READ 2,MAX,AINT,ID

	READ3,FF12,FF23,FF34,FF45,FF56

	READ4,FF67,FF78,FF89,FF910

	PRINT 20,MAX

	PRINT 21,AINT

	PRINT23,FF12,FF23,FF34,FF45,FF56,FF67,FF78,FF89,FF910

	PRINT12

	PRINT1

	A=1.0

	M=0

6	IF(15-M)10,10,11

10	PRINT12

	PRINT1

	M=0

11	FTWO=A*FF12

	FTHR=FTWO*FF23

	FFOU=FTHR*FF34

	FFIV=FFOU*FF45

	FSIX=FFIV*FF56

	FSEV=FSIX*FF67

	FEIG=FSEV*FF78

	FNIN=FEIG*FF89

	FTEN=FNIN*FF910

	PRINT5,A,FTWO,FTHR,FFOU,FFIV,FSIX,FSEV,FEIG,FNIN,FTEN

	N=A

	IF(MAX-N)7,7,8

8	A=A+AINT

	M=M+1

	GO TO 6

7	IF(ID-1)50,55,50

55	GO TO 60

50	CALL EXIT

70	FORMAT(1H1)

61	FORMAT(80A1)

1	FORMAT(107H READING    FILM 2     FILM 3     FILM 4     FILM 5

	1FILM 6     FILM 7     FILM 8     FILM 9     FILM 10,/)

2	FORMAT(I5,F5.0,I1)

3	FORMAT(5F10.4)

4	FORMAT(4F10.4)

20	FORMAT(26H MAXIMUM INTENSITY VALUE =,I5/)

21	FORMAT(42H THE INTERVAL BETWEEN INTENSITY VALUES IS ,F5.0/)

23	FORMAT(22H THE FILM FACTORS ARE ,9F10.4/)

12	FORMAT(1H0)

5	FORMAT(F6.0,4XF9.0,8(2XF9.0))

