(* Program *)

els = ReadList["elsloss.txt", {Number, Number}];
inels = ReadList["inelsloss.txt", {Number, Number}];
cx = ReadList["cxloss.txt", {Number, Number}];
jr = ReadList["jr.txt", {Number, Number}];
ijr = ReadList["ijr.txt", {Number, Number}];
jdote = ReadList["jdote.txt", {Number, Number}];
ijdote = ReadList["ijdote.txt", {Number, Number}];
jzdotez = ReadList["jzdotez.txt", {Number, Number}];
ijzdotez = ReadList["ijzdotez.txt", {Number, Number}];
fe1 = ReadList["fe1.txt", {Number, Number}];
fe2 = ReadList["fe2.txt", {Number, Number}];
iz = ReadList["iz.txt", {Number, Number}];
exc = ReadList["exc.txt", {Number, Number}];

n = 151;
n1 = 201;
n2 = 51;
delr = 6.66667*10^-5 ;
height = 0.2;
r1 = 0.01;
t1 = 0.000128;
t2 = 0.000129033;
delt = t2-t1;
q = 1.6*10^-19;
area1 = 2*Pi*r1*height;
nc2p = 5*10^5;
eion = 15.76;
eexc = 11.62;
Ez = 73;

volnorm = Pi*delr*delr*height;
Array[vol,n];
vol[1] = 0.25*volnorm;
sumvol = vol[1];
For[i = 2, i < n,i++, 
    ravg = (2*i - 1)/2;
    vol[i] = volnorm*ravg*ravg - sumvol;
    sumvol = sumvol + vol[i];
];
vol[n] = volnorm*n*n - sumvol;

sumPel = Sum[vol[j]*els[[j,2]], {j,1,n}]
sumPie = Sum[vol[j]*inels[[j,2]], {j,1,n}]
sumPcx = Sum[vol[j]*cx[[j,2]], {j,1,n}]
sumjdote = Sum[vol[j]*jdote[[j,2]], {j,1,n}]
sumijdote = Sum[vol[j]*ijdote[[j,2]], {j,1,n}]
sumjzdotez = Sum[vol[j]*jzdotez[[j,2]], {j,1,n}]
sumijzdotez = Sum[vol[j]*ijzdotez[[j,2]], {j,1,n}]
sumiz = Sum[vol[j]*iz[[j,2]], {j,1,n}]
sumexc = Sum[vol[j]*exc[[j,2]], {j,1,n}]

sume1 = Sum[fe1[[j,1]]*fe1[[j,2]], {j,1,n1}]
sumfe1 = Sum[fe1[[j,2]], {j,1,n1}]
sume2 = Sum[fe2[[j,1]]*fe2[[j,2]], {j,1,n2}]
sumfe2 = Sum[fe2[[j,2]], {j,1,n2}]

pcreated = sumiz*delt/nc2p;
pelost = sumfe1
pilost = sumfe2
pelost2 = jr[[n,2]]*area1*delt/(q*nc2p)
pilost2 = ijr[[n,2]]*area1*delt/(q*nc2p)

Pe1loss = nc2p*q*sume1/delt
Pe2loss = nc2p*q*sume2/delt

izloss = sumiz*eion*q
excloss = sumexc*eexc*q
ieloss = izloss + excloss
ielossb = sumPie

aveIz = sumjzdotez/(Ez*height)

CollLoss = sumPel + sumPcx + sumPie
WallLoss = Pe1loss+Pe2loss
PowInput = sumjzdotez + sumijzdotez

Print["done"];