Start incremental block phixs

globals.h

@@ -88,6 +88,7 @@ struct EnergyLevel {
   int ndowntrans{0};
   PhotoionTarget *phixstargets{};  // pointer to table of target states and probabilities
   float *photoion_xs{};            // Pointer to a lookup-table providing photoionisation cross-sections for this level.
+  int phixsstart{-1};              // index to start of photoionisation cross-sections table in global::allphixs
   int nphixstargets{0};            // length of phixstargets array:
   float stat_weight{0.};           // Statistical weight of this level.
 

input.cc

@@ -87,6 +87,8 @@ constexpr std::array<std::string_view, 24> inputlinecomments = {
 
 CellCachePhixsTargets *chphixstargetsblock{};
 
+std::vector<float> tmpallphixs;
+
 void read_phixs_data_table(std::fstream &phixsfile, const int nphixspoints_inputtable, const int element,
                            const int lowerion, const int lowerlevel, const int upperion, int upperlevel_in,
                            size_t *mem_usage_phixs, const int phixs_file_version) {
@@ -175,6 +177,8 @@ void read_phixs_data_table(std::fstream &phixsfile, const int nphixspoints_input
   }
 
   *mem_usage_phixs += globals::NPHIXSPOINTS * sizeof(float);
+  globals::elements[element].ions[lowerion].levels[lowerlevel].phixsstart = tmpallphixs.size();
+  tmpallphixs.resize(tmpallphixs.size() + globals::NPHIXSPOINTS);
   globals::elements[element].ions[lowerion].levels[lowerlevel].photoion_xs =
       static_cast<float *>(malloc(globals::NPHIXSPOINTS * sizeof(float)));
   assert_always(globals::elements[element].ions[lowerion].levels[lowerlevel].photoion_xs != nullptr);