Budget DC Coil Primary Circuit Simulation

.MODEL D1N4007 D(IS=2.55E-9 RS=0.042 N=1.75 TT=5.76E-6 CJO=1.85E-11 VJ=0.75 M=0.333 BV=1000 IBV=9.86E-5)

Xnst      1   0   2              NST15KV
Rnstmeas  1   2                  1GOhms
Xrect     1   2   3   4          HVRECTIFIER1
Csmooth   3   4                  400nF IC=0
Rsmooth   3   4                  25.9MEGOhms

*            +HV GND -HV
.SUBCKT NST15KV 4 2 5
Vtrans1      1  2   SIN(0 10607 60)
Vtrans2      3  2   SIN(0 -10607 60)
R1           1  4   2
R2           3  5   2
.ENDS

*                AC AC DC+ DC-
.SUBCKT HVRECTIFIER1 1 3 5 6
X1 1 2 DiodeString25
X2 3 2 DiodeString25
X3 4 3 DiodeString25
X4 4 1 DiodeString25
R1 2 5 2
R2 4 6 2
.ENDS

.SUBCKT DiodeString25 1 26
D1  1  2  D1N4007
D2  2  3  D1N4007
D3  3  4  D1N4007
D4  4  5  D1N4007
D5  5  6  D1N4007
D6  6  7  D1N4007
D7  7  8  D1N4007
D8  8  9  D1N4007
D9  9  10 D1N4007
D10 10 11 D1N4007
D11 11 12 D1N4007
D12 12 13 D1N4007
D13 13 14 D1N4007
D14 14 15 D1N4007
D15 15 16 D1N4007
D16 16 17 D1N4007
D17 17 18 D1N4007
D18 18 19 D1N4007
D19 19 20 D1N4007
D20 20 21 D1N4007
D21 21 22 D1N4007
D22 22 23 D1N4007
D23 23 24 D1N4007
D24 24 25 D1N4007
D25 25 26 D1N4007
.ENDS

.options vmin=-1e8 vmax=1e8
.print tran vout(Rnstmeas) vout(Csmooth)
.op
.tran .0001 .05 UIC > primary.out
.end