!    This module contains information on the grid employed by SP_Ace 
!    to move through the parameter space and some absorption lines that 
!    behave differently from the all the others contained in the line list.
!
!    It is part of the program SP_Ace, which derives stellar parameters, 
!    such as gravity, temperature, and element abundances from optical
!    stellar spectra, assuming Local Thermodynamic Equilibrium (LTE) 
!    and 1D stellar atmosphere models.
!
!    Copyright (C) 2016 Corrado Boeche
!    
!    This program is free software: you can redistribute it and/or modify
!    it under the terms of the GNU General Public License as published by
!    the Free Software Foundation, either version 3 of the License, or
!    (at your option) any later version.

!
!    This program is distributed in the hope that it will be useful,
!    but WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!    GNU General Public License for more details.
!    
!    You should have received a copy of the GNU General Public License
!    along with this program.  If not, see <http://www.gnu.org/licenses/>.


MODULE data_lib
        USE num_type
        IMPLICIT NONE

        REAL(DP), DIMENSION(16), PARAMETER :: temp_grid=&
          &(/4000.,4200.,4400.,4600.,4800.,5000.,5200.,5400.,5600. &
          & ,5800.,6000.,6200.,6400.,6600.,6800.,7000./)
        REAL(DP), DIMENSION(10), PARAMETER :: logg_grid=&
          &(/1.0,1.4,1.8,2.2,2.6,3.0,3.4,3.8,4.2,4.6/)
        REAL(DP), DIMENSION(11), PARAMETER :: met_grid=&
          &(/-2.0,-1.8,-1.6,-1.4,-1.2,-1.0,-0.8,-0.6,-0.4,-0.2,1e-4/)
        REAL(DP), DIMENSION(20), PARAMETER :: temp_gridL=&
          &(/3600.,3800.,4000.,4200.,4400.,4600.,4800.,5000.,5200.,5400.,5600. &
          & ,5800.,6000.,6200.,6400.,6600.,6800.,7000.,7200.,7400./)
        REAL(DP), DIMENSION(14), PARAMETER :: logg_gridL=&
          &(/0.2,0.6,1.0,1.4,1.8,2.2,2.6,3.0,3.4,3.8,4.2,4.6,5.0,5.4/)
        REAL(DP), DIMENSION(15), PARAMETER :: met_gridL=&
          &(/-2.4,-2.2,-2.0,-1.8,-1.6,-1.4,-1.2,-1.0,-0.8,-0.6,-0.4,-0.2,0.0,0.2,0.4/)
        REAL(DP), DIMENSION(6), PARAMETER :: aFe_gridL=&
          &(/-0.4,-0.2,0.0,0.2,0.4,0.6/)
        REAL(DP), DIMENSION(89), PARAMETER :: w_str_ll=&
          (/5229.845,5229.866,5269.537,5341.024,5345.796,5353.374,5371.489,5397.128,5404.117,5404.151,&
          &5405.775,5409.808,5415.199,5429.696,5434.524,5446.916,5476.900,5481.396,5487.740,5487.745,&
          &5497.516,5501.465,5506.779,5586.756,5602.842,5602.945,5615.644,5624.542,5675.410,5675.417,&
          &5791.017,5889.951,5895.924,5948.541,5949.282,6016.632,6122.217,6141.731,6162.173,6230.722,&
          &6256.361,6318.018,6318.109,6343.308,6361.786,6439.075,6546.238,6562.797,8401.400,8413.322,&
          &8426.506,8434.950,8434.957,8437.959,8446.359,8446.758,8467.258,8496.990,8502.487,8514.072,&
          &8536.164,8538.015,8538.175,8545.387,8556.777,8598.396,8648.465,8665.022,8680.282,8680.411,&
          &8712.682,8717.815,8728.010,8736.016,8742.446,8750.476,8752.007,8763.966,8772.865,8773.897,&
          &8790.389,8806.756,8824.220,8862.787,8866.931,8892.720,8899.231,8912.068,8923.555/)
        REAL(DP), DIMENSION(89), PARAMETER :: dd_str_ll=&
          (/4.0,4.0,0.9,1.5,0.9,1.4,0.9,0.9,2.0,2.0,0.9,4.0,1.1,0.9,0.8,0.8,0.9,5.0,2.0,2.0,&
          &0.9,0.9,0.85,0.9,1.5,1.5,0.9,1.3,2.0,2.0,3.0,0.9,0.9,1.5,1.5,2.0,0.9,5.0,0.9,1.1,&
          &1.7,2.0,5.0,8.0,6.0,0.9,1.5,1.5,3.0,20.0,0.8,3.0,0.8,20.0,1.5,1.5,20.0,3.0,20.0,&
          &0.8,3.0,3.0,3.0,20.0,2.0,10.0,1.5,15.0,2.5,2.5,2.0,2.0,2.0,1.5,1.5,10.0,1.5,1.3,&
          &2.0,1.6,2.0,1.2,0.85,5.0,1.2,1.5,1.5,1.3,1.5/)
         REAL(DP), DIMENSION(7), PARAMETER :: w_rej_op=(/5226.7,5267.0,5328.5,5341.0,5371.5,5429.5,6318.0/)
         REAL(DP), DIMENSION(7), PARAMETER :: r_rej_op=(/2.0,6.0,3.5,0.5,1.5,1.5,1.1/)
         REAL(DP), DIMENSION(11), PARAMETER :: w_rej_nlte=&
         &(/5269.537,5328.039,6562.8,8413.322,8437.959,8467.258,8498.023,8542.091,8598.396,8662.141,8862.787/)!,8750.476/)
         REAL(DP), DIMENSION(11), PARAMETER :: r_rej_nlte=(/0.5,0.5,20.,5.,5.,5.,10.,10.,10.,10.,10./)!,10./)
         REAL(DP), DIMENSION(10), PARAMETER :: w_rej_unknown=&
         &(/5394.670,5407.455,5420.355,&
         &5470.580,5503.070,5637.405,5991.870,6019.970,6042.040,6089.565/)
         REAL(DP), DIMENSION(10), PARAMETER :: r_rej_unknown=&
          &(/0.2,0.3,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2/)
         REAL(DP), DIMENSION(10), PARAMETER :: w_rej_bad=&
         &(/5212.219,5298.5,5369.8,5403.0,5415.2,5890.0,5896.0,6343.3,6361.8,8806.756/)
         REAL(DP), DIMENSION(10), PARAMETER :: r_rej_bad=&
         &(/0.2,1.5,0.2,9.0,0.2,7.0,7.0,2.0,2.0,1.0/)

        CHARACTER(LEN=6),DIMENSION(95), PARAMETER :: ELE_symb=&
        &(/' [H/H]','[He/H]','[Li/H]','[Be/H]',' [B/H]',' [C/H]',&
        &' [N/H]',' [O/H]',' [F/H]','[Ne/H]','[Na/H]','[Mg/H]','[Al/H]',&
        &'[Si/H]',' [P/H]',' [S/H]','[Cl/H]','[Ar/H]',' [K/H]','[Ca/H]',&
        &'[Sc/H]','[Ti/H]',' [V/H]','[Cr/H]','[Mn/H]','[Fe/H]','[Co/H]',&
        &'[Ni/H]','[Cu/H]','[Zn/H]','[Ga/H]','[Ge/H]','[As/H]','[Se/H]',&
        &'[Br/H]','[Kr/H]','[Rb/H]','[Sr/H]',' [Y/H]','[Zr/H]','[Nb/H]',&
        &'[Mo/H]','[Tc/H]','[Ru/H]','[Rh/H]','[Pd/H]','[Ag/H]','[Cd/H]',&
        &'[In/H]','[Sn/H]','[Sb/H]','[Te/H]',' [I/H]','[Xe/H]','[Cs/H]',&
        &'[Ba/H]','[La/H]','[Ce/H]','[Pr/H]','[Nd/H]','[Pm/H]','[Sm/H]',&
        &'[Eu/H]','[Gd/H]','[Tb/H]','[Dy/H]','[Ho/H]','[Er/H]','[Tm/H]',&
        &'[Yb/H]','[Lu/H]','[Hf/H]','[Ta/H]','[Wl/H]','[Re/H]','[Os/H]',&
        &'[Ir/H]','[Pt/H]','[Au/H]','[Hg/H]','[Tl/H]','[Pb/H]','[Bi/H]',&
        &'[Po/H]','[At/H]','[Rn/H]','[Fr/H]','[Ra/H]','[Ac/H]','[Th/H]',&
        &'[Pa/H]',' [U/H]',' [m/H]',' [a/H]','[CN/H]'/)



END MODULE data_lib