!
!     Description of the calculation parameters:
!
! *=======+=========+=========+=========+=========+===========================*
! *                           GENERAL PARAMETERS :                            *
! *=======+=========+=========+=========+=========+===========================*
!
!       * Q_MIN    : minimum value of q (in units of k_F)
!       * Q_MAX    : maximum value of q (in units of k_F)
!       * N_Q      : number of q-points
!
!       * E_MIN    : minimum value of energy (in units of E_F)
!       * E_MAX    : maximum value of energy (in units of E_F)
!       * N_E      : number of E-points
!
!       * R_MIN    : minimum value of distance r (in units of 1/k_F)
!       * R_MAX    : maximum value of distance r (in units of 1/k_F))
!       * N_R      : number of r-points
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * RS       : average distance between 2 electrons (in units of a0)
!                                            ~
!       * MSOM     : m*/m (for semiconductors)
!                                            ~
!       * MAT_TYP  : type of material
!                       MAT_TYPE  = 'SCHRO' standard solid
!                       MAT_TYPE  = 'DIRAC' massless Fermions 
!                       MAT_TYPE  = 'WATER' liquid water
!                                            ~
!       * EPS_B    : material's dielectric constant
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * T        : temperature (in SI)
!                                            ~
!       * E        : external electric field (in SI)
!                                            ~
!       * H        : external magnetic field (in SI)
!                                            ~
!       * FLD      : strength of the magnetic field
!                       FLD = 'NO' no field        
!                       FLD = 'WF' weak field        
!                       FLD = 'IF' intermediate field
!                       FLD = 'LF' large field
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * DIM      : dimension of the system
!                       DIM  = '3D' 
!                       DIM  = '2D' 
!                       DIM  = '1D' 
!                       DIM  = 'Q0'  quasi-0D
!                       DIM  = 'Q1'  quasi-1D
!                       DIM  = 'Q2'  quasi-2D
!                       DIM  = 'BL'  bilayer
!                       DIM  = 'ML'  multilayer
!                       DIM  = 'ML'  multilayer
!
!                                            ~
!
!       * R0       : wire radius
!                                            ~
!       * L        : length of quantum well
!                                            ~
!       * OM0      : frequency of the confinement potential (SI)
!                                            ~
!       * CONFIN   : type of confinement
!                       CONFIN = 'NO-CONF' no confinement 
!                       CONFIN = 'DSEPLAY' layer within a stacking of layers
!                       CONFIN = 'CC-1111' cylindrical within subband 1
!                       CONFIN = 'CC-1122' cylindrical between subbands 1 and 2
!                       CONFIN = 'CC-1221' cylindrical between subbands 1 and 2
!                       CONFIN = 'CC-2222' cylindrical within subband 2
!                       CONFIN = 'HC-1111' harmonic within subband 1
!                       CONFIN = 'HC-1122' harmonic between subbands 1 and 2
!                       CONFIN = 'HC-1221' harmonic between subbands 1 and 2
!                       CONFIN = 'HC-2222' harmonic within subband 2
!                       CONFIN = 'INVLAYE' inversion layer in semiconductor
!                       CONFIN = 'IQWE_LB' square well with an infinite barrier
!                       CONFIN = 'PC1_QWI' parabolic 
!                       CONFIN = 'PC2_QWI' parabolic
!                       CONFIN = 'SOFTCOR' soft-core potential
!                       CONFIN = 'SWC_QWI' square well with an infinite barrier
! 
!                                            ~
!
!       * DL       : interlayer distance 
!                                            ~
!       * D1       : distance between the two layers in the unit cell
!                                            ~
!       * N_DEP    : electron concentration in depletion layer (SI)
!                                            ~
!       * N_INV    : electron concentration in inversion layer (SI)
! 
!                                            ~
!
!       * H_TYPE   : heterostructure type
!                       H_TYPE  =  'SSL1' semiconductor superlattice of type I
!                       H_TYPE  =  'SSL2' semiconductor superlattice of type II
!                       H_TYPE  =  'BILA' bilayer
!                       H_TYPE  =  'MLA1' multilayer with with one layer  / unit cell
!                       H_TYPE  =  'MLA2' multilayer with with two layers / unit cell
!                                            ~
!       * EPS_1    : background/layer dielectric constant
!                                            ~
!       * EPS_2    : interlayer dielectric constant
!
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * UNIT     : system unit  
!                       UNIT  = 'SIU' international system
!                       UNIT  = 'CGS' CGS system
!                       UNIT  = 'ATU' atomic units
! 
!       * UNIK     : K unit  
!                       UNIK  = 'SI' international system
!                       UNIK  = 'AU' atomic units
!
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * SC_TYPE  : type of screeening
!                       SC_TYPE  = 'NO' no screening
!                       SC_TYPE  = 'DH' Debye-Hückel
!                       SC_TYPE  = 'KL' Kleinman
!                       SC_TYPE  = 'ST' Streitenberger
!                       SC_TYPE  = 'TF' Thomas-Fermi
!                       SC_TYPE  = 'UI' Utsumi-Ichimaru
!                       SC_TYPE  = 'YT' Yasuhara-Takada
!
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * PL_TYPE  : type of plasma considered
!                      PL_TYPE  = 'OCP' --> one-component plasma (~ electron gas)
!                      PL_TYPE  = 'DCP' --> two-component plasma
! 
!       * ZION     : atomic number of (first) ion
! 
!       * ZION2    : atomic number of second ion
!
! *=======+=========+=========+=========+=========+===========================*
! *                        DIELECTRIC FUNCTION :                              *
! *=======+=========+=========+=========+=========+===========================*
!
!
!       * ESTDY    : static vs dynamic dielectric function
!                      STDY = ' STATIC'
!                      STDY = 'DYNAMIC'
!
!       * EPS_T    : type of dielcectric function
!                       EPS_T  = 'LONG' longitudinal --> D_FUNCL
!                       EPS_T  = 'TRAN' transverse   --> D_FUNCT
!
!       * D_FUNC   : model of dielectric function :
!
!                    1) Static:
!
!                       ---> longitudinal:
!
!                      D_FUNCL = 'LRPA' random phase approximation          (3D,2D,1D)
!                      D_FUNCL = 'THFE' Thomas-Fermi approximation          (3D,2D,1D)
!
!                    2) Dynamic:
!
!                       ---> transverse:
!
!                      D_FUNCT = 'RPA1' random phase approximation          (3D,2D)
!                      D_FUNCT = 'RPA2' random phase approximation          (3D)
!                      D_FUNCT = 'LVLA' linearized Vlasov                   (3D)
!                      D_FUNCT = 'MER1' Mermin                              (3D)
!                      D_FUNCT = 'BLTZ' Boltzmann                           (3D)
!
!                       ---> longitudinal:
!
!                      D_FUNCL = 'ATAS' Atwal-Ashcroft                      <-- T-dependent
!                      D_FUNCL = 'BLZ1' Boltzmann
!                      D_FUNCL = 'BLZ2' damped Boltzmann
!                      D_FUNCL = 'DACA' Arista-Brandt                       <-- T-dependent
!                      D_FUNCL = 'HEAP' Hertel-Appel
!                      D_FUNCL = 'HAFO' Hartree-Fock        
!                      D_FUNCL = 'HUCO' Hu-O'Connell                        <-- damping (3D,2D,Q1)
!                      D_FUNCL = 'HYDR' hydrodynamic                        <-- damping
!                      D_FUNCL = 'KLEI' Kleinman                            <-- T-dependent
!                      D_FUNCL = 'KLKD' Klimontovich-Kraeft                 <-- T-dependent
!                      D_FUNCL = 'KLKN' Klimontovich-Kraeft                 <-- T-dependent
!                      D_FUNCL = 'LAND' Landau parameters-based
!                      D_FUNCL = 'LVL1' linearized Vlasov (weak coupling)   <-- T-dependent
!                      D_FUNCL = 'LVL2' linearized Vlasov (strong coupling) <-- T-dependent
!                      D_FUNCL = 'MER1' Mermin 1                            <-- damping
!                      D_FUNCL = 'MER2' Mermin 2                            <-- T-dependent
!                      D_FUNCL = 'MSAP' mean spherical approximation
!                      D_FUNCL = 'NEVA' Nevanlinna                          <-- T-dependent
!                      D_FUNCL = 'PLPO' plasmon pole
!                      D_FUNCL = 'RDF1' Altshuler et al                     <-- damping
!                      D_FUNCL = 'RDF2' Altshuler et al                     <-- damping
!                      D_FUNCL = 'RPA1' RPA
!                      D_FUNCL = 'RPA2' RPA                                 <-- T-dependent
!                      D_FUNCL = 'RPA3' random phase approximation          <-- magnetic field (2D)
!                      D_FUNCL = 'UTIC' Utsumi-Ichimaru                     <-- T-dependent
!                      D_FUNCL = 'VLFP' Vlasov-Fokker-Planck                <-- damping
!
!       * NEV_TYPE : type of Nevalinna function used (only for D_FUNCL = 'NEVA')
!                      NEV_TYPE  = 'NONE' --> no function
!                      NEV_TYPE  = 'STA1' --> static value h(q)
!                      NEV_TYPE  = 'STA2' --> static value h(q)
!                      NEV_TYPE  = 'CLCO' --> Classical Coulomb OCP
!                      NEV_TYPE  = 'AMTA' --> Adamjan-Meyer-Tkachenko
!                      NEV_TYPE  = 'PEEL' --> Perel'-Eliashberg function
!                      NEV_TYPE  = 'PE76' --> Perel'-Eliashberg by Arkhipov et al
!
!                                            ~
!
!       * I_T      :  way to incorporate temperature dependence
!                      I_T  = 0  --> no temperature (T = 0)
!                      I_T  = 1  --> analytical models used
!                      I_T  = 2  --> temperature convolution
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * PL_DISP   : method used to compute the plasmon dispersion (3D real case)
!                       PL_DISP = 'HYDRODY' hydrodynamic model
!                       PL_DISP = 'RPA_MOD' RPA model
!                       PL_DISP = 'TWA_MOD' Toigo-Woodruff model
!                       PL_DISP = 'GOA_MOD' Gorobchenko model
!                       PL_DISP = 'HUBBARD' Hubbard model
!                       PL_DISP = 'ELASTIC' elastic model
!                       PL_DISP = 'SGBBN_M' SGBBN model            
!                       PL_DISP = 'UTI_MOD' Utsumi-Ichimaru model
!       * PL_DISP   : method used to compute the plasmon dispersion (2D real case)
!                       PL_DISP = 'HYDRODY' hydrodynamic model
!                       PL_DISP = 'RPA_MOD' RPA model
!       * PL_DISP   : method used to compute the plasmon dispersion (1D real case)
!                       PL_DISP = 'HYDRODY' hydrodynamic model
!                       PL_DISP = 'RPA_MOD' RPA model
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * GSTDY    : static vs dynamic local-field corrections
!                      GSTDY = ' STATIC'  --> GQ_TYPE
!                      GSTDY = 'DYNAMIC'  --> GQO_TYPE 
!
!       * GQ_TYPE  : local-field correction type (3D) static
!                       GQ_TYPE = 'NONE' no local field correction
!                       GQ_TYPE = 'ALDA' adiabatic local density
!                       GQ_TYPE = 'ALFL' Alvarellos-Flores
!                       GQ_TYPE = 'BEBR' Bedell-Brown
!                       GQ_TYPE = 'CDOP' TDDFT Corradini et al correction
!                       GQ_TYPE = 'GEV2' Geldart-Vosko 2
!                       GQ_TYPE = 'GEVO' Geldart-Vosko correction
!                       GQ_TYPE = 'GOCA' Gold-Calmels
!  temperature-dep. --> GQ_TYPE = 'HNCA' hypernetted chain
!                       GQ_TYPE = 'HORA' Holas-Rahman
!                       GQ_TYPE = 'HUBB' Hubbard correction (only exchange)
!                       GQ_TYPE = 'ICUT' Ichimaru-Utsumi correction
!                       GQ_TYPE = 'IKPA' Iwamoto-Krotscheck-Pines
!                       GQ_TYPE = 'IWA1' Iwamoto G_{-1}
!                       GQ_TYPE = 'IWA2' Iwamoto G_{3} approx.
!  temperature-dep. --> GQ_TYPE = 'IWA3' Iwamoto G_{-1}  
!                       GQ_TYPE = 'IWA4' Iwamoto G_{3} exact
!                       GQ_TYPE = 'JGDG' Jung-Garcia-Gonzalez-Dobson-Godby
!                       GQ_TYPE = 'KLLA' Kleinman-Langreth correction
!                       GQ_TYPE = 'LDAC' LDA correction
!                       GQ_TYPE = 'MCSC' Moroni-Ceperley-Senatore correction
!                       GQ_TYPE = 'NAGY' Nagy correction
!                       GQ_TYPE = 'PAVA' Pavas-Vashishta correction
!                       GQ_TYPE = 'PGGA' Petersilka-Gossmann-Gross
!                       GQ_TYPE = 'RICE' Rice correction
!                       GQ_TYPE = 'SHAW' Shaw correction
!                       GQ_TYPE = 'SLAT' Slater correction
!                       GQ_TYPE = 'STLS' Singwi et al correction
!  temperature-dep. --> GQ_TYPE = 'TKAC' Tkachenko correction
!                       GQ_TYPE = 'TRMA' Tripathy-Mandal
!                       GQ_TYPE = 'VASI' Vashishta-Singwi correction
!                       GQ_TYPE = 'UTI1' Utsumi-Ichimaru correction (only exchange)
!
!       * GQ_TYPE  : local-field correction type (2D) static
!                       GQ_TYPE = 'NONE' no local field correction
!                       GQ_TYPE = 'BUTO' Bulutay-Tomak
!                       GQ_TYPE = 'DPGT' Davoudi-Giuliani-Giuliani-Tosi
!                       GQ_TYPE = 'GOCA' Gold-Calmels
!                       GQ_TYPE = 'HUBB' Hubbard correction (only exchange)
!                       GQ_TYPE = 'IWA1' Iwamoto G_{-1}
!                       GQ_TYPE = 'IWA2' Iwamoto G_{3}
!                       GQ_TYPE = 'SAIC' Sato-Ichimaru correction
!
!       * GQ_TYPE  : local-field correction type (1D) static
!                       GQ_TYPE = 'NONE' no local field correction
!                       GQ_TYPE = 'GOCA' Gold-Calmels
!                       GQ_TYPE = 'HUBB' Hubbard correction (only exchange)
! 
!       * IQ_TYPE  : type of approximation for I(q)
!                       IQ_TYPE  = 'GKM' Gorobchenko-Kohn-Maksimov
!                       IQ_TYPE  = 'HKA' Hong-Kim
!                       IQ_TYPE  = 'IKP' Iwamoto-Krotscheck-Pines parametrization
!                       IQ_TYPE  = 'KU1' Kugler 1
!                       IQ_TYPE  = 'KU2' Kugler 1
!
!                                            ~
!
!       * LANDAU  : model chosen for the calculation of the Landau parameters (3D)
!                       LANDAU = 'NONE' Landau's theory not used
!                       LANDAU = 'CHEN' Chen's approach
!                       LANDAU = 'RASC' Rayleigh-Schrödinger expansion
!                       LANDAU = 'ANBR' Anderson-Brinkman model
!                       LANDAU = 'GUTZ' Gutzwiller model
!                       LANDAU = 'IWPI' Iwamoto-Pines model (hard-sphere)
!                       LANDAU = 'GCYO' Giuliani-Vignale parametrization of 
!                                        Yasuhara-Ousaka approach
!                       LANDAU = 'SBOH' slave-boson one-band Hubbard model
!
!       * LANDAU  : model chosen for the calculation of the Landau parameters (2D)
!                       LANDAU = 'NONE' Landau's theory not used
!                       LANDAU = 'ERZA' Engelbrecht-Randeria-Zhang approach
!                       LANDAU = 'GVYO' Giuliani-Vignale parametrization of 
!                                        Yasuhara-Ousaka approach
!                       LANDAU = 'KCMP' Kwoon-Ceperley-Martin parametrization
!       * GQO_TYPE : local-field correction type (3D)
!                       GQO_TYPE = 'NONE' no local field correction
!                       GQO_TYPE = 'ALFL' Alvarellos-Flores correction
!                       GQO_TYPE = 'BACA' Barriga-Carrasco correction
!                       GQO_TYPE = 'BBSA' Bachlechner-Böhm-Schinner
!                       GQO_TYPE = 'COPI' Constantin-Pitarke
!                       GQO_TYPE = 'DABR' Dabrowski
!                       GQO_TYPE = 'FWRA' Forstmann-Wierling-Röpke
!                       GQO_TYPE = 'HOK1' Hong-Kim correction
!                       GQO_TYPE = 'HOK2' Hong-Kim correction
!                       GQO_TYPE = 'JEWS' Jewsbury approximation
!                       GQO_TYPE = 'KUG1' Kugler q --> 0 approximation
!                       GQO_TYPE = 'KUG2' Kugler approximation
!                       GQO_TYPE = 'MDGA' Mithen-Daligault-Gregori
!                       GQO_TYPE = 'NLGA' Nagy-Laszlo-Giber approximation
!                       GQO_TYPE = 'RIA1' Richardson-Ashcroft G_s
!                       GQO_TYPE = 'RIA2' Richardson-Ashcroft G_n
!                       GQO_TYPE = 'RIA3' Richardson-Ashcroft G_a
!                       GQO_TYPE = 'SHMU' Shah-Mukhopadhyay 
!                       GQO_TYPE = 'TOWO' Toigo-Woodruff
!                       GQO_TYPE = 'UTI2' Utsumi-Ichimaru approximation
!                       GQO_TYPE = 'VISC' viscosity approximation
! 
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * DAMPING : type of damping used
!                       DAMPING = 'NONE' no damping
!                       DAMPING = 'LFTM' lifetime
!                       DAMPING = 'RELA' relaxation time
!                       DAMPING = 'DECA' decay rate
!                       DAMPING = 'DIFF' diffusion coefficient
!                       DAMPING = 'VISC' viscosity
!
!       * LT_TYPE    : approximation used for lifetime (3D)
!                       LT_TYPE = 'DAVI' Davies formula
!                       LT_TYPE = 'GALI' Galitskii formula
!                       LT_TYPE = 'GIQU' Giuliani-Quinn formula
!                       LT_TYPE = 'GIVI' Giuliani-Vignale formula 
!                       LT_TYPE = 'INPE' Inogamov-Petrov formula
!                       LT_TYPE = 'LUBR' Lugovskoy-Bray formula
!                       LT_TYPE = 'NAEC' Nagy-Echenique formula
!                       LT_TYPE = 'QIVI' Qian-Vignale formula
!                       LT_TYPE = 'QUFE' Quinn-Ferrell formula
!
!       * LT_TYPE    : approximation used for lifetime (2D)
!                       LT_TYPE = 'GIQ1' Giuliani-Quinn formula for e-h loss
!                       LT_TYPE = 'GIQ2' Giuliani-Quinn formula for plasmon loss
!                       LT_TYPE = 'GIVI' Giuliani-Vignale formula
!                       LT_TYPE = 'HAWR' Hawrylak formula
!                       LT_TYPE = 'MELA' Menashe-Laikhtman formula
!                       LT_TYPE = 'QIVI' Qian-Vignale formula
!
!       * RT_TYPE  : relaxation time
!                       RT_TYPE = ' NO'   --> no relaxation time
!                       RT_TYPE = 'E-E'   --> electron-electron interaction
!                       RT_TYPE = 'E-P'   --> electron-phonon interaction
!                       RT_TYPE = 'E-I'   --> electron-phonon impurity
!                       RT_TYPE = 'ALL'   --> all three taken into account
!
!                                            ~
! 
!       * DR_TYPE  : decay rate in 3D
!                       DR_TYPE = 'UTIC'   --> Utsumi-Ichimaru approximation
!                       DR_TYPE = 'VLAS'   --> Vlasov approximation
!
!       * DC_TYPE  : diffusion coefficient in 3D
!                       DC_TYPE = 'ASHO'   --> Ashurst-Hoover
!
!       * VI_TYPE  : viscosity in 3D
!                       VI_TYPE = 'AMPP' Angilella et al hard-sphere fluid     --> T-dependent
!                       VI_TYPE = 'DRBA' Daligault-Rasmussen-Baalrud (plasmas) --> T-dependent
!                       VI_TYPE = 'KHRA' Khrapak for Yukawa fluid  --> T-dependent
!                       VI_TYPE = 'LHPO' Longuet-Higgins-Pope      --> T-dependent
!                       VI_TYPE = 'SCHA' Schäfer                   --> T-dependent
!                       VI_TYPE = 'SCHD' Schäfer (dynamic)         --> T-dependent
!                       VI_TYPE = 'SHTE' Shternin                  --> T-dependent
!                       VI_TYPE = 'STEI' Steinberg low-temperature --> T-dependent
!
!       * VI_TYPE  : viscosity in 2D
!                       VI_TYPE = 'SCHA' Schäfer                     --> T-dependent
!         graphene <--  VI_TYPE = 'KISC' Kiselev-Schmalian (dynamic) --> T-dependent
!         graphene <--  VI_TYPE = 'MSFA' Müller-Schmalian-Fritz      --> T-dependent  
!
!                                            ~
!
!       * EE_TYPE  : e-e relaxation time in 3D
!                       RT_TYPE = 'ALAR'   --> Al'tshuler-Aronov (e-e + impurities)
!                       RT_TYPE = 'BACA'   --> Barriga-Carrasco approximation (e-e)
!                       RT_TYPE = 'FSTB'   --> Fann et al approximation       (e-e)
!                       RT_TYPE = 'QIVI'   --> Qian-Vignale                   (e-e)
!                       RT_TYPE = 'RASM'   --> Rammer-Smith                   (e-e)
!                       RT_TYPE = 'UTIC'   --> Utsumi-Ichimaru approximation  (e-e)
!                       RT_TYPE = 'TAIC'   --> Tanaka-Ichimaru approximation  (e-e)
!
!
!       * EE_TYPE  : relaxation time in 2D
!                       EE_TYPE = 'FUAB'   --> Fukuyama-Abrahams (disordered metals)
!          graphene <-- EE_TYPE = 'LUFO'   --> Lucas-Fong      (e-e)
!                       EE_TYPE = 'QIVI'   --> Qian-Vignale    (e-e)
!                       EE_TYPE = 'RASM'   --> Rammer-Smith    (e-e)
!  heterostructures <-- EE_TYPE = 'REWI'   --> Reizer-Wilkins  (e-e)
!                       EI_TYPE = 'SHAS'   --> Sharma-Ashraf   (e-e + impurities)
!                       EE_TYPE = 'ZHDA'   --> Zhang-Das Sarma (e-e)
!
!       * EE_TYPE  : relaxation time in 1D
!                       EI_TYPE = 'SHAS'   --> Sharma-Ashraf   (e-e + impurities)
!
!       * EP_TYPE  : e-phonon relaxation time in 3D
!                       EP_TYPE = 'STEL'   --> Steinberg low-temperature  
!                       EP_TYPE = 'STEH'   --> Steinberg High-temperature
!
!       * EI_TYPE  : e-impurit relaxation time in 3D
!                       EI_TYPE = 'HEAP'   --> Hertel-Appel approximation
!
!                                            ~
!
!       * IP_TYPE  : ion plasma relaxation time in 3D
!                       IP_TYPE = 'SEMO'   --> Selchow-Morawetz approximation
!                       IP_TYPE = 'SPIT'   --> Spitzer approximation 
!
!       * PD_TYPE  : method used to compute the plasmon damping (3D)
!                       PD_TYPE = 'NONE'   --> no plasmon damping
!                       PD_TYPE = 'CALL'   --> Callen approximation
!                       PD_TYPE = 'DGKA'   --> DuBois-Gilinsky-Kivelson approximation
!                       PD_TYPE = 'FEWA'   --> Fetter and Walecka approximation
!                       PD_TYPE = 'JEWS'   --> Jewsbury approximation
!                       PD_TYPE = 'LITI'   --> Giuliani-Quinn lifetime approximation
!                       PD_TYPE = 'MOPE'   --> Molinari-Peerani approximation
!                       PD_TYPE = 'NPSA'   --> Ninham-Powel-Swanson approximation
!                       PD_TYPE = 'SGAA'   --> Segui-Gervasoni-Arista approximation
!
!                                            ~
!
!       * QD_TYPE  : method used to compute q-dependent relaxation time
!                       QD_TYPE = 'NONE'   --> no q-dependence
!                       QD_TYPE = 'GAUS'   --> Gaussian
!                       QD_TYPE = 'LORE'   --> Lorentzian
!
!       * ZETA     : Value of Tanaka-Ichimaru parameter
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * INT_POT  : type of interaction potential (3D)
!                       INT_POT = 'COULO' Coulomb interaction
!                       INT_POT=  'YUKAW' Yukawa interaction
!                       INT_POT=  'RPAPO' RPA interaction
!                       INT_POT = 'OVER1' Overhauser interaction
!                       INT_POT = 'OVER2' modified Overhauser interaction
!                       INT_POT = 'DEUTS' Deutsch interaction
!                       INT_POT = 'PHOLE' particle-hole interaction
!                       INT_POT = 'KELBG' Kelbg interaction
!
!        * S       :  
!
!        * EPS     :  
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!       * EK       :  electron kinetic energy (eV)
!       * EP_C     :  electron-phonon coupling
!       * DEBYE_T  :  material Debye temperature
!
!                                            ~
!
!       * NA       :  number of atoms per unit volume
!       * MA       :  mass of the atoms
!       * RA       :  radius of the atoms
!
! *-------+---------+---------+---------+---------+---------------------------*
!
!        * NI      :  impurity concentration
!        * EI_C    :  strength of impurity scattering
!
! *-------+---------+---------+---------+---------+---------------------------*
!

!        * CF_TYPE  : type of classical fluid calculation         --> 3D
!                       CF_TYPE = 'SHS' smooth hard spheres
!                       CF_TYPE = 'RH1' rough hard spheres (Pidduck)
!                       CF_TYPE = 'RH2' rough hard spheres (Condiff-Lu-Dahler)
!                       CF_TYPE = 'RH3' rough hard spheres (McCoy-Sandler-Dahler)
!                       CF_TYPE = 'DCE' dilute Chapman-Enskog
!                       CF_TYPE = 'HCE' heavy (i.e. dense) Chapman-Enskog
!                       CF_TYPE = 'LJF' Lennard-Jones fluid

!                       CF_TYPE = 'DHD' dense hard disks          --> 2D 

!        * PF_TYPE  : type of packing fraction                    --> 2D    
!                       PF_TYPE = 'HDM'   --> hard disk model
!
!                                                                 --> 3D
!                       PF_TYPE = 'HSM'   --> hard sphere model
!                       PF_TYPE = 'RCP'   --> random closed-packed
!                       PF_TYPE = 'FCC'   --> FCC closed-packed
!                       PF_TYPE = 'FRE'   --> freezing
!                       PF_TYPE = 'MEL'   --> melting
!
!        * SL_TYPE  : type of scattering length calculation
!                       SL_TYPE = 'HSP'   --> hard sphere potential
!                       SL_TYPE = 'ASW'   --> attractive square well (without bound state)
!                       SL_TYPE = 'RSW'   --> repulsive square well
!                       SL_TYPE = 'DSP'   --> delta-shell potential
!                       SL_TYPE = 'AYP'   --> attractive Yukawa potential
!                       SL_TYPE = 'CCO'   --> Coulomb cut-off potential
!                       SL_TYPE = 'HUL'   --> Hulthén potential
!
! *=======+=========+=========+=========+=========+===========================*
! *                        STRUCTURE FACTOR :                                 *
! *=======+=========+=========+=========+=========+===========================*
!
!
!       * SSTDY    : static vs dynamic local-field corrections
!                      SSTDY = ' STATIC'  --> SQ_TYPE
!                      SSTDY = 'DYNAMIC'  --> SQO_TYPE 
!
!       * SQ_TYPE   : structure factor approximation (3D) --> static
!                       SQ_TYPE  = 'DEH' Debye-Hückel approximation
!                       SQ_TYPE  = 'GEA' generalized approximation
!                       SQ_TYPE  = 'GOR' Gorobchenko approximation
!                       SQ_TYPE  = 'GR2' computed from g(r) (GR_TO_SQ.f code)
!                       SQ_TYPE  = 'GSB' Gori-Giorgi-Sacchetti-Bachelet approximation
!                       SQ_TYPE  = 'HFA' Hartree-Fock approximation (only exchange)
!                       SQ_TYPE  = 'HUB' Hubbard approximation
!                       SQ_TYPE  = 'ICH' Ichimaru approximation
!                       SQ_TYPE  = 'MSA' mean spherical approximation
!                       SQ_TYPE  = 'PKA' Pietiläinen-Kallio
!                       SQ_TYPE  = 'RPA' RPA approximation
!                       SQ_TYPE  = 'SHA' Shaw approximation
!                       SQ_TYPE  = 'SIN' Singh
!                       SQ_TYPE  = 'SPA' Singh-Pathak
!                       SQ_TYPE  = 'TWA' Toigo-Woodruff approximation
!! 
!       * SQO_TYPE   : structure factor approximation (3D) --> dynamic
!                       SQO_TYPE  = 'ABA' Arista-Brandt approximation
!                       SQO_TYPE  = 'HFA' Hartree-Fock approximation
!                       SQO_TYPE  = 'HYD' hyrodynamic approximation
!                       SQO_TYPE  = 'IGA' ideal gas approximation
!                       SQO_TYPE  = 'ITA' Ichimaru-Tanaka approximation
!                       SQO_TYPE  = 'MFA' Hansen-McDonald-Pollock approximation
!                       SQO_TYPE  = 'NIC' Nakano-Ichimaru approximation
!                       SQO_TYPE  = 'UTI' Utsumi-Ichimaru approximation (3D)
!                       SQO_TYPE  = 'VLA' linearized Vlasov approximation
! 
! *=======+=========+=========+=========+=========+===========================*
! *                        PAIR CORRELATION FUNCTION :                        *
! *=======+=========+=========+=========+=========+===========================*
!
!       * GR_TYPE     : pair correlation function approximation g(r) (3D)
!                       GR_TYPE  = 'CDF' from chain diagram formula of PDF (long distance)
!                       GR_TYPE  = 'DHA' Debye-Hückel approximation 
!                       GR_TYPE  = 'DWA' DeWitt approximation 
!                       GR_TYPE  = 'FBA' Frieman-Book approximation 
!                       GR_TYPE  = 'HFA' Hartree-Fock approximation (only exchange)
!                       GR_TYPE  = 'HUB' Hubbard approximation
!                       GR_TYPE  = 'LLA' Lee-Long approximation
!                       GR_TYPE  = 'ORB' Ortiz-Ballone approximation
!                       GR_TYPE  = 'PDF' from pair distribution function 
!                       GR_TYPE  = 'SHA' Shaw approximation
!                       GR_TYPE  = 'SQ2' computed from S(q) (SQ_TO_GR.f code)
!                       GR_TYPE  = 'WIG' Wigner approximation
! 
!       * GR0_MODE : g(0) (3D)
!                       GR0_MODE  = 'CAGO' --> Calmels-Gold
!                       GR0_MODE  = 'DPGT' --> Davoudi-Polini-Giuliani-Tosi
!                       GR0_MODE  = 'HASA' --> Holas-Aravind-Singwi (small r_s)
!                       GR0_MODE  = 'ICHI' --> Ichimaru
!                       GR0_MODE  = 'KIMB' --> Kimball
!                       GR0_MODE  = 'OVE1' --> Overhauser 1
!                       GR0_MODE  = 'OVE2' --> Overhauser 2
!                       GR0_MODE  = 'QIAN' --> Qian
!       * GR0_MODE : g(0) (2D)
!                       GR0_MODE  = 'CAGO' --> Calmels-Gold
!                       GR0_MODE  = 'HAFO' --> Hartree-Fock 
!                       GR0_MODE  = 'MOMA' --> Moreno-Marinescu 
!                       GR0_MODE  = 'NSOA' --> Nagano-Singwi-Ohnishi
!                       GR0_MODE  = 'QIAN' --> Qian
! 
! *=======+=========+=========+=========+=========+===========================*
! *                        PAIR DISTRIBUTION FUNCTION :                       *
! *=======+=========+=========+=========+=========+===========================*
! 
!       * RH_TYPE     : pair distribution function approximation (3D)
!                       RH_TYPE  = 'CDI' chain diagram improved
!                       RH_TYPE  = 'CEG' classical electron gas
!                       RH_TYPE  = 'DEB' Debye electron gas 
!                       RH_TYPE  = 'FUA' correct to order 2 in epsilon
!                       RH_TYPE  = 'SDC' short-distance correlations
!                       RH_TYPE  = 'WDA' watermelon diagrams summed
! 
 *=======+=========+=========+=========+=========+============================*
 *                         ENERGY CALCULATIONS :                              *
 *=======+=========+=========+=========+=========+============================*
!
!       * EC_TYPE  : type of correlation energy functional (3D)
!                       EC_TYPE = 'GEBR_W'   --> Gell-Mann and Brueckner
!                       EC_TYPE = 'CAMA_W'   --> Carr and Maradudin
!                       EC_TYPE = 'EHTY_S'   --> Endo-Horiuchi-Takada-Yasuhara
!                       EC_TYPE = 'HELU_W'   --> Hedin and Lundqvist
!                       EC_TYPE = 'VBHE_W'   --> von Barth and Hedin
!                       EC_TYPE = 'PEZU_W'   --> Perdew and Zunger
!                       EC_TYPE = 'WIGN_S'   --> Wigner
!                       EC_TYPE = 'NOPI_S'   --> Nozières and Pines
!                       EC_TYPE = 'LIRO_S'   --> Lindgren and Rosen
!                       EC_TYPE = 'PEZU_S'   --> Perdew and Zunger
!                       EC_TYPE = 'REHI_S'   --> Rebei and Hitchon
!                       EC_TYPE = 'GGSB_G'   --> Gori-Giorgi-Sacchetti-Bachelet
!                       EC_TYPE = 'PRKO_G'   --> Proynov and Kong
!                       EC_TYPE = 'VWNU_G'   --> Vosko, Wilk and Nusair
!                       EC_TYPE = 'PEWA_G'   --> Perdew and Wang
!                       EC_TYPE = 'HUBB_G'   --> Hubbard
!                       EC_TYPE = 'CHAC_G'   --> Chachiyo
!                       EC_TYPE = 'ISKO_T'   --> Isihara and Kojima
!       * EC_TYPE  : type of correlation energy functional (2D)
!                       EC_TYPE = 'TACE_G'   --> Tanatar-Ceperley
!                       EC_TYPE = 'CPPA_G'   --> Seidl-Perdew_Levy
!                       EC_TYPE = 'AMGB_G'   --> Attaccalite-Moroni-Gori-Giorgi-Bachelet
!                       EC_TYPE = 'SEID_G'   --> Seidl
!                       EC_TYPE = 'LOOS_W'   --> Loos
!                       EC_TYPE = 'WIGN_S'   --> Wigner
!                       EC_TYPE = 'ISTO_T'   --> Isihara-Toyoda
!       * EC_TYPE  : type of correlation energy functional (1D)
!                       EC_TYPE = 'LOOS_W'   --> Loos
!                       EC_TYPE = 'WIGN_S'   --> Wigner
! 
! 
!       * FXC_TYPE : type of XC free energy functional           --> 3D
!                       FXC_TYPE = 'NO'   --> None
!                       FXC_TYPE = 'EB'   --> Ebeling et al
!                       FXC_TYPE = 'IC'   --> Ichimaru et al
!                       FXC_TYPE = 'KS'   --> Karasiev et al
!                       FXC_TYPE = 'VS'   --> Vashishta and Singwi
!                       FXC_TYPE = 'PD'   --> Perrot and Dharma-Wardana
!                       FXC_TYPE = 'EK'   --> Ebeling-Kraeft-Kremp-Röpke
 
!       * EXC_TYPE : type of exchange-correlation energy functional   --> 3D
!                       EXC_TYPE = 'NO'   --> None
!                       EXC_TYPE = 'GT'   --> Goedeker-Tetter-Hutter
!                       EXC_TYPE = 'ST'   --> 
!                       EXC_TYPE = 'BD'   --> Brown-DuBois-Holzmann-Ceperley
!
!                                            ~
!
! 
 *=======+=========+=========+=========+=========+============================*
 *                         SPIN POLARIZATION :                                *
 *=======+=========+=========+=========+=========+============================*
! 
!       * IMODE    : choice of spin parameters
!                       IMODE = 1 : no spin polarization 
!                       IMODE = 2 : fully spin-polarized
!             
!       * XI       : spin polarization : (n+ - n-) / n
!
! *=======+=========+=========+=========+=========+===========================*
! *                        THERMODYNAMIC PROPERTIES :                         *
! *=======+=========+=========+=========+=========+===========================*
! 
!       * TH_PROP   : type of calculation  --> thermodynamic properties
!                      TH_PROP  = 'CLAS' : classical approximation 
!                      TH_PROP  = 'QUAN' : quantum approximation 
!
!       * GP_TYPE  : grand partition function type (3D)
!                       GP_TYPE = 'IK0' Isihara-Kojima formulation
!                       GP_TYPE = 'RH0' Rebei-Hitchon formulation
!                       GP_TYPE = 'IKM' Isihara-Kojima with magnetic field
!
!       * GP_TYPE  : grand partition function type (2D)
!                       GP_TYPE = 'I20' Isihara-Kojima formulation
!                       GP_TYPE = 'I2M' Isihara-Kojima with magnetic field
! 
! *=======+=========+=========+=========+=========+===========================*
! *                        INCOMING ION BEAM :                                *
! *=======+=========+=========+=========+=========+===========================*
!
!       * Z_BEAM    :  charge of ions in incoming beam
!       * EK_BEAM   :  kinetic energy of incoming beam ions (eV)
!
! *=======+=========+=========+=========+=========+===========================*
! *                        OUTPUT CALCULATIONS :                              *
! *=======+=========+=========+=========+=========+===========================*
! 
!       * I_DF      : switch for dielectric function printing
!                      I_DF  = 0 : dielectric function not printed
!                      I_DF  = 1 : dielectric function printed in file 'diel_func.dat'
! 
!       * I_PZ      : switch for polarization function printing
!                      I_PZ  = 0 : polarization function not printed
!                      I_PZ  = 1 : polarization function printed in file 'pola_func.dat'
! 
!       * I_SU      : switch for susceptibility function printing
!                      I_SU  = 0 : susceptibility function not printed
!                      I_SU  = 1 : susceptibility function printed in file 'susc_func.dat'
! 
!       * I_CD      : switch for electrical conductivity printing
!                      I_CD  = 0 : electrical conductivity not printed
!                      I_CD  = 1 : electrical conductivity printed in file 'cond_func.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_PD      : switch for plasmon dispersion printing
!                      I_PD  = 0 :  plasmon dispersion not calculated
!                      I_PD  = 1 :  plasmon dispersion printed in file 'plas_disp.dat'
! 
!       * I_EH      : switch for electron-hole dispersion printing
!                      I_EH  = 0 :  electron-hole dispersion not calculated
!                      I_EH  = 1 :  electron-hole dispersion printed in file 'elec_hole.dat'
! 
!       * I_E2      : switch for two-electron-hole dispersion printing
!                      I_E2  = 0 :  two-electron-hole dispersion not calculated
!                      I_E2  = 1 :  two-electron-hole dispersion printed in file 'elec_hol2.dat'
! 
!       * I_CK      : switch for k-space e-e interaction potential printing
!                      I_CK  = 0 :  potential not calculated
!                      I_CK  = 1 :  potential printed in file 'int_pot_k.dat'
! 
!       * I_CR      : switch for real-space e-e interaction potential printing
!                      I_CR  = 0 :  potential not calculated
!                      I_CR  = 1 :  potential printed in file 'int_pot_r.dat'
! 
!       * I_SA      : switch for scattering amplitude printing
!                      I_SA  = 0 :  scattering amplitude not calculated
!                      I_SA  = 1 :  scattering amplitude printed in file 'scat_ampl.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_LF      : switch for local-field corrections G(q,omega) printing
!                      I_LF  = 0 : local-field corrections not calculated
!                      I_LF  = 1 : local-field corrections printed in file 'loca_fiel.dat'
! 
!       * I_IQ      : switch for G(q,inf) printing
!                      I_IQ  = 0 : G(q,inf) not calculated
!                      I_IQ  = 1 : G(q,inf) printed in file 'ginf_fiel.dat'
! 
!       * I_SF      : switch for structure factor S(q,omega) printing
!                      I_SF  = 0 : structure factor not calculated
!                      I_SF  = 1 : structure factor printed in file 'stru_fact.dat'
! 
!       * I_PC      : switch for pair correlation function g(r) printing
!                      I_PC  = 0 :  pair correlation function not calculated
!                      I_PC  = 1 :  pair correlation function printed in file 'pair_corr.dat'
! 
!                                            ~
! 
!       * I_P2      : switch for pair distribution rho2(r) printing
!                      I_P2  = 0 :  pair distribution function not calculated
!                      I_P2  = 1 :  pair distribution function printed in file 'pair_dist.dat'
! 
!       * I_VX      : switch for  vertex function Gamma(q,omega) printing
!                      I_VX  = 0 :  vertex function not calculated
!                      I_VX  = 1 :  vertex function printed in file 'vertex_fu.dat'
! 
!       * I_DC      : switch for plasmon damping coefficient Im[eps]/q^2 printing
!                      I_DC  = 0 :  plasmon damping not calculated
!                      I_DC  = 1 :  plasmon damping printed in file 'plas_damp.dat'
! 
!       * I_MD      : switch for momentum distribution printing
!                      I_MD  = 0 :  momentum distribution not calculated
!                      I_MD  = 1 :  momentum distribution printed in file 'mome_dist.dat'
! 
!                                            ~
! 
!       * I_LD      : switch for Landau parameters printing
!                      I_LD  = 0 : Landau parameters not calculated
!                      I_LD  = 1 : Landau parameters printed in file 'landau_pa.dat'
! 
!       * I_DP      : switch for damping printing
!                      I_DP  = 0 : damping not calculated
!                      I_DP  = 1 : damping printed in file 'damp_file.dat'
! 
!       * I_LT      : switch for plasmon lifetime printing
!                      I_LT  = 0 : plasmon lifetime not calculated
!                      I_LT  = 1 : plasmon lifetime printed in file 'life_time.dat'
! 
!       * I_BR      : switch for plasmon broadening printing
!                      I_BR  = 0 : plasmon broadening not calculated
!                      I_BR  = 1 : plasmon broadening printed in file 'broadenin.dat'
! 
!                                            ~
! 
!       * I_PE      : switch for plasmon energy printing
!                      I_PE  = 0 : plasmon energy not calculated
!                      I_PE  = 1 : plasmon energy printed in file 'plas_ener.dat'
! 
!       * I_QC      : switch for plasmon q-bounds printing
!                      I_QC  = 0 : plasmon q-bounds not calculated
!                      I_QC  = 1 : plasmon q-bounds printed in file 'qc_bounds.dat'
! 
!       * I_RL      : switch for relaxation time printing
!                      I_RL  = 0 : relaxation time not calculated
!                      I_RL  = 1 : relaxation time printed in file 'rela_time.dat'
! 
!       * I_TF      : switch for Thomas-Fermi wave vector printing
!                      I_TF  = 0 : Thomas-Fermi wave vector not calculated
!                      I_TF  = 1 : Thomas-Fermi wave vector printed in file 'thomas-fe.dat'
! 
!                                            ~
! 
!       * I_DY      : switch for Debye wave vector printing
!                      I_DY  = 0 : Debye wave vector not calculated
!                      I_DY  = 1 : Debye wave vector printed in file 'debye_wav.dat'
! 
!       * I_ME      : switch for moments of epsilon(q,omega) printing
!                      I_ME  = 0 : moments of epsilon not calculated
!                      I_ME  = 1 : moments of epsilon printed in file 'moments_e.dat'
! 
!       * I_MS      : switch for moments of S(q,omega) printing
!                      I_MS  = 0 : moments of structure factor not calculated
!                      I_MS  = 1 : moments of structure factor printed in file 'moments_s.dat'
! 
!       * I_ML      : switch for moments of loss function printing
!                      I_ML  = 0 : moments of loss function not calculated
!                      I_ML  = 1 : moments of loss function printed in file 'moments_l.dat'
! 
!                                            ~
! 
!       * I_MC      : switch for moments of conductivity printing
!                      I_MC  = 0 : moments of conductivity not calculated
!                      I_MC  = 1 : moments of conductivity printed in file 'moments_c.dat'
! 
!       * I_DE      : switch for derivative of Re[ dielectric function ] printing
!                      I_DE  = 0 : derivative not calculated
!                      I_DE  = 1 : derivative printed in file 'deri_epsi.dat'
! 
!       * I_ZE      : switch for Re[ dielectric function ] = 0 printing
!                      I_ZE  = 0 : function not calculated
!                      I_ZE  = 1 : function printed in file 'ree0_file.dat'
! 
!       * I_SR      : switch for sum rules for epsilon printing
!                      I_SR  = 0 : sum rules not calculated
!                      I_ST  = 1 : sum rules printed in file 'sum_rules.dat'
! 
!                                            ~
! 
!       * I_CW      : switch for confinement wave function printing
!                      I_CW  = 0 : confinement wave function not calculated
!                      I_CW  = 1 : confinement wave function printed in file 'confin_wf.dat'
! 
!       * I_CF      : switch for confinement potential printing
!                      I_CF  = 0 : confinement potential not calculated
!                      I_CF  = 1 : confinement potential printed in file 'confin_pt.dat'
! 
!       * I_EM      : switch for effective mass printing
!                      I_EM  = 0 : effective mass not calculated
!                      I_EM  = 1 : effective mass printed in file 'effe_mass.dat'
! 
!       * I_MF      : switch for mean free path printing
!                      I_MF  = 0 : mean free path not calculated
!                      I_MF  = 1 : mean free path printed in file 'mean_path.dat'
! 
!                                            ~
! 
!       * I_SP      : switch for spectral function printing
!                      I_SP  = 0 : spectral function not calculated
!                      I_SP  = 1 : spectral function printed in file 'spec_func.dat'
! 
!       * I_SE      : switch for self-energy printing
!                      I_SE  = 0 : self-energy not calculated
!                      I_SE  = 1 : self-energy printed in file 'self_ener.dat'
! 
!       * I_NV      : switch for Nevanlinaa function printing
!                      I_NV  = 0 : Nevanlinaa function not calculated
!                      I_NV  = 1 : Nevanlinaa function printed in file 'nevanlina.dat'
! 
!       * I_ES      : switch for Eliashberg function printing
!                      I_ES  = 0 : Eliashberg function not calculated
!                      I_ES  = 1 : Eliashberg function printed in file 'elia_func.dat'
! 
!                                            ~
! 
!       * I_GR      : switch for Grüneisen parameter printing
!                      I_GR  = 0 : Grüneisen parameter not calculated
!                      I_GR  = 1 : Grüneisen parameter printed in file 'grune_par.dat'
! 
!       * I_FD      : switch for Fermi-Dirac distribution printing
!                      I_FD  = 0 : Fermi-Dirac distribution not calculated
!                      I_FD  = 1 : Fermi-Dirac distribution printed in file 'fermi_dir.dat'
! 
!       * I_BE      : switch for Bose-Einstein distribution printing
!                      I_BE  = 0 : Bose-Einstein distribution not calculated
!                      I_BE  = 1 : Bose-Einstein distribution printed in file 'bose_eins.dat'
! 
!       * I_MX      : switch for Maxwell distribution printing
!                      I_MX  = 0 : Maxwell distribution not calculated
!                      I_MX  = 1 : Maxwell distribution printed in file 'maxwell_d.dat' 
! 
!                                            ~
! 
!       * I_SC      : switch for scale parameters  printing
!                      I_SC  = 0 : scale parameters not calculated
!                      I_SC  = 1 : scale parameters printed in file 'scale_par.dat'
! 
!       * I_DS      : switch for density of states printing
!                      I_DS  = 0 : density of states not calculated
!                      I_DS  = 1 : density of states printed in file 'dens_stat.dat'
! 
!       * I_SB      : switch for  subband energies printing
!                      I_SB  = 0 : subband energies not calculated
!                      I_SB  = 1 : subband energies printed in file 'subb_ener.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_GP      : switch for grand partition function printing
!                      I_GP  = 0 : grand partition function not calculated
!                      I_GP  = 1 : grand partition function printed in file 'gran_part.dat'
! 
!       * I_PR      : switch for electronic pressure printing
!                      I_PR  = 0 : electronic pressure not calculated
!                      I_PR  = 1 : electronic pressure printed in file 'epressure.dat'
! 
!       * I_CO      : switch for compressibility printing
!                      I_CO  = 0 : compressibility not calculated
!                      I_CO  = 1 : compressibility printed in file 'comp_file.dat'
! 
!       * I_CP      : switch for chemical potential printing
!                      I_CP  = 0 : chemical potential not calculated
!                      I_CP  = 1 : chemical potential printed in file 'chem_pote.dat' 
! 
!                                            ~
! 
!       * I_BM      : switch for bulk modulus printing
!                      I_BM  = 0 : bulk modulus not calculated
!                      I_BM  = 1 : bulk modulus printed in file 'bulk_modu.dat'
! 
!       * I_SH      : switch for shear modulus printing
!                      I_SH  = 0 : shear modulus not calculated
!                      I_SH  = 1 : shear modulus printed in file 'shear_mod.dat'
! 
!       * I_S0      : switch for zero sound velocity printing
!                      I_S0  = 0 : zero sound velocity not calculated
!                      I_S0  = 1 : zero sound velocity printed in file 'zero_soun.dat' 
! 
!       * I_S1      : switch for first sound velocity printing
!                      I_S1  = 0 : first sound velocity not calculated
!                      I_S1  = 1 : first sound velocity printed in file 'firs_soun.dat' 
! 
!                                            ~
! 
!       * I_DT      : switch for Debye temperature printing
!                      I_DT  = 0 : Debye temperature not calculated
!                      I_DT  = 1 : Debye temperature printed in file 'Debye_tmp.dat'
! 
!       * I_PS      : switch for Pauli paramagnetic susceptibility printing
!                      I_PS  = 0 : Pauli paramagnetic susceptibility not calculated
!                      I_PS  = 1 : Pauli paramagnetic susceptibility printed in file 'para_susc.dat'
! 
!       * I_IE      : switch for internal energy printing
!                      I_IE  = 0 : internal energy not calculated
!                      I_IE  = 1 : internal energy printed in file 'inter_ene.dat'
! 
!       * I_EI      : switch for excess internal energy printing
!                      I_EI  = 0 : excess internal energy not calculated
!                      I_EI  = 1 : excess internal energy printed in file 'exces_ene.dat'
! 
!                                            ~
! 
!       * I_FH      : switch for Helmholtz free energy printing
!                      I_FH  = 0 : Helmholtz free energy not calculated
!                      I_FH  = 1 : Helmholtz free energy printed in file 'helm_free.dat'
! 
!       * I_EY      : switch for entropy printing
!                      I_EY  = 0 : entropy not calculated
!                      I_EY  = 1 : entropy printed in file 'entropy_f.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_EF      : switch for Fermi energy printing
!                      I_EF  = 0 : Fermi energy not calculated
!                      I_EF  = 1 : Fermi energy printed in file 'fermi_ene.dat'
! 
!       * I_KF      : switch for Fermi momentum printing
!                      I_KF  = 0 : Fermi momentum not calculated
!                      I_KF  = 1 : Fermi momentum printed in file 'fermi_vec.dat'
! 
!       * I_VF      : switch for Fermi velocity printing
!                      I_VF  = 0 : Fermi velocity not calculated
!                      I_VF  = 1 : Fermi velocity printed in file 'fermi_vel.dat'
! 
!       * I_TE      : switch for Fermi temperature printing
!                      I_TE  = 0 : Fermi temperature not calculated
!                      I_TE  = 1 : Fermi temperature printed in file 'fermi_tmp.dat'
! 
!                                            ~
! 
!       * I_DL      : switch for density of states at Fermi level printing
!                      I_DL  = 0 subband energies: density of states at Fermi level not calculated
!                      I_DL  = 1 : density of states at Fermi level printed in file 'fermi_dos.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_TW      : switch for thermal De Broglie wavelength printing
!                      I_TW  = 0 : thermal De Broglie wavelength not calculated
!                      I_TW  = 1 : thermal De Broglie wavelength printed in file 'thermal_w.dat'
! 
!       * I_VT      : switch for thermal velocity printing
!                      I_VT  = 0 : thermal velocity not calculated
!                      I_VT  = 1 : thermal velocity printed in file 'thermal_v.dat'
! 
!       * I_TC      : switch for thermal conductivity printing
!                      I_TC  = 0 : thermal conductivity not calculated
!                      I_TC  = 1 : thermal conductivity printed in file 'thermal_c.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_EG      : switch for ground state energy printing
!                      I_EG  = 0 : ground state energy not calculated
!                      I_EG  = 1 : ground state energy printed in file 'ground_st.dat'
! 
!       * I_EX      : switch for exchange energy printing
!                      I_EX  = 0 : exchange energy not calculated
!                      I_EX  = 1 : exchange energy printed in file 'ex_energy.dat'
! 
!       * I_XC      : switch for exchange correlation energy printing
!                      I_XC  = 0 : exchange correlation energy not calculated
!                      I_XC  = 1 : exchange correlation energy printed in file 'xc_energy.dat'
! 
!       * I_EC      : switch for correlation energy printing
!                      I_EC  = 0 : correlation energy not calculated
!                      I_EC  = 1 : correlation energy printed in file 'corr_ener.dat'
! 
!                                            ~
! 
!       * I_HF      : switch for Hartree-Fock energy printing
!                      I_HF  = 0 : Hartree-Fock energy not calculated
!                      I_HF  = 1 : Hartree-Fock energy printed in file 'hf_energy.dat'
! 
!       * I_EK      : switch for kinetic energy printing
!                      I_EK  = 0 : kinetic energy not calculated
!                      I_EK  = 1 : kinetic energy printed in file 'kine_ener.dat'
! 
!       * I_EP      : switch for potential energy printing
!                      I_EP  = 0 : potential energy not calculated
!                      I_EP  = 1 : potential energy printed in file 'pote_ener.dat
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_VI      : switch for shear viscosity printing
!                      I_VI  = 0 : shear viscosity not calculated
!                      I_VI  = 1 : shear viscosity printed in file 'visc_coef.dat'
! 
!       * I_DI      : switch for diffusion coefficient printing
!                      I_DI  = 0 : diffusion coefficient not calculated
!                      I_DI  = 1 : diffusion coefficient printed in file 'diff_coef.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_FP      : switch for fluctuation potential calculation
!                      I_FP  = 0 :  fluctuation potential not calculated
!                      I_FP  = 1 :  fluctuation potential printed in file 'fluct_pot.dat'
! 
!       * I_EL      : switch for loss function calculation
!                      I_EL  = 0 :  loss function not calculated
!                      I_EL  = 1 :  loss function printed in file 'ener_loss.dat'
! 
!       * I_PO      : switch for stopping power calculation
!                      I_PO  = 0 :  stopping power not calculated
!                      I_PO  = 1 :  stopping power printed in file 'stop_powe.dat'
! 
!       * I_RF      : switch for refractive index calculation
!                      I_RF  = 0 :  refractive index not calculated
!                      I_RF  = 1 :  refractive index printed in file 'refrac_in.dat' 
! 
!                                            ~
! 
!       * I_VC      : switch for dynamic screened Coulomb potential V(q,omega) calculation
!                      I_RF  = 0 : dynamic screened Coulomb potential not calculated
!                      I_RF  = 1 : dynamic screened Coulomb potential printed in file 'dyna_coul.dat'
! 
! *-------+---------+---------+---------+---------+---------------------------*
! 
!       * I_FN      : switch for appending the calculation type string 
!                        to the output filename
!                      I_FN  = 0 :  standard output filename
!                      I_FN  = 1 :  parameter added to filename
! 
!       * I_WR      : switch for writing physical properties into the log file 
!                      I_WR  = 0 :  does not write
!                      I_WR  = 1 :  writes
! 







 
