bad h-factor units in dragon data structure manual

H-factor units MACROLIB are not in J/cm. H-factor are in eV/cm.

NJOY heat units are in ev-barns. LIBTR2 transcript these values in the H-FACTOR vector from the MATXS file.

      IF(HVPS(2:).EQ.'HEAT') THEN
         VECT(:NGRO+1)=0.0
         IF(.NOT.LSUBM1) CALL LCMGET(KPLIB,HVECT(I),VECT)
         DO 261 IK=1,NK
         IF(XS(L5+IK).EQ.0.0) GO TO 261
         JJ=IA(NEX1+IR)+IK-1
         TERPZ=1.0
         IF(.NOT.LSUBM1) TERPZ=TERP(NGRO*(ISUBM-1)+JJ)
         VECT(JJ)=VECT(JJ)+TERPZ*XS(L5+IK)
  261    CONTINUE
         LDEP(1)=.TRUE.
         CALL LCMPUT(KPLIB,'H-FACTOR',NGRO,2,VECT)

Then LIBDEN.f multiplies these values by the isotopes' densities [in cm barns] as follows

         CALL LCMLEN(JPLIB,'H-FACTOR'//NORD(IXSPER),ILONG,ITYLCM)
         IF(ILONG.GT.0) THEN
            LH=.TRUE.
            CALL LCMGET(JPLIB,'H-FACTOR'//NORD(IXSPER),GA1) !eV.barn
            DO 290 LLL=1,NGROUP
            GAF(IBM,LLL,7)=GAF(IBM,LLL,7)+GA1(LLL)*DENISO   !eV/cm
  290       CONTINUE
         ENDIF

So final units are eV/cm and not J/cm as stated in IGE351.pdf.

Ahmed-Naceur / dragon5

bad h-factor units in dragon data structure manual #31