projects / pdp8.git / blame
| Commit | Line | Data |
|---|---|---|
| 81e70d48 PH |
1 | C -+-+-+-+-+ \ e H A E C F L . F T \ e -+-+-+-+-+\r |
| 2 | C\r | |
| 3 | C HERE WE PRINT THE SUMMED INTENSITY OF THE POWER SPECTRUM STORED ON UNIT 5\r | |
| 4 | C\r | |
| 5 | SUBROUTINE CFLACH\r | |
| 6 | INCLUDE HAEBUF.FI\r | |
| 7 | INCLUDE HAEGSA.FI\r | |
| 8 | INCLUDE HAEPTI.FI\r | |
| 9 | INCLUDE HABRK.FI\r | |
| 10 | INCLUDE HAEI85.FI\r | |
| 11 | INCLUDE HAECSZ.FI\r | |
| 12 | INCLUDE HAETTY.FI\r | |
| 13 | C\r | |
| 14 | \fC\r | |
| 15 | REAL COND\r | |
| 16 | INTEGER NCOND,TOASCI\r | |
| 17 | EXTERNAL COND,NCOND,TOASCI\r | |
| 18 | C\r | |
| 19 | \f REAL FFTX(85), @ PART OF THE FFT BUFFER\r | |
| 20 | * SUMM, @ SUMM OF ALL INTENSITIES\r | |
| 21 | * PERCNT(10),CNTPER(10) @ PERCENTAGE OF THE WHOLE POWER, FREQUENCY OF THIS VALUE (COMPUTED)\r | |
| 22 | INTEGER I,J,INCR2,PNT,\r | |
| 23 | * RANGE @ POINTER TO THE PERCNT VECTOR\r | |
| 24 | LOGICAL F1 @ SAMPLE BUFFER IS EMPTY FLAG\r | |
| 25 | C\r | |
| 26 | DATA PERCNT /10.,20.,40.,50.,60.,70.,80.,90.,95.,99./\r | |
| 27 | C\r | |
| 28 | C\r | |
| 29 | C OPEN THE INPUT DATA FILE UNIT 5\r | |
| 30 | C READ THE FILE HEADER FIRST BLOCK AND\r | |
| 31 | C THE DATA HEADER, THE SECOND BLOCK OF UNIT 5\r | |
| 32 | C INSERT THE COMMON CEGESA FROM THE DATA FILE\r | |
| 33 | C\r | |
| 34 | IF (REC5.LT.1) DEFINE FILE5(MAXBL5,85,U,REC5)\r | |
| 35 | REC5=1\r | |
| 36 | F1=SAMCNT.EQ.0 .OR. .NOT. COMP @ HEADER OF FILE 5 DOES NOT MATCH THE HEADER OF THE ACTUAL SAMPLE BUFFER\r | |
| 37 | READ (5'REC5) (RCRD0(I),I=1,85) @ READ THE FILE HEADER\r | |
| 38 | IF (F1) SAMCNT=0 @ CLEAR THE ACTUAL SAMPLE BUFFER\r | |
| 39 | IF (F1) COMP=.FALSE. @ AND THE COMPUTATIONS FLAG TOO\r | |
| 40 | IF (REDVAL.NE.6H2DPWLD) GOTO 20 @ RECORD IS EMPTY\r | |
| 41 | REDVAL=0 @ FOR SAVETY ONLY\r | |
| 42 | READ (5'REC5) LCEGSA(1),(LCEGSA(I),I=2,LCEGSA(1)) @ READ THE DATA HEADER BLOCK ( COMMON CEGESA)\r | |
| 43 | SETNUM=SETNU5\r | |
| 44 | CHANEL=CHANE5\r | |
| 45 | C\r | |
| 46 | INCR2=INCR/2 @ NUMBER OF DATA POINTS IN POWER SPECTRUM\r | |
| 47 | PLTEND=NCOND(PLTEND,PLTEND,SAMRAT/2) @ DEFAULT END OF FREQUENCY SCALE IN HERTZ!\r | |
| 48 | C\r | |
| 49 | IF BREAK(11) RETURN @ USER ABORTS PROGRAM\r | |
| 50 | WRITE (8,2) PERCNT @ WRITE A HEADER \r | |
| 51 | C\r | |
| 52 | C\r | |
| 53 | C COMPUTE THE WHOLE SUMM OF ALL INTENSITIES\r | |
| 54 | C\r | |
| 55 | XMAXI=0 @ THIS HOLDS THE SUMM \r | |
| 56 | DO 40 I=1,INCR/2\r | |
| 57 | J=MOD(I-1,85)\r | |
| 58 | IF (J.EQ.0) READ (5'REC5) FFTX\r | |
| 59 | 40 XMAXI=FFTX(J+1)+XMAXI\r | |
| 60 | C\r | |
| 61 | IF (OPTION.EQ.12) GOTO 30 @ PF* TAKES THE DEFAULT VALUES AND DOESNT ASK QUESTIONS\r | |
| 62 | IF BREAK(11) RETURN @ USER SWITCHED OFF THE TASK\r | |
| 63 | 30 CONTINUE\r | |
| 64 | C\r | |
| 65 | C\r | |
| 66 | L1=.FALSE.\r | |
| 67 | REC5=3 @ POWER SPECTRUM IS STORED ON REC 3 TO 17\r | |
| 68 | \r | |
| 69 | SUMM=0\r | |
| 70 | RANGE=1 @ SET THE POINTER FOR THE INTENSITY RANGE SCAN\r | |
| 71 | CALL MOVE (-10,99.,CNTPER)\r | |
| 72 | DO 10 I=PLTBEG*SPAN,PLTEND*SPAN,1\r | |
| 73 | J=MOD(I,85)\r | |
| 74 | IF (J.EQ.0) READ (5'REC5) FFTX\r | |
| 75 | SUMM=FFTX(J+1)+SUMM\r | |
| 76 | IF (SUMM.LT.XMAXI*PERCNT(RANGE)*.01) GOTO 10\r | |
| 77 | CNTPER(RANGE)=FLOAT(I)/FLOAT(SPAN) @ THE UNITS ARE HERTZ (HZ!)\r | |
| 78 | RANGE=RANGE+1\r | |
| 79 | 10 CONTINUE\r | |
| 80 | SS=TOASCI(MOD(PLTBEG,60))\r | |
| 81 | MM=PLTBEG/60 @ COMPUTE THE TIME OF THE POWER SPECTRUM\r | |
| 82 | WRITE (8,3) SETNUM,CHANEL,(LABEL(I),I=1,9),\r | |
| 83 | * MM,SS,SPAN,CNTPER,XMAXI\r | |
| 84 | RETURN\r | |
| 85 | \fC\r | |
| 86 | C\r | |
| 87 | 20 CONTINUE @ FILE 5 IS EMPTY\r | |
| 88 | COMPUT=0\r | |
| 89 | SAMCNT=0\r | |
| 90 | WRITE (TTO,1) @ NO DATA TO PLOT\r | |
| 91 | RETURN\r | |
| 92 | C\r | |
| 93 | C\r | |
| 94 | 1 FORMAT (' KEINE DATEN IN DER DATEI 5')\r | |
| 95 | 2 FORMAT (1H1,T72,10F6.1)\r | |
| 96 | 3 FORMAT (I4,I2,':'9A6,I4,1H:,A2,1H-,I2,10F6.1,1P,E11.2)\r | |
| 97 | END\r | |
| 98 | \1a\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0 |