/* Copyright 2019 Axel Huebl, Maxence Thevenet
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
// This script was automatically generated!
// Edit dev/Source/Utils/write_atomic_data_cpp.py instead!
#ifndef WARPX_IONIZATION_TABLE_H_
#define WARPX_IONIZATION_TABLE_H_

#include <AMReX_AmrCore.H>
#include <map>

std::map<std::string, int> ion_map_ids = {
    {"H", 0},
    {"He", 1},
    {"Li", 2},
    {"Be", 3},
    {"B", 4},
    {"C", 5},
    {"N", 6},
    {"O", 7},
    {"F", 8},
    {"Ne", 9},
    {"Na", 10},
    {"Mg", 11},
    {"Al", 12},
    {"Si", 13},
    {"P", 14},
    {"S", 15},
    {"Cl", 16},
    {"Ar", 17},
    {"Kr", 18},
    {"Rb", 19},
    {"Xe", 20},
    {"Rn", 21} };

const int nelements = 22;

const int ion_atomic_numbers[nelements] = {
    1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
    11, 12, 13, 14, 15, 16, 17, 18, 36, 37,
    54, 86};

const int ion_energy_offsets[nelements] = {
    0, 1, 3, 6, 10, 15, 21, 28, 36, 45,
    55, 66, 78, 91, 105, 120, 136, 153, 171, 207,
    244, 298};

const int energies_tab_length = 384;

const amrex::Real table_ionization_energies[energies_tab_length]{
    // H
    13.59843449,
    // He
    24.58738880, 54.4177650,
    // Li
    5.39171495, 75.6400964, 122.4543581,
    // Be
    9.322699, 18.21115, 153.896203, 217.7185843,
    // B
    8.298019, 25.15483, 37.93058, 259.3715, 340.226020,
    // C
    11.2602880, 24.383154, 47.88778, 64.49352, 392.090515, 489.993194,
    // N
    14.53413, 29.60125, 47.4453, 77.4735, 97.8901, 552.06732, 667.046116,
    // O
    13.618055, 35.12112, 54.93554, 77.41350, 113.8990, 138.1189, 739.32682,
    871.40988,
    // F
    17.42282, 34.97081, 62.70798, 87.175, 114.249, 157.16311, 185.1868,
    953.89804, 1103.11747,
    // Ne
    21.564540, 40.96297, 63.4233, 97.1900, 126.247, 157.934, 207.271,
    239.0970, 1195.80783, 1362.19915,
    // Na
    5.1390769, 47.28636, 71.6200, 98.936, 138.404, 172.23, 208.504,
    264.192, 299.856, 1465.13449, 1648.70218,
    // Mg
    7.646236, 15.035271, 80.1436, 109.2654, 141.33, 186.76, 225.02,
    265.924, 327.99, 367.489, 1761.80487, 1962.66365,
    // Al
    5.985769, 18.82855, 28.447642, 119.9924, 153.8252, 190.49, 241.76,
    284.64, 330.21, 398.65, 442.005, 2085.97700, 2304.14005,
    // Si
    8.15168, 16.34585, 33.49300, 45.14179, 166.767, 205.279, 246.57,
    303.59, 351.28, 401.38, 476.273, 523.415, 2437.65813, 2673.17753,
    // P
    10.486686, 19.76949, 30.20264, 51.44387, 65.02511, 220.430, 263.57,
    309.60, 372.31, 424.40, 479.44, 560.62, 611.741, 2816.90876,
    3069.8415,
    // S
    10.36001, 23.33788, 34.86, 47.222, 72.5945, 88.0529, 280.954,
    328.794, 379.84, 447.7, 504.55, 564.41, 651.96, 706.994,
    3223.7807, 3494.1879,
    // Cl
    12.967632, 23.81364, 39.80, 53.24, 67.68, 96.94, 114.2013,
    348.306, 400.851, 456.7, 530.0, 591.58, 656.30, 750.23,
    809.198, 3658.3437, 3946.2909,
    // Ar
    15.7596117, 27.62967, 40.735, 59.58, 74.84, 91.290, 124.41,
    143.4567, 422.60, 479.76, 540.4, 619.0, 685.5, 755.13,
    855.5, 918.375, 4120.6656, 4426.2228,
    // Kr
    13.9996053, 24.35984, 35.838, 50.85, 64.69, 78.49, 109.13,
    125.802, 233.0, 268, 308, 350, 391, 446,
    492, 540, 591, 640, 785, 831.6, 882.8,
    945, 999.0, 1042, 1155.0, 1205.23, 2928.9, 3072,
    3228, 3380, 3584, 3752.0, 3971, 4109.083, 17296.420,
    17936.209,
    // Rb
    4.1771280, 27.28954, 39.247, 52.20, 68.44, 82.9, 98.67,
    132.79, 150.628, 277.12, 313.1, 356.0, 400, 443,
    502, 550, 601, 654, 706.0, 857, 905.3,
    958.9, 1024, 1080, 1125, 1242.5, 1294.57, 3133.3,
    3281, 3443, 3600, 3815, 3988, 4214, 4356.865,
    18305.884, 18965.516,
    // Xe
    12.1298436, 20.975, 31.05, 42.20, 54.1, 66.703, 91.6,
    105.9778, 179.84, 202.0, 229.02, 255.0, 281, 314,
    343, 374, 404, 434, 549, 582, 616,
    650, 700, 736, 818, 857.0, 1493, 1571,
    1653, 1742, 1826, 1919, 2023, 2113, 2209,
    2300, 2556, 2637, 2726, 2811, 2975, 3068,
    3243, 3333.8, 7660, 7889, 8144, 8382, 8971,
    9243, 9581, 9810.37, 40271.724, 41299.71,
    // Rn
    10.74850, 21.4, 29.4, 36.9, 52.9, 64.0, 88.0,
    102.0, 154.0, 173.9, 195.0, 218.0, 240, 264,
    293, 317, 342, 367, 488, 520, 550,
    580, 640, 680, 760, 800, 850, 920,
    980, 1050, 1110, 1180, 1250, 1310, 1390,
    1460, 1520, 1590, 1660, 1720, 2033, 2094,
    2158, 2227, 2293, 2357, 2467, 2535, 2606,
    2674, 2944, 3010, 3082, 3149, 3433, 3510,
    3699, 3777, 6169, 6318, 6476, 6646, 6807,
    6964, 7283, 7450, 7630, 7800, 8260, 8410,
    8570, 8710, 9610, 9780, 10120, 10290, 21770,
    22160, 22600, 22990, 26310, 26830, 27490, 27903.1,
    110842.0, 112843.7
};

#endif // #ifndef WARPX_IONIZATION_TABLE_H_