A large commit.

[pdp8.git] / sw / rescue / lab8e_goettingen / disk2_11 / rkb / paroff / haespl.ft

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