C -+-+-+-+-+  H A E C F L . F T  -+-+-+-+-+ C C HERE WE PRINT THE SUMMED INTENSITY OF THE POWER SPECTRUM STORED ON UNIT 5 C SUBROUTINE CFLACH INCLUDE HAEBUF.FI INCLUDE HAEGSA.FI INCLUDE HAEPTI.FI INCLUDE HABRK.FI INCLUDE HAEI85.FI INCLUDE HAECSZ.FI INCLUDE HAETTY.FI C C REAL COND INTEGER NCOND,TOASCI EXTERNAL COND,NCOND,TOASCI C REAL FFTX(85), @ PART OF THE FFT BUFFER * SUMM, @ SUMM OF ALL INTENSITIES * PERCNT(10),CNTPER(10) @ PERCENTAGE OF THE WHOLE POWER, FREQUENCY OF THIS VALUE (COMPUTED) INTEGER I,J,INCR2,PNT, * RANGE @ POINTER TO THE PERCNT VECTOR LOGICAL F1 @ SAMPLE BUFFER IS EMPTY FLAG C DATA PERCNT /10.,20.,40.,50.,60.,70.,80.,90.,95.,99./ C C C OPEN THE INPUT DATA FILE UNIT 5 C READ THE FILE HEADER FIRST BLOCK AND C THE DATA HEADER, THE SECOND BLOCK OF UNIT 5 C INSERT THE COMMON CEGESA FROM THE DATA FILE C IF (REC5.LT.1) DEFINE FILE5(MAXBL5,85,U,REC5) REC5=1 F1=SAMCNT.EQ.0 .OR. .NOT. COMP @ HEADER OF FILE 5 DOES NOT MATCH THE HEADER OF THE ACTUAL SAMPLE BUFFER READ (5'REC5) (RCRD0(I),I=1,85) @ READ THE FILE HEADER IF (F1) SAMCNT=0 @ CLEAR THE ACTUAL SAMPLE BUFFER IF (F1) COMP=.FALSE. @ AND THE COMPUTATIONS FLAG TOO IF (REDVAL.NE.6H2DPWLD) GOTO 20 @ RECORD IS EMPTY REDVAL=0 @ FOR SAVETY ONLY READ (5'REC5) LCEGSA(1),(LCEGSA(I),I=2,LCEGSA(1)) @ READ THE DATA HEADER BLOCK ( COMMON CEGESA) SETNUM=SETNU5 CHANEL=CHANE5 C INCR2=INCR/2 @ NUMBER OF DATA POINTS IN POWER SPECTRUM PLTEND=NCOND(PLTEND,PLTEND,SAMRAT/2) @ DEFAULT END OF FREQUENCY SCALE IN HERTZ! C IF BREAK(11) RETURN @ USER ABORTS PROGRAM WRITE (8,2) PERCNT @ WRITE A HEADER C C C COMPUTE THE WHOLE SUMM OF ALL INTENSITIES C XMAXI=0 @ THIS HOLDS THE SUMM DO 40 I=1,INCR/2 J=MOD(I-1,85) IF (J.EQ.0) READ (5'REC5) FFTX 40 XMAXI=FFTX(J+1)+XMAXI C IF (OPTION.EQ.12) GOTO 30 @ PF* TAKES THE DEFAULT VALUES AND DOESNT ASK QUESTIONS IF BREAK(11) RETURN @ USER SWITCHED OFF THE TASK 30 CONTINUE C C L1=.FALSE. REC5=3 @ POWER SPECTRUM IS STORED ON REC 3 TO 17 SUMM=0 RANGE=1 @ SET THE POINTER FOR THE INTENSITY RANGE SCAN CALL MOVE (-10,99.,CNTPER) DO 10 I=PLTBEG*SPAN,PLTEND*SPAN,1 J=MOD(I,85) IF (J.EQ.0) READ (5'REC5) FFTX SUMM=FFTX(J+1)+SUMM IF (SUMM.LT.XMAXI*PERCNT(RANGE)*.01) GOTO 10 CNTPER(RANGE)=FLOAT(I)/FLOAT(SPAN) @ THE UNITS ARE HERTZ (HZ!) RANGE=RANGE+1 10 CONTINUE SS=TOASCI(MOD(PLTBEG,60)) MM=PLTBEG/60 @ COMPUTE THE TIME OF THE POWER SPECTRUM WRITE (8,3) SETNUM,CHANEL,(LABEL(I),I=1,9), * MM,SS,SPAN,CNTPER,XMAXI RETURN C C 20 CONTINUE @ FILE 5 IS EMPTY COMPUT=0 SAMCNT=0 WRITE (TTO,1) @ NO DATA TO PLOT RETURN C C 1 FORMAT (' KEINE DATEN IN DER DATEI 5') 2 FORMAT (1H1,T72,10F6.1) 3 FORMAT (I4,I2,':'9A6,I4,1H:,A2,1H-,I2,10F6.1,1P,E11.2) END