Skip to content

Bug in LAW 3 Level Inelastic Scatter #22

Open

Description

@HunterBelanger

Currently, this documentation states that for level inelastic scattering, LDAT(1) = (1+A)*abs(Q)/A, implying that LDAT(1) should always be positive. This is typically true as Q values are generally negative, but this is not the case when looking at excited states, which may have positive Q values. One such example is Am244m1. Looking at line 10361 in acefc.f90, NJOY actually writes LDAT(1) = (1+A)*(-Q)/A. This is an important distinction for error checking in ACE files, as what NJOY actually writes can indeed be negative for excited states.

I have a fix prepared for this, but I am waiting to submit the PR, as there are several things that are waiting to be merged, and the PDF will is likely causing enough merge conflicts already.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions