C .TITLE CONP C LOGICAL B1,B2,B3,B4,B5,B6,B7,B8,B9,B10,F INTEGER S DIMENSION FILE(2) COMMON/EYCEP/EYE(3),CEN(3),UP(3) COMMON/POINTS/INIT,PEC(4,160) COMMON/PULDAT/N COMMON/PAINT/I2STRT,I2END,PIC(4,200) COMMON/OIL/IVECT(800) COMMON/LIMITS/XMIN,XMAX,YMIN,YMAX COMMON/SUBS/NN,LL,LI,DICT,ICS,ICL DATA FILE(1),FILE(2)/4HCOND,4H SRC/ DATA ENFOR,EN2,EN3/5H(1X,2,5HHN=,I,2H4)/ DATA EMFOR,EM2,EM3/5H(1X,3,5HHSET,,3HI2)/ DATA ELFOR,EL2,EL3/5H(1X,3,5HHLI=,,3HI2)/ DATA IDS,IDL/2,1/ CALL CLEAR(.TRUE.) CALL DISPLY (0,4000,IGET) CALL BUTTON(1) F=.FALSE. XMIN=-5.0 XMAX=5.0 YMIN=-5.0 YMAX=5.0 CALL SEEK (1,FILE) READ(1,501)((PEC(I,J),I=1,4),J=1,28) 501 FORMAT (4F5.1) DO 503 J=1,4 503 PIC(J,12)=PEC(J,12) 504 FORMAT(1X,2HN=,I4,5H S=,I2,5H LI=,I2) LL=3 IBX=1 EYE(1)=0.0 EYE(2)=0.0 EYE(3)=10.0 CEN(1)=0.0 CEN(2)=0.0 CEN(3)=0.0 UP(1)=0.0 UP(2)=10.0 UP(3)=0.0 DIST=10.0 I2STRT=12 I2END=12 IBX=1 MODE=1 NN=0 LI=5 505 FORMAT(1X,5HILD1=,I3) 506 FORMAT(1X,5HILD2=,I3) 507 FORMAT(1X,5HILA1=,I3) 508 FORMAT(1X,5HILA2=,I3) 509 FORMAT(1X,5HILD3=,I3) 510 FORMAT(1X,5HILA2=,I3) 511 FORMAT(1X,6F5.1) 512 FORMAT(1X,5HIARC=,I3) RD=3.1415926/192.0 CALL VECTOR(1,PIC,I2STRT,I2END,IVECT,1,LAST,LI,0,0) CALL DISPLY(6,IDL,IVECT,1,LAST) CALL DISPLY(6,IDS,IVECT,1,LAST) CALL DISPLY(3,IDL) IH=IDS IDS=IDL IDL=IH I2STRT=11 DICT=DIST ICS=IDS ICL=IDL 1 CALL TEXT(0,0,0,0,0) S=0 IF(LL.EQ.3) GO TO 3 WRITE(4,504) NN,S,LI GO TO 4 3 CALL TEXT(465,0,7,1,EMFOR,S) CALL TEXT(535,0,7,1,ELFOR,LI) CALL TEXT(500,1000,7,1,ENFOR,NN) 4 CALL BUTT(B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) CALL BUTTON(3,B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) 2 CALL BUTTON(2,B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) IF(B1) GO TO 15 IF(B2) GO TO 25 IF(B3) GO TO 35 IF(B4) GO TO 45 IF(B5) GO TO 50 IF(B6) GO TO 60 IF(B7) GO TO 70 IF(B8) GO TO 80 IF(B9) GO TO 90 IF(B10) GO TO 100 GO TO 2 15 CALL TEXT(0,0,0,0,0) S=1 IF(LL.EQ.3) GO TO 16 WRITE(4,504) NN,S,LI GO TO 17 16 CALL TEXT(465,0,7,1,EMFOR,S) CALL TEXT(535,0,7,1,ELFOR,LI) CALL TEXT(500,1000,7,1,ENFOR,NN) 17 CALL BUTT(B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) CALL BUTTON(3,B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) 11 CALL BUTTON(2,B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) IF(B1) GO TO 110 IF(B2) GO TO 120 IF(B3) GO TO 130 IF(B4) GO TO 140 IF(B5) GO TO 1 IF(B6) GO TO 160 IF(B7) GO TO 170 IF(B8) GO TO 180 IF(B9) GO TO 1 IF(B10) GO TO 1 GO TO 11 25 CALL TEXT(0,0,0,0,0) S=2 IF(LL.EQ.3) GO TO 23 WRITE(4,504) NN,S,LI GO TO 24 23 CALL TEXT(465,0,7,1,EMFOR,S) CALL TEXT(535,0,7,1,ELFOR,LI) CALL TEXT(500,1000,7,1,ENFOR,NN) 24 CALL BUTT(B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) CALL BUTTON(3,B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) 22 CALL BUTTON(2,B1,B2,B3,B4,B5,B6,B7,B8,B9,B10) IF(B1) GO TO 210 IF(B2) GO TO 220 IF(B3) GO TO 230 IF(B4) GO TO 240 IF(B5) GO TO 1 IF(B6) GO TO 260 IF(B7) GO TO 270 IF(B8) GO TO 280 IF(B9) GO TO 290 IF(B10) GO TO 294 GO TO 22 35 LL=3 GO TO 1 45 LL=4 GO TO 1 50 IBX=2 GO TO 1 60 DO 61 K=1,24 CALL PULSE(N,N) 61 NN=NN+1 GO TO 1 70 NN=0 GO TO 1 80 DO 81 I=14,28 DO 82 J=1,4 82 PIC(J,I)=PEC(J,I) I2END=I CALL ROTP(8) 81 CONTINUE CALL ROTP(96) GO TO 1 90 GO TO 1 100 GO TO 1 110 DO 111 I=11,13 DO 111 J=1,4 111 PIC(J,I)=PEC(J,I) I2END=28 CALL EXP(29,100) CALL ROTP(48) II=100 CALL PROJ(0,20.0,11,12,13,1,II) ILD1=II CALL PROJ(0,20.0,11,12,13,3,II) I2END=II ILD2=II CALL ROTP(96) GO TO 15 120 CALL EXP(29,200) DO 122 J=1,4 122 PIC(J,11)=PEC(J,11) DO 121 I=12,13 DO 121 J=1,4 II=I-11 121 PIC(J,I)=PEC(J,II) CALL ROTP(48) II=100 CALL PROJ(0,20.0,11,12,13,1,II) I2END=II CALL ROTP(24) ILD1=II CALL PROJ(0,2.0,11,12,13,3,II) ILD2=II CALL ROTP(48) IF(IBX.EQ.2) IBX=3 / GO TO 1 130 CS=COS(-6.0*RD) SN=SIN(-6.0*RD) IF(IBX.EQ.3) IBX=2 IK=24 GO TO 131 140 CS=COS(+4.0*RD) SN=SIN(+4.0*RD) IK=12 IF(IBX.EQ.2)IBX=3 ICX=2 131 CALL EXP(13,ILD2) DO 132 I=14,17 DO 132 J=1,4 132 PIC(J,I)=PEC(J,I) II=44 PIC(1,II)=PIC(1,12) PIC(2,II)=PIC(2,12) CALL DASP(20.0,1,11,II) BS=PIC(4,12) II=II+1 PIC(1,II)=PIC(3,12) PIC(2,II)=BS CALL DASP(20.0,3,11,II) ILD1=II I2END=II CALL ROTP(48) DO 133 J=1,4 PIC(J,13)=PIC(J,11) 133 PIC(J,18)=PIC(J,12) ILD2=ILD1+1 IA=19 PIC(1,19)=PIC(3,12) PIC(2,19)=BS PIC(4,19)=BS DO 134 K=1,IK IA=IA+1 IS=ILD1+1 CALL EXP(IS,ILD2) PIC(1,IA)=PIC(3,11) PIC(2,IA)=PIC(4,11) I2END=11 CALL RORP(CS,SN,-2.0,2.0) PIC(3,IA)=PIC(3,11) PIC(4,IA)=PIC(4,11) PIC(3,12)=PIC(3,11) PIC(3,19)=PIC(3,11) II=ILD1+1 PIC(1,II)=PIC(3,12) PIC(2,II)=PIC(4,12) CALL DASP(20.0,3,11,II) I2END=II ILD2=II CALL ROTP(4) 134 CONTINUE CALL ROTP(96) DO 135 J=1,4 PIC(J,11)=PIC(J,13) 135 PIC(J,12)=PIC(J,18) CALL EXP(18,IA) CALL EXP(IS,ILD2) GO TO(15,140,15),IBX 160 DO 164 J=1,4 PIC(J,10)=PIC(J,11) PIC(J,29)=PIC(J,12) 164 PIC(J,30)=PIC(J,13) XVS=PIC(3,12) XPS=PIC(3,13) BS=PIC(4,12) XV=PIC(1,12) YV=PIC(2,12) XP=PIC(1,13) YP=PIC(2,13) IK=384/6 CS=COS(-6.0*RD) SN=SIN(-6.0*RD) ILD3=ILD2 IA=36 DO 161 K=1,IK IA=IA+1 IS=ILD1+1 CALL EXP(IS,ILD3) PIC(1,IA)=PIC(3,11) PIC(2,IA)=PIC(4,11) I2END=11 I2SRT=11 CALL RORP(CS,SN,-2.0,2.0) I2STRT=10 I2END=ILD3 PIC(3,IA)=PIC(3,11) PIC(4,IA)=PIC(4,11) X=PIC(3,11) Y=PIC(4,11) PIC(3,12)=X PIC(3,13)=Y PIC(4,13)=BS PIC(4,12)=BS PIC(2,12)=YV PIC(2,13)=YP PIC(1,13)=XP PIC(1,12)=XV II=IS DM=20.0 CALL PROJ(2,DM,11,12,13,3,II) ILD2=II DN=10.0 611 IF(K.GT.7)GO TO 165 PIC(1,31)=XVS PIC(2,31)=BS PIC(3,31)=X PIC(4,31)=BS PIC(1,32)=XPS PIC(2,32)=BS PIC(3,32)=Y PIC(4,32)=BS CALL EXP(33,36) GO TO 600 165 IF(K.GT.8)GO TO 166 DO 601 J=1,4 601 PIC(J,33)=PIC(J,13) PIC(3,31)=X PIC(3,32)=Y GO TO 600 166 IF(K.GT.16)GO TO 167 II=ILD2+1 PIC(2,II)=BS PIC(1,II)=PIC(3,13) CALL DASP(DN,1,13,II) ILD2=II CALL EXP(13,13) PIC(3,31)=X GO TO 600 167 IF(K.GT.23)GO TO 168 CALL EXP(30,30) PIC(3,31)=X PIC(2,36)=BS PIC(1,36)=XPS PIC(3,36)=Y PIC(4,36)=BS GO TO 600 168 IF(K.GT.24)GO TO 169 DO 602 J=1,4 602 PIC(J,34)=PIC(J,12) PIC(3,31)=X PIC(3,36)=Y GO TO 600 169 IF(K.GT.39)GO TO 603 PIC(3,36)=Y II=ILD2+1 PIC(1,II)=PIC(3,12) PIC(2,II)=BS CALL DASP(DN,1,12,II) ILD2=II CALL EXP(12,12) GO TO 600 603 IF(K.GT.40)GO TO 604 DO 605 J=1,4 605 PIC(J,30)=PIC(J,13) PIC(3,36)=Y II=ILD2+1 PIC(1,II)=PIC(3,12) PIC(2,II)=BS CALL DASP(DN,1,12,II) ILD2=II CALL EXP(12,12) GO TO 600 604 IF(K.GT.47)GO TO 606 II=ILD2+1 PIC(1,II)=PIC(3,12) PIC(2,II)=BS CALL DASP(DN,1,12,II) ILD2=II CALL EXP(12,12) GO TO 600 606 IF(K.GT.48)GO TO 607 CALL EXP(12,12) GO TO 600 607 IF(K.GT.55)GO TO 608 PIC(1,35)=XVS PIC(2,35)=BS PIC(3,35)=X PIC(4,35)=BS CALL EXP(29,29) GO TO 618 608 IF(K.GT.56)GO TO 618 DO 610 J=1,4 610 PIC(J,29)=PIC(J,12) PIC(3,35)=X GO TO 618 600 II=ILD2+1 PIC(1,II)=XVS+0.1 GO TO 619 618 II=ILD2+1 PIC(1,II)=PIC(3,35)+0.1 619 PIC(2,II)=BS CALL DASP(DN,3,32,II) ILD2=II I2END=ILD2 CALL ROTP(4) 161 CONTINUE I2STRT=11 IAR2=IA CALL ROTP(96) GO TO (15,170),IBX 170 CS=COS(-6.0*RD) SN=SIN(-6.0*RD) GO TO 182 180 CS=COS(+6.0*RD) SN=SIN(+6.0*RD) ICX=2 182 BS=PIC(4,12) DO 171 K=1,24 I2END=11 CALL RORP(CS,SN,-2.0,2.0) I2END=ILD2 X=PIC(3,11) Y=PIC(4,11) PIC(3,12)=X PIC(3,13)=Y PIC(4,13)=BS PIC(4,12)=BS PIC(1,13)=XP PIC(2,13)=YP PIC(1,12)=XV PIC(2,12)=YV IS=ILD1+1 CALL EXP(IS,ILD2) II=ILD1 CALL PROJ(2,20.0,11,12,13,3,II) I2END=II DN=10.0 II=II+1 PIC(1,II)=PIC(3,35)+0.1 PIC(2,II)=BS CALL DASP(DN,3,32,II) I2END=II ILD2=II IF(K.GT.7) GO TO 173 GO TO 179 173 IF(K.GT.8) GO TO 174 CALL ROTP(48) GO TO 179 174 IF(ICX.EQ.2) GO TO 172 II=ILD2+1 PIC(1,II)=PIC(3,13) PIC(2,II)=BS CALL DASP(DN,1,13,II) ILD2=II CALL EXP(13,13) GO TO 179 172 II=ILD2+1 PIC(1,II)=PIC(3,12) PIC(2,II)=BS CALL DASP(DN,1,12,II) ILD2=II CALL EXP(12,12) 179 I2END=ILD2 CALL ROTP(3) 171 CONTINUE DO 175 J=1,4 175 PIC(J,11)=PEC(J,11) PIC(1,12)=XV PIC(2,12)=YV PIC(3,12)=PIC(3,11) PIC(4,12)=BS PIC(1,13)=XP PIC(2,13)=YP PIC(3,13)=PIC(4,11) PIC(4,13)=BS IS=ILD1+1 CALL EXP(IS,ILD2) II=ILD1 CALL PROJ(0,20.0,11,12,13,3,II) I2END=II ILD2=II CALL ROTP(48) IF(ICX.EQ.2) ICX=1 GO TO (15,130,140),IBX 220 GO TO 25 290 GO TO 25 294 GO TO 25 260 DO 264 J=1,4 PIC(J,29)=PIC(J,11) 264 PIC(J,30)=PIC(J,13) XHS=PIC(3,11) BHS=PIC(4,11) XH=PIC(1,11) YH=PIC(2,11) XP=PIC(1,13) YP=PIC(2,13) YPS=PIC(4,13) IK=384/6 CS=COS(+6.0*RD) SN=SIN(+6.0*RD) ILD3=ILD2 IA=36 DO 261 K=1,IK IA=IA+1 IS=ILD1+1 CALL EXP(IS,ILD3) PIC(1,IA)=PIC(3,12) PIC(2,IA)=PIC(4,12) I2END=12 I2STRT=12 CALL RORP(CS,SN,-2.0,-2.0) I2END=ILD3 I2STRT=11 PIC(3,IA)=PIC(3,12) PIC(4,IA)=PIC(4,12) X=PIC(3,12) Y=PIC(4,12) PIC(3,11)=X PIC(4,13)=Y PIC(3,13)=BHS PIC(4,11)=BHS PIC(2,11)=YH PIC(2,13)=YP PIC(1,13)=XP PIC(1,11)=XH II=IS DM=20.0 CALL PROJ(34,DM,11,12,13,3,II) ILD2=II DN=10.0 II=II+1 PIC(1,II)=BHS PIC(2,II)=BHS+0.1 IB=13 IF(K.GT.1) IB=32 CALL DASP(DN,3,IB,II) ILD2=II IF(K.GT.7)GO TO 265 PIC(1,31)=XHS PIC(2,31)=BHS PIC(3,31)=X PIC(4,31)=BHS PIC(1,32)=BHS PIC(2,32)=YPS PIC(3,32)=BHS PIC(4,32)=Y CALL EXP(33,36) GO TO 700 265 IF(K.GT.8)GO TO 266 DO 701 J=1,4 701 PIC(J,33)=PIC(J,13) PIC(3,31)=X PIC(4,32)=Y GO TO 700 266 IF(K.GT.16)GO TO 267 II=ILD2+1 PIC(2,II)=Y PIC(1,II)=PIC(3,13) CALL DASP(DN,1,13,II) ILD2=II CALL EXP(13,13) PIC(3,31)=X GO TO 700 267 IF(K.GT.23)GO TO 268 CALL EXP(30,30) PIC(3,31)=X PIC(2,36)=YPS PIC(1,36)=BHS PIC(3,36)=BHS PIC(4,36)=Y GO TO 700 268 IF(K.GT.24)GO TO 269 DO 702 J=1,4 702 PIC(J,34)=PIC(J,11) PIC(3,31)=X PIC(4,36)=Y GO TO 700 269 IF(K.GT.39)GO TO 703 PIC(4,36)=Y II=ILD2+1 PIC(1,II)=PIC(3,11) PIC(2,II)=BHS CALL DASP(DN,1,11,II) ILD2=II CALL EXP(11,11) GO TO 700 703 IF(K.GT.40)GO TO 704 DO 705 J=1,4 705 PIC(J,30)=PIC(J,13) PIC(4,36)=Y II=ILD2+1 PIC(1,II)=PIC(3,11) PIC(2,II)=BHS CALL DASP(DN,1,11,II) ILD2=II CALL EXP(11,11) GO TO 700 704 IF(K.GT.47)GO TO 706 II=ILD2+1 PIC(1,II)=PIC(3,11) PIC(2,II)=BHS CALL DASP(DN,1,11,II) ILD2=II CALL EXP(11,11) GO TO 700 706 IF(K.GT.48)GO TO 707 CALL EXP(11,11) GO TO 700 707 IF(K.GT.55)GO TO 708 PIC(1,35)=XHS PIC(2,35)=BHS PIC(3,35)=X PIC(4,35)=BHS CALL EXP(29,29) GO TO 700 708 IF(K.GT.56)GO TO 700 DO 710 J=1,4 710 PIC(J,29)=PIC(J,11) PIC(3,35)=X 700 II=ILD2+1 IB=29 IF(K.GT.48) IB=35 PIC(1,II)=BHS+0.1 PIC(2,II)=BHS CALL DASP(DN,3,IB,II) ILD2=II I2END=ILD2 CALL ROTP(4) 261 CONTINUE IAR2=IA CALL ROTP(96) GO TO (25,170),IBX ICX=2 270 CS=COS(6.0*RD) SN=SIN(6.0*RD) GO TO 282 280 CS=COS(-6.0*RD) SN=SIN(-6.0*RD) IF(IBX.EQ.2) IBX=3 282 DO 275 J=1,4 PIC(J,1)=PIC(J,11) PIC(J,2)=PIC(J,12) 275 PIC(J,3)=PIC(J,13) DN=20.0 DO 271 K=1,24 I2END=12 I2STRT=12 CALL RORP(CS,SN,-2.0,-2.0) I2STRT=11 X=PIC(3,12) Y=PIC(4,12) DO 276 J=1,4 PIC(J,11)=PIC(J,1) 276 PIC(J,13)=PIC(J,3) PIC(3,11)=X PIC(4,13)=Y IS=ILD1+1 CALL EXP(IS,ILD2) II=ILD1 CALL PROJ(34,20.0,11,12,13,3,II) II=II+1 PIC(1,II)=BHS+0.1 PIC(2,II)=BHS CALL DASP(DN,3,35,II) II=II+1 PIC(1,II)=BHS PIC(2,II)=BHS+0.1 CALL DASP (DN,3,32,II) ILD2=II DN=10.0 DN=10.0 IF(K.GT.7) GO TO 273 GO TO 279 273 IF(K.GT.8) GO TO 274 CALL ROTP (24) GO TO 279 274 IF(ICX.EQ.2) GO TO 272 DN=10.0 II=ILD2+1 PIC(1,II)=XP PIC(2,II)=YP PIC(3,13)=BHS PIC(4,13)=Y CALL DASP (DN,3,13,II) CALL EXP(13,13) ILD2=II GO TO 279 272 II=ILD2+1 PIC(1,II)=XH PIC(2,II)=YH CALL DASP(DN,3,11,II) CALL EXP(11,11) ILD2=II 279 I2END=ILD2 CALL ROTP(3) 271 CONTINUE CALL ROTP(48) DO 277J=1,4 PIC(J,11)=PIC(J,1) PIC(J,12)=PIC(J,2) 277 PIC(J,13)=PIC(J,3) IS=ILD1+1 CALL EXP(IS,ILD2) II=ILD1 CALL PROJ(34,20.0,11,12,13,3,II) II=II+1 PIC(1,II)=BHS+0.1 PIC(2,II)=BHS CALL DASP(DN,3,35,II) II=II+1 PIC(1,II)=BHS PIC(2,II)=BHS+0.1 CALL DASP (DN,3,32,II) ILD2=II I2END=ILD2 CALL ROTP(48) GO TO(25,280,230),IBX 210 DO 211 J=1,4 PIC(J,11)=PEC(J,3) PIC(J,12)=PEC(J,1) 211 PIC(J,13)=PEC(J,4) CALL ROTP(24) II=100 CALL PROJ(0,20.0,11,12,13,1,II) ILD1=II CALL PROJ(0,20.0,11,12,13,3,II) ILD2=II I2END=II CALL ROTP(96) GO TO (25,260),IBX 230 CS=COS(+6.0*RD) SN=SIN(+6.0*RD) IF(IBX.EQ.3) IBX=2 IK=24 CALL EXP(13,ILD2) DO 232 I=14,17 DO 232 J=1,4 232 PIC(J,I)=PEC(J,I) GO TO 231 240 CS=COS(-3.0*RD) SN=SIN(-3.0*RD) IK=16 IF(IBX.EQ.2) IBX=3 231 II=45 DN=10.0 PIC(1,II)=PIC(1,12) PIC(2,II)=PIC(2,12) CALL DASP(DN,1,11,II) II=II+1 PIC(1,II)=PIC(3,12) PIC(2,II)=PIC(4,12) CALL DASP(DN,3,11,II) ILD1=II I2END=II CALL ROTP(48) DO 233 J=1,4 PIC(J,13)=PIC(J,12) 233 PIC(J,18)=PIC(J,11) IA=19 IS=ILD1+1 ILD2=IS PIC(1,19)=XHS PIC(2,19)=BHS PIC(4,19)=BHS DO 234 K=1,IK IA=IA+1 PIC(1,IA)=PIC(3,12) PIC(2,IA)=PIC(4,12) I2STRT=12 I2END=12 CALL RORP(CS,SN,-2.,-2.) I2STRT=11 PIC(3,IA)=PIC(3,12) PIC(4,IA)=PIC(4,12) PIC(3,11)=PIC(3,12) PIC(3,19)=PIC(3,12) CALL EXP(IS,ILD2) II=ILD2+1 PIC(1,II)=PIC(3,12) PIC(2,II)=PIC(4,12) CALL DASP(DN,3,11,II) ILD2=II I2END=II CALL ROTP(4) 234 CONTINUE CALL ROTP(96) DO 235 J=1,4 PIC(J,11)=PIC(J,18) 235 PIC(J,12)=PIC(J,13) CALL EXP(18,IA) CALL EXP(IS,ILD2) GO TO (25,240,25),IBX END