File: CONT2.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1

(Source file text) 

---

	SUBROUTINE CONT(IMP,R,T,IANZ,AVERAG,TO1,DELTA))

C	SPECIFIC HEAT/CONT./CALCULATION

	DIMENSION RRR(25)
	DIMENSION X(10)
	DIMENSION IANZ(100)
	DIMENSION AVERAG(100)
	DIMENSION R(100)
	DIMENSION T(100)
	DIMENSION TEMP(100)

	WRITE(4,100)
	READ(4,110)M1
	WRITE(4,120)
	READ(4,110) N1
	WRITE(4,130)
	READ(4,110) IMAX

	REWIND 5
	K=1
900	READ(5,200) IEND,RRR(K)
	K=K+1
	IF(IEND) 800,800,900
800	KM=K-1
	DO 1000 I=1,IMP
        IK=1
714	IF(R(I)-RRR(IK)) 711,712,712
712	IF(R(I)-RRR(IK+1)) 710,710,711
711	IK=IK+1
	IF(IK.GT.KM) GOTO 730
	GOTO 714

710	M1=IK-1
	N1=IK+2
	IF(M1-1) 721,721,722
721	M1=1
722	IF(N1-KM) 723,723,724
723	N1=KM
724	CALL RFIT(IMAX,M1,N1,X)
C	WRITE(4,150)
	CALL KELVIN(IMAX,X,R(I),TEMP(I))
	WRITE(4,170) R(I),TEMP(I)
	GOTO 1000
730	WRITE(4,210)
1000	CONTINUE

C	WRITE(3,180)
	L= IMP-1
	WRITE(4,180)
	DO 1100  I=1,L
	SLOPE=(TEMP(I)-TEMP(I+1))/(T(I)-T(I+1))
	SL1=1./SLOPE
	TSTR=(TEMP(I)+TEMP(I+1))/2.
	WRITE(4,190) TSTR,SLOPE,SL1
	CALL AV(IANZ,AVERAG,TO1,DELTA,TSTR,SL1)
C	WRITE(3,190) TSTR,SLOPE,SL1
1100	CONTINUE

100	FORMAT(/1X,'ERSTES PUNKTEPAAR FUER T-FIT:',T60,$)
110	FORMAT(I3)
120	FORMAT(/1X,'LETZTES PUNKTEPAAR FUER T-FIT:',T60,$)
130	FORMAT(/1X,'ANZAHL TERME (T-FIT):',T60,$)
140	FORMAT(1X,I3)
150	FORMAT(1X,T31,'T(K)',T18,'R(KOHM)')
160	FORMAT(1X,F8.4,E16.5)
170	FORMAT(1X,F20.4,E20.5)
180	FORMAT(1X,'       T(K)      TEMP.DERIVATIVE   TIME DERIV.')
190	FORMAT(1X,3E16.5)
200	FORMAT(1X,I5,E16.9)
210	FORMAT(1X,'KEIN FIT MOEGLICH,SR CONT')
	RETURN
	END

---