diff options
| author | HEBERT Alain <alain.hebert@polymtl.ca> | 2025-12-19 11:17:25 -0500 |
|---|---|---|
| committer | HEBERT Alain <alain.hebert@polymtl.ca> | 2025-12-19 11:17:25 -0500 |
| commit | 0086fa5c672a20b4fcefea39fca31313b9c08d68 (patch) | |
| tree | 837b47785e7876046a80ea6391c5d3b445b2d5ba /Trivac/src/BIVSPS.f | |
| parent | 338c590068f0f867d89c78598af32f121df1936a (diff) | |
| parent | 682f80f4357f30780939ac84b8e2937233c647da (diff) | |
Merge branch '20-correct-a-ntot1-issue-in-spn-bivac-trivac-solutions' into 'main'
#20: Correct a NTOT1 issue in SPN Bivac/Trivac solutions
See merge request dragon/5.1!34
Diffstat (limited to 'Trivac/src/BIVSPS.f')
| -rwxr-xr-x | Trivac/src/BIVSPS.f | 9 |
1 files changed, 6 insertions, 3 deletions
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 |