diff options
| author | HEBERT Alain <alain.hebert@polymtl.ca> | 2025-12-28 15:55:41 -0500 |
|---|---|---|
| committer | HEBERT Alain <alain.hebert@polymtl.ca> | 2025-12-28 15:55:41 -0500 |
| commit | 754ef58dfd2880f95dd9765d035389f391917492 (patch) | |
| tree | d7056a5fcb559893c91df8d7533fa5fdb03d8480 /Dragon/src/BREKOE.f | |
| parent | ec64ba52445d2d06deba1216471ccf3d289c78a3 (diff) | |
| parent | 744b40856a035580b786378cae13d453edd26689 (diff) | |
Merge branch '19-depreciate-use-of-version-4-and-5-0-draglibs' into 'main'
Resolve "Depreciate use of Version 4 and 5.0 Draglibs"
See merge request dragon/5.1!40
Diffstat (limited to 'Dragon/src/BREKOE.f')
| -rw-r--r-- | Dragon/src/BREKOE.f | 15 |
1 files changed, 11 insertions, 4 deletions
diff --git a/Dragon/src/BREKOE.f b/Dragon/src/BREKOE.f index d98ca1f..855bec2 100644 --- a/Dragon/src/BREKOE.f +++ b/Dragon/src/BREKOE.f @@ -1,6 +1,6 @@ *DECK BREKOE - SUBROUTINE BREKOE(IPMAC1,NC,NG,NL,NMIX1,ISPH,B2,ENER,DC1,TOT1, - 1 SCAT1,JXM,FHETXM,IPRINT) + SUBROUTINE BREKOE(IPMAC1,NC,NG,NL,NMIX1,ISPH,IH,B2,ENER,DC1,TOT1, + 1 SCAT1,HFACT1,JXM,FHETXM,IPRINT) * *----------------------------------------------------------------------- * @@ -24,10 +24,12 @@ * scattering in LAB). * NMIX1 number of mixtures in the nodal calculation. * ISPH SPH flag (=0: use discontinuity factors; =1: use SPH factors). +* IH H-FACTOR flag (=0: not used; =1: recovered). * B2 buckling. * ENER energy limits. * TOT1 total cross sections. * SCAT1 scattering P0 cross sections. +* HFACT1 H-FACTOR values. * JXM left boundary currents. * FHETXM left boundary fluxes. * IPRINT edition flag. @@ -39,9 +41,10 @@ * SUBROUTINE ARGUMENTS *---- TYPE(C_PTR) IPMAC1 - INTEGER NC,NG,NL,NMIX1,ISPH,IPRINT + INTEGER NC,NG,NL,NMIX1,ISPH,IH,IPRINT REAL B2(NC),ENER(NG+1),DC1(NMIX1,NG,NC),TOT1(NMIX1,NG,NL,NC), - 1 SCAT1(NMIX1,NG,NG,NL,NC),JXM(NMIX1,NG,NC),FHETXM(NMIX1,NG,NL,NC) + 1 SCAT1(NMIX1,NG,NG,NL,NC),HFACT1(NMIX1,NG,NC),JXM(NMIX1,NG,NC), + 2 FHETXM(NMIX1,NG,NL,NC) *---- * LOCAL VARIABLES *---- @@ -192,6 +195,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 |