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/BRERT.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/BRERT.f')
| -rw-r--r-- | Dragon/src/BRERT.f | 19 |
1 files changed, 13 insertions, 6 deletions
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 |