X-Git-Url: http://gitweb.hachti.de/?a=blobdiff_plain;f=sw%2Frescue%2Flab8e_goettingen%2Fdisk2_11%2Frkb%2Fparoff%2Fhaenya.ft;fp=sw%2Frescue%2Flab8e_goettingen%2Fdisk2_11%2Frkb%2Fparoff%2Fhaenya.ft;h=8c031c6ad00fe69d1e5a5b2fdfd4734748c417ee;hb=81e70d488b71bf995c459ca3a02c025993460ffa;hp=0000000000000000000000000000000000000000;hpb=07ec0278333ed187ac242dedcff13c56cf1b0b91;p=pdp8.git diff --git a/sw/rescue/lab8e_goettingen/disk2_11/rkb/paroff/haenya.ft b/sw/rescue/lab8e_goettingen/disk2_11/rkb/paroff/haenya.ft new file mode 100644 index 0000000..8c031c6 --- /dev/null +++ b/sw/rescue/lab8e_goettingen/disk2_11/rkb/paroff/haenya.ft @@ -0,0 +1,489 @@ +C -+-+-+-+-+  H A E N Y A . F T  -+-+-+-+-+ +C +C * UP 1-OCT-82 +C * UP 3-NOV-82 +C * UP 26-NOV-82 +C * UP 12-JAN-83 +C * UP 25-APR-83 +C * UP 26-APR-83 CREATES HAENYA FROM HAESPL +C * UP 28-APR-83 +C * UP 2-MAY-83 +C * UP 3-MAY-83 +C * UP 5-MAY-83 +C * UP 10-MAY-83 +C * UP 20-NOV-83 +C * UP 7-DEC-83 DEBUGGING +C * UP 8-DEC-83 DEBUGGING +C * UP 15-DEC-83 REARANGEMENT OF DATA CHANNELS ON OUTPUT PLOT +C * UP 19-DEC-83 DEBUGGING +C * UP 25-JAN-84 INSERTS NORMALISATION FACTORS INTO PLOT +C * UP 26-APR-84 INSERTS SPECIAL CH2 HANDLING ( TURN TABLE) +C * UP 26-JAN-84 OPTION '+' ASKS FOR CALIBRATION & NORMALISATION VECTOR +C * UP 28-MAY-84 SACCADE DIRECTION SEPARATION INTO L/R +C * UP 30-NOV-84 SOME NEW ENHANCEMENTS FOR I/O AND PLOTTER +C * UP 10-DEC-84 ONLE ONE NORMALISATION FACTOR FOR EACH CHANNEL +C * UP 27-SEP-85 INSERTS STRIGG CHANNELS +C * UP 19-NOV-85 REMOVES SOME BUGS +C * UP 2-DEC-85 CORRECT HANDLING FOR TRIGGER CHANNELS + +C PLOT THE DATA POINTS AND THE SUMM OF CH0 AND CH1 +C + SUBROUTINE NYADD + INCLUDE HAEPTI.FI + INCLUDE HAEBUF.FI + INCLUDE HAEGSA.FI + INCLUDE HAETTY.FI + INCLUDE HABRK.FI + INCLUDE HAECSZ.FI + INCLUDE HAEHAS.FI + INCLUDE HAECNY.FI + INCLUDE HAECTR.FI + INCLUDE HAPPEN.FI + REAL XYSAM,COND,NYXYSM + LOGICAL BETW + INTEGER NCOND,TOASCI,IXTA,FLD + EXTERNAL XYSAM,NCOND,COND,TOASCI,BETW,NYXYSM,IXTA,FLD +C + C + LOGICAL TIMABZ,TIMORD, @ ABZISSE / ORDINATE IS TIME SCALE + * ABISMI, @ IF THE ABZISSA SCALE EXCEEDS 5 MINUTES THEN .TRUE. ELSE .FALSE. + * BOTH,L1,ADPLT,TRGPLT @ PLOT MIN&MAX OF SOME Y POINTS, WE ARE PLOTTING A/D CHANNELS, WE ARE PLOTTING STRIG CHANNELS + INTEGER I,K,K1,J,PENPOS,MSG, + * CO, @ CHANNEL OF ORDINATE + * PLBEG,PLEND, @ START/END OF TIME SCALE + * NUMCN, @ NUMBER OF CHANNELS TO SCAN (NO MEANING IF TIMABZ IS SET .FALSE.) + * EXPAND @ IF THE TIME SCALE .LE. 5 SECONDS THEN WE MARK EACH 1/10 SECOND ELSE EACH SECOND + * ,TS1 @ HOLDS END OF TIME SCALE + * ,REARAN(8) @ REARANGE VECTOR + * ,NORMSG (1) @ MESSAGE FOR CALIBRATION FACTORS + REAL TEMSET, @ DIFFERENCE PLOTTER 0,0 AND PICTURE 0,0 (NO MEANING IF TIMABZ IS SET TO .FALSE.) + * LEFT,RIGHT,UP,DOWN, @ PLOT SCALES ABZISSE/ ORDINATE + * TE1,TE2,TE3, @ TEMPS + * X,Y,YFMIN,YFMAX, @ USED TO COMPRESS PLOT OUTPUT + * OLDX @ USED TO COMPRESS PLOT OUTPUT + * ,OMAX,OMIN @ MAXIMUM AND MINIMUM OF ORDINATE DATA + INTEGER TEXT0(4) @ DIMENSION OF THE TIME AXIS (MIN OR SEC) + * ,TEXT1(2) @ MSG FOR FORMAT 3 OUTPUT + DATA CA,CO,NUMCN,STEP,NUMTRI,ALLCHS /0,8,2,3,0,2/ + * ,TEXT1 /'A/D','TRG'/ + * ,TEXT0 /'TIME [MIN] TIME [SEC]'/ + * ,REARAN /2,5,4,1,0,3,6,7/ @ PERMUTATION: NUMBER OF PLOT ^ CHANNEL + * ,NORMSG /'GRAD'/ +C SIGNAL NAME CHANNEL NUMBER COMMENT +C CH0 A/D 0 NYSTAGMUS +C CH1 A/D 1 HEAD +C CH2 A/D 2 TABLE +C TRG 1 SPIKE 3 NEURON SPIKES +C ADD 0 4 COMPUTED FROM NYSTAGMUS+HEAD +C ADD 1 5 COMPUTED NYSTAGMUS+HEAD+TABLE +C +C +C STATEMENT FUNCTIONS ARE: +C + BOX(AP1)=AMAX1(AMIN1(AP1,UP-DOWN),-.01) + XF(IP1)=XYSAM(IP1,CA)-LEFT @ RETURNS THE VALUE OF THE ABZISSA WITH CORRECT OFFSET + YF(IP1)=NYXYSM(IP1,CO)-DOWN @ RETURNS THE VALUES OF THE ORDINATE WITH THE CORRECT OFFSET + DEFALT(I1,I2)=NCOND(I1.EQ.1H ,I2,IXTA(FLD(0,6,I1))-48) @ RETURNS DEFAULT VALUES +C + IF (SAMCNT.EQ.0) RETURN @ NO SAMPLE, NOTHING TO DO +C +C COMPUTE THE DEFAULT VALUES FOR THE PLOT +C + IF (OPTION.EQ.11) GOTO 15 @ SP- SKIPS QUESTIONS AND TAKES THE LAST VALUES USED + NUMCN=CHNLS-COSTRT+2 @ ALL CHANNELS TO PLOT + NYAFIL=0 @ NO FILTER FOR TABLE POSITION + IF (OPTION.EQ.12) GOTO 62 + CALL ASKHIM (2) @ SP* DOES NOT ASK SCILLY QUESTIONS + WRITE (TTO,1) CA,COSTRT,STEP @ READ CHANNEL ## (COSTRT SET TO 0 ALLWAYS)! + READ (TTI,2) I,J,K @ AND STEP INCREMENTS + IF BREAK(11) RETURN + CA=DEFALT(I,CA) + COSTRT=0 @ CAUSED BY REARANGE VECTOR, I DON'T KNOW HOW TO DO BETTER???? + STEP=NCOND(K,K,STEP) @ CA,CO, STEP DEFAULTS FOR BLANK INPUT +62 CONTINUE + TIMABZ=CA.GT.7 @ ABZISSE IS THE TIME SCALE + TIMORD=COSTRT.GT.7 @ ORDINATE IS THE TIME SCALE + IF (.NOT.TIMABZ) GOTO 15 @ ABZISSA IS NO TIME SCALE: ONLY TWO CHANNELS (CO,CO+1) TO PLOT + NUMCN=MIN0(NUMCN,CHNLS+1-COSTRT) @ COMPUTE CORRECT NUMBER OF CHANNELS + IF (OPTION.EQ.12) GOTO 61 + WRITE (TTO,3) TEXT1(1),NUMCN @ READ THE NUMBER OF CHANNELS TO SCAN + READ (TTI,2) I + I=DEFALT(I,NUMCN) @ FOR BLANK INPUT DEFAULT IS NUMCN + NUMCN=MIN0(CHNLS+STRIGS-COSTRT+1,NCOND(I,I,NUMCN)) @ NO INPUT MEANS CHNLS-CO + ALLCHS=NUMCN+NUMTRI @ NUMBER OF A/D AND STRIGG CHANNELS TO PLOT +C +C ASK FOR NUMBER OF STRIGG CHANNELS TO PLOT +C + IF (STRIGS.EQ.0) GOTO 11 @ NO STRIGS, NO QUESTIONS + NUMTRI=STRIGS @ INSERT DEFAULT VALUE + WRITE (TTO,3) TEXT1(2),NUMTRI @ READ THE NUMBER OF STRIGG CHANNELS TO SCAN + READ (TTI,2) I + I=DEFALT(I,NUMTRI) @ BLANK INPUTS DEFAULTS TO NUMTRI + NUMTRI=MIN0(STRIGS,MAX0(0,I)) @ INSERT CORRECT VALUE + ALLCHS=MIN0(NUMCN+NUMTRI,8) @ ONLY 8 CHANNELS ALLOWED +11 CONTINUE +61 CONTINUE +C +C ASK FOR SCALE OF THE ORDINATE (MINIMUM AND MAXIMUM) +C + DO 80 K=COSTRT+1,NUMCN+COSTRT+NUMTRI + K1=K+1 + ADPLT=K.LT.CHNLS @ QUESTIONS FOR A/D CHANNEL TO PLOT + TRGPLT=BETW(CHNLS,K,CHNLS+STRIGS) @ FOR SCHMITT TRIGGER CHNNELS + I=NCOND(ADPLT,'A/D',TRGPLT,'TRG','ADD') @ COMPUTE CORRECT TEXT + WRITE (TTO,5) I,K,ORDMAX(K1) @ READ THE MAXIMUM + READ (TTI,7) TE1 + ORDMAX(K1)=COND(TE1,TE1,ORDMAX(K1)) + IF (BREAK(11)) RETURN + WRITE (TTO,6) I,K,ORDMIN(K1) @ READ THE MINIMUM + READ (TTI,7) TE1 + ORDMIN(K1)=COND(TE1,TE1,ORDMIN(K1)) + IF (BREAK(11)) RETURN + IF (TRGPLT) GOTO 80 @ NO NORMALISATION FOR TRIGGER CHANNELS + WRITE (TTO,8) NORMSG(1),NORM(K1) @ READ THE CALIBRATION FACTOR + READ (TTI,7) TE1 + NORM(K1)=COND(TE1,TE1,NORM(K1)) @ +1 HERE SINCE REARAN'S CH # IS FORM 0 TO 7 AND NORM FROM 1 TO 8! +80 CONTINUE + C +C AND HERE WE ASK FOR TABLE POSITION HANDLING (REAL,FILTER,POLYNOM) +C NYAFIL: DETERMINE DESIRED FILTER +C NO FILTER: 0 +C 11 POINT FILTER:1 +C LINE: 2 +C 3 POINT FILTER: 3 +C +C + WRITE (TTO,104) NYAFIL,TABDIR @ WELCHES FILTER FUER TISCHPOSITION? + READ (TTI,105) NYAFIL,TABDIR + IF (NYAFIL.LE.0) GOTO 82 @ 0: KEINE BEHANDLUNG + IF (NYAFIL.EQ.2) GOTO 125 @ 2: GERADE DURCH STUETZPUNKTE + NYAFIL=MIN0(NYAFIL,3) @ 1: 11 POL FILTER, 3: 3 POL FILTER + GOTO 82 +C +125 CONTINUE @ LESE STUETZPUNKTE FUER GERADE EIN + CALL MOVE (-16,-10,SMOTAB) + DO 110 I=1,8 +111 WRITE (TTO,101) I + READ (TTI,102) SMOTAB(1,I) @ UNTERER STUETZPUNKT + IF (I.EQ.1) GOTO 112 + IF (SMOTAB(1,I).LT.SMOTAB(2,I-1)) GOTO 120 +112 IF (SMOTAB(1,I).LT.0) GOTO 111 @ FEHLERHAFTE EINGABE +113 WRITE (TTO,103) + READ (TTI,102) SMOTAB(2,I) @ OBERER STUETZPUNKT + IF (SMOTAB(2,I).LE.SMOTAB(1,I)) GOTO 113 +110 CONTINUE + GOTO 122 +120 DO 121 J=I,8 + SMOTAB(1,J)=-10 +121 SMOTAB(2,J)=-10 +122 CONTINUE + DO 130 I=1,8 + SMOTAB(1,I)=SMOTAB(1,I)*SAMRAT+.5 @ CONVERT TO INDICES AND ROUND + SMOTAB(2,I)=SMOTAB(2,I)*SAMRAT+.5 + IF (SMOTAB(1,I).LT.0) GOTO 140 + SMOTBL(1,I)=(XYSAM(INT(SMOTAB(2,I)),2)- + * XYSAM(INT(SMOTAB(1,I)),2)) + * / (SMOTAB(2,I)-SMOTAB(1,I)) @ STEIGUNG +130 SMOTBL(2,I)=XYSAM(INT(SMOTAB(1,I)),2)- + * SMOTBL(1,I)*(SMOTAB(1,I)-.5) +140 CONTINUE +101 FORMAT (' STUETZPUNKT PAAR NR.:',I2,/ + * ' UNTERER STUETZPUNKT IN SEC (F6.0)? ',$) +102 FORMAT (F6.0) +103 FORMAT (' OBERER STUETZPUNKT ? ',$) +104 FORMAT (' TISCH: FILTER 0/1/2/3, DREHRICHTUNG L/R/* (I1,A1) ', + * I1,A1,' ?',$) +105 FORMAT (I1,A1) +C +C +82 CONTINUE +15 CONTINUE + C +C COMPUTE THE SCALING LIMITS +C + S0=BEGIN*SAMRAT @ HERE WE START THE OUTPUT PLOT + S1=ENDS*SAMRAT @ AND HERE WE END THE PLOT + TS1=S1 @ END OF TIME SCALE +C PLBEG HOLDS THE FIRST SECOND OF THE TIME SCALE (LEFTMOST ) +C PLEND HOLDS THE LAST SECOND OF THE TIME SCALE ( RIGHTMOST DOT) +C TS1 HOLDS FACTOR FOR TIME SCALE (PLEND*SAMRAT) +C S0 HOLDS THE STARTING INDEX OF DATA TO PLOT ( PLBEG*SAMRAT) +C S1 HOLDS THE INDEX OF LAST DATA POINT (MIN0(SAMCNT,TS1)) + PLBEG=BEGIN + PLEND=ENDS +C + IF (S1.LE.SAMCNT) GOTO 17 @ IF THE INTERVALL TO PLOT EXCEEDS THE + TE1=SAMCNT-1/SAMRAT @ DATA POINTS THEN + TE2=S1/SAMRAT + WRITE (TTO,4) TE1,TE2 @ WE PRINT THE OVERFLOW MSG AND + S1=SAMCNT @ INSERT THE CORRECT PLEND VALUE +17 CONTINUE @ INTERVALL OK, COMPUTE UPPER AND LOWER MARGIN OF PLOT DATA + RIGHT=COND(TIMABZ,FLOAT(TS1),512.) @ ABZISSE MARGIN EITHER TIME OR 10 BIT CONVERTER + LEFT =COND(TIMABZ,FLOAT(S0),-512.) + UP =COND(TIMORD,FLOAT(TS1),512.) @ ORDINATE MARGIN + DOWN =COND(TIMORD,FLOAT(S0),-512.) + IF BREAK(11) RETURN @ USER GETS RID OF PRGRM + CALL STPLT @ START THE PLOTTER ( CALL PLOTS) + CALL XYPLOT (XOFSET,YOFSET,-PENUP) @ DRIVE PEN TO PICTURE ZERO +C +C MAIN PLOT LOOP +C +C LOOP TO PLOT EACH CHANNEL WITH CORRECT OFFSET, FACTOR AND SCALE +C TEMSET HOLDS THE ORIGIN FOR EACH CHANNEL +C TEMSET=0 PLOT ONE CHANNEL/TIME OR X/Y +C TEMSET=YLEN/NUMCN FOR MORE THAN ONE CHANNEL +C THE TEMSET ORIGIN OFFSET IS RESET AT THE END OF THE LOOP ( ST. 25) +C +C + TEMSET=0 + DO 20 K=1,NCOND(TIMABZ,ALLCHS,1) @ SCAN NUMCN CHANNELS WITH INDEX K FROM 1 TO NUMCN + CO=REARAN(K) @ REARANGEMENT OF GRAPHICS + K1=CO+1 @ INDEX TO ORDMIN/MAX SINCE CO STARTS WITH 0 +CXCX IF (CO.LT.COSTRT) GOTO 20 @ COMPUTE ONLY DESIRED CHNNELS +C +C HERE WE COMPUTE THE ACTUAL CHANNEL NUMBER: +C EITHER A/D CO TO NUMCN-1 (0 TO CHNLS-1) OR +C TRIGG EVENT FROM CHNLS TO CHNLS+STRIGS-1 +C +C ADPLT IS SET TRUE FOR A/D CHANNELS TO SCAN +C STEPSIZE IS THE INCREMENT FOR THE SCAN; +C FOR A/D CHANNELS THE MAXIMUM VALUE OF THE SPECIFIED: STEP OR +C IF LOWER SAMPLE FREQUENCY THAN MAXRAT THE PACK-VALUE. +C FOR TRIGG EVENT CHANNELS STEP MUST ALLWAYS BE 1! +C + TRGPLT=BETW(CHNLS,CO,CHNLS+STRIGS) @ WE ARE PLOTTING SCHMITT TRIGGER CHANNELS + ADPLT=CO.LT.CHNLS @ WE ARE PLOTTING A/D CHANNELS + STEPS=NCOND(ADPLT,MAX0(STEP,PACK(CO+1)), + * TRGPLT,1,MAX0(STEP,PACK(1),PACK(2),PACK(3))) @ COMPUTE STEPSIZE FOR SCAN +C +C HERE WE DECIDE : X OVER Y PLOT +C OR X OVER TIME +C + IF ( .NOT.TIMABZ) GOTO 51 +C +C X OVER TIME +C +C IF THERE IS MORE THAN ONE CHANNEL WE HAVE TO PREPARE A LOT OF SMALL +C Y - AXIS FOR EACH CHANNEL TO SCAN. SO WE SET A TEMPORARY Y OFFSET CALLED TEMSET + TEMSET=FLOAT(ALLCHS-K)*(YLEN-ALLCHS*.5+.5)/ + * FLOAT(ALLCHS)+(ALLCHS-K)*.5 @ TEMPORARY Y OFFSET + CALL XYPLOT (0,TEMSET,-PENUP) @ SWITCH TO CORRECT PART OF PICTURE +C +C COMPUTE THE FACTOR TO GET A PRETTY PLOT +C + UP=ORDMAX(K1) + DOWN=ORDMIN(K1) @ INSERT THE DEFAULT SCALING LIMITS + IF (UP.LE.DOWN) GOTO 51 @ AUTOMATIC SCALING DESIRED + OMAX=UP + OMIN=DOWN + GOTO 51 @ SCALING ALREADY DONE + C HERE WE DO AUTOMATIC SCALING +C +21 CONTINUE + OMAX=-99999. + OMIN=-OMAX + DO 50 J=S0,S1,STEPS @ FETCH MIN & MAX OF THE DATA + IF (BREAK(11)) GOTO 25 @ IF THE USER GETS RID OF PLOT THEN WE EXIT VIA 25 + TE1=NYXYSM(J,CO) @ DATA POINT --> TE1 + OMAX=AMAX1(OMAX,TE1) @ MAXIMUM --> OMAX +50 OMIN=AMIN1(OMIN,TE1) @ MINIMUM --> OMIN +C + IF (UP .GT. DOWN-1.) GOTO 52 @ IF OMAX.EQ.OMIN THEN WE MAKE A FIXED SCALE +C FIXED SCALE IF POSSIBLE: + IF (OMIN.GE.OMAX-(DOWN-UP)) GOTO 53 @ MAKE A FIXED SCALE AT 53 + DOWN=OMIN @ DIFFERENCE BETWEEN MIN AND MAX EXCEEDS THE FIX SCALE LENGTH, SO WE TAKE DEFAULT VALUES + UP=OMAX @ ELSE DOWN=UP-(OMIN-OMAX) (FIXED DIFFERENCE SCALE) + GOTO 51 @ WORK DONE +C FIXED SCALE POSSIBLE, MAKE IT AND CENTER THE PLOT +53 TE1=OMIN-((DOWN-UP)-(OMAX-OMIN))*.5 + UP=TE1+DOWN-UP @ UPPER LIMIT + DOWN=TE1 @ WE SHIFT THE LINE INTO THE CENTER + GOTO 51 @ FIXED DIFFERENCE SCALE DONE +C AUTOMATIC SCALE BUT WITH INCREMENTAL STEPPED LIMITS +C +52 CONTINUE + DOWN=COND(OMIN.LT.-500.,OMIN,OMIN.LT.-250.,-500., + * OMIN.LT.0,-250.,OMIN.LT.250.,0, 250.) @ MAKE 250 STEPS INCREMENT + UP=COND(OMAX.GT.500.,OMAX,OMAX.GT.250.,500., + * OMAX.GT.0,250.,OMAX.GT.-250.,0, -250.) @ 250 STEPS INCREMENT +C + UP=COND(DOWN.GE.UP,DOWN+255.,UP) @ SOMETIMES IT WILL HAPPEN: UP=DOWN +51 CONTINUE @ HERE WE COMPUTE THE FACTOR FOR ONE LINE + WRITE(0,4711) RIGHT,LEFT,ALLCHS,UP,DOWN +4711 FORMAT(5I7) + IF (TIMABZ) CALL FACTOR (XLEN/(RIGHT-LEFT), + * (YLEN-ALLCHS*.5+.5)/FLOAT(ALLCHS)/(UP-DOWN)) + IF (.NOT.TIMABZ) CALL FACTOR (AMIN1(XLEN,YLEN)/(RIGHT-LEFT), + * AMIN1(XLEN,YLEN)/(UP-DOWN)) + IF BREAK(11) GOTO 25 @ USER GETS RID OF PRGRM + C HERE WE COMPRESS ( TO SAVE PLOTTER TIME) AND PLOT THE DATA POINTS +C + CALL WHERE (TE1,TE1,TE1,TE2) @ FETCH THE X FACTOR --> TE1 + YFMIN=9999. + YFMAX=-YFMIN @ PREP TEMPS TO FETCH THE MIN/MAX OF THE ORDINATE POINTS + BOTH=.FALSE. +C + OLDX=XF(S0)-0.4/TE1*XLEN + PENPOS=PENUP + L1=.FALSE. @ L1 HOLDS: YF(I) INSIDE LOWER/UPPER LIMIT OF PLOT +C +C IF THERE IS A X MOVEMENT LESS THAN .3 MM THEN WE DO NOT PLOT THE +C POINT. WE FETCH THE MIN/MAX OF Y AND PLOT A LINE FROM MIN TO MAX AT THIS X LOCATION. +C THIS IS TO COMPRESS USELESS OUTPUT TO THE PLOTTER. +C + DO 10 I=S0,S1,STEPS @ HERE WE PLOT THE DATA POINTS ONE BY ONE + X=XF(I) @ TEMP TO SAVE COMPUTING TIME +75 Y=YF(I) @ TEMP + IF BREAK(11) GOTO 25 @ USER GETS RID OF PROGRAM + IF(.NOT.TRGPLT) + * PENPOS=NCOND(BETW(-.1,Y,UP-DOWN).OR.L1,PENPOS,PENUP) @ WE DO NOT PLOT ANY POINTS EXCEEDING THE UPPER/LOWER LIMIT + IF (ABS(OLDX-X)*TE1.GE.0.35/XLEN) GOTO 71 @ PLOT THIS POINT + IF(TRGPLT) GOTO 74 @ TRIGGER EVENTS WITHOUT ANY COMPRESSION! + YFMIN=AMIN1(Y,YFMIN) @ NO X MOVEMENT, GET THE MIN Y VALUE + YFMAX=AMAX1(Y,YFMAX) @ FETCH THE MAX Y VALUE + BOTH=.TRUE. @ TWO Y POINTS TO PLOT!! (MIN & MAX) + GOTO 10 @ WITH THE NEXT REAL X MOVEMENT +C +74 CALL XYPLOT(X,BOX(-DOWN),PENPOS) @ DELTA SPIKE FOR TRIGGER CHANNEL +71 IF (BOTH) GOTO 72 @ IF THERE ARE TWO POINTS TO PLOT THEN GO 72 ELSE + CALL XYPLOT (X,BOX(Y),PENPOS) @ WE PLOT THIS ONE AND PRESS IT INTO U/L LIMITS + IF(TRGPLT) + * CALL XYPLOT(X,BOX(-DOWN),PENPOS) @ DELTA SPIKE FOR TRIGGER CHANNEL + GO TO 73 +72 CALL XYPLOT (X,BOX(YFMIN),PENPOS) @ UPPER AND LOWER Y VALUE TO PLOT + CALL XYPLOT (X,BOX(YFMAX),PENPOS) + YFMIN=9999. @ SET UP NEW MIN/MAX FETCH CYCLE + YFMAX=-YFMIN + BOTH =.FALSE. @ CLEAR UPPER/LOWER FLAG +73 CONTINUE + OLDX=X @ SAVE LAST X VALUE + L1=BETW(-.1,Y,UP-DOWN) @ L1:= Y INSIDE U/L LIMITS! + PENPOS=PENDWN + IF (NXTBIT.GE.0) GOTO 10 @ NO MORE SPIKES IN THIS TRIGG TIME WORD, SKIP TO 10 + X=X+(XF(I+1)-XF(I))/TTRGW @ COMPUTE NEW X FOR MULTIPLE SPIKE EVENTS + GOTO 75 @ REPEAT FOR NEXT SPIKE IN THIS TIME WORD +10 CONTINUE + IF (.NOT. TIMABZ) GOTO 25 @ WORK DONE FOR X-Y GRAPHICS +C +C MAKE A SCALE FOR THE ORDINATE ( INPUT VOLTAGE ) +C + CALL WHERE (TE1,TE1,TE1,TE2) +C HERE WE WRITE THE CALIBRATION AND NORMALISATION FACTOR ONTO THE RIGHT SIDE + CALL FACTOR (1.,1.) @ RESET TO THE NORMAL FACTOR FOR SYMBOL + CALL NUMBER (XLEN+YN,TE2*(UP-DOWN)/2.-1.2*YN + * ,YN*.8,NORM(CO+1),90.,2) +C + CALL FACTOR (1.,TE2) + CALL XYPLOT (-.5,0,-PENUP) + PENPOS=PENUP + DO 55 I=0,10 + TE1=DOWN+I*(UP-DOWN)/10. + CALL XYPLOT (0,TE1-DOWN,PENPOS) + PENPOS=PENDWN + J=MOD(I,5) + TE3=COND(J.EQ.0,-.5,-.3) + CALL XYPLOT (TE3,TE1-DOWN,PENDWN) + CALL XYPLOT (0,TE1-DOWN,PENDWN) + IF (MOD(I,10).NE.0) GOTO 55 + CALL FACTOR (1.,1.) @ WRITE LOWER AND UPPER LIMIT ONTO THE + CALL NUMBER (-XOFSET+2.*YZ,(TE1-DOWN)*TE2-YN*.38,YN*.75,TE1,0, + * -1) @ AXIS + CALL FACTOR (1.,TE2) + CALL XYPLOT (0,TE1-DOWN,PENUP) +55 CONTINUE + CALL XYPLOT (.5,0,-PENUP) + CALL FACTOR (1.,1.) + CALL SYMBOL (-XOFSET+YZ, + * AMAX1(0.,((YLEN-ALLCHS*.5)/ALLCHS-4.*YZ))*.5,YZ, + * NCOND(ADPLT,'CH ',TRGPLT,'TRG','ADD'),90.,3) + CALL NUMBER (999.,999.,YZ, + * NCOND(ADPLT,CO,TRGPLT,CO-CHNLS+1,CO-STRIGS-CHNLS),90.,-1) @ AD CHANNEL # : 0 .. CHNLS, TRG CHAN#: 1 .. STRIGS + IF(.NOT.TRGPLT) + * CALL SYMBOL (-XOFSET+YZ+1.2*YN,YLEN/NUMCN*.33-YZ,YN*.8 + * ,'[DEG]',90.,5) @ A/D CHANNELS, SO WE PLOT DEGREES +C +C + CO=CO+1 @ HERE WE INCREMENT THE CHANNEL NUMBER +25 CALL FACTOR (1.,1.) + CALL XYPLOT (0.,-TEMSET,-PENUP) @ RESET CORRECT ORIGIN OF DIAGRAM + IF BREAK(11) GOTO 40 +20 CONTINUE + C +C IF ABZISSA IS TIME SCALE THEN DRAW A TIME SCALE +C + IF (.NOT.TIMABZ) GOTO 40 + CALL FACTOR (XLEN/(PLEND-PLBEG),1.) + CALL XYPLOT (0.,-.5,-PENUP) +C HERE WE PLOT THE LITTLE BAR INDICATING A SECOND OR MINUTE + ABISMI=PLEND-PLBEG.GT.300 @ ABZISSE EXCEEDS 5 MINUTES SO WE DRAW A MINUTE SCALING + EXPAND=NCOND(PLEND-PLBEG.LE.5,10,1) @ IF THE TIME SCALE IS .LE. 5 SECONDS THEN WE MARK EACH 1/10 SEC +C @ REMEMBER : IF (ABISMI) EXPAND=.FALSE. FOR CORRECT ACTION + DO 30 I=PLBEG*EXPAND,PLEND*EXPAND @ 1 SEC INCREMENTS + TE1=I-PLBEG*EXPAND + CALL XYPLOT (TE1/EXPAND,0,PENDWN) + IF (.NOT.ABISMI) TE2=COND(MOD(I,10).EQ.0,-.7, MOD(I,5).EQ.0,-.5, + * -.3) + IF (ABISMI) TE2=COND(MOD(I,600).EQ.0,-.7, + * MOD(I,300).EQ.0,-.5, MOD(I,60).EQ.0,-.3) + CALL XYPLOT (TE1/EXPAND,TE2,PENDWN) @ HERE WE PLOT THE SMALL BAR + CALL XYPLOT (TE1/EXPAND,0,PENDWN) +C HERE WE COMPUTE WHERE TO WRITE A NUMBER ONTO THE TIME AXIS + TE2=PLEND-PLBEG @ COMPUTE THE SIZE OF THE TIME AXIS + IF (.NOT.ABISMI) GOTO 31 @ IF WE PLOT A SECONDS AXIS THEN GOTO 31 + IF (MOD(I, + * NCOND(TE2.GT.5400,1800,TE2.GE.1800,600,TE2.GT.900,300, + * TE2.GT.300,120,60)).GT.0) GOTO 30 @ COMPUTE THE MINUTE WHERE TO PLOT THE NUMBER + GOTO 32 +31 IF (MOD(I,NCOND(TE2.GT.180,30,TE2.GT.60,20,TE2.GT.30,10, + * TE2.GT.10,5,TE2.GT.5,2,10)).GT.0) @ WE ASSUME: TE2.LE.5 THEN EXPAND:=10 ELSE EXPAND:=1 !!! + * GOTO 30 +C +C HERE WE PLOT THE NUMBER ONTO THE AXIS +C +32 CONTINUE + CALL WHERE (TE2,TE2,TE2,TE3) @ GET THE X FACTOR + CALL FACTOR (1.,1.) @ RESET TO THE NORMAL FACTOR FOR SYMBOL + CALL NUMBER (TE1/EXPAND*TE2-INT(ALOG10(AMAX0(COND(ABISMI, + * I/60.,I/EXPAND),1))+.01) + * *YN/2.-YN*.3, + * -1.3,YN,COND(ABISMI,I/60.,I/EXPAND),0,-1) + CALL FACTOR (TE2,TE3) + CALL XYPLOT (TE1/EXPAND,0,PENUP) + CALL XYPLOT (TE1/EXPAND,0,PENDWN) + I=NCOND(ABISMI,MIN0(I+59,PLEND-1),I) @ MINUTE SCALE HAS 60 SEC INCREMENT + +30 CONTINUE + CALL FACTOR (1.,1.) + CALL XYPLOT (0,.5,-PENUP) + CALL SYMBOL (XLEN/1.5,-YOFSET,YZ,TEXT0(NCOND(ABISMI,1,3)),0,10) +C +C +40 CALL FACTOR (1.,1.) + IF (BREAK(11)) GOTO 60 @ USER SWITCHED OFF THE PLOTTER + CALL LABPLT (0,-YOFSET,PLEND) + IF (BREAK(11)) GOTO 60 + CALL SYMBOL (YZ,YLEN+YZ,YZ,LABEL,0,MIN0(42,LABCNT*6)) @ PLOT 42 CHARACTERS INTO THE FIRST LINE OF LABEL + IF (LABCNT.GT.7) + * CALL SYMBOL (YZ,YLEN,YZ,LABEL(8),0,LABCNT*6-42) @ AND THE REMAINING INTO THE NEXT LINE +60 CONTINUE + CALL XYPLOT (38.,25.,PENUP) + CALL EXPLT + RETURN +1 FORMAT (' KANAL ABZISSE/ORDINATE (ZEITBASIS=8), SCHRITTWEITE', + * ' (2I1,I3) (', + * 2I1,I3,')? ',$) +2 FORMAT (2A1,I3) +3 FORMAT (' ANZAHL DER DARZUSTELLENDEN ',A3,' KANAELE (I1) (',I1, + * ') ? ',$) +4 FORMAT (' NUR FUER',F6.1,' SEK. MESSDATEN JEDOCH',F6.1, + * ' SEC. ANGEFORDERT. ') +5 FORMAT (' ORDINATEN MAXIMUM FUER ',A3,' KANAL', + * I2,' (',F5.0,')?',$) +6 FORMAT (' " " MINIMUM " ',A3,' " " ', + * I2,' (',F5.0,')?',$) +7 FORMAT (F5.0) +8 FORMAT (' FAKTOR FUER ',A5,14X,'(',F5.3,')?',$) +9 FORMAT (' ABZISSEN/ ORDINATEN KANAL ODER SCHRITTWEITE ', + * ' FEHLERHAFT: ',3I4) + END + \ No newline at end of file