A large commit.
[pdp8.git] / sw / rescue / lab8e_goettingen / disk2_11 / rkb / paroff / haespl.ft
diff --git a/sw/rescue/lab8e_goettingen/disk2_11/rkb/paroff/haespl.ft b/sw/rescue/lab8e_goettingen/disk2_11/rkb/paroff/haespl.ft
new file mode 100644 (file)
index 0000000..47bdd09
--- /dev/null
@@ -0,0 +1,411 @@
+C      -+-+-+-+-+ \ e  H A E S P L . F T  \ e -+-+-+-+-+\r
+C\r
+C      * UP  1-OCT-82\r
+C      * UP  3-NOV-82\r
+C      * UP 26-NOV-82\r
+C      * UP 12-JAN-83\r
+C      * UP 25-APR-83\r
+C      * UP 26-APR-83\r
+C      * UP 28-APR-83\r
+C      * UP  5-MAY-83\r
+C      * UP 10-MAY-83\r
+C      * UP 11-NOV-83 TIME SCALE OF PLOT MAY EXCEED SAMPLE TIME \r
+C      * UP  8-DEC-83 DEBUGGING\r
+C      * UP 30-NOV-84 SOME NEW ENHANCEMENTS FOR I/O AND PLOTTER\r
+C      * UP 10-DEC-84\r
+C      * UP  3-JUL-85 OUTPUT SCHMITT TRIGGER CHANNELS\r
+C      * UP 15-AUG-85 \r
+C      * UP 22-AUG-85 REMOVES SOME BUGS CONCERNING TRIGG OUTPUT \r
+C      * UP  3-SEP-85 IF SW 8 UP: INSERTS TRIGG FREQUENCY PLOT\r
+C      * UP  1-NOV-85 DEBUGGS ASZISSA SCALE ROUTINE\r
+C      * UP 19-NOV-85 REMOVES SOME BUGS\r
+C\r
+C      PLOT THE DATA POINTS\r
+C\r
+       SUBROUTINE SPLOT\r
+       INCLUDE HAEPTI.FI\r
+       INCLUDE HAEBUF.FI\r
+       INCLUDE HAEGSA.FI\r
+       INCLUDE HAETTY.FI\r
+       INCLUDE HABRK.FI\r
+       INCLUDE HAECSZ.FI\r
+       INCLUDE HAEHAS.FI\r
+       INCLUDE HAECSP.FI\r
+       INCLUDE HAECTR.FI\r
+       INCLUDE HAPPEN.FI\r
+       REAL XYSAM,COND\r
+       LOGICAL BETW\r
+       INTEGER NCOND,TOASCI,IXTA,FLD\r
+       EXTERNAL XYSAM,NCOND,COND,TOASCI,IXTA,FLD,BETW\r
+C\r
+\fC\r
+       REAL ORDMIN(8),ORDMAX(8)\r
+       EQUIVALENCE (SRDMIN(1),ORDMIN(1)),(SRDMAX(1),ORDMAX(1)) @ EQUIVALENCE INTO THE HAECSP COMMON BLOCK\r
+       LOGICAL TIMABZ,TIMORD,          @ ABZISSE / ORDINATE IS TIME SCALE\r
+     *  ABISMI,                                @ IF THE ABZISSA SCALE EXCEEDS 5 MINUTES THEN .TRUE. ELSE .FALSE.\r
+     *  BOTH,L1,ADPLT                          @ PLOT MIN&MAX OF SOME Y POINTS, WE ARE PLOTTING A/D CHANNELS\r
+       INTEGER I,K,K1,J,PENPOS,MSG,    @ MSG HOLDS TEMP MESSAGE FOR FORMAT 5 AND 6\r
+     *  STEP,STEPS,                            @ INCREMENT IN SAMBUF SCAN\r
+     *  CA,CO,                         @ CHANNEL OF ABZISSE, ORDINATE\r
+     *  PLBEG,PLEND,                   @ START/END OF TIME SCALE\r
+     *  NUMCN,                         @ NUMBER OF CHANNELS TO SCAN (NO MEANING IF TIMABZ IS SET .FALSE.)\r
+     *   NUMTRI,ALLCHS,                        @ NUMBER OF TRIGGER INPUTS TO PLOTT, TOTAL # OF CHANNELS TO SCAN (A/D +STRIGG)\r
+     *  EXPAND                         @ IF THE TIME SCALE .LE. 5 SECONDS THEN WE MARK EACH 1/10 SECOND ELSE EACH SECOND\r
+     *  ,TS1                           @ HOLDS PLEND*SAMRAT FOR SCALEING THE ABZISSA\r
+       REAL TEMSET,                    @ DIFFERENCE PLOTTER 0,0 AND PICTURE 0,0 (NO MEANING IF TIMABZ IS SET TO .FALSE.)\r
+     *  LEFT,RIGHT,UP,DOWN,            @ PLOT SCALES ABZISSE/ ORDINATE\r
+     *  TE1,TE2,TE3,                   @ TEMPS\r
+     *  X,Y,YFMIN,YFMAX,               @ USED TO COMPRESS PLOT OUTPUT\r
+     *  OLDX                           @ USED TO COMPRESS PLOT OUTPUT\r
+     *   ,OMAX,OMIN                    @ MAXIMUM AND MINIMUM OF ORDINATE DATA\r
+       INTEGER TEXT0(4)                @ DIMENSION OF THE TIME AXIS (MIN OR SEC)\r
+     *  ,TEXT1(2)                      @ MSG FOR FORMAT 3 OUTPUT\r
+       DATA CA,CO,NUMCN,STEP,NUMTRI,ALLCHS /0,8,2,3,0,2/,\r
+     *  TEXT0 /'TIME [MIN]  TIME [SEC]'/\r
+     *  ,TEXT1 /'A/D','TRG'/\r
+C\r
+C      STATEMENT FUNCTIONS ARE:\r
+C\r
+       BOX(AP1)=AMAX1(AMIN1(AP1,UP-DOWN),-.01)\r
+       XF(IP1)=XYSAM(IP1,CA)-LEFT      @ RETURNS THE VALUE OF THE ABZISSA WITH CORRECT OFFSET\r
+       YF(IP1)=XYSAM(IP1,CO)-DOWN      @ RETURNS THE VALUES OF THE ORDINATE WITH THE CORRECT OFFSET\r
+       DEFALT(I1,I2)=NCOND(I1.EQ.1H ,I2,IXTA(FLD(0,6,I1))-48) @ RETURNS DEFAULT VALUES\r
+\fC\r
+       IF (SAMCNT.EQ.0) RETURN         @ NO SAMPLE, NOTHING TO DO\r
+C      \r
+C      COMPUTE THE DEFAULT VALUES FOR THE PLOT\r
+C\r
+       IF (OPTION.EQ.11) GOTO 15       @ SP- SKIPS QUESTIONS AND TAKES THE LAST VALUES USED\r
+       IF (OPTION.NE.12) CALL ASKHIM (2) @ SP* DOES NOT ASK SCILLY QUESTIONS\r
+14     CO=0                            @ THE ORDINATE IS CHANNEL 0\r
+       CA=8                            @ THE ABZISSA IS THE TIME SCALE\r
+       NUMCN=CHNLS-CO                  @ ALL CHANNELS TO PLOT\r
+       STEP=1                          @ WE PLOT EACH POINT OF THE DATA\r
+       WRITE (TTO,1) CA,CO,STEP        @ READ CHANNEL ##\r
+       READ (TTI,2) I,J,K              @   AND STEP INCREMENTS\r
+       IF BREAK(11) RETURN\r
+       CA=DEFALT(I,CA)\r
+       CO=DEFALT(J,CO)\r
+       STEP=NCOND(K,K,STEP)            @ CA,CO, STEP DEFAULT VALUES FOR BLANK INPUT\r
+       IF (.NOT.(CA.LT.0 .OR. CO.LT.0 .OR. STEP.LE.0)) GOTO 12\r
+13     WRITE (TTO,9)CA,CO,STEP\r
+       GOTO 14\r
+12     IF (CA.GE.CHNLS .AND. CA.LT.8) GOTO 13 @ TEST FOR VALID CHANNEL NUMBER\r
+       IF (CO.GE.CHNLS .AND. CO.LT.8) GOTO 13\r
+C\r
+       TIMABZ=CA.GT.7                  @ ABZISSE IS THE TIME SCALE\r
+       TIMORD=CO.GT.7                  @ ORDINATE IS THE TIME SCALE\r
+       IF (.NOT.TIMABZ) GOTO 15        @ ABZISSA IS NO TIME SCALE: ONLY TWO CHANNELS (CO,CO+1) TO PLOT\r
+       NUMCN=MIN0(NUMCN,CHNLS-CO)      @ COMPUTE CORRECT NUMBER OF CHANNELS\r
+       WRITE (TTO,3) TEXT1(1),NUMCN    @ READ THE NUMBER OF CHANNELS TO SCAN\r
+       READ (TTI,2) I\r
+       CA=DEFALT(I,CA)\r
+       CO=DEFALT(J,CO)\r
+       I=DEFALT(I,NUMCN)               @ FOR BLANK INPUT DEFAULT IS NUMCN\r
+       NUMCN=MIN0(CHNLS-CO,NCOND(I,I,NUMCN)) @ NO INPUT MEANS CHNLS-CO\r
+       NUMTRI=STRIGS\r
+       ALLCHS=NUMCN+NUMTRI             @ NUMBER OF A/D AND STRIGG CHANNELS TO PLOT\r
+C\r
+C      ASK FOR NUMBER OF STRIGG CHANNELS TO PLOT\r
+C\r
+       IF (STRIGS.EQ.0) GOTO 11        @ NO STRIGS, NO QUESTIONS\r
+       NUMTRI=STRIGS                   @ INSERT DEFAULT VALUE\r
+       WRITE (TTO,3) TEXT1(2),NUMTRI   @ READ THE NUMBER OF STRIGG CHANNELS TO SCAN\r
+       READ (TTI,2) I\r
+       I=DEFALT(I,NUMTRI)              @ BLANK INPUTS DEFAULTS TO NUMTRI\r
+       NUMTRI=MIN0(STRIGS,MAX0(0,I))   @ INSERT CORRECT VALUE\r
+       ALLCHS=MIN0(NUMCN+NUMTRI,7)     @ ONLY 8 CHANNELS ALLOWED\r
+11     CONTINUE\r
+C\r
+C      ASK FOR SCALE OF THE ORDINATE (MINIMUM AND MAXIMUM)\r
+C\r
+       K=CO-1\r
+       DO 80 I=1,ALLCHS\r
+       K=K+1\r
+       MSG=NCOND(K.GE.CO+NUMCN,TEXT1(2),TEXT1(1)) @ EITHER A/D OR TRG MESSAGE TO PRINT\r
+       J=NCOND(K.GE.CO+NUMCN,1+K-CO-NUMCN,K) @ NUMBER OF A/D OR TRIGGER CHANNEL\r
+       WRITE (TTO,5) MSG,J,ORDMAX(K+1) @ READ THE MAXIMUM\r
+       READ  (TTI,7) TE1\r
+       ORDMAX(K+1)=COND(TE1,TE1,ORDMAX(K+1))\r
+       IF (BREAK(11)) RETURN\r
+       WRITE (TTO,6) MSG,J,ORDMIN(K+1) @ READ THE MINIMUM\r
+       READ  (TTI,7) TE1\r
+       ORDMIN(K+1)=COND(TE1,TE1,ORDMIN(K+1))\r
+       IF (BREAK(11)) RETURN\r
+80     CONTINUE\r
+15     CONTINUE\r
+\fC\r
+C      COMPUTE THE SCALING LIMITS\r
+C\r
+       S0=BEGIN*SAMRAT                 @ HERE WE START THE OUTPUT PLOT\r
+       S1=ENDS*SAMRAT                  @ AND HERE WE END THE PLOT\r
+       TS1=S1\r
+C      PLBEG   HOLDS THE FIRST SECOND OF THE TIME SCALE (LEFTMOST )\r
+C      PLEND   HOLDS THE LAST SECOND OF THE TIME SCALE ( RIGHTMOST DOT)\r
+C      TS1     HOLDS FACTOR FOR TIME SCALE (PLEND*SAMRAT)\r
+C      S0      HOLDS THE STARTING INDEX OF DATA TO PLOT ( PLBEG*SAMRAT)\r
+C      S1      HOLDS THE INDEX OF LAST DATA POINT (MIN0(SAMCNT,TS1))\r
+       PLBEG=BEGIN\r
+       PLEND=ENDS\r
+C\r
+       IF (S1.LE.SAMCNT) GOTO 17       @ IF THE INTERVALL TO PLOT EXCEEDS THE\r
+       TE1=SAMCNT/SAMRAT               @  DATA POINTS THEN\r
+       TE2=S1/SAMRAT\r
+       WRITE (TTO,4) TE1,TE2           @  WE PRINT THE OVERFLOW MSG AND\r
+       S1=SAMCNT-1                     @ TAKE CARE FOR TRUNCATION ERRORS\r
+17     CONTINUE                        @ INTERVALL OK, COMPUTE UPPER AND LOWER MARGIN OF PLOT DATA\r
+       RIGHT=COND(TIMABZ,FLOAT(TS1),512.) @ ABZISSE MARGIN EITHER TIME OR 10 BIT CONVERTER\r
+       LEFT =COND(TIMABZ,FLOAT(S0),-512.)\r
+       UP   =COND(TIMORD,FLOAT(TS1),512.) @ ORDINATE MARGIN\r
+       DOWN =COND(TIMORD,FLOAT(S0),-512.)\r
+       IF BREAK(11) RETURN             @ USER GETS RID OF PRGRM\r
+       CALL STPLT                      @ START THE PLOTTER ( CALL PLOTS)\r
+       CALL XYPLOT (XOFSET,YOFSET,-PENUP) @ DRIVE PEN TO PICTURE ZERO\r
+C\r
+C      LOOP TO PLOT EACH CHANNEL WITH CORRECT OFFSET, FACTOR AND SCALE\r
+C      TEMSET HOLDS THE ORIGIN FOR EACH CHANNEL\r
+C      TEMSET=0  PLOT ONE CHANNEL/TIME OR X/Y\r
+C      TEMSET=YLEN/(NUMCN+NUMTRI) FOR MORE THAN ONE CHANNEL\r
+C      THE TEMSET ORIGIN OFFSET IS RESET AT THE END OF THE LOOP ( ST. 25)\r
+C\r
+C\r
+       TEMSET=0                        @ CO HOLDS THE FIRST CHANNEL TO PLOT (INCREMENTED AT END OF LOOP)\r
+       DO 20 K=1,NCOND(TIMABZ,ALLCHS,1)        @ SCAN ALLCHS CHANNELS WITH INDEX K FROM 1 TO ALLCHS\r
+C\r
+C      HERE WE COMPUTE THE ACTUAL CHANNEL NUMBER:\r
+C      EITHER A/D  CO TO NUMCN-1  (0 TO CHNLS-1)  OR\r
+C       TRIGG EVENT  FROM CHNLS TO CHNLS+STRIGS\r
+C\r
+C      ADPLT IS SET TRUE FOR A/D CHANNELS TO SCAN\r
+C      STEPSIZE IS THE INCREMENT FOR THE SCAN;\r
+C      FOR A/D CHANNELS THE MAXIMUM VALUE OF THE SPECIFIED: STEP OR\r
+C      IF LOWER SAMPLE FREQUENCY THAN MAXRAT THE PACK-VALUE.\r
+C      FOR TRIGG EVENT CHANNELS STEP MUST ALLWAYS BE 1!\r
+C\r
+       IF (K.GT.NUMCN)CO=CHNLS+K-NUMCN-1 @ CO RUNS FROM CO TO CO+NUMCN FOR AD, FOR TRIGG'S FROM CHNLS TO CHNLS+TRIGS\r
+       ADPLT=CO.LT.CHNLS               @ WE ARE PLOTTING A/D CHANNELS\r
+       STEPS=NCOND(ADPLT,MAX0(STEP,PACK(CO+1)),1) @ COMPUTE STEPSIZE FOR SCAN\r
+C\r
+C      HERE WE DECIDE : X OVER Y PLOT \r
+C           OR          X OVER TIME\r
+C\r
+       IF (.NOT. TIMABZ) GOTO 51\r
+C\r
+C      X OVER TIME\r
+C      \r
+C      IF THERE IS MORE THAN ONE CHANNEL WE HAVE TO PREPARE A LOT OF SMALL\r
+C      Y - AXIS FOR EACH CHANNEL TO SCAN. SO WE SET A TEMPORARY Y OFFSET CALLED TEMSET\r
+       TEMSET=FLOAT(ALLCHS-K)*(YLEN-ALLCHS*.5+.5)/\r
+     *   FLOAT(ALLCHS)+(ALLCHS-K)*.5   @ TEMPORARY Y OFFSET\r
+       CALL XYPLOT (0,TEMSET,-PENUP)   @ SWITCH TO CORRECT PART OF PICTURE\r
+C\r
+C      COMPUTE THE FACTOR TO GET A PRETTY PLOT\r
+C\r
+       K1=CO+1                         @ INDEX TO ORDMIN/MAX SINCE CO STARTS WITH 0\r
+       UP=ORDMAX(K1)\r
+       DOWN=ORDMIN(K1)                 @ INSERT THE DEFAULT SCALING LIMITS\r
+       IF (ORDMAX(K1).LE.ORDMIN(K1)) GOTO 21 @ AUTOMATIC SCALING DESIRED\r
+       OMAX=UP\r
+       OMIN=DOWN\r
+       GOTO 51                         @ SCALING ALREADY DONE\r
+\fC\r
+C      HERE WE DO AUTOMATIC SCALING\r
+C\r
+21     CONTINUE\r
+       OMAX=-99999.\r
+       OMIN=-OMAX\r
+       DO 50 J=S0,S1,STEPS              @ FETCH MIN & MAX OF THE DATA\r
+       IF (BREAK(11)) GOTO 25          @ IF THE USER GETS RID OF PLOT THEN WE EXIT VIA 25\r
+       TE1=XYSAM(J,CO)                 @ DATA POINT --> TE1\r
+       OMAX=AMAX1(OMAX,TE1)            @ MAXIMUM --> OMAX\r
+50     OMIN=AMIN1(OMIN,TE1)            @ MINIMUM --> OMIN\r
+C\r
+       IF (UP .GT. DOWN-1.) GOTO 52    @ IF UP.LE.DOWN THEN WE MAKE A FIXED SCALE\r
+C      FIXED SCALE IF POSSIBLE:\r
+       IF (OMIN.GE.OMAX-(DOWN-UP)) GOTO 53 @ MAKE A FIXED SCALE AT 53\r
+       DOWN=OMIN                       @ DIFFERENCE BETWEEN MIN AND MAX EXCEEDS THE FIX SCALE LENGTH, SO WE TAKE DEFAULT VALUES\r
+       UP=OMAX                         @ ELSE DOWN=UP-(OMIN-OMAX) (FIXED DIFFERENCE SCALE)\r
+       GOTO 51                         @ WORK DONE\r
+C      FIXED SCALE POSSIBLE, MAKE IT AND CENTER THE PLOT\r
+53     TE1=OMIN-((DOWN-UP)-(OMAX-OMIN))*.5\r
+       UP=TE1+DOWN-UP                  @ UPPER LIMIT\r
+       DOWN=TE1                        @ WE SHIFT THE LINE INTO THE CENTER\r
+       GOTO 51                         @ FIXED DIFFERENCE SCALE DONE\r
+C      AUTOMATIC SCALE BUT WITH INCREMENTAL STEPPED LIMITS\r
+52     CONTINUE\r
+       DOWN=COND(OMIN.LT.-512.,OMIN,OMIN.LT.-255.,-512.,\r
+     *   OMIN.LT.0,-255.,OMIN.LT.256.,0, 255.) @ MAKE 255 STEPS INCREMENT\r
+       UP=COND(OMAX.GT.512.,OMAX,OMAX.GT.255.,512.,\r
+     *   OMAX.GT.0,255.,OMAX.GT.-255.,0, -255.)        @ 255 STEPS INCREMENT\r
+C\r
+       UP=COND(DOWN.GE.UP,DOWN+255.,UP) @ SOMETIMES IT WILL HAPPEN: UP=DOWN\r
+51     CONTINUE                        @ HERE WE COMPUTE THE FACTOR FOR ONE LINE\r
+       IF (TIMABZ) CALL FACTOR (XLEN/(RIGHT-LEFT),\r
+     *   (YLEN-ALLCHS*.5+.5)/FLOAT(ALLCHS)/(UP-DOWN))\r
+       IF (.NOT.TIMABZ) CALL FACTOR (AMIN1(XLEN,YLEN)/(RIGHT-LEFT),\r
+     *   AMIN1(XLEN,YLEN)/(UP-DOWN))\r
+       IF BREAK(11) GOTO 25            @ USER GETS RID OF PRGRM\r
+\fC\r
+C      HERE WE COMPRESS ( TO SAVE PLOTTER TIME) AND PLOT THE DATA POINTS\r
+C\r
+       CALL WHERE (TE1,TE1,TE1,TE2)    @ FETCH THE X FACTOR --> TE1\r
+       YFMIN=9999.\r
+       YFMAX=-YFMIN                    @ PREP TEMPS TO FETCH THE MIN/MAX OF THE ORDINATE POINTS\r
+       BOTH=.FALSE.\r
+C\r
+       OLDX=XF(S0)-0.4/TE1*XLEN\r
+       PENPOS=PENUP\r
+       L1=.FALSE.                      @ L1 HOLDS: YF(I) INSIDE LOWER/UPPER LIMIT OF PLOT\r
+C\r
+C      IF THERE IS A X MOVEMENT LESS THAN .3 MM THEN WE DO NOT PLOT THE\r
+C      POINT. WE FETCH THE MIN/MAX OF Y AND PLOT A LINE FROM MIN TO MAX AT THIS X LOCATION.\r
+C      THIS IS TO COMPRESS USELESS OUTPUT TO THE PLOTTER.\r
+C\r
+C\r
+C\r
+       DO 10 I=S0,S1,STEPS              @ HERE WE PLOT THE DATA POINTS ONE BY ONE \r
+       X=XF(I)                         @ TEMP TO SAVE COMPUTING TIME\r
+75     Y=YF(I)                         @  IF MORE THAN ONE SPIKE SET IN THE TRIGGER TIME WORD, WE REENTER HERE\r
+       IF BREAK(11) GOTO 25            @ USER GETS RID OF PROGRAM\r
+       IF(ADPLT)PENPOS=NCOND(BETW(-.1,Y,UP-DOWN+.5).OR.L1,PENPOS,PENUP) @ WE DO NOT PLOT ANY POINTS EXCEEDING THE UPPER/LOWER LIMIT\r
+       IF (.NOT.ADPLT) GOTO 74         @ TRIGGER EVENTS WITHOUT ANY COMPRESSION!\r
+       IF (ABS(OLDX-X)*TE1.GE.0.25/XLEN) GOTO 71 @ PLOT THIS POINT\r
+       YFMIN=AMIN1(Y,YFMIN)            @ NO X MOVEMENT, GET THE MIN Y VALUE\r
+       YFMAX=AMAX1(Y,YFMAX)            @ FETCH THE MAX Y VALUE\r
+       BOTH=.TRUE.                     @ TWO Y POINTS TO PLOT!! (MIN & MAX)\r
+       GOTO 10                         @  WITH THE NEXT REAL X MOVEMENT\r
+C\r
+74      CALL XYPLOT(X,BOX(-DOWN),PENPOS) @ DELTA SPIKE FOR TRIGGER CHANNEL\r
+71     IF (BOTH) GOTO 72               @ IF THERE ARE TWO POINTS TO PLOT THEN GO 72 ELSE\r
+       CALL XYPLOT (X,BOX(Y),PENPOS)    @ WE PLOT THIS ONE AND PRESS IT INTO U/L LIMITS\r
+       IF (.NOT.ADPLT)\r
+     *   CALL XYPLOT(X,BOX(-DOWN),PENPOS) @ DELTA SPIKE FOR TRIGGER CHANNEL\r
+       GO TO 73\r
+72     CALL XYPLOT (X,BOX(YFMIN),PENPOS) @ UPPER AND LOWER Y VALUE TO PLOT\r
+       CALL XYPLOT (X,BOX(YFMAX),PENPOS)\r
+       YFMIN=9999.                     @ SET UP NEW MIN/MAX FETCH CYCLE\r
+       YFMAX=-YFMIN\r
+       BOTH =.FALSE.                   @ CLEAR UPPER/LOWER FLAG\r
+73     CONTINUE\r
+       OLDX=X                          @ SAVE LAST X VALUE\r
+       L1=BETW(-.1,Y,UP-DOWN)          @ L1:= Y INSIDE U/L LIMITS!\r
+       PENPOS=PENDWN\r
+       IF (NXTBIT.GE.0) GOTO 10        @ NO MORE SPIKES IN THIS TRIGG TIME WORD, SKIP TO 10\r
+       X=X+(XF(I+1)-XF(I))/TTRGW       @ COMPUTE NEW X FOR MULTIPLE SPIKE EVENTS\r
+       GOTO 75                         @ REPEAT FOR NEXT SPIKE IN THIS TIME WORD\r
+10     CONTINUE\r
+       IF (.NOT. TIMABZ) GOTO 25       @ WORK DONE FOR X-Y GRAPHICS\r
+C\r
+C      MAKE A SCALE FOR THE ORDINATE ( INPUT VOLTAGE )\r
+C\r
+       CALL WHERE (TE1,TE1,TE1,TE2)\r
+       CALL FACTOR (1.,TE2)\r
+       CALL XYPLOT (-.5,0,-PENUP)\r
+       PENPOS=PENUP\r
+       DO 55 I=0,10\r
+       TE1=DOWN+I*(UP-DOWN)/10.\r
+       CALL XYPLOT (0,TE1-DOWN,PENPOS)\r
+       PENPOS=PENDWN\r
+       J=MOD(I,5)\r
+       TE3=COND(J.EQ.0,-.5,-.3)\r
+       CALL XYPLOT (TE3,TE1-DOWN,PENDWN)\r
+       CALL XYPLOT (0,TE1-DOWN,PENDWN)\r
+       IF (MOD(I,10).NE.0) GOTO 55\r
+       CALL FACTOR (1.,1.)             @ WRITE LOWER AND UPPER LIMIT ONTO THE\r
+       CALL NUMBER (-XOFSET+2.*YZ,(TE1-DOWN)*TE2-YN*.38,YN*.75,TE1,0,\r
+     *   -1) @ AXIS\r
+       CALL FACTOR (1.,TE2)\r
+       CALL XYPLOT (0,TE1-DOWN,PENUP)\r
+55     CONTINUE\r
+       CALL XYPLOT (.5,0,-PENUP)\r
+       CALL FACTOR (1.,1.)\r
+       CALL SYMBOL (-XOFSET+YZ,\r
+     *  AMAX1(0.,((YLEN-ALLCHS*.5)/ALLCHS-4.*YZ))*.5,YZ,\r
+     *   NCOND(ADPLT,'CH ','TRG'),90.,3)\r
+       CALL NUMBER (999.,999.,YZ,\r
+     *   NCOND(ADPLT,CO,CO-CHNLS+1),90.,-1) @ AD CHANNEL # : 0 .. CHNLS, TRG CHAN#: 1 .. STRIGS\r
+C\r
+C\r
+       CO=CO+1                 @ HERE WE INCREMENT THE CHANNEL NUMBER\r
+25     CALL FACTOR (1.,1.)\r
+       CALL XYPLOT (0.,-TEMSET,-PENUP) @ RESET CORRECT ORIGIN OF DIAGRAM\r
+       IF BREAK(11) GOTO 40\r
+20     CONTINUE\r
+\fC\r
+C      IF ABZISSA IS TIME SCALE THEN DRAW A TIME SCALE\r
+C\r
+       IF (.NOT.TIMABZ) GOTO 40\r
+       CALL FACTOR (XLEN/(PLEND-PLBEG),1.)\r
+       CALL XYPLOT (0.,-.5,-PENUP)\r
+C      HERE WE PLOT THE LITTLE BAR INDICATING A SECOND OR MINUTE\r
+       ABISMI=PLEND-PLBEG.GT.300       @ ABZISSE EXCEEDS 5 MINUTES SO WE DRAW A MINUTE SCALING\r
+       EXPAND=NCOND(PLEND-PLBEG.LE.5,10,1) @ IF THE TIME SCALE IS .LE. 5 SECONDS THEN WE MARK EACH 1/10 SEC\r
+C      @ REMEMBER : IF (ABISMI) EXPAND=.FALSE. FOR CORRECT ACTION \r
+       DO 30 I=PLBEG*EXPAND,PLEND*EXPAND @ 1 SEC INCREMENTS\r
+       TE1=I-PLBEG*EXPAND\r
+       CALL XYPLOT (TE1/EXPAND,0,PENDWN)\r
+       IF (.NOT.ABISMI) TE2=COND(MOD(I,10).EQ.0,-.7, MOD(I,5).EQ.0,-.5,\r
+     *   -.3)\r
+       IF (ABISMI) TE2=COND(MOD(I,600).EQ.0,-.7,\r
+     *    MOD(I,300).EQ.0,-.5, MOD(I,60).EQ.0,-.3)\r
+       CALL XYPLOT (TE1/EXPAND,TE2,PENDWN) @ HERE WE PLOT THE SMALL BAR\r
+       CALL XYPLOT (TE1/EXPAND,0,PENDWN)\r
+C      HERE WE COMPUTE WHERE TO WRITE A NUMBER ONTO THE TIME AXIS\r
+       TE2=PLEND-PLBEG                 @  COMPUTE THE SIZE OF THE TIME AXIS\r
+       IF (.NOT.ABISMI) GOTO 31        @ IF WE PLOT A SECONDS AXIS THEN GOTO 31\r
+       IF (MOD(I,\r
+     *   NCOND(TE2.GT.5400,1800,TE2.GE.1800,600,TE2.GT.900,300,\r
+     *   TE2.GT.290,120,60)).GT.0) GOTO 29 @ COMPUTE THE MINUTE WHERE TO PLOT THE NUMBER\r
+       GOTO 32                         @ SKIP TO PLOT THE NUMBER\r
+31     IF (MOD(I,NCOND(TE2.GT.180,30,TE2.GT.60,20,TE2.GT.30,10,\r
+     *   TE2.GT.10,5,TE2.GT.5,2,10)).GT.0) @ WE ASSUME: TE2.LE.5 THEN EXPAND:=10 ELSE EXPAND:=1 !!!\r
+     *   GOTO 30\r
+C\r
+C      HERE WE PLOT THE NUMBER ONTO THE AXIS\r
+C\r
+32     CONTINUE\r
+       CALL WHERE (TE2,TE2,TE2,TE3)    @ GET THE X FACTOR\r
+       CALL FACTOR (1.,1.)             @ RESET TO THE NORMAL FACTOR FOR SYMBOL\r
+       CALL NUMBER (TE1/EXPAND*TE2-INT(ALOG10(AMAX0(COND(ABISMI,\r
+     *    I/60.,I/EXPAND),1))+.01)\r
+     *   *YN/2.-YN*.3,\r
+     *   -1.3,YN,COND(ABISMI,I/60.,I/EXPAND),0,-1)\r
+       CALL FACTOR (TE2,TE3)\r
+       CALL XYPLOT (TE1/EXPAND,0,PENUP)\r
+       CALL XYPLOT (TE1/EXPAND,0,PENDWN)\r
+29     I=NCOND(ABISMI,MIN0(I+59,\r
+     *    NCOND(I.EQ.PLEND*EXPAND,I,PLEND*EXPAND-1)),I) @ MINUTE SCALE HAS 60 SEC INCREMENT\r
+C      SECOND NCOND ABOVE PREVENTS ENDLESS LOOP\r
+C\r
+30     CONTINUE\r
+       CALL FACTOR (1.,1.)\r
+       CALL XYPLOT (0,.5,-PENUP)\r
+        CALL SYMBOL (XLEN/1.5,-YOFSET,YZ,TEXT0(NCOND(ABISMI,1,3)),0,10)\r
+C\r
+C\r
+40     CALL FACTOR (1.,1.)\r
+       IF (BREAK(11)) GOTO 60          @ USER SWITCHED OFF THE PLOTTER\r
+       CALL LABPLT (0,-YOFSET,PLEND)\r
+       IF (BREAK(11)) GOTO 60\r
+       CALL SYMBOL (YZ,YLEN+YZ,YZ,LABEL,0,MIN0(42,LABCNT*6)) @ PLOT 42 CHARACTERS INTO THE FIRST LINE OF LABEL\r
+       IF (LABCNT.GT.7)\r
+     *  CALL SYMBOL (YZ,YLEN,YZ,LABEL(8),0,LABCNT*6-42) @ AND THE REMAINING INTO THE NEXT LINE\r
+60     CONTINUE\r
+       CALL XYPLOT (38.,25.,PENUP)\r
+       CALL EXPLT\r
+       RETURN\r
+1      FORMAT (' KANAL ABZISSE/ORDINATE (ZEITBASIS=8), SCHRITTWEITE',\r
+     *  ' (2I1,I3) (',\r
+     *  2I1,I3,')? ',$)\r
+2      FORMAT (2A1,I3)\r
+3      FORMAT (' ANZAHL DER DARZUSTELLENDEN ',A3,' KANAELE (I1) (',I1,\r
+     *  ') ? ',$)\r
+4      FORMAT (' NUR FUER',F6.1,' SEK. MESSDATEN JEDOCH',F6.1,\r
+     *  ' SEC. ANGEFORDERT. ')\r
+5      FORMAT (' ORDINATEN MAXIMUM FUER ',A3,' KANAL',I2,' (',F5.0,')?',$)\r
+6      FORMAT ('  "   "    MINIMUM FUER ',A3,' KANAL',I2,' (',F5.0,')?',$)\r
+7      FORMAT (F5.0)\r
+9      FORMAT (' ABZISSEN/ ORDINATEN KANAL ODER SCHRITTWEITE ',\r
+     *  ' FEHLERHAFT: ',3I4)\r
+       END\r
+\1a\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0
\ No newline at end of file