From 744b40856a035580b786378cae13d453edd26689 Mon Sep 17 00:00:00 2001
From: HEBERT Alain <alain.hebert@polymtl.ca>
Date: Sun, 28 Dec 2025 15:55:41 -0500
Subject: Resolve "Depreciate use of Version 4 and 5.0 Draglibs"

---
 Dragon/src/BRERT.f | 19 +++++++++++++------
 1 file changed, 13 insertions(+), 6 deletions(-)

(limited to 'Dragon/src/BRERT.f')

diff --git a/Dragon/src/BRERT.f b/Dragon/src/BRERT.f
index 19afde9..91e2b8b 100644
--- a/Dragon/src/BRERT.f
+++ b/Dragon/src/BRERT.f
@@ -1,7 +1,7 @@
 *DECK BRERT
       SUBROUTINE BRERT(IPMAC1,IELEM,ICOL,NG,NL,LX1,NMIX1,IMIX,ICODE,
-     1 ISPH,IDIFF,ZKEFF,B2,ENER,XXX1,VOL1,FLX1,DC1,TOT1,CHI1,SIGF1,
-     2 SCAT1,JXM,JXP,FHETXM,FHETXP,ADF1,NGET,ADFREF,IMPX)
+     1 ISPH,IH,IDIFF,ZKEFF,B2,ENER,XXX1,VOL1,FLX1,DC1,TOT1,CHI1,SIGF1,
+     2 SCAT1,HFACT1,JXM,JXP,FHETXM,FHETXP,ADF1,NGET,ADFREF,IMPX)
 *
 *-----------------------------------------------------------------------
 *
@@ -30,6 +30,7 @@
 * IMIX    mix index of each node.
 * ICODE   physical albedo index on each side of the domain.
 * ISPH    SPH flag (=0: use discontinuity factors; =1: use SPH factors).
+* IH      H-FACTOR flag (=0: not used; =1: recovered).
 * IDIFF   PN calculation option (=0: diffusion theory; =1: SPN theory
 *         with 'NTOT1'; =2: SPN theory with 1/(3*D)).
 * ZKEFF   effective multiplication factor.
@@ -43,6 +44,7 @@
 * CHI1    fission spectra.
 * SIGF1   nu*fission cross sections.
 * SCAT1   scattering P0 cross sections.
+* HFACT1  H-FACTOR values.
 * JXM     left boundary currents.
 * JXP     right boundary currents.
 * FHETXM  left boundary fluxes.
@@ -61,12 +63,13 @@
 *  SUBROUTINE ARGUMENTS
 *----
       TYPE(C_PTR) IPMAC1
-      INTEGER IELEM,ICOL,NG,NL,LX1,NMIX1,IMIX(LX1),ICODE(2),ISPH,IDIFF,
-     1 NGET,IMPX
+      INTEGER IELEM,ICOL,NG,NL,LX1,NMIX1,IMIX(LX1),ICODE(2),ISPH,IH,
+     1 IDIFF,NGET,IMPX
       REAL ZKEFF,B2,ENER(NG+1),XXX1(LX1+1),VOL1(NMIX1),FLX1(NMIX1,NG),
      1 DC1(NMIX1,NG),TOT1(NMIX1,NG,NL),CHI1(NMIX1,NG),SIGF1(NMIX1,NG),
-     2 SCAT1(NMIX1,NG,NG,NL),JXM(NMIX1,NG),JXP(NMIX1,NG),
-     3 FHETXM(NMIX1,NG,NL),FHETXP(NMIX1,NG,NL),ADF1(NMIX1,NG),ADFREF(NG)
+     2 SCAT1(NMIX1,NG,NG,NL),HFACT1(NMIX1,NG),JXM(NMIX1,NG),
+     3 JXP(NMIX1,NG),FHETXM(NMIX1,NG,NL),FHETXP(NMIX1,NG,NL),
+     4 ADF1(NMIX1,NG),ADFREF(NG)
 *----
 *  LOCAL VARIABLES
 *----
@@ -298,6 +301,9 @@
                SCAT1(IBM,IGR,JGR,IL)=SCAT1(IBM,IGR,JGR,IL)*FDXM(IBM,IGR)
               ENDDO
             ENDDO
+            IF(IH.EQ.1) THEN
+              HFACT1(IBM,IGR)=HFACT1(IBM,IGR)/FDXM(IBM,IGR)
+            ENDIF
           ENDDO
         ENDDO
         IF(ICODE(2).NE.0) THEN
@@ -409,6 +415,7 @@
           CALL LCMPUT(KPMAC1,'IJJS'//CM,NMIX1,1,IJJ)
           CALL LCMPUT(KPMAC1,'IPOS'//CM,NMIX1,1,IPOS)
         ENDDO
+        IF(IH.EQ.1) CALL LCMPUT(KPMAC1,'H-FACTOR',NMIX1,2,HFACT1(:,IGR))
       ENDDO
 *----
 *  SCRATCH STORAGE DEALLOCATION
-- 
cgit v1.2.3