File: ABFALL.FT of Disk: Disks/MyPDP/m8-1-rka1-rkb1
(Source file text)
C DECAYMEASUREMENT AND CALCULATION DIMENSION R(100) DIMENSION T(100) DIMENSION A(5) LOGICAL GOON IVAL=9999 IERR=50 IRASC=1 MODESC=2 MODDV=1+8 IRADV=0+4 WRITE(4,1700) READ(4,900) CONST WRITE(4,300) READ(4,900) TSC WRITE(4,200) READ(4,900) STEP 9 K=1 TARGET=0 WRITE(4,1800) READ(4,1200) GOON IF(.NOT.GOON) GOTO 6 WRITE(4,1700) READ(4,900) CONST 6 WRITE(4,800) READ(4,1200) GOON IF(.NOT.GOON) GOTO 7 WRITE(4,700) READ(4,900) STEP 7 WRITE(4,600) READ(4,1200) GOON IF(.NOT.GOON) GOTO 13 WRITE(4,300) READ(4,900) TSC 13 WRITE(4,2400) READ(4,2200) IMP WRITE(4,2000) READ(4,1200) GOON IF(.NOT.GOON) GOTO 19 21 WRITE(4,2100) READ(4,2200) IMAX CALL BASET(IMAX,RO,MODDV,IRADV,IERR,CONST) WRITE(4,2300) RO WRITE(4,2500) READ(4,1200) GOON IF (.NOT.GOON) GOTO 14 22 WRITE(4,2100) READ(4,2200) IMAX1 CALL BASET(IMAX1,RR,MODDV,IRADV,IERR,CONST) WRITE(4,2600) RR WRITE(4,2700) READ(4,2200)L IF (L.EQ.1) GOTO 21 IF (L.EQ.2) GOTO 22 14 WRITE(4,400) READ(4,1200) GOON IF(.NOT.GOON) GOTO 20 8 WRITE(4,3700) IOUT=DINP(1,1,IERRDI,INPUT) IF(IOUT.EQ.0)GOTO 8 OUT=SCINP(MODESC,IRASC,IERR,TARGET) CALL DVINP(MODDV,IRADV,IERR,IVAL) XZERO=OUT TIME=0 RT=FLOAT(IVAL)*CONST T(K)=TIME R(K)=RT IF (K.EQ.1) GOTO 100 120 K=K+1 TARGET=XZERO+STEP*FLOAT(K) IERR=50 OUT=SCINP(MODESC,IRASC,IERR,TARGET) CALL DVINP(MODDV,IRADV,IERR,IVAL) XEFF= OUT+TARGET TIME=(XEFF-XZERO)*TSC RT=FLOAT(IVAL)*CONST R(K)=RT T(K)=TIME 100 WRITE(4,1900) TIME,RT IF (K-IMP) 120,110,110 110 WRITE(4,1100) READ(4,1200) GOON IF(.NOT.GOON) GOTO 19 C EVALUATION CONTINUE DO 140 I=1,5 140 A(I)=0 DO 130 J=1,IMP A(1)=A(1)+T(J) Y=ALOG(ABS(RO-R(J))) A(5)=A(5)+Y*Y A(2)=A(2)+T(J)*T(J) A(3)=A(3)+Y*T(J) A(4)=A(4)+Y 130 CONTINUE ANZ=1./FLOAT(IMP) DET=A(2)-A(1)*A(1)*ANZ TAU1= (A(3)-ANZ*A(1)*A(4))/DET B=A(3)-ANZ*A(1)*A(4) C=A(2)-A(1)*A(1)*ANZ D=A(5)-A(4)*A(4)*ANZ RQUA=B*B/(C*D) TAU=1./TAU1 WRITE(4,1300) WRITE(4,2300)RO WRITE(4,2600)RR WRITE(4,1400)TAU WRITE(4,1500) RQUA 19 WRITE(4,1600) READ(4,1200) GOON IF (GOON) GOTO 9 3700 FORMAT(1H+,T70,'') 200 FORMAT(/1X,'STEP= ',T60,$) 300 FORMAT(/1X,'TIME SCALE (SEC/VOLTS): ',T60,$) 400 FORMAT(1X,'START OF MEASUREMENT (T/F):',T60,$) 500 FORMAT(/1X,'NEW TIME SCALE (T/F):',T60,$ 600 FORMAT(/1X,'NEW TIME SCALE WANTED? (T/F):',T60,$) 700 FORMAT(/1X,'NEW STEP :',T60,$) 800 FORMAT(/1X,'STEP CHANGE WANTED ?(T/F):'T60,$) 900 FORMAT(E16.8) 1100 FORMAT(/1X,'CALCULATION WANTED?(T/F)',T60,$) 1200 FORMAT(L1) 1300 FORMAT(//,' RESULT:') 1400 FORMAT(/1X,'DECAYCONSTANT:',T60,E14.5,' SEC') 1500 FORMAT(/1X,'KORRELATION: SQR(R)=',T60,E16.8) 1600 FORMAT(/1X,'ONE MORE POINT? (T/F):',T60,$) 1700 FORMAT(/1X,'DVM-CONSTANT:',T60,$) 1800 FORMAT(/1X,'NEW DVM-CONSTANT WANTED? (T/F):',T60,$) 1900 FORMAT(1X,F8.4,E16.5) 2000 FORMAT(/1X,'START READING BASETEMPERATURE? (T/F):',T60,$) 2100 FORMAT(/1X,'AVERAGE OVER IMAX POINTS:IMAX=:',T60,$) 2200 FORMAT(I3) 2300 FORMAT(/1X,'BASETEMPERATURE-RESISTANCE: RO=',T60,E16.8,' KOHMS') 2400 FORMAT(/1X,'NUMBER OF POINTS (DECAY):',T60,$) 2500 FORMAT(/1X,'MEASUREMENT OF RAISED TEMPERATURE?(T/F):',T60,$) 2600 FORMAT(/1X,'RES.AT RAISED TEMPERATURE: RR=',T60,E16.8,' KOHMS') 2700 FORMAT(/1X,'ONE MORE READING OF RO(1),RR(2) OR GO ON?(O)',T60,$) 2800 FORMAT(/1X,'RO=',T60,$) 20 STOP END