SAVDOS as defualt. sigm_fbz. --pdos mode(on the way half). Cleanup

ctrlgen2.ctrl.* now have
SAVDOS as default.

lmf (when RDSIG=12), read sigm and write sigm_fbz.

lmf --pdos read sigm_fbz.

Clean up codes a little.
PDOS mode. Usage is followings.

How to get PDOS (lda or QSGW in the same manner but need sigm.* sigm_fbz.*).
-------
 1. lmf (when sigm is read, generate sigm_fbz.*).
 2. lmf --pdos si (it read sigm_fbz.si).
 3. lmdos --pdos si
 4. pldos dos.si -lst='1; 2:4; 5:9; 10:16; 17:25'
    (to know which is which channel, see "Channels in dos file generated by LMDOS:")
    In this exapmle, 5 panels which for s, p, d ,f,g channels are
    shown.(separator is ;.  I is possible to say -lst='5,6,8;
    7,9'. Then you have two panels with dos of 5+6+8 and dos of 7+9.

-----------------------------
Cr  For s,p,d,f the indices 1-16 correspond to:
Cr  index   l    m     polynomial
Cr     1    0    0        1
-----------------------------
Cr     2    1   -1        y
Cr     3    1    0        z
Cr     4    1    1        x
-----------------------------
Cr     5    2   -2        xy
Cr     6    2   -1        yz
Cr     7    2    0        3z^2-1
Cr     8    2    1        xz
Cr     9    2    2        x^2-y^2
-----------------------------
Cr     10   3   -3        y(3x^2-y^2)
Cr     11   3   -2        xyz
Cr     12   3   -1        y(5z^2-1)
Cr     13   3    0        z(5z^2-3)
Cr     14   3    1        x(5z^2-1)
Cr     15   3    2        z(x^2-y^2)
Cr     16   3    3        x(x^2-3y^2)

------------------
Development plan.

In future, we may need to make dos devide into atoms such as
dos.1.si
dos.2.si
dos.3.si
(site id is by lmchk).
Instead of pldos, gnuplot will be better. 
We like to unuse pldos.