41                     nx = ny = nz
-----------------------------------
5.d0                    freq (omega=2*pi*freq)
-200.                    qual_fac (if <0., quality_factor=infinity)
-20.                      qual_fac2(if <0., artificial_quality_factor=infinity)
1                       idautodamp(if=1,auto, else,choose above)
3                       ncc: choose the coef c(x,y)
1.0d-4                  tol: tolerance
1.0d-2                  tolCGNR: tolerance for CGNR
5000                    itmax
-----------------------------------
1                       level
5                       npcol: number of columns in print out(<=5)
8                       inorm: method of norm(2:L2-norm, else:L8-norm)
0                       iddirac(0:ananytic solv, 1:point source)
1.0d0                   srcmag(the magnitute of the source $f$)
-----------------------------------
0                       idchdmn: change domain(0:no, else:changed by below)
0.0d0                   ax
1.0d0                   bx
0.0d0                   ay
1.0d0                   by
0.0d0                   az
1.0d0                   bz
