clear;

i_anode = pinch(getdata("CHARGE","ac_anode.dI"));
i_cathode = pinch(getdata("CHARGE","ac_cathode.dI"));
f = pinch(getdata("CHARGE","ac_cathode.f"));
Vdc = pinch(getdata("CHARGE","ac_cathode.V_cathode"));

kern = size(i_anode);

Nv = kern(1);
Nf = kern(2);

v = 0.001;

y_anode = i_anode/v;
y_cathode = i_cathode/v;

w = 2*pi*f;

A = 5e-6 * getnamed('CHARGE','norm length') * 1e12;  # um^2

C_anode = -imag(y_anode)/w*1e15/A;  # fF/um^2
C_cathode = imag(y_cathode)/w*1e15/A;  # fF/um^2

plot(f/1e9,C_anode,C_cathode,"Freq. (GHz)","Capacitance (fF/um^2)","","log10x");