C MAIN COMMON IBUFF(255),DUMI(85),MANU,IBLK,CP,CV,CF COMMON/BLK1/VOLN(100),VOLX(100),FLON(110),FLOX(110),PRSN(100) COMMON/BLK2/ACCN(100),ACCX(100),PRSX(100) COMMON/BLK3/INEG,IPOS,IDIM,ITMP,IDIM1,IPAS,ISTOP,IREC COMMON/BLK4/HOLD(30),K,BUFF(1400),X,Y,APOS,J COMMON/BLK5/AMAT(5,5),COEF(4),COV(5,5),COR(5,5),A21(4) COMMON/BLK6/SD(5),A12(4),CEPT,DET,EPS,SX(5),SXY(5),SY,RATE IBLK=0 IPAS=0 10 WRITE(4,20) 20 FORMAT(' MANU,RATE,CALFLO',$) READ(4,30)MANU,RATE,CF 30 FORMAT(I2,1X,I3,1X,E13.4) CP=0.1487 CV=0.0111 40 CALL BRTH CALL DSPL IF(IPOS.EQ.0)GO TO 100 CALL ACCEL(1,FLON,ACCN,ITMP) CALL ACCEL(2,FLOX,ACCX,IDIM) CALL REGRS(PRSN,VOLN,FLON,ACCN,ITMP) WRITE(3,90)AMAT 90 FORMAT(' AMAT',/(5E13.4)) CALL SIMUL(4,EPS,0) WRITE(3,90)AMAT WRITE(3,80)(PRSN(I),VOLN(I),FLON(I),ACCN(I),I=1,ITMP) 80 FORMAT(' ',4E13.4) CALL PRIN CALL REGRS(PRSX,VOLX,FLOX,ACCX,IDIM) WRITE(3,90)AMAT CALL SIMUL(4,EPS,0) WRITE(3,90)AMAT WRITE(3,80)(PRSX(I),VOLX(I),FLOX(I),ACCX(I),I=1,IDIM) CALL PRIN 100 IF(IBLK.LT.16)GO TO 40 STOP END