projects / pdp8.git / blame
| Commit | Line | Data |
|---|---|---|
| 81e70d48 PH |
1 | C -+-+-+-+-+ \ e H A E I M A . F T \ e -+-+-+-+-+\r |
| 2 | C\r | |
| 3 | C\r | |
| 4 | REAL FUNCTION IMAXNT (X,LFREQU,HFREQU,INDEX)\r | |
| 5 | REAL X(1024),LFREQU,HFREQU\r | |
| 6 | INTEGER INDEX @ OUTPUT PARAMETER\r | |
| 7 | C\r | |
| 8 | INCLUDE HAEBUF.FI\r | |
| 9 | INCLUDE HAETTY.FI\r | |
| 10 | INCLUDE HAEGSA.FI\r | |
| 11 | C\r | |
| 12 | REAL MAXI,STEP\r | |
| 13 | INTEGER LCHNL,HCHNL,I\r | |
| 14 | C\r | |
| 15 | C STATEMENT FUNCTION, RETURNS THE CHANNEL NUMBER FOR THE FREQUENCY RP\r | |
| 16 | ICHNL(RP)=INT(RP/FLOAT(SAMRAT)*FLOAT(INCR)+0.5)\r | |
| 17 | C\r | |
| 18 | STEP=FLOAT(SAMRAT)/FLOAT(INCR) @ HZ PER CHANNEL\r | |
| 19 | LCHNL=ICHNL(LFREQU)+1 @ LOWER FREQUENCY INCLUDED\r | |
| 20 | HCHNL=ICHNL(HFREQU) @ HIGH FREQUENCY EXCLUDED\r | |
| 21 | HCHNL=MIN0(HCHNL,INCR/2+1) @ MAKE A VALID NUMBER\r | |
| 22 | IF (LCHNL.GE.HCHNL) GOTO 20 @ HIGH/LOW MIXED\r | |
| 23 | C COMPUTE THE MAXIMUM OF THE INTENSITY IN THIS BAND LFREQU,HFREQU\r | |
| 24 | MAXI=0\r | |
| 25 | DO 10 I=LCHNL,HCHNL\r | |
| 26 | IF (MAXI.GT.X(I)) GOTO 10\r | |
| 27 | IMAXNT=FLOAT(I)*STEP\r | |
| 28 | INDEX=I\r | |
| 29 | MAXI=X(I)\r | |
| 30 | 10 CONTINUE\r | |
| 31 | RETURN\r | |
| 32 | 20 WRITE (TTO,1) LFREQU,HFREQU\r | |
| 33 | RETURN\r | |
| 34 | 1 FORMAT (' IMAXNT: FREQUENZPARAMETER DURCHEINANDER',1P,2E12.4) \r | |
| 35 | END\r | |
| 36 | \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 |