C	COMTST
C	PDP-9 ADVANCED SOFTWARE
C	FORTRAN IV TEST 1/2/68
	DIMENSION ADIFE(6),AFAC(3),APIPE(6),IMRCD(3),PP(27)
	1,ACPRI(3)
	DATA ASPVA,APIPE(1),APIPE(2),APIPE(4),APIPE(5)/
	1.60,12.09,12.09,5.0,5.0/,IMRCD(1),IMRCD(2)
	2/2*30/
	DATA ADIFE(1),ADIFE(2),ADIFE(4),ADIFE(5)/47.,47.,508.,3857048./
	DATA PP(1),PP(2),PP(3),PP(4)/4*2.0/,PP(5),PP(6),PP(7),PP(8)/
	14*3.0/,PP(9),PP(10),PP(11),PP(12)/4*4.0/,PP(13),PP(14),PP(15),
	2PP(16)/4*5.0/,PP(17),PP(18),PP(19),PP(20)/4*6.0/,PP(21),
	3PP(22),PP(23),PP(24)/4*7.0/,PP(25),PP(26),PP(27)/3*8.0/
	DATA AF,SC,AMEAS,FSUBB,ASUVA/37.96,3.1416,9.02,10.00,100.98/
	AGAST=38
	INORU=2
	WRITE (4,1)
1	FORMAT(16H FORTRAN IV TEST /20H COMTST ** BEGIN RUN)
205	BSPVA=(1./ASPVA)**.5
	DO 550 JCB=1,INORU
	AVEDE=FLOAT(IMRCD(JCB))
	BE=APIPE(JCB+3)/APIPE(JCB)
	IF(BE-.75)471,472,472
472	AK=.731
	GO TO 16
471	AG=.075
	DO 100 IE=1,27
	AG=AF+.025
	IF(AG-BE)100,100,110
100	CONTINUE
110	TOTA=PP(IE)
	TOTB=PP(IE-1)
	SC=.025-(AG-BE)
	IF(TOTA-TOTB)120,120,130
120	AK=TOTA
	GO TO 16
130	AK=TOTB+(SC*(TOTA-TOTB))/.025
16	FRD=830.-5000.*BE+9000.*BE**2-4200.*BE**3+(530./APIPE(JCB)**.5)
	AMEAB=AMEAS+14.4
	FR=1.+((FRD/(12835.*AK))/((AMEAB*AVEDE)**.5))
	XSUB2=AVEDE/(27.7*AMEAB)
	YTOTA=(XSUB2+1.)**.5
	YTOTB=.35*BE**4.+41.
	YTOTC=XSUB2/(1.3*YTOTA)
	YSUB2=YTOTA-YTOTB*YTOTC
	ACPRI(JCB)=YSUB2*FR*1.0177*FSUBB
550	AFAC(JCB)=ADIFE(JCB)*BSPVA
	AFTF=(520./(460.+AGAST))**.5
	AFPV=(1.+(ASUVA*AMEAS)/((AGAST+460.)**3.825))**.5
	FLOW=0
	DO 38 I=1,INORU
	AMWP=(ADIFE(1)*AMEAS)**.5/1000.
38	FLOW=FLOW+(AFAC(I)*AFPV*AMWP)
	WRITE(4,14)FLOW
14	FORMAT(F10.1)
	WRITE (4,2)
2	FORMAT (23H COMTST -- RUN COMPLETE)
	STOP
	END