calculation of head
 &inputph
  trans=.false.
  l_head=.true.
  tr2_ph=1.d-4,
  prefix='si',
  omega_gauss=20.0
  n_gauss=97
  grid_type=5
  second_grid_i=1
  second_grid_n=10
  niter_ph=1
  nsteps_lanczos=30
  outdir='./'
 /
0.0 0.0 0.0