#Select connectivity matrix 
C = C;
#Select x, y, z coordinates 
nvtx = 18905 ;
x = x;
x = reshape(x, [18905, 1]); 
y = y;
y = reshape(y, [18905, 1]); 
z = z;
z = reshape(z, [18905, 1]); 
DSB_test = unstructureddataset("DSB_test", x, y, z, C) ; 
#Add default parameter 
DSB_test.addparameter("a", 0) ;
#Add attributes 
N = reshape(N, [18905,1,1]) ;
DSB_test.addattribute("N",N) ; 
#Create a structure 
extractstructure(DSB_test) ; 
set("name", "silicon"); 
set("material", "Si (Silicon)"); 
#Create doping imports 
N_ds = unstructureddataset("N", x, y, z, C); 
N_ds.addparameter('v',0); 
N_ds.addattribute("N",1e+06*N) ; 
addimportdope; 
set('x',0); 
set('y',0); 
set('z',0); 
importdataset(N_ds); 
set('name','doping'); 
set('dopant type','n');