C -+-+-+-+-+  H A E P L O . F T  -+-+-+-+-+ C C PLOT THE POWER SPECTRUM STORED ON RECORD 3 TO STSCAN ON SCRATCH FILE 5 C SUBROUTINE PLOTIT INCLUDE HAEBUF.FI INCLUDE HAEGSA.FI INCLUDE HAEPTI.FI INCLUDE HABRK.FI INCLUDE HAEI85.FI INCLUDE HAECSZ.FI INCLUDE HAPPEN.FI INCLUDE HAETTY.FI C C REAL COND INTEGER NCOND EXTERNAL COND,NCOND C REAL FFTX(85), @ PART OF THE FFT BUFFER * ORD,ABZ,TEMSET, * TE1,TE2,TE3, @ TEMPORARY STORAGE * INTFAC,ORDMAX @ FACTOR TO SCALE THE INTENSITY OF EACH BAND, MAXIMUM OF THE ORDINATE INTEGER I,J,K,INCR2,PNT,BAND, * RELABS @ ABS INTENSITY =1, REL INTENSITY =2 LOGICAL F1, @ SAMPLE BUFFER IS EMPTY FLAG * L1 @ TEMPORARY STORAGE INTEGER PENPOS @ TEMPORARY PEN STATUS C C C C OPEN THE INPUT DATA FILE UNIT 5 C READ THE FILE HEADER FIRST BLOCK AND C THE DATA HEADER, THE SECOND BLOCK OF UNIT 5 C INSERT THE COMMON CEGESA FROM THE DATA FILE C IF (REC5.LT.1) DEFINE FILE5(MAXBL5,85,U,REC5) REC5=1 F1=SAMCNT.EQ.0 .OR. .NOT. COMP @ HEADER OF FILE 5 DOES NOT MATCH THE HEADER OF THE ACTUAL SAMPLE BUFFER READ (5'REC5) (RCRD0(I),I=1,85) @ READ THE FILE HEADER IF (F1) SAMCNT=0 @ CLEAR THE ACTUAL SAMPLE BUFFER IF (F1) COMP=.FALSE. @ AND THE COMPUTATIONS FLAG TOO IF (REDVAL.NE.6H2DPWLD) GOTO 20 @ RECORD IS EMPTY REDVAL=0 @ FOR SAVETY ONLY READ (5'REC5) LCEGSA(1),(LCEGSA(I),I=2,LCEGSA(1)) @ READ THE DATA HEADER BLOCK ( COMMON CEGESA) SETNUM=SETNU5 CHANEL=CHANE5 C HERE WE PLOT THE LAST POWER SPECTRUM STORED ON UNIT 5 REC 3 TO STSCAN C INCR2=INCR/2 @ NUMBER OF DATA POINTS IN POWER SPECTRUM PLTEND=NCOND(PLTEND,PLTEND,INCR2/SPAN) @ DEFAULT END OF FREQUENCY SCALE IF (OPTION.EQ.12) GOTO 30 @ PL* TAKES THE DEFAULT VALUES AND DOESNT ASK QUESTIONS C READ THE FREQUENCY AND INTENSITY RANGE OF THE PLOT 32 WRITE (TTO,4) PLTEND,PLTBEG READ (TTI,6) I,K,J @ ASK FOR THE FREQUENCY WINDOW TO PLOT IF BREAK(11) RETURN @ USER ABORTS PROGRAM PLTEND=NCOND(I,I,PLTEND) @ NO INPUT LOADS DEFAULT VALUES PLTBEG=NCOND(J,J,PLTBEG) IF (K.NE.1H ) PLTEND=I IF (K.NE.1H ) PLTBEG=J IF (PLTEND-PLTBEG .LE.0) GOTO 32 @ UPPER AND LOWER LIMITS MIXED, USER HAS TO RETYPE 33 WRITE (TTO,5) XMAXI @ ASK FOR THE SIZE OF INTENSITY SCALE READ (TTI,7) TE1 IF BREAK(11) RETURN @ USER SWITCHED OFF THE TASK XMAXI=COND(TE1,TE1,XMAXI) @ NO INPUT MEANS THE ACTUAL MAXIMUM IF (XMAXI.LE.0) GOTO 33 @ NEGATIVE POWER SPECTRUM INTENSITY IS NOT VALID 30 CONTINUE C CXCX XMAXI=ALOG(XMAXI) @ CXCX DISABLES LOG INTENSITY VERSION OF HAEPLO CALL STPLT @ START THE PLOTTER (CALL PLOTS) CALL XYPLOT (XOFSET,YOFSET,-PENUP) CALL FACTOR (XLEN/(SPAN*(PLTEND-PLTBEG)),YLEN/XMAXI) @ SIZE OF DIAGRAM C L1=.FALSE. PENPOS=PENUP DO 10 I=PLTBEG*SPAN,PLTEND*SPAN,1 J=MOD(I,85) IF (J.EQ.0) READ (5'REC5) FFTX ABZ=I ORD=AMAX1(0,FFTX(J+1)) CXCX ORD=AMIN1(ALOG(AMAX1(.1,FFTX(J+1))),XMAXI) L1=ORD.GT.XMAXI .AND. L1 @ TRUNCATE SUCH POINTS WHICH ARE TOO LARGE FOR THIS INTENSITY SCALE PENPOS=NCOND(L1,PENUP,PENPOS) L1=ORD.GT.XMAXI ORD=AMIN1(ORD,XMAXI) @ PLOT TO THE BOUNDER OF THIS DIAGRAM IF BREAK(11) GOTO 15 @ USER SWITCHED OFF PLOTTING CALL XYPLOT (ABZ,ORD,PENPOS) PENPOS=PENDWN 10 CONTINUE C C MAKE A SCALE C C FREQUENCY SCALE FROM PLTBEG TO PLTEND HERTZ C CALL FACTOR (XLEN/(PLTEND-PLTBEG),1.) CALL XYPLOT (0,-.5,-PENUP) IF BREAK(11) GOTO 15 @ USER SWITCHED OFF THE PLOTTER DO 12 I=PLTBEG,PLTEND,NCOND(PLTEND-PLTBEG.GE.100,10,1) @ 1 HERTZ STEP SIZE TE1=FLOAT(I-PLTBEG) CALL XYPLOT (TE1,0,PENDWN) J=MOD(I,NCOND(PLTEND-PLTBEG.GT.100,50,PLTEND.GT.30,10, * PLTEND.GT.15,5, * PLTEND.GT.6,2,1)) TE2=COND(J.EQ.0,-.7, J.EQ.5,-.5, -.3) CALL XYPLOT (TE1,TE2,PENDWN) IF (J.NE.0) GOTO 12 CALL WHERE (TE2,TE2,TE2,TE3) @ GET THE X FACTOR CALL FACTOR (1.,1.) CALL NUMBER (TE1*TE2-INT(ALOG10(AMAX1(1.,TE1)) * +.01)*YN/2.-YN*.3, * -1.3,YN,TE1,0,-1) CALL FACTOR (TE2,TE3) CALL XYPLOT (TE1,0,PENUP) 12 CALL XYPLOT (TE1,0,PENDWN) CALL XYPLOT (0,.5,-PENUP) CALL FACTOR (1.,1.) CALL SYMBOL (XLEN/2.,-YOFSET,YZ, * 'FREQUENCY [HZ]',0,14) C C INTENSITY SCALE LINEAR FROM 0 TO XMAXI C IF BREAK(11) GOTO 15 CALL FACTOR (1.,YLEN/XMAXI) CALL XYPLOT (-.5,0,-PENUP) TE2=XMAXI/10. DO 14 I=0,10 TE1=TE2*I CALL XYPLOT (0,TE1,PENDWN) J=MOD(I,10) TE3=COND(J.EQ.0,-.7, J.EQ.5,-.5, -.3) CALL XYPLOT (TE3,TE1,PENDWN) 14 CALL XYPLOT (0,TE1,PENDWN) C INSERT THE INTENSITY SCALE MAXIMUM CALL WHERE (TE2,TE2,TE2,TE3) CALL FACTOR (1.,1.) CALL ENUMBR (-1.1,TE1*TE3-4.*YN,YN,XMAXI,90.,1) CALL FACTOR (TE2,TE3) CALL XYPLOT (.5,0,-PENUP) CALL FACTOR (1.,1.) CALL SYMBOL (-XOFSET,YLEN/5.,YZ, * 'INTENSITY [ARB.UNITS]',90.,21) C C CALL FACTOR (1.,1.) I=BEGIN @ SAVE THE ORIGINAL BEGIN BEGIN=BEGFFT @ AND REPLACE WITH THE BEGIN OF THIS FFT IF BREAK(11) GOTO 15 @ USER SWITCHES OFF THE PLOTTER BEGIN=I @ RESTORE THE ORIGINAL BEGIN OF THE SCAN IF BREAK(11) GOTO 15 @ USER SWITCHES OFF THE PLOTTER CALL SYMBOL (YZ,YLEN+YZ,YZ,LABEL,0,MIN0(42,LABCNT*6)) IF (LABCNT.GT.7) * CALL SYMBOL (3.*YZ,YLENYZ,YZ,LABEL(8),0,LABCNT*6-42) C 15 CONTINUE @ ENTRY HERE IF USER SWITCHED OFF PLOTTING CALL FACTOR (1.,1.) IF (.NOT.BREAK(11)) CALL LABPLT (0,-YOFSET,PZEND) CALL XYPLOT(38.,25.,PENUP) CALL EXPLT RETURN C C 20 CONTINUE @ FILE 5 IS EMPTY COMPUT=0 SAMCNT=0 WRITE (TTO,1) @ NO DATA TO PLOT RETURN C C 1 FORMAT (' KEINE DATEN IN DER DATEI 5') 2 FORMAT (' SUMME DER LEISTUNGSSPEKTREN UEBER ',F4.1,' -->', * F5.1,' GEZEICHNET.') 4 FORMAT (' FREQUENZSKALA ENDE, ANFANG? (',I2,1H,,I2, * ') ? ',$) 5 FORMAT (' ORDINATEN MAXIMUM? (',1P,E8.1,') ? ',$) 6 FORMAT (I2,A1,I2) 7 FORMAT (E8.1) END