File: EX3FNV.FT of Tape: Various/ETH/eth11-3
(Source file text)
DIMENSION INPUT(256),OUTPUT(7,6),EMU(6),SIGMA(6),SIZE(6) DIMENSION ITABLE(51),VTYPE(3,2) DATA VTYPE/' ',' VA','LLEY',' T','HRES','HOLD'/ DATA EMU/20.,75.,200.,500.,600.,900./ DATA SIGMA/20.,10.,75.,30.,35.,100./ DATA SIZE/500.,100.,200.,800.,700.,300./ DATA IT,IC,X/151,0,0./ DATA ITABLE/0,3,50,1,47*0/ DATA J,INPTR,IDIMO,NPEAKS/0,0,6,0/ DO 3 I=1,1024 A=0. X=X+1. DO 2 JJ=1,6 2 A=A+SIZE(JJ)*EXP(-.5*((X-EMU(JJ))/SIGMA(JJ))**2) IF(A.LE.IT) GO TO 4 IF(IC.EQ.0) GO TO 6 J=J+1 INPUT(J)=IC ASSIGN 21 TO IRET GO TO 30 21 J=J+1 INPUT(J)=0 ASSIGN 22 TO IRET GO TO 30 22 IC=0 6 J=J+1 IF(I.NE.1) GO TO 7 INPUT(1)=0 INPUT(2)=0 J=3 7 INPUT(J)=A ASSIGN 3 TO IRET GO TO 30 4 IF(J.EQ.0.OR.IC.NE.0) GO TO 20 J=J+1 INPUT(J)=0 ASSIGN 20 TO IRET GO TO 30 20 IC=IC+1 3 CONTINUE IF(J.EQ.0) GO TO 10 ASSIGN 10 TO IRET 30 IF(J.NE.256.AND.I.LT.1024) GO TO IRET CALL NVELOP(ITABLE,INPUT,J,INPTR,OUTPUT,IDIMO,NPEAKS) J=0 GO TO IRET 10 TYPE 1002 1002 FORMAT(1H1,'PEAK NO.',8X,'TYPE',8X,'AREA',4X,'CENTROID',5X, 1 'P WIDTH',/,15X,'P TIME',4X,'P HEIGHT',3X,'LEAD TIME') DO 11 L=1,NPEAKS KK=1+OUTPUT(1,L) 11 TYPE 1003,L,(VTYPE(K,KK),K=1,3),(OUTPUT(I,L),I=2,7) 1003 FORMAT(I9,3A4,3F12.0,/,9X,3F12.0) END