From 682f80f4357f30780939ac84b8e2937233c647da Mon Sep 17 00:00:00 2001 From: Alain Hebert Date: Fri, 19 Dec 2025 10:57:23 -0500 Subject: #20: Correct a NTOT1 issue in SPN Bivac/Trivac solutions --- Trivac/src/TRIRCA.f | 14 +++++++++----- 1 file changed, 9 insertions(+), 5 deletions(-) (limited to 'Trivac/src/TRIRCA.f') diff --git a/Trivac/src/TRIRCA.f b/Trivac/src/TRIRCA.f index 34cf544..8202f29 100755 --- a/Trivac/src/TRIRCA.f +++ b/Trivac/src/TRIRCA.f @@ -1,5 +1,6 @@ *DECK TRIRCA - SUBROUTINE TRIRCA(IPMACR,IPMACP,NGRP,NBMIX,NANI,LDIFF,IL,IPR,RCAT) + SUBROUTINE TRIRCA(IPMACR,IPMACP,NGRP,NBMIX,NANI,NW,LDIFF,IL,IPR, + 1 RCAT) * *----------------------------------------------------------------------- * @@ -22,6 +23,7 @@ * NGRP number of energy groups. * NBMIX total number of material mixtures in the macrolib. * NANI maximum scattering order recovered from tracking and macrolib. +* NW maximum Legendre order (0 or 1) for the total cross sections. * LDIFF flag set to .true. to use 1/3D as 'NTOT1' cross sections. * IL scattering Legendre order. * IPR type of assembly: @@ -41,7 +43,7 @@ * SUBROUTINE ARGUMENTS *---- TYPE(C_PTR) IPMACR,IPMACP - INTEGER NGRP,NBMIX,NANI,IL,IPR + INTEGER NGRP,NBMIX,NANI,NW,IL,IPR LOGICAL LDIFF DOUBLE PRECISION RCAT(NGRP,NGRP,NBMIX) *---- @@ -59,7 +61,7 @@ *---- ALLOCATE(IJJ(NBMIX),NJJ(NBMIX),IPOS(NBMIX)) ALLOCATE(SGD(NBMIX,3),WORK(NBMIX*NGRP)) -* +* JPMACR=LCMGID(IPMACR,'GROUP') JPMACP=LCMGID(IPMACP,'GROUP') WRITE(CM,'(I2.2)') IL-1 @@ -128,9 +130,11 @@ ENDIF GO TO 100 ELSE - IF(LENGT.EQ.NBMIX) THEN + IF(NW.EQ.0) THEN + CALL LCMGET(KPMACP,'NTOT0',SGD(1,2)) + ELSE IF(LENGT.EQ.NBMIX) THEN CALL LCMGET(KPMACP,TEXT12,SGD(1,2)) - ELSE IF(LENGT1.EQ.NBMIX) THEN + ELSE IF((NW.GE.1).AND.(LENGT1.EQ.NBMIX)) THEN CALL LCMGET(KPMACP,'NTOT1',SGD(1,2)) ELSE CALL LCMGET(KPMACP,'NTOT0',SGD(1,2)) -- cgit v1.2.3