summaryrefslogtreecommitdiff
path: root/Dragon/src/BRELLB.f
diff options
context:
space:
mode:
authorHEBERT Alain <alain.hebert@polymtl.ca>2025-12-28 15:55:41 -0500
committerHEBERT Alain <alain.hebert@polymtl.ca>2025-12-28 15:55:41 -0500
commit744b40856a035580b786378cae13d453edd26689 (patch)
treed7056a5fcb559893c91df8d7533fa5fdb03d8480 /Dragon/src/BRELLB.f
parentec64ba52445d2d06deba1216471ccf3d289c78a3 (diff)
Resolve "Depreciate use of Version 4 and 5.0 Draglibs"
Diffstat (limited to 'Dragon/src/BRELLB.f')
-rw-r--r--Dragon/src/BRELLB.f14
1 files changed, 11 insertions, 3 deletions
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