Filling out the F71Stream.

radixio/f71stream.i.hh

@@ -40,10 +40,10 @@ bool F71Stream::load_next(f71case_type& inventory)
   mStream.readInt();                                   // NOBLND =
   mStream.readInt();                                   // NDSET  =
   mStream.readInt();                                   // NTYPE  =
-  mStream.readInt();                                   // NGRP   =
+  int ngrp = mStream.readInt();                        // NGRP   =
   mStream.readInt();                                   // NELEM  =
   mStream.readInt();                                   // NVERT  =
-  mStream.readInt();                                   // NG     =
+  int ng = mStream.readInt();                          // NG     =
   mStream.readInt();                                   // MMN    =
   mStream.readInt();                                   // MOUT   =
   mStream.readInt();                                   // INDEX  =
@@ -53,10 +53,10 @@ bool F71Stream::load_next(f71case_type& inventory)
   int num_light_element_groups   = mStream.readInt();  // NENLE  =
   int num_actinide_groups        = mStream.readInt();  // NENAC  =
   int num_fission_product_groups = mStream.readInt();  // NENFP  =
-  mStream.readInt();                                   // L1     =
-  mStream.readInt();                                   // L2     =
-  mStream.readInt();                                   // L3     =
-  mStream.readInt();                                   // L4     =
+  int l1                         = mStream.readInt();  // L1     =
+  int l2                         = mStream.readInt();  // L2     =
+  int l3                         = mStream.readInt();  // L3     =
+  int l4                         = mStream.readInt();  // L4     =
 
   mStream.readFloat();    // TMO    =
   mStream.readFloat();    // FRACPW =
@@ -78,13 +78,62 @@ bool F71Stream::load_next(f71case_type& inventory)
   mStream >> abundances;
   std::string title = trim_string(mStream.readString(80));
   std::string basis = trim_string(mStream.readString(40));
+
   std::vector<float> gamma_le_bounds((size_t(num_light_element_groups)));
   mStream >> gamma_le_bounds;
+
   std::vector<float> gamma_ae_bounds((size_t(num_actinide_groups)));
   mStream >> gamma_ae_bounds;
+
   std::vector<float> gamma_fpe_bounds((size_t(num_fission_product_groups)));
   mStream >> gamma_fpe_bounds;
 
+  std::vector<float> spectrum_le((size_t(l1)));
+  mStream >> spectrum_le;
+
+  std::vector<float> spectrum_ac((size_t(l2)));
+  mStream >> spectrum_ac;
+
+  std::vector<float> spectrum_fp((size_t(l3)));
+  mStream >> spectrum_fp;
+
+  std::vector<float> energy_bounds(ng + 1);
+  mStream >> energy_bounds;
+
+  std::vector<float> dsav(ng);
+  mStream >> dsav;
+
+  std::vector<float> esav(ng);
+  mStream >> esav;
+
+  // photon spectrum title
+  std::string itsav = trim_string(mStream.getString(80));
+
+  // values of total neutron spectrum in neutron/s
+  std::vector<float> spnnuc(l4);
+  mStream >> spnnuc;
+
+  std::vector<float> alpnuc(l4);
+  mStream >> alpnuc;
+
+  // values of total neutron spectrum in neutrons/s
+  std::vector<float> spneut(ngrp);
+  mStream >> spneut;
+
+  // values of (alpha,n) spectrum in neutrons/s
+  std::vector<float> specan(ngrp);
+  mStream >> specan;
+
+  // values of spontaneous fission neutron spectrum in neutrons/s
+  std::vector<float> specsp(ngrp);
+  mStream >> specsp;
+
+  std::vector<float> eneuts(ngrp + 1);
+  mStream >> eneuts;
+
+  // read end of record
+  mStream.readFooter();
+
   mStream.peek();  // peek into the next bit. This ensures EOF is seen in good()
                    // call
 }