sw/rescue/lab8e_goettingen/disk2_11/rkb/paroff/haeint.ft

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