C -+-+-+-+-+  H A E I M A . F T  -+-+-+-+-+ C C REAL FUNCTION IMAXNT (X,LFREQU,HFREQU,INDEX) REAL X(1024),LFREQU,HFREQU INTEGER INDEX @ OUTPUT PARAMETER C INCLUDE HAEBUF.FI INCLUDE HAETTY.FI INCLUDE HAEGSA.FI C REAL MAXI,STEP INTEGER LCHNL,HCHNL,I C C STATEMENT FUNCTION, RETURNS THE CHANNEL NUMBER FOR THE FREQUENCY RP ICHNL(RP)=INT(RP/FLOAT(SAMRAT)*FLOAT(INCR)+0.5) C STEP=FLOAT(SAMRAT)/FLOAT(INCR) @ HZ PER CHANNEL LCHNL=ICHNL(LFREQU)+1 @ LOWER FREQUENCY INCLUDED HCHNL=ICHNL(HFREQU) @ HIGH FREQUENCY EXCLUDED HCHNL=MIN0(HCHNL,INCR/2+1) @ MAKE A VALID NUMBER IF (LCHNL.GE.HCHNL) GOTO 20 @ HIGH/LOW MIXED C COMPUTE THE MAXIMUM OF THE INTENSITY IN THIS BAND LFREQU,HFREQU MAXI=0 DO 10 I=LCHNL,HCHNL IF (MAXI.GT.X(I)) GOTO 10 IMAXNT=FLOAT(I)*STEP INDEX=I MAXI=X(I) 10 CONTINUE RETURN 20 WRITE (TTO,1) LFREQU,HFREQU RETURN 1 FORMAT (' IMAXNT: FREQUENZPARAMETER DURCHEINANDER',1P,2E12.4) END