A large commit.
[pdp8.git] / sw / rescue / lab8e_goettingen / disk2_11 / rkb / paroff / haeint.ft
CommitLineData
81e70d48
PH
1C -+-+-+-+-+ \ e H A E I N T . F T \ e -+-+-+-+-+\r
2C\r
3C\r
4 REAL FUNCTION INTENS (X,LFREQU,HFREQU)\r
5 REAL X(1024),LFREQU,HFREQU\r
6C\r
7 INCLUDE HAEBUF.FI\r
8 INCLUDE HAEGSA.FI\r
9 INCLUDE HAETTY.FI\r
10C\r
11 INTEGER LCHNL,HCHNL,I\r
12C\r
13C STATEMENT FUNCTION, RETURNS THE CHANNEL NUMBER FOR THE FREQUENCY RP\r
14 ICHNL(RP)=INT(RP/FLOAT(SAMRAT)*FLOAT(INCR)+0.5)\r
15C\r
16 LCHNL=ICHNL(LFREQU)+1 @ LOWER FREQUENCY INCLUDED\r
17 HCHNL=ICHNL(HFREQU) @ HIGH FREQUENCY EXCLUDED\r
18 HCHNL=MIN0(HCHNL,INCR/2+1) @ MAKE A VALID NUMBER\r
19 IF (LCHNL.GE.HCHNL) GOTO 20 @ HIGH/LOW MIXED\r
20C COMPUTE THE INTENSITY (SUMM THE X(I))\r
21 INTENS=0\r
22 DO 10 I=LCHNL,HCHNL\r
2310 INTENS=X(I)+INTENS\r
24 RETURN\r
2520 WRITE (TTO,1) LFREQU,HFREQU\r
26 RETURN\r
271 FORMAT (' INTENS: FREQUENZPARAMETER DURCHEINANDER',1P,2E12.4) \r
28 END\r
29\1a\0\0\0\0\0