SUBROUTINE LORGEN (J, TMASS, TP, DCTGT, BETLAB, BPLAB,
1 DCBM, EIN, PIN, DCIN, EOUT, POUT, DCOUT ) 0030
CLORGEN LORENTZ TRANSFORMATION GENERAL 0020
C J = O GIVES TRANSFORM FROM LAB INTO C OF M 0040
C J NOT ZERO GIVES TRANSFORM FROM C OF M INTO LAB 0050
C ERROR RETURN IS POUT ZERO OR NEGATIVE 0060
C BETLAB, BPLAB, DCTRA ARE TOTAL ENERGY, MOMENTUM AND DCS OF 0070
C 1 INCIDENT PART 0080
C EIN, PIN, DCIN ARE SAME FOR TRANSFORMED PART 0090
C POUT, DCOUT ARE SAME (RETURNED) AFTER TRANSFORMATION 0100
C TMASS IS TARGET MASS * ZERO FORM DECAY 0110
C NO COMMON REQUIRED 0120
DIMENSION DCTOT(3), DCTGT(3), DCBM(3) 0130
BETA = 0 0150
DO 50 K = 1,3 0160
DCTOT(K) = DCTGT(K)*TP+DCBM(K)*BPLAB 0170
50 BETA = BETA + DCTOT(K)**2 0180
BETA = SQRT(BETA) 0190
DO 60 K = 1,3 0200
60 DCTOT(K) = DCTOT(K)/BETA 0210
BETA = BETA/(BETLAB + SQRT(TP**2 + TMASS**2)) 0220
IF (J) 63, 65, 63 0230
63 BETA = -BETA 0240
65 CONTINUE 0250
CALL LORTRA(BETA, DCTOT(1),EIN,PIN,DCIN,EOUT,POUT,DCOUT) 0260
RETURN 0270
END 0280
SUBROUTINE LORTRA(BETA,DCTRA,EIN,PIN,DCIN,EOUT,POUT,DCOUT)
CLORTRA LORENTZ TRANSFORMATION SUBROUTINE 0020
DIMENSION DCTRA(3), DCIN(3), DCOUT(3), TRAN(4,4),PUNT(4),PTRAN(4) 0030
GAMMA = 1.0 / SQRT( 1.0 - BETA**2 ) 0050
5 CONTINUE 0060
DO 10 L = 1,4 0070
DO 9 M = 1,4 0080
9 TRAN(L,M) = 0.0 0090
10 TRAN(L,L) = 1.0 0100
DO 15 L = 1,3 0110
DO 12 M = L,3 0120
12 TRAN (L,M) = TRAN (L,M) + (GAMMA - 1.0) * DCTRA(L) * DCTRA(M) 0130
TRAN (L,4) = -BETA * DCTRA(L) * GAMMA 0140
LL = L + 1 0150
DO 15 N = LL, 4 0160
TRAN (N,L) = TRAN(L,N) 0170
15 CONTINUE 0180
TRAN(4, 4) = GAMMA 0190
DO 20 K = 1,3 0200
20 PUNT(K) = PIN * DCIN(K) 0210
PUNT(4) = EIN 0220
DO 25 L = 1,4 0230
PTRAN(L) = 0.0 0240
DO 25 M = 1,4 0250
25 PTRAN(L) = PTRAN(L) + TRAN(L,M) * PUNT (M) 0260
POUT = SQRT( PTRAN(1) **2 + PTRAN(2)**2 + PTRAN(3)**2) 0270
DO 30 K = 1,3 0280
30 DCOUT(K) = PTRAN(K) / POUT 0290
EOUT = PTRAN(4) 0300
RETURN 0310
END 0320
SUBROUTINE LTON(N,A,L,LENGTH)
COMMON /LET/ NBLANK,LET(30),NPLUS,NUM(10)
DIMENSION A(2),L(2)
EQUIVALENCE (KA,AA)
LFLAG=+1
NN=N
DO 10 I=1,NN
10 L(I)=0
NZ=1
DO 70 I=1,NN
AA=A(I)
DO 30 M=1,30
IF(KA.NE.LET(M)) GO TO 30
KA=M
GO TO 60
30 CONTINUE
DO 40 M=1,10
IF(KA.NE.NUM(M)) GO TO 40
KA=M-1
GO TO 60
40 CONTINUE
IF(KA.NE.NPLUS) GO TO 50
LFLAG=-1
GO TO 70
50 IF(KA.EQ.NBLANK) GO TO 70
KA=0
60 L(NZ)=LFLAG*KA
NZ=NZ+1
LFLAG=+1
70 CONTINUE
LENGTH=NZ-1
RETURN
END
SUBROUTINE NTOL(N,A,L)
COMMON /LET/ NBLANK,LET(30),NPLUS,NUM(10)
DIMENSION A(2),L(2)
EQUIVALENCE (KA,AA)
DO 200 I=1,N
LL=L(I)
IF(LL)115,105,110
105 KA=NBLANK
GO TO 120
110 IF(LL.LE.30) GO TO 125
115 KA=0
120 A(I)=AA
GO TO 300
125 KA=LET(LL)
A(I)=AA
200 CONTINUE
300 RETURN
END
FUNCTION NOTABL (LA, LB)
CNOTABL FORTRAN FUNCTION TO FIND OTABLE INDEX GIVEN END LABELS 0020
C ************************* COMMON COMMON ************************** 0030
COMMON MAP(2000),PARS(1000),MISC(27),KLIST(500),MTABLE(2)
DIMENSION LTABLE(9,20,2),ITABLE(6,20)
DIMENSION PARA(1000),NPARA(1000),SNAME(1000),NAME(1000),TABLE(100)
DIMENSION HEAD(11), NBRNCH(10),HTABLE(7,100)
EQUIVALENCE (MAP,ZMAP),(TABLE(1),PARA(101))
EQUIVALENCE (OTABLE,JTABLE,MAP(701)), (RTABLE,LTABLE,MAP(1051)),
1 (ITABLE,MAP(1411)), (VAL,IVAL,MAP(1531)),
2 (WGT,MAP(1631))
EQUIVALENCE (NCFLAG,MAP(1869)), (WEIGHT,MAP(1978)),
1 (NTAPE,MAP(1988)), (EINC,MAP(1998)),
2 (PINC,MAP(1999)), (BINC,MAP(2000))
EQUIVALENCE (PARA,NPARA,PARS),(SNAME,NAME,MAP(1))
EQUIVALENCE (PI, MISC), (RADIAN, MISC(2)), (NIT, MISC(3)),
1 (NOT, MISC(4)), (HEAD, MISC(5)), (NBRNCH, MISC(16)),
2 (NPAGE, MISC(26)), (NORD, MISC(27))
EQUIVALENCE (LTAPE,NBRNCH(9)), (LINK,NBRNCH(10)), (HTABLE,MAP)
C ************************* END OF C, D, E STATEMENTS ************** 0260
C 0270
IF (LA.LT.1 .OR. LA.GT.20) GO TO 15
IF (LA.EQ.1 .AND. LB.EQ.1) GO TO 18
DO 10 L = 1,8 0280
IF (LTABLE(L,LA,2) - LB) 10, 20, 10 0290
10 CONTINUE 0300
C ERROR RETURN ON FAILURE TO FIND TRACK 0310
15 NOTABL = 0
GO TO 30 0330
C BEAM TRACK
18 NOTABL = 1
GO TO 30
20 NOTABL = ITABLE(2,LA) + L 0340
30 RETURN 0350
END 0360
SUBROUTINE OFIX (NFLG)
COFIX SUBROUTINE TO SET UP OTABLE INDICES 0020
C ************************* COMMON COMMON ************************** 0030
COMMON MAP(2000),PARS(1000),MISC(27),KLIST(500),MTABLE(2)
DIMENSION ZMAP(2000) NBOD0110
DIMENSION REMARK(500)
DIMENSION OTABLE(7,50), JTABLE(7,50), RTABLE(9,20,2), NBOD0120
1 LTABLE(9,20,2), ITABLE(6,20), VAL(100), IVAL(100), NBOD0130
2 WGT(100) NBOD0140
DIMENSION PARA(1000),NPARA(1000),SNAME(1000),NAME(1000)
DIMENSION HEAD(11), NBRNCH(10) NBOD0170
DIMENSION HTABLE(7,100) NBOD0180
EQUIVALENCE (MAP,ZMAP) NBOD0340
EQUIVALENCE (REMARK,MAP(1001))
EQUIVALENCE (OTABLE,JTABLE,MAP(701)), (RTABLE,LTABLE,MAP(1051)), NBOD0350
1 (ITABLE,MAP(1411)), (VAL,IVAL,MAP(1531)), NBOD0360
2 (WGT,MAP(1631)) NBOD0370
EQUIVALENCE (NCFLAG,MAP(1869)), (WEIGHT,MAP(1978)), NBOD0380
1 (NTAPE,MAP(1988)), (EINC,MAP(1998)), NBOD0390
2 (PINC,MAP(1999)), (BINC,MAP(2000)) NBOD0400
EQUIVALENCE (PARA,NPARA,PARS),(SNAME,NAME,MAP(1))
EQUIVALENCE (PI, MISC), (RADIAN, MISC(2)), (NIT, MISC(3)), NBOD0490
1 (NOT, MISC(4)), (HEAD, MISC(5)), (NBRNCH, MISC(16)), NBOD0500
2 (NPAGE, MISC(26)), (NORD, MISC(27)) NBOD0510
EQUIVALENCE (HTABLE,MAP)
C ************************* END OF C, D, E STATEMENTS ************** 0260
C 0270
IF (NFLG) 440, 400, 440
400 NFLG = NFLG+1
C OTABLE(I,1) IS RESERVED FOR THE INCIDENT BEAM TRACK
JTABLE(6,1) = 1
JTABLE(7,1) = 1
II = 2
KK = 1
C ALLOT A POSITION IN OTABLE TO EACH REMAINING TRACK
DO 430 K=1,20
NPROD = ITABLE(1,K)
IF (NPROD) 430, 430, 410
410 DO 420 I=1,NPROD
JTABLE(6,II) = K
JTABLE(7,II) = LTABLE(I,K,2)
420 II = II+1
ITABLE(2,K) = KK
KK = KK+NPROD
430 CONTINUE
440 RETURN
END
SUBROUTINE PAREAD(NIT,NOT,NBRNCH,HEAD,PARA,LPARA,SNAME,RMAP,LMAP)
DIMENSION NBRNCH(10), HEAD(10), REMARK(9)
DIMENSION PARA(2), SNAME(2), RMAP(10,2) PARE0060
DIMENSION IADDRS(10), ILAST(10)
EQUIVALENCE (TEMP, ITEMP), (TNAME, NAME)
EQUIVALENCE (ATABLE,ITABLE), (BLANKS,IBLNKS)
DATA BLANKS /' '/
DATA IH9999/'H9999'/, IOBBBB/'O '/
DATA ILARGE /16777216/
C PARE0120
NTABLE = 0 PARE0130
20 READ (NIT,1) NBRNCH, HEAD
1 FORMAT (10I2,10A4)
IF(NBRNCH(1)) 25,20,20
25 IF(LPARA.GT.0) GO TO 28
C
C FIRST CALL - CLEAR PARA
LPARA=IABS(LPARA)
26 DO 27 I=1,LPARA
27 PARA(I) = 0.
C
28 DO 10 I=1,LPARA
10 SNAME(I)=BLANKS
29 LIMAP = LMAP/10 JULY 12
IF(LMAP) 30,30,11
11 CONTINUE
DO 12 J = 1,LIMAP PARE0340
DO 12 I=1,10 OCT 1
12 RMAP(I,J)=BLANKS OCT 1
30 NBRNCH(1) = IABS(NBRNCH(1))
WRITE (NOT,2) NBRNCH, HEAD
2 FORMAT (1H1 10I3, 10X, 10A4 )
II = 0 PARE0460
C BEGIN READ-IN LOOP OCT 1
DO 110 I = 1,LPARA PARE0470
READ (NIT,3)NN, ATABLE, TEMP, TNAME, REMARK PARE0490
3 FORMAT (I3, A1, 1X, F12.0, 3X,A4,2X9A4)
IF (NN) 32, 120, 32 PARE0510
32 IF (NN - LPARA) 36, 36, 34 PARE0520
34 WRITE (NOT,4)NN, TEMP, TNAME, REMARK PARE0530
4 FORMAT(26H0THIS CARD EXCEEDS LIMITS I3, 2X, F12.5, 2X, A4, 3X,
1 9A4)
GO TO 110 PARE0560
36 IF (ITABLE - IBLNKS) 380, 38, 380
380 ILONG = TEMP PARE0580
IBEG = NN PARE0590
IEND = IBEG + ILONG - 1 PARE0600
NTABLE = NTABLE + 1 PARE0610
IADDRS(NTABLE)=NN
ILAST(NTABLE)=IEND
IF (IEND - LPARA) 400, 400, 500 PARE0650
400 READ (NIT,REMARK)(PARA(MN), MN = IBEG,IEND) PARE0660
IF (NTABLE - 1) 390, 385, 390 PARE0670
385 WRITE (NOT,386) PARE0680
386 FORMAT (15H0TABLES READ IN )
390 WRITE (NOT,391)TNAME, ILONG, IBEG, IEND PARE0700
391 FORMAT (1H0 A4, 1H( I3, 9H) = PARA( I3, 1H- I3, 1H) /)
WRITE (NOT,REMARK)(PARA(MN), MN = IBEG, IEND) PARE0720
GO TO 110 PARE0730
C PARE0740
C CHECK NAME FOR BLANKS AND DETERMINE FIXED OR FLOATING PARE0750
C PARE0760
38 II=II+1 OCT 1
IF (NAME.GT.IH9999 .AND. NAME.LT.IOBBBB) ITEMP = TEMP
70 PARA(NN) = TEMP PARE0860
C PARE0870
C PLACE SYMBOLIC NAMES IN SNAME AND REMARKS IN RMAP PARE0880
C PARE0890
SNAME(NN) = TNAME PARE0900
IF (II - LIMAP) 80, 80, 110 PARE0910
80 RMAP(1,II) = NN PARE0920
DO 90 JJ = 2,10 PARE0930
RMAP(JJ,II) = REMARK(JJ - 1) PARE0940
90 CONTINUE PARE0950
110 CONTINUE PARE0960
C PARE0970
C PRINT OUT PARAMETERS PARE0980
C PARE0990
120 WRITE (NOT,5) PARE1000
5 FORMAT (25H0PARAMETERS CURRENTLY ARE /)
C BEGIN PRINT-OUT LOOP OCT 1
DO 180 M = 1,LPARA
IF(NTABLE) 211,211,111
111 DO 200 KCOR=1,NTABLE
IF(M-IADDRS(KCOR)) 200,180,191
191 IF(M-ILAST(KCOR))180,180,200
200 CONTINUE
211 TEMP = PARA(M)
IF (IABS(ITEMP).LT.ILARGE) GO TO 190
185 IFLAG=0 OCT 1
GO TO 210 OCT 1
190 IFLAG=1 OCT 1
210 IF (PARA(M)) 130, 888, 130
888 IF (SNAME(M).EQ.BLANKS) GO TO 180
130 FM = M PARE1090
DO 140 II = 1,LIMAP PARE1100
IJ = II PARE1110
IF (RMAP(1,II) - FM) 140, 150, 140 PARE1120
140 CONTINUE PARE1130
GO TO 160 PARE1140
150 IF (IFLAG)152,152,151 OCT 1
151 WRITE(NOT,88)M, SNAME(M), ITEMP, (RMAP(K,IJ),K=2,10)
88 FORMAT(1H I4, 3X,A4,2XI12, 6X, 9A4)
GO TO 180 OCT 1
152 WRITE(NOT,8) M,SNAME(M),TEMP, (RMAP(K,IJ),K=2,10)
8 FORMAT (1H I4, 3X,A4,2XF12.5, 6X, 9A4)
GO TO 180 PARE1170
160 IF(IFLAG) 162,162,161 OCT 1
161 WRITE (NOT,88)M, SNAME(M), ITEMP
GO TO 180 OCT 1
162 WRITE (NOT,8) M, SNAME(M), TEMP
180 CONTINUE PARE1200
GO TO 1000 PARE1210
500 WRITE (NOT,501)TNAME, ILONG, IBEG PARE1220
501 FORMAT (12H0 THE TABLE A4, 1H( I3, 20H) BEGINNING AT PARA( I3, PARE1230
1 64H) EXCEEDS DIMENSION LIMITS. FURTHER EXECUTION TERMINATED BY EXPARE1240
2IT) PARE1250
CALL EXIT PARE1260
1000 RETURN PARE1270
END PARE1280
SUBROUTINE SCALW
COMMON MAP(2000), PARS(1000)
EQUIVALENCE (WSCALE,MAP(1972))
WSCALE = 1.0
IF (PARS(93) .GT. 0.) WSCALE = PARS(93)
RETURN
END
FUNCTION SPACE (KTABLE, LIM)
CSPACE SPACES OUT MTABLE ENTRIES FOR WEIGHTED HISTOGRAMS 0020
C 0030
DIMENSION KTABLE(7,100) 0040
C 0060
SPACE = 0.0 0070
MNO = 1 0080
C 0090
DO 60 K = 1, 100 0100
KC = KTABLE(1,K) 0110
KN = KTABLE(5,K) 0120
IF (KC) 10, 80, 10 0130
10 KTABLE(2,K) = MNO 0140
IF (KC - 1000) 20, 30, 30 0150
20 MNOINC = 3 * (IABS(KN) + 2) 0160
GO TO 60 0170
30 IF (KC - 2000) 40, 50, 50 0180
40 MNOINC = 0 0190
NX = IABS(KN) + 1 0200
GO TO 60 0210
50 MNOINC= NX* (IABS(KN) + 1) + 5 0220
60 MNO = MNO + MNOINC 0230
C 0240
80 IF (MNO - LIM) 100, 100, 90 0250
90 SPACE = 1.0 0260
100 RETURN 0270
END 0280
FUNCTION SQMASS (KB, KE)
C ****************** COMMON COMMON *************************** 0020
COMMON MAP(2000),PARS(1000),MISC(27),KLIST(500),MTABLE(2)
DIMENSION OTABLE(7,50), VM(3)
EQUIVALENCE (OTABLE,MAP(701))
C ************************* END OF C, D, E STATEMENTS **************
C 0250
EMS = 0.0 0260
DO 5 N = 1,3 0270
5 VM(N) = 0.0 0280
DO 30 K = KB,KE 0290
KL = KLIST(K) 0300
EMS = EMS + OTABLE(5,KL) 0310
DO 10 L = 1,3 0320
10 VM(L) = VM(L) + OTABLE(L,KL) * OTABLE(4,KL) 0330
30 CONTINUE 0340
C 0350
SQMASS = EMS**2 - VM(1)**2 - VM(2)**2 - VM(3)**2 0360
RETURN 0370
END 0380
SUBROUTINE TRHIST(NOT,NPRINT,NAME,NHST,HST,ERROR,M,A,B)
C
C TRHIST*1 - SUBROUTINE TO PRINT HISTOGRAMS PREPARED BY HISTO
C NPRINT = 0 FOR INTEGER HISTOGRAM (HISTO) TRHI0030
C NPRINT = 1 FOR IDEOGRAM PRINT (HISTOI) TRHI0040
C NPRINT = 2 FOR WEIGHTED HISTOGRAM (HISTOW) TRHI0050
C
DIMENSION NHST(10),HST(10),ERROR(10),HIST(60),NAME(3)
IF (NOT) 5000, 5000, 102 TRHI0090
102 CONTINUE TRHI0100
NSUM = 0 TRHI0110
NSCALE = 0 TRHI0120
SUM = 0.0 TRHI0130
SCALE = 0.0 TRHI0140
MORE = M+1 TRHI0150
LIM = M+2 TRHI0160
AA = A TRHI0170
DO 105 I = 1, LIM TRHI0180
IF (NHST(I)) 105, 105, 106 TRHI0190
105 CONTINUE TRHI0200
WRITE (NOT,9105)NAME TRHI0210
C CHANGE 3A4 TO 3A5 FOR PDP-10
9105 FORMAT ( 11H1HISTOGRAM 3A5, 13H IS UNFILLED )
GO TO 5000 TRHI0230
106 CONTINUE TRHI0240
IF (NPRINT - 1) 180, 110, 110 TRHI0250
110 DO 120 I = 1,M
IF (HST(I) - SCALE) 120, 120, 115 TRHI0270
115 SCALE = HST(I) TRHI0280
120 SUM = SUM + HST(I) TRHI0290
SUM = SUM + HIST(MORE) + HIST(LIM)
180 DO 190 I = 1,M
IF (NHST(I) - NSCALE) 190, 190, 185 TRHI0310
185 NSCALE = NHST(I) TRHI0320
190 NSUM = NSUM + NHST(I) TRHI0330
NSUM = NSUM + NHST(MORE) + NHST(LIM)
SCALEN = NSCALE TRHI0340
200 EM = M TRHI0350
D = (B-A)/EM TRHI0360
530 WRITE (NOT,9010)NAME,NSUM,SUM,NSCALE TRHI0370
C CHANGE 3A4 TO 3A5 FOR PDP-10
9010 FORMAT (34H1 THE FOLLOWING IS A HISTOGRAM OF 3A5, 6H WITH I5, TRHI0380
1 23H UNWEIGHTED EVENTS AND F8.2, 34H WEIGHTED EVENTS. SCALE FACTRHI0390
2TOR = I4/// 102H0 INTERVAL NOT WEIGHTED TRHI0400
3 ERROR HISTOGRAM BASED ON SCALE FACTOR = 60 / TRHI0410
4 23X 9HWEIGHTED /)
IF (NPRINT - 1) 535, 540, 545 TRHI0430
535 NUM = (NHST(LIM) * 60 ) / NSCALE TRHI0440
WTD = 0.0 TRHI0450
GO TO 542 TRHI0460
540 NUM = (HST(LIM) * 60.0 ) / SCALE TRHI0470
WTD = HST(LIM) TRHI0480
542 ERR = NHST(LIM) TRHI0490
ERR = SQRT(ERR) TRHI0500
GO TO 550 TRHI0510
545 NUM = (HST(LIM) * 60.0 ) / SCALE TRHI0520
WTD = HST(LIM) TRHI0530
ERR = SQRT(ERROR(LIM)) TRHI0540
550 CALL LOADX(HIST, NUM) TRHI0550
WRITE (NOT,9020)AA, NHST(LIM), WTD, ERR, HIST TRHI0560
9020 FORMAT (13H0LESS THAN F8.3,4XI5,5XF8.3,5XF8.3,4X60A1 /)
DO 300 I = 1,M TRHI0580
AB =AA + D TRHI0590
220 IF (NPRINT-1) 230, 235, 240 TRHI0600
230 NUM = (FLOAT(NHST(I)) * 60.0) / SCALEN TRHI0610
WTD = 0.0 TRHI0620
GO TO 237 TRHI0630
235 NUM = ( HST(I) * 60.0 ) / SCALE TRHI0640
WTD = HST(I) TRHI0650
237 ERR = NHST(I) TRHI0660
ERR = SQRT(ERR) TRHI0670
GO TO 245 TRHI0680
240 NUM = ( HST(I) * 60.0 ) / SCALE TRHI0690
WTD = HST(I) TRHI0700
ERR = SQRT(ERROR(I)) TRHI0710
245 CALL LOADX ( HIST, NUM) TRHI0720
WRITE (NOT,9030)AA, AB, NHST(I), WTD, ERR, HIST TRHI0730
9030 FORMAT (1H F8.3,4H TO F8.3,4XI5,5XF8.3,5XF8.3,4X60A1 ) 11/14/64
AA= AB TRHI0750
300 CONTINUE TRHI0760
IF (NPRINT - 1) 335, 340, 345 TRHI0770
335 NUM = (NHST(MORE)* 60 ) / NSCALE TRHI0780
WTD = 0.0 TRHI0790
GO TO 342 TRHI0800
340 NUM = (HST(MORE)* 60.0 ) / SCALE TRHI0810
WTD = HST(MORE) TRHI0820
342 ERR = NHST(MORE) TRHI0830
ERR = SQRT(ERR) TRHI0840
GO TO 350 TRHI0850
345 NUM = (HST(MORE)* 60.0 ) / SCALE TRHI0860
WTD = HST(MORE) TRHI0870
ERR = SQRT(ERROR(MORE)) TRHI0880
350 CALL LOADX(HIST, NUM) TRHI0890
WRITE (NOT,9040)AA, NHST(MORE),WTD, ERR, HIST TRHI0900
9040 FORMAT(13H0GREATER THAN F8.3,4XI5,5XF8.3,5XF8.3,4X60A1 )
5000 RETURN TRHI0920
END
SUBROUTINE LOADX (XSS,NUM)
C
C
DIMENSION XSS(60)
DATA X /'X'/, BLANK /' '/
NA = MIN0(NUM,60) 11/14/64
DO 110 J = 1, NA 11/14/64
110 XSS(J)=X
IF (NA-60) 112, 120, 120 11/14/64
112 NA = NA + 1 11/14/64
DO 114 J = NA,60 11/14/64
114 XSS(J)=BLANK
120 CONTINUE 11/14/64
RETURN 11/14/64
END