from NanoTCAD_ViDES import *
from module_UTB import *

Lx=1;
Ly=2;
a=0.5;

nx=int(Lx/a)+1;
ny=int(Ly/a)+1;

device=UTB(nx,ny);
device.a=0.5;

device.Eupper=2;
device.Elower=0;

device.compute_charge_T()
a=[device.E,device.T]
savetxt("T.out",transpose(a));