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

(Source file text) 

---

C	BERECHNET ABFALLKONSTANTE
	COMMON X
	WRITE(4,599) 
599	FORMAT(X,/,' ANZAHL TERME (I3)')
	READ(4,600) IM
600	FORMAT(X,I3)
	CONST = 1.E-09
	CALL RFIT(IM)
	WRITE(4,2010)
	READ(4,2011) XSCALE
	WRITE(4,2012)
	READ(4,2011) YSCALE
	WRITE(4,2013)
	READ(4,2011) RB
	WRITE(4,2014)
	READ(4,2011) RR
	CALL KELVIN(IM,X,RB,TO)
C	TO: BASE TEMPERATURE
	CALL KELVIN(IM,X,RR,TR)
C	TR: RAISED TEMPERATURE
	DELTA= TR-TO
	REAL LN
	WRITE(4,2030)
	READ(4,2011) ATGEW
	WRITE(4,2031)
	READ(4,2011) SM
	RATIO=ATGEW/SM

	DIMENSION A(5)
	REWIND 6
	DO 102 I=1,5
102	A(I)=0
100	READ(6,2015) K,XA,Y
	IF(K) 101,103,101
101	T = XA*XSCALE
	RATT=Y*YSCALE+RR
	CALL KELVIN(IM,X,RATT,TAT1)
	A(1)=A(1) +T
C	SUMME TI
	DATT=TAT1-TO
C	DELTA AT T
	LN=ALOG(DATT)
	A(5) = A(5) + LN*LN
	A(2) = A(2) + T*T
	A(3) = A(3) + LN*T
	A(4) = A(4) + LN 
	J=K
	GOTO 100

C	BERECHNE TAU,A,B,DELTA

103	ANZ = 1./FLOAT(J)
	DET= A(2) - A(1)*A(1)*ANZ
	TAU1 = (A(3) -ANZ*A(1)*A(4))/DET
	DFIT = EXP(A(4)*ANZ - TAU1*A(1)*ANZ)
	B = A(3) - ANZ*A(1)*A(4)
	C= A(2) - A(1)*A(1)*ANZ
	D=A(5) - A(4)*A(4)*ANZ
	RQOU = B*B/(D*C)
	WRITE(4,2016)
	READ(4,2011) WID
	WRITE(4,2017)
	READ(4,2011) STROM
	CO=CONST*STROM*STROM*WID/DELTA
	TAU=-1./TAU1
	SPEZ= TAU*CO*RATIO
	TEMP=(TR+TO)/2.

C	OUTPUT

	WRITE(4,2023) TEMP
	WRITE(4,2020) DELTA
	WRITE(4,2021) DFIT
	WRITE(4,2022) CO
	WRITE(4,2024) TAU
	WRITE(4,2025) RQOU
	WRITE(4,2026) SPEZ


2010	FORMAT(X,/,' X - SCALE: (SEC/MM) ')
2012	FORMAT(X,/,' YSCALE: (KOHM/MM) ')
2013	FORMAT(X,/,' RESISTANCE AT BASE TEMPERATURE (KOHM)')
2014	FORMAT(X,/,' RESISTANCE AT RAISED TEMPERATURE (KOHM)')
2016	FORMAT(X,/,' HEATER (KOHM) =')
2017	FORMAT(X,/,' HEATERCURRENT (MICROAMPS) =')
2011	FORMAT(X,E16.6)
2015	FORMAT(X,I3,2E16.6)
2023	FORMAT(X,/,30H MITTLERE TEMPERATUR         = E14.4,'  K')
2021	FORMAT(X,/,30H DELTA AUS FIT FUER TAU      = E14.4,'  K')
2025	FORMAT(X,/,30H STANDARTABWEICHUNG          = E14.4)
2020	FORMAT(X,/,30H DELTA-T AUS TEMP.-FIT       = E14.4,'  K')
2022	FORMAT(X,/,30H THERMAL CONDUCTIVITY        = E14.4,'  W/K')
2024	FORMAT(X,/,30H TIME CONSTANT               = E14.4,'  SEC')
2026	FORMAT(X,/,30H SPEZ.HEAT                   = E14.4,'  J/(MOL*K)')
2030	FORMAT(X,/,' MOLEKULARGEWICHT (G) =')
2031	FORMAT(X,/,' MASSE DER PROBE (G) =')
	STOP
	END

---