****************************************************
.SUBCKT D3a-Triode 1 2 3; A G C;
*      Extract V3.000
* Model created: 10-Jun-2020
*
*
* Curves traced and model developed by Ale Moglia
* (c) 2020 Bartola Valves
* www.bartola.co.uk/valves email: valves@bartola.co.uk
*
*
X1 1 2 3 TriodeK MU= 77.92 EX=1.515 KG1=  25.7 KP= 587.2 KVB= 2940. RGI=2000
+ CCG=6.7P  CGP=3.3P CCP=1.0P  ;
.ENDS
 
****************************************************
.SUBCKT TriodeK 1 2 3; A G C
*
* NOTE: LOG(x) is base e LOG or natural logarithm.
* For some Spice versions, e.g. MicroCap, this has to be changed to LN(x).
*
E1 7 0 VALUE=
+{V(1,3)/KP*LOG(1+EXP(KP*(1/MU+V(2,3)/SQRT(KVB+V(1,3)*V(1,3)))))}
RE1 7 0 1G
G1 1 3 VALUE={0.5*(PWR(V(7),EX)+PWRS(V(7),EX))/KG1}
RCP 1 3 1G    ; TO AVOID FLOATING NODES IN MU-FOLLOWER
C1 2 3 {CCG}  ; CATHODE-GRID
C2 2 1 {CGP}  ; GRID-PLATE
C3 1 3 {CCP}  ; CATHODE-PLATE
D3 5 3 DX     ; FOR GRID CURRENT
R1 2 5 {RGI}  ; FOR GRID CURRENT
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N)
.ENDS TriodeK