lmf-MPI for sigm. it sometimes stop. Maybe because HAMindex
SYMOPS. lmfgw-MPI is not yet---> job divider for sugw.F and how to
divide jobs. Modification on lmf2gw is also required (or something
else. How to reform output, supplied to GW calculations?).

NOTE: extension is proc dependent because
  lmv7.F:  print *,' xxxx file extension: ',procid,ext(1:i+1)
BUGFIX. We removed qq mechanism.

This is simple and numerically correct.
Next, we need to revive local field effect again.
This was clear for the case of Si100H sample.

----------------
We were using |E_\mu(q=0)> for |E_\mu(Q0P)>.
The basis of the Coulomb interaction. Here we discuss about \mu=1
basis for largest eigenvelue of the Coulomb matrix (1/q^2).
(E_\mu are defined in C.Freidlich et at PRB81 125102.)

These E_\mu=1 basis should be almost look like exp(i q \bfr).
But "how the basis |E_\mu(q=0)> is dependent on q" can cause problem.
Instead of analysing it in detail, we move to safe algolism (not use
q=0 basis for Q0P points).
bugfix for hx0fp0.m.F
bug fix for alp/=1.
tests for W(0) are included. (but commented out).
temporay version. test LFC included for dielectric. (note screeenfac
in switches.F). Full inversion to obtain eps(with LFC).

funciton in hx0fp0_sc mode.
TEST code for block inversion to get eps(k=0) with LFC. But
this gives not so good agreement with EPS*.dat
this works OK. testgw.py pass but slight difference because
we do full matrix inversion for epsion with LFC (llw and llwI in hx0fp0_sc).
gwsc1shot_diag is given. Z=1 diagonal only one-shot GW.
minor changes. gwsc1shot_mpi_diag. W0W0I file hx0fp0.sc.m.F
W0W0I modification to WVR.0001 WVI.0001 (1011 mode for hx0fp0)
selectqp (NEED QPNT in GWinput) for 100* mode of hsfp0_sc
New getwklm.
nqxx ->nq for gwsc1shot_mpi_diag
WQFAC is included. test version.
bugfix wbzz added.
convergence test. A Special test code for Si100H slab case.
lxklm=6 version. WQFAC removed.

WQFAC(averaged largest Coulomb in microcell) did not work well
->commented out. Not smooth as functions of Nz (=n3 of n1n2n3).

WKlm (spherical weight Klm in Computer Phys. Comm 176(1007)1-13)
is detemined by auxially function. It is for small NxNx (e.g. 333
supercel for Si100h model).
lxklm=0 only version. LFC@Gamma=F (GWinput) as default.

This corresponds to isotropic case in Computer Phys. Comm 176(1007)1-13.
See Fig.1 in it.
I think the anisotroipc treatment in the paper can be problematic.
lxklm=6 (for inversion of epsilon), but set spherical weight at gamma

wklm(10:)=0 (mkqg.F). Default LFC@Gamma=T (GWinput).