.SUBCKT MEC201 PLUS MINUS c=0.001 soc=1 ns=1 tscale=60
E_Erate N8 0 VOL='I(V_batt)/C'
E_Ens N5 N1 N4 N1 'NS'
C_CT_S N7 N6 '(Cts/C)/Tscale'
E_Eloss SOC N3 PWL(1) rate 0
+0.33,-0.1 1.65,0.08 3.3,0.16 16.5,0.3 25,0.4 40,0.51
R_Rself SOC0 0 10MEG
C_CT_L PLUS N7 '(Ctl/C)/Tscale'
V_batt MINUS N1 0Vdc
E_Esoc SOC 0 VOL='V(SOC0,0)' MIN=-0.1 MAX=1
R_RT_L PLUS N7 'NS*Rtl/C'
R_RS N6 N5 'NS*Rs/C'
R_R1 N8 RATE 1
X_U1 N3 N4 N1 ETABLE
C_C_capacity SOC0 0 '(3600*C)/Tscale'
.IC V(SOC0)='SOC'
R_RT_S N7 N6 'NS*Rts/C'
C_C1 RATE 0 '3/Tscale'
G_G1 SOC0 0 CUR='I(V_batt)'
.Param rs=0.027 rtl=0.00095 ctl=15000 rts=0.00035 cts=1800
.ENDS MEC201
.SUBCKT ETABLE 1 2 3
R1 1 0 10MEG
E1 2 3 PWL(1) 1 0
+0.0010,0.100 0.0024,1.600 0.0037,2.700 0.0052,3.220
+0.0061,3.360 0.0075,3.480 0.0098,3.600 0.0130,3.700
+0.0185,3.785 0.0290,3.840 0.0400,3.860 0.0600,3.885
+0.1000,3.893 0.2000,3.900 0.3000,3.905 0.4000,3.912
+0.5000,3.925 0.6000,3.945 0.7000,3.973 0.8000,4.000
+0.9000,4.030 1.0000,4.060 1.1000,4.094
.ENDS