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/BRELLB.f | 14 +++++++++++---
 1 file changed, 11 insertions(+), 3 deletions(-)

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

diff --git a/Dragon/src/BRELLB.f b/Dragon/src/BRELLB.f
index 205a994..e8ebcd1 100644
--- a/Dragon/src/BRELLB.f
+++ b/Dragon/src/BRELLB.f
@@ -1,5 +1,6 @@
 *DECK BRELLB
-      SUBROUTINE BRELLB(IPMAC1,NC,NG,NL,NMIX1,ENER,JXM,FHETXM,IPRINT)
+      SUBROUTINE BRELLB(IPMAC1,NC,NG,NL,NMIX1,IH,ENER,HFACT1,JXM,
+     1 FHETXM,IPRINT)
 *
 *-----------------------------------------------------------------------
 *
@@ -22,7 +23,9 @@
 * NL      Legendre order of TOT1 and SCAT1 arrays (=1 for isotropic
 *         scattering in LAB).
 * NMIX1   number of mixtures in the nodal calculation.
+* IH      H-FACTOR flag (=0: not used; =1: recovered).
 * ENER    energy limits.
+* HFACT1  H-FACTOR values.
 * JXM     left boundary currents.
 * FHETXM  left boundary fluxes.
 * IPRINT  edition flag.
@@ -38,8 +41,9 @@
 *  SUBROUTINE ARGUMENTS
 *----
       TYPE(C_PTR) IPMAC1
-      INTEGER NC,NG,NL,NMIX1,IPRINT
-      REAL ENER(NG+1),JXM(NMIX1,NG,NC),FHETXM(NMIX1,NG,NL,NC)
+      INTEGER NC,NG,NL,NMIX1,IH,IPRINT
+      REAL ENER(NG+1),HFACT1(NMIX1,NG,NC),JXM(NMIX1,NG,NC),
+     1 FHETXM(NMIX1,NG,NL,NC)
 *----
 *  LOCAL VARIABLES
 *----
@@ -145,6 +149,10 @@
         CALL LCMPUT(KPMAC1,'NJJS00',NMIX1,1,NJJ)
         CALL LCMPUT(KPMAC1,'IJJS00',NMIX1,1,IJJ)
         CALL LCMPUT(KPMAC1,'IPOS00',NMIX1,1,IPOS)
+        IF(IH.EQ.1) THEN
+          WORK(1)=0.5*(HFACT1(IBM,IGR,1)+HFACT1(IBM,IGR,2))
+          CALL LCMPUT(KPMAC1,'H-FACTOR',NMIX1,2,WORK)     
+        ENDIF
       ENDDO
 *----
 *  SCRATCH STORAGE DEALLOCATION
-- 
cgit v1.2.3