Fix incident neutron energy subtraction

openmc/data/fission_energy.py

@@ -603,12 +603,12 @@ class FissionEnergyRelease(object):
                                  data=self.prompt_neutrons.coef)
 
             q_prompt = (self.fragments + self.prompt_neutrons +
-                        self.prompt_photons + Polynomial((-1.0, 0.0)))
+                        self.prompt_photons + Polynomial((0.0, -1.0)))
             group.create_dataset('q_prompt', data=q_prompt.coef)
             q_recoverable = (self.fragments + self.prompt_neutrons +
                              self.delayed_neutrons + self.prompt_photons +
                              self.delayed_photons + self.betas +
-                             Polynomial((-1.0, 0.0)))
+                             Polynomial((0.0, -1.0)))
             group.create_dataset('q_recoverable', data=q_recoverable.coef)
         elif self.form == 'Sher-Beck':
             group.attrs['format'] = np.string_('Sher-Beck')

src/endf.F90

@@ -60,6 +60,10 @@ contains
       string = "events"
     case (SCORE_INVERSE_VELOCITY)
       string = "inverse-velocity"
+    case (SCORE_FISS_Q_PROMPT)
+      string = "fission-q-prompt"
+    case (SCORE_FISS_Q_RECOV)
+      string = "fission-q-recoverable"
 
     ! Normal ENDF-based reactions
     case (TOTAL_XS)

tests/test_tallies/results_true.dat

@@ -1 +1 @@
-f8ad60a94994e6b126b64b8a3e10ac3647314fff37558f9e32f174e6baf944aa2e617ad15f3aa6d6b3fc4e4ead944afb3ca2ea606b40e50bf96e7ad86c86a12b
+a6e5480c66e6510687bf281983b3f387da45005bb08d8cd0aa629e53c5a42a1b2fa8d76345ad2df497d85f2b273fbc5edc36105a71a73c1107479ca0af8329f6