Get rid of n_nuclides dataset in summary/statepoint

docs/source/usersguide/output/statepoint.rst

@@ -206,10 +206,6 @@ if (run_mode == MODE_EIGENVALUE)
 
             Value for each filter bin of this type.
 
-    **/tallies/tally i/n_nuclides** (*int*)
-
-        Number of nuclide bins. If none are specified, this is just one.
-
     **/tallies/tally i/nuclides** (*char[][]*)
 
         Values of specified nuclide bins.

docs/source/usersguide/output/summary.rst

@@ -219,8 +219,6 @@ do i = 1, n_tallies
 
     end do
 
-    **/tallies/tally <uid>/n_nuclides** (*int*)
-
     **/tallies/tally <uid>/nuclides** (*char[][]*)
 
     **/tallies/tally <uid>/n_score_bins** (*int*)

openmc/statepoint.py

@@ -408,8 +408,6 @@ class StatePoint(object):
                     tally.add_filter(filter)
 
                 # Read Nuclide bins
-                n_nuclides = self._f['{0}{1}/n_nuclides'.format(base, tally_key)].value
-
                 nuclide_names = self._f['{0}{1}/nuclides'.format(base, tally_key)].value
 
                 # Add all Nuclides to the Tally
@@ -430,7 +428,7 @@ class StatePoint(object):
                 # Compute and set the filter strides
                 for i in range(n_filters):
                     filter = tally.filters[i]
-                    filter.stride = n_score_bins * n_nuclides
+                    filter.stride = n_score_bins * len(nuclide_names)
 
                     for j in range(i+1, n_filters):
                         filter.stride *= tally.filters[j].num_bins

src/state_point.F90

@@ -276,8 +276,6 @@ contains
             call close_group(filter_group)
           end do FILTER_LOOP
 
-          call write_dataset(tally_group, "n_nuclides", tally%n_nuclide_bins)
-
           ! Set up nuclide bin array and then write
           allocate(str_array(tally%n_nuclide_bins))
           NUCLIDE_LOOP: do j = 1, tally%n_nuclide_bins

src/summary.F90

@@ -551,9 +551,6 @@ contains
         call close_group(filter_group)
       end do FILTER_LOOP
 
-      ! Write number of nuclide bins
-      call write_dataset(tally_group, "n_nuclides", t%n_nuclide_bins)
-
       ! Create temporary array for nuclide bins
       allocate(str_array(t%n_nuclide_bins))
       NUCLIDE_LOOP: do j = 1, t%n_nuclide_bins