Removed surface_current attribute on Tally in favor of type attribute.

src/input_xml.F90

@@ -735,6 +735,9 @@ contains
        t % n_filter_bins = 0
        t % stride = 0
 
+       ! Set tally type to volume by default
+       t % type = TALLY_VOLUME
+
        ! Copy material id
        t % id = tally_(i) % id
 
@@ -892,7 +895,7 @@ contains
                 t % score_bins(j) % scalar = SCORE_NU_FISSION
              case ('current')
                 t % score_bins(j) % scalar = SCORE_CURRENT
-                t % surface_current = .true.
+                t % type = TALLY_SURFACE_CURRENT
 
                 ! Check to make sure that current is the only desired response
                 ! for this tally

src/tally.F90

@@ -61,7 +61,7 @@ contains
        filter_bins = 1
 
        ! handle surface current tallies separately
-       if (t % surface_current) then
+       if (t % type == TALLY_SURFACE_CURRENT) then
           m => meshes(t % mesh)
 
           t % stride(SURF_FILTER_SURFACE) = filter_bins
@@ -273,7 +273,7 @@ contains
        t => tallies(i)
 
        ! Surface current tallies are treated separately
-       if (t % surface_current) cycle
+       if (t % type == TALLY_SURFACE_CURRENT) cycle
 
        ! =======================================================================
        ! DETERMINE SCORING BIN COMBINATION
@@ -552,7 +552,7 @@ contains
        t => tallies(i)
 
        ! Skip non-surface-current tallies
-       if (.not. t % surface_current) cycle
+       if (t % type /= TALLY_SURFACE_CURRENT) cycle
 
        ! Determine indices for starting and ending location
        m => meshes(t % mesh)
@@ -969,7 +969,7 @@ contains
        call header("TALLY " // trim(to_str(t % id)), unit=UNIT_TALLY, level=3)
 
        ! Handle surface current tallies separately
-       if (t % surface_current) then
+       if (t % type == TALLY_SURFACE_CURRENT) then
           call write_surface_current(t)
           cycle
        end if

src/tally_header.F90

@@ -62,7 +62,6 @@ module tally_header
      integer :: type      ! volume, surface current
      integer :: estimator ! collision, track-length
      real(8) :: volume    ! volume of region
-     logical :: surface_current = .false.
 
      ! Tally bin specifications