From 682f80f4357f30780939ac84b8e2937233c647da Mon Sep 17 00:00:00 2001
From: Alain Hebert <alainhebert@Alains-MacBook-Pro.local>
Date: Fri, 19 Dec 2025 10:57:23 -0500
Subject: #20: Correct a NTOT1 issue in SPN Bivac/Trivac solutions

---
 Trivac/src/BIVSPS.f | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

(limited to 'Trivac/src/BIVSPS.f')

diff --git a/Trivac/src/BIVSPS.f b/Trivac/src/BIVSPS.f
index d9e1ff4..e03ce67 100755
--- a/Trivac/src/BIVSPS.f
+++ b/Trivac/src/BIVSPS.f
@@ -1,6 +1,6 @@
 *DECK BIVSPS
-      SUBROUTINE BIVSPS(IPTRK,IPMACR,IPSYS,IMPX,NGRP,NEL,NLF,NANI,NBFIS,
-     1 NALBP,LDIFF,MAT,VOL,NBMIX)
+      SUBROUTINE BIVSPS(IPTRK,IPMACR,IPSYS,IMPX,NGRP,NEL,NLF,NANI,NW,
+     1 NBFIS,NALBP,LDIFF,MAT,VOL,NBMIX)
 *
 *-----------------------------------------------------------------------
 *
@@ -27,6 +27,7 @@
 * NEL     total number of finite elements.
 * NLF     number of Legendre orders for the flux (even number).
 * NANI    number of Legendre orders for the scattering cross sections.
+* NW      maximum Legendre order (0 or 1) for the total cross sections.
 * NBFIS   number of fissionable isotopes.
 * NALBP   number of physical albedos per energy group.
 * LDIFF   flag set to .true. to use 1/3D as 'NTOT1' cross sections.
@@ -41,7 +42,7 @@
 *  SUBROUTINE ARGUMENTS
 *----
       TYPE(C_PTR) IPTRK,IPMACR,IPSYS
-      INTEGER IMPX,NGRP,NEL,NLF,NANI,NBFIS,NALBP,MAT(NEL),NBMIX
+      INTEGER IMPX,NGRP,NEL,NLF,NANI,NW,NBFIS,NALBP,MAT(NEL),NBMIX
       REAL VOL(NEL)
       LOGICAL LDIFF
 *----
@@ -126,6 +127,8 @@
             DO 5 IBM=1,NBMIX
             SGD(IBM,2)=1.0/(3.0*SGD(IBM,2))
     5       CONTINUE
+         ELSE IF(NW.EQ.0) THEN
+            CALL LCMGET(KPMACR,'NTOT0',SGD(1,2))
          ELSE IF(LENGT.EQ.NBMIX) THEN
             CALL LCMGET(KPMACR,TEXT12,SGD(1,2))
          ELSE IF(LENGT1.EQ.NBMIX) THEN
-- 
cgit v1.2.3