diff options
Diffstat (limited to 'Dragon/src/BRENEM.f')
| -rw-r--r-- | Dragon/src/BRENEM.f | 15 |
1 files changed, 10 insertions, 5 deletions
diff --git a/Dragon/src/BRENEM.f b/Dragon/src/BRENEM.f index 862ea80..defd511 100644 --- a/Dragon/src/BRENEM.f +++ b/Dragon/src/BRENEM.f @@ -1,7 +1,7 @@ *DECK BRENEM - SUBROUTINE BRENEM(IPMAC1,NG,LX1,NMIX1,ITRIAL,IMIX,ICODE,ISPH, - 1 ZKEFF,B2,ENER,VOL1,FLX1,DC1,TOT1,CHI1,SIGF1,SCAT1,JXM,JXP,FHETXM, - 2 FHETXP,ADF1,NGET,ADFREF,IPRINT) + SUBROUTINE BRENEM(IPMAC1,NG,LX1,NMIX1,ITRIAL,IMIX,ICODE,ISPH,IH, + 1 ZKEFF,B2,ENER,VOL1,FLX1,DC1,TOT1,CHI1,SIGF1,SCAT1,HFACT1,JXM, + 2 JXP,FHETXM,FHETXP,ADF1,NGET,ADFREF,IPRINT) * *----------------------------------------------------------------------- * @@ -37,6 +37,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. @@ -59,8 +60,8 @@ 1 IPRINT REAL ZKEFF,B2,ENER(NG+1),VOL1(NMIX1),FLX1(NMIX1,NG),DC1(NMIX1,NG), 1 TOT1(NMIX1,NG),CHI1(NMIX1,NG),SIGF1(NMIX1,NG), - 2 SCAT1(NMIX1,NG,NG),JXM(NMIX1,NG),JXP(NMIX1,NG),FHETXM(NMIX1,NG), - 3 FHETXP(NMIX1,NG),ADF1(NMIX1,NG),ADFREF(NG) + 2 SCAT1(NMIX1,NG,NG),HFACT1(NMIX1,NG),JXM(NMIX1,NG),JXP(NMIX1,NG), + 3 FHETXM(NMIX1,NG),FHETXP(NMIX1,NG),ADF1(NMIX1,NG),ADFREF(NG) *---- * LOCAL VARIABLES *---- @@ -241,6 +242,9 @@ DO JGR=1,NG SCAT1(IBM,IGR,JGR)=SCAT1(IBM,IGR,JGR)/FDXM(IBM,JGR) ENDDO + IF(IH.EQ.1) THEN + HFACT1(IBM,IGR)=HFACT1(IBM,IGR)/FDXM(IBM,IGR) + ENDIF ENDDO ENDDO IF(ICODE(2).NE.0) THEN @@ -347,6 +351,7 @@ 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) CALL LCMPUT(KPMAC1,'H-FACTOR',NMIX1,2,HFACT1(:,IGR)) ENDDO *---- * SCRATCH STORAGE DEALLOCATION |