search

Ahmed-Naceur / dragon5

Three-Dimensional Discrete Ordinates Neutron, Photon, Electron and Positron Transport Code
7 stars 1 forks source link

error in c-factor treatement in LIB: #34

Open Ahmed-Naceur opened 2 years ago

Ahmed-Naceur commented 2 years ago

In LIBDEN.f, there is an error while transforming microscopic charge deposition cross sections to macroscopic ones.

         IF(ILONG.GT.0) THEN
            LC=.TRUE.
            CALL LCMGET(JPLIB,'C-FACTOR'//NORD(IXSPER),GA1) !e-barns
            DO 295 LLL=1,NGROUP
            GAF(IBM,LLL,13)=GAF(IBM,LLL,13)+GA1(LLL)*DENISO !e/cm
  295       CONTINUE
         ENDIF

GAF(IBM,LLL,13) has not been adequately initiated to zero. It contains negative values. This is now corrected as:

         DO 255 IP=1,14 
         DO 250 LLL=1,NGROUP
         GAF(IBM,LLL,IP)=0.0
  250    CONTINUE
  255    CONTINUE