File: TOTAL.FT of Disk: Disks/MyPDP/m8-backup-rka1-rkb1

(Source file text) 

---

C	WAERMELEITUNG IM NORMALEN ZUSTAND MIT BERUECK-
C	SICHTIGUNG VON ELEKTRON-PHONON, ELEKTRON-DEFEKT,
C	PHONON-ELEKTRON UND PHONON-DEFEKT STREUUNG.
	WRITE(4,1000)
	READ(4,1100)TE
1000	FORMAT(1X,'ANFANGSTEMPERATUR',T40,$)
1100	FORMAT(E16.8)
	WRITE(4,1200)
	READ(4,1300)DETE
1200	FORMAT(1X,'TEMPERATURSCHRITTE?',T40,$)
1300	FORMAT(E16.8)
	WRITE(4,1400)	
	READ(4,1500)ALF
1400	FORMAT(1X,'ALPHA=?',T40,$)
1500	FORMAT(E16.8)
	WRITE(4,1550)
	READ(4,1600)BETA
1550	FORMAT(1X,'BETA=?',T40,$)
1600	FORMAT(E16.8)
	WRITE(4,1700)
	READ(4,1800)YK
1700	FORMAT(1X,'BOUNDARY SCATTERING=?',T40,$)
1800	FORMAT(E16.8)
	WRITE(4,1900)
	READ(4,2000)CY
1900	FORMAT(1X,'PONON-ELEKTRONST.=?',T40,$)
2000	FORMAT(E16.8)
	WRITE(4,2100)	
	READ(4,2200)GY
2100	FORMAT(1X,'PUNKTDEFEKTSTR.=?',T40,$)
2200	FORMAT(E16.8)
	WRITE(3,2300)ALF
2300	FORMAT(1X,'ALPHA=',T40,E16.8)
	WRITE(3,2400)BETA
2400	FORMAT(1X,'BETA=',T40,E16.8)
	WRITE(3,2500)YK
2500	FORMAT(1X,'BOUNDARY.PAR=',T40,E16.8)
	WRITE(3,2600)CY
2600	FORMAT(1X,'PHONON-ELEKT.=',T40,E16.8)
	WRITE(3,2700)GY
2700	FORMAT(1X,'PUNKTDEF.PARAM.=',T40,E16.8)
	WRITE(3,2800)
2800	FORMAT(1H0,'TEMPERATUR',T21,'ELEKTR.WAERMEL.',T41,
     $	'PHONON.WAERM.',T61,'TOTALE WAERMEL.',/)
100	TE=TE+DETE
	ELEKT=((TE**3*ALF+BETA)/TE)**(-1.)
	AH=0.
	DX=0.1
	PHON=0.
200	AH=AH+DX
	PHON=PHON+AH**4.*EXP(AH)*TE**3./((EXP(AH)-1.)**2.*
     $	(YK+CY*TE*AH+TE**4.*GY*AH**4.))
	IF(AH.LE.10.)GOTO200
	PHON=PHON*DX
	TOTAL=PHON+ELEKT
	WRITE(3,2900)TE
2900	FORMAT(E16.8)
	WRITE(3,3000)ELEKT
3000	FORMAT(1H+,T20,E16.8)
	WRITE(3,3100)PHON
3100	FORMAT(1H+,T40,E16.8)
	WRITE(3,3200)TOTAL
3200	FORMAT(1H+,T61,E16.8)
	IF(TE.LE.50)GOTO100
	STOP
	END

---