!
! 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
!
Sylvain Tricot
04e936c881