 minrefine = 1
 maxrefine = 3
 length = 1.0 
 nLength = 32
 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
