length = 0.5e-6 
nLength = 32*2^3
dx = length/nLength 
--length = 2.2
nLength_bnd = 2*nLength+4
--nLength = 200
--dx = length_bnd/(2^level)
--dx_ini = length/nLength
level = math.ceil(math.log(nLength_bnd)/math.log(2))
length_bnd = (2^level)*dx
--level = 9
--dx = length_bnd/2^level
--nLength = math.ceil((length)/dx)
dx_half = dx*0.5
zpos = dx_half
sph_pos = {0.2,0.2,zpos}
radius = 0.05
Dia = radius*2.0
dx_eps = length_bnd/2^20

dt =dx*dx
ref_pot = 10e-3--25e-3
permit = 6.95e-10
moleDens0 = 1e-1
valence_sqr = 1
charge = 1.60217657e-19
faraday = 96485.3365
gasConst = 8.3144621
k_b = 1.3805e-23
temp = 273
N_A = 6.02e23 
