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

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