C PROGRAMM SVS - SIMULATION VON SCHAUERENTLADUNGEN
C NACH G.W.MILLS : "THE MECHANISMS OF THE SHOWERING ARC"
C FORTRAN II
	WRITE(1,100)
100	FORMAT('SIMULATION VON SCHAUERENTLADUNGEN'/33('-'))
	READ(1,99)UBATT
99	FORMAT('      VERSORGUNGS-SPANNUNG [V] : 'F5.1)
	READ(1,101)C
101	FORMAT('   SCHALTKREIS-KAPAZITAET [PF] : 'F9.1)
	READ(1,102)FI0
102	FORMAT('        GESCHALTETER STROM [A] : 'F9.4)
	READ(1,103)UMAX
103	FORMAT('     UMAX (GESAMT-ENERGIE) [V] : 'F9.1)
	READ(1,104)V
104	FORMAT('   TRENN-GESCHWINDIGKEIT [M/S] : 'F9.3)
	READ(1,105)FBK
105	FORMAT('SCHMELZBRUECKEN-LAENGE [MMM/A] : 'F8.3)
	READ(1,106)Z
106	FORMAT('          KABEL-IMPEDANZ [OHM] : 'F7.1)
	READ(1,110)R
110	FORMAT('           HF-WIDERSTAND [OHM] : 'F8.3)
	READ(1,107)VM
107	FORMAT('    MINIMALE BOGENSPANNUNG [V] : 'F7.1)
	READ(1,108)FIM
108	FORMAT('      MINIMALER BOGENSTROM [A] : 'F7.1)
	READ(1,109)FB
109	FORMAT(' KRITISCHE FELDSTAERKE [V/MMM] : 'F9.3)
	READ(1,120)DT
120	FORMAT('      ITERATIONS-SCHRITT [MMS] : 'F9.3)
	C=C*1.E-12
	EA0=.5*C*UMAX**2
	DT=DT*1.E-6
	FBK=FBK*1.E-6
	DJ0=FBK*FI0
C INITIATION
	TJ1=0.
	VJ1=0.
	FIJ1=FI0
	EUJ1=0.
	VA1=0.
	VA2=0.
C ITERATION
30	TJ=TJ1+DT
	VJ=VJ1+(FIJ1*DT)/C
	EUJ=EUJ1+(VJ**2-VJ1**2)*C*.5
	DJ=V*TJ+DJ0
C FIELD BREAKDOWN ?
	DB=VJ/FB
	IF(DJ-DB)31,31,31
31	CALL BOSPA(VJ,VA1,VA2,DJ,VM,FIM,K)
	IF(K-1)33,33,33
C AIR BREAKDOWN ?
33	CONTINUE
	WRITE(1,200)TJ,VJ,DJ,VA1,VA2,K
200	FORMAT(5E12.7,'  'I4)
	TJ1=TJ
	VJ1=VJ
	EUJ1=EUJ
	GOTO 30
	END