|
This is the issue-tracking page for Vers. 99 of the NJOY Nuclear Data Processing System. Issues are entered by the developer (ryxm@lanl.gov) in response to postings to the NJOY list (njoy@nea.fr), items reported directly by email, observations by NJOY users at Los Alamos, problems noted by the developer, or new ideas from the developer. The summary table includes an issue number, issue status, and issue title. Click on the issue number to jump to a more detailed description of the issue and its resolution, if any. Issue status can be S (submitted), P (pending, meaning there are some comments posted on the problem), or R (resolved, meaning that there is an NJOY update that handles the issue). |
Number |
Status |
Title |
---|---|---|
R |
Incorrect initialization in HEATR |
|
R |
Allow more Legendre terms in HEATR |
|
R |
Vertical segments in File 6 distributions |
|
R |
Incorrect law=61 distributions |
|
R |
More compact charged-particle energy-angle distributions |
|
R |
Possible infinite loop in PLOTR |
|
R |
MT=0 in ACER message |
|
R |
Interpolation for JENDL-3.2 U-235 nubar |
|
R |
Problem with redundant MT=10 |
|
R |
Error in MCNP/sabcol reported for nonthermal problem |
|
R |
Law=5 used in JEF U-233 fission spectrum |
|
R |
Need special skip6 routine in more places |
|
R |
Error in up11 for law=7 and newfor=0 |
|
R |
Missing initialization in acer/ptlegc |
|
R |
Reaction self-shielded cross section not printed correctly |
|
R |
Reich-Moore channel-spin format added |
|
R |
Allow more groups for photons |
|
R |
Improve sigfig truncation for multigroup matrix elements |
|
R |
Incorrect threshold for particle production cross sections in ACER |
|
R |
DTFR doesn't find correct material and temperature on PENDF file |
|
R |
Constant nlmax used before being initialized |
|
R |
Problem in DLAYXS with 69 groups |
|
R |
Unresolved with no resolved for JENDL materials |
|
R |
PDF same as CDF for law-7 charged-particle emission |
|
R |
Problems processing recent EFF-3.0 Be-9 evaluation |
|
S |
Make MT=875-890 available to MCNP |
|
R |
Problems with multigroup runs for 150 MeV data |
|
S |
ACER has problem opening some scratch files under Forte6 |
|
S |
Problem with 181-point tabulated angulur distributions in MF=4 |
|
R |
GROUPR error when running with 300 gamma groups |
|
R |
Failure for Pb-208 from ENDF/B-VI Release 6 |
|
R |
RECONR fails with 100 GeV photoatomic data |
|
R |
Charged-particle heating zero for some parts of Be-9 ACE file. |
|
R |
ACE interpolation-law printout |
|
R |
Excess timers in PURR listing. |
|
R |
Heating zero at first energy point in ACE file. |
|
R |
Allow for larger ACE files |
|
R |
Length of some ACE thermal files too large by one word |
|
R |
ACE parameter LANDH wrong for Be-9 |
|
R |
Problem with ACE law3/33 for isomeric targets |
|
R |
Error message for JENDL photonuclear evaluations |
|
R |
Two-body recoils for photonuclear evaluations |
|
R |
Tabulated fission spectra for photonuclear evaluations |
|
R |
Isotropic discrete recoils |
|
R |
Levels of the (n,2n) reaction in RECONR |
|
R |
Levels of the (n,2n) reaction in HEATR |
|
R |
Histogram interpolation for law=7 |
|
R |
Conversion of law=7 to law=1 |
|
R |
ACER fails for charged-particle set with tabulated angular distribution. |
|
R |
Levels of the (n,2n) reaction in ACER |
|
R |
Diffs in energy grids from BROADR |
|
R |
Bad vertical scale for 3-D plots. |
|
R |
Reactions with nonzero cross section at threshold |
|
R |
Missing threshold for charged-particle reaction |
|
R |
Common variable coeff not set |
|
R |
Photoatomic heating value given in eV-barns instead of MeV-barns |
|
R |
Problem with heating in probability tables |
|
R |
ACER doesn't include delayed neutron data as allowed for MCNP4c |
|
R |
Unitialized variables introduced in up57 |
|
R |
Undefined variables introduced in up63 |
|
R |
Fix radionuclide production for EAF2001 |
|
R |
Bad MT name string in GROUPR |
|
R |
Bad Kalbach a factor for photonuclear cases |
|
R |
Add photonuclear capability for MATXS |
|
R |
Confusion of pi versus pi |
|
R |
Temperature field in RECONR |
|
R |
MCNP Fluorescence Data |
|
R |
Delayed neutron consistency checks and plots |
|
R |
Cross sections or yields that are zero at all energies |
|
R |
8-group delayed neutron data in ACER |
|
R |
8-group delayed neutron data in GROUPR |
|
R |
Editing error in LWR EPRI weight function |
|
R |
Upscatter in getsed |
|
R |
Recent typographical errors |
|
R |
Increased energy resolution in PLOTR output |
|
R |
Plot differences or ratios |
|
R |
Tighten up tolerances for integral thinning in BROADR |
|
R |
Cross section grids in BROADR |
|
R |
Tighten up integral thinning tolerances in RECONR |
|
R |
Cross section grids in RECONR |
|
R |
RECONR printouts |
|
R |
Photoatomic edge discontinuities |
|
R |
Particle production thresholds |
|
R |
More energy points in ACE files |
|
R |
Plots with proton production present |
|
R |
ENDF-6 covariance processing in MODER |
|
R |
Photoatomic heating |
|
R |
Negative percent differences |
|
R |
Unresolved probability tables strange for ENDF66 |
|
R |
Photon production eV-barns |
|
R |
Coherent thermal scattering for aluminum, iron, lead |
|
R |
ACE files for incident charged particles |
|
R |
MATXS self shielding and group bounds |
|
R |
Increase arrays for some evaluations |
|
R |
PURR zeroes two flags |
Number |
|
---|---|
Status |
Resolved |
Submitted |
23 Jul 2000 |
By |
Margarete Mattes (mattes@ike.uni-stuttgart.de) |
Title |
Incorrect initialization in HEATR |
Details |
To process Ni-58 of JEFF-3.0, I had to make a small correction.... |
Comments |
The parameter ir was being initialized incorrectly. |
Impact |
Unknown beyond the case cited. |
Resolution |
Fixed in up15 (31 Jul 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
31 Jul 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Allow more Legendre terms in HEATR |
Details |
Saw a crash in heatr/h6ddx while testing the proposed JEFF-3T evaluation for Fe-56. |
Comments |
It turned out that this new evaluation used more Legendre terms than previous evaluations, so we increased the size of the arrays to handle up to 15 orders. |
Impact |
Only affects Fe-56 from the new JEFF-3T set. |
Resolution |
Fixed in up15 (31 Jul 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
21 Jul 2000 |
By |
Zhang Jingshang (zhangjs@iris.ciae.ac.cn) |
Title |
Vertical segments in File 6 distributions |
Details |
He observed divides by zero in HEATR for Zr-nat from CENDL. |
Comments |
When we looked at the evaluation, we found vertical segments in the distributions for some energies. These vertical segments were causing the divide by zero in the interpolation routine. We fixed the code to find such vertical segments, move the upper energy up slightly, and print a message to the user. Evaluators should be careful not to include such features in their files. |
Impact |
This only affects a few CENDL evaluations currently. |
Resolution |
Fixed in up15 for HEATR and up16 for GROUPR (21 Jul 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
31 Jul 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Incorrect law=61 distributions |
Details |
Originally noticed strange anisotropic distributions at threshold where the scattering cross section is zero. |
Comments |
It turned out that this was caused by an incorrect index in one array, and we originally thought that the patch would only be cosmetic in that it would change the shape of distributions with zero cross section. However, we now realize that it affects all the law=61 distributions. The corrected coding was tested against analytic calculations of the scattering probability at the forward and backward limits. |
Impact |
Affects all distributions that use law=61. This means MCNP4C libraries for such important ENDF/B-VI materials as Fe, Ni, Cr, Cu, and Pb. This is an important patch. |
Resolution |
Fixed in up17 (31 Jul 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
31 Jul 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
More compact charged-particle energy-angle distributions |
Details |
We want to use the compact law=4 format for isotropic charged-particle distributions and the more general law=61 for anisotropic cp distributions. This problem was noticed when processing a proposed evaluation for Fe-56 from JEFF-3. |
Comments |
The problem is that we cannot tell which is which until reading past the first few energies, so the patch is pretty complex! |
Impact |
Basically affects appearance and reduces the size of files. Numerical results when using MCNP should be the same either way. |
Resolution |
Fixed in up17 (31 Jul 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
3 Aug 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Possible infinite loop in PLOTR |
Details |
There is a possibility for an infinite loop. |
Comments |
Simple fix. |
Impact |
Probably not very significant. |
Resolution |
Fixed in up18 (3 Aug 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
3 Aug 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
MT=0 in ACER message |
Details |
The ACER message about energy-dependent yields gives MT=0 for the reaction type. |
Comments |
Just a small patch required. |
Impact |
This is a trivial cosmetic patch and will not affect any results. |
Resolution |
Fixed in up18 (3 Aug 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
3 Aug 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Interpolation for JENDL-3.2 U-235 nubar |
Details |
Problems seen in nubar. |
Comments |
Must set the nonlin flag. |
Impact |
Probably only this one material |
Resolution |
Fixed in up18 (3 Aug 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
3 Aug 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Problem with redundant MT=10 |
Details |
While studying the JEF-2.2 evaluation for lead to find the problem summarized above, we also noted that it used the reaction MT=10. |
Comments |
MT=10 was originally proposed to represent "continuum neutron emission," but it never really became a part of the ENDF scheme and NJOY never has handled it. Unfortunately, it is leading to a double counting of these reactions in the total cross section for ACER and affecting results from HEATR. |
Impact |
So far, it seems that only Pb from JEF is affected. |
Resolution |
Fixed in up18 and up19 (3 Aug 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
11 Aug 2000 |
By |
Waclaw Gudowski (wacek@neutron.kth.se) |
Title |
Error in MCNP/sabcol reported for nonthermal problem |
Details |
MCNP crashes with error in sabcol when using data from ENDF/B-VI Release 6. |
Comments |
There turned out to be two problems here. When processing angular distributions using ltt=3 with newfor=0 (MCNP4B compatibility), there was an extra tosend that caused the code to skip over the first reaction after the elastic MF=4. This lead to a bad tyr=0 value, which somehow causes MCNP to issue a confusing error message about sabcol, even when no S(alpha,beta) data are being used. Second, the code was finding the wrong value for the lct parameter (lab or cm) for sections of File 4 that were fully isotropic. This also can result in an incorrect value for the ACE tyr parameter, which leads to the strange message from sabcol. This problem was introduced when the NJOY97 coding was being converted to block structuring. |
Impact |
The first part of the problem affects runs using ENDF/B-VI Release 6 (which contains the 150-MeV extensions) with newfor=0. It can be avoided by using Release 5, which is identical to Release 6 for MCNP4B purposes. The second part is more serious, and it affects many evaluations for the minor materials where evaluators haven't gone to the trouble to construct realistic angular distributions. |
Resolution |
Fixed in up20 (16 Aug 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
11 Aug 2000 |
By |
Waclaw Gudowski (wacek@neutron.kth.se) |
Title |
Law=5 used in JEF U-233 fission spectrum |
Details |
There are some funny results for U-233.... |
Comments |
The unique feature of U-233 from JEF-2.2 is the use of law=5 to describe part of the secondary energy distribution for fission. This law is normally used only for delayed neutrons, and not for prompt neutron emission. Unfortunately, NJOY does not really process this representation (and never has). The numbers written into the file are just junk. MCNP does have a formal representation for law5 that uses equally probable bins. Equally probable bins are not very good for representing fission; the cummulative law4 is much better. The best solution here is to convert the law=5 part of MF=5/MT=18 to law=1 tabulated format. We are putting an error message in NJOY to alert users to this problem. |
Impact |
Only affects U-233 from JEF-2.2. |
Resolution |
Error message installed in up20 (16 Aug 2000), but that doesn't really solve the problem! |
Number |
|
---|---|
Status |
Resolved |
Submitted |
16 Aug 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Need special skip6 routine in more places |
Details |
The reading of MF6 is sometimes getting out of position. |
Comments |
There are some additional places where we need to use the special version of skip6 called skip6a. For law=7 sections, the tab2 is changed to a tab1 in order to pass on the integrated angular distribution. |
Impact |
This only affects processing of materials that use law=7 with newfor=1 (MCNP4C libraries). For ENDF/B-VI, this is Be-9. |
Resolution |
Fixed in up20 (16 Aug 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
16 Aug 00 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Error in up11 for law=7 and newfor=0 |
Details |
An error was made in the part of up11 that handles law=7. |
Comments |
This problem showed up when processing Be-9 from ENDF/B-VI. |
Impact |
Only affects ACER runs using newfor=0, that is files for MCBP4B, for evaluations that use law=7 for energy-angle distributions. For ENDF, this is only Be-9. |
Resolution |
Fixed in up20 (16 Sep 00). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
16 Sep 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Missing initialization in acer/ptlegc |
Details |
Missing initiazation for dco |
Comments |
Some platforms may have problems for incident charged particles in ACE files if this parameter is not initially set to zero. |
Impact |
Only affects ACE files for incident charged particles (e.g., the APT incident-proton 150-MeV data) on some platforms. |
Resolution |
Fixed in up20 (16 Sep 00). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
28 Sep 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Reaction self-shielded cross section not printed correctly |
Details |
The self-shielded cross sections for the reactions appear to be all the same! |
Comments |
The self-shielded cross sections are not being printed out correctly for the reactions, but the total is OK. The gendf file is OK, so libraries made with NJOY99 are OK. |
Impact |
Only affects the listings, not the data that goes on to production libraries. |
Resolution |
Fixed in up21 (28 Sep 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10 Jun 1999 |
By |
Nancy Larson (nml@ornl.gov) |
Title |
Reich-Moore channel-spin format added |
Details |
There has been a problem in the ENDF Reich-Moore resonance format for a long time for materials with spin greater than zero in that it doesn't specify the widths for different channel spins explicitly. R. Q. Wright (ORNL) suggested using a minus sign on the J value as a flag to use the other channel spin, and this format modification was approved by CSEWG. Nancy Larson provided an NJOY patch to implement this option for NJOY97. John White provided test results from SAMMY to verify the operation of the patch in NJOY. |
Comment |
We made a few changes to Nancy's contribution to make it fit into the block structured coding used in NJOY99. We tested the results against the SAMMY numbers, and Red Cullen and Bob MacFarlane also did tests of NJOY against RECENT running at LLNL. |
Impact |
At present, this change does not impact any existing processing. Several evaluation under development will use it (e.g., Al-27 from ORNL). |
Resolution |
Fixed in up22 (28 Sep 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
28 Sep 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Allow more groups for photons |
Details |
LLNL requested that we allow for up to 400 groups for photoatomic photons and for photon production from particle reactions. |
Comments |
We needed to change some dimensions and sizes in gaminr, dtfr, and groupr. |
Impact |
No impact on existing processing. The change only has an impact if you want more than 150 gamma groups. |
Resolution |
Fixed in up23, up24, and up25 (28 Sep 2000 and 11 Oct 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
11 Oct 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Improve sigfig truncation for multigroup matrix elements |
Details |
There are too many differences in very small matrix elements when going from platform to platform. |
Comments |
The previous algorithm was only operating on the in-group element. We now work on all of them to try to make the results for elastic and discrete inelastic matrices the same on different machines. The basic idea is that these numbers are obtained by subtraction of numbers on the order of unity, so any results less than about 1e-7 are just random numbers and can be removed. |
Impact |
The affect of this change on real calculations is entirely negligible, but it does help to make it easier to do diffs between different platforms. |
Resolution |
Fixed in up25 (11 Oct 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Oct 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Incorrect threshold for particle production cross sections in ACER |
Details |
The current coding sometimes gets the threshold for charged particle production off by one point. |
Comments |
It is necessary to supply a slight shift in the test for the threshold. |
Impact |
Very minor, but it looks better after being fixed. |
Resolution |
Fixed in up26 (12 Oct 2000) |
Number |
|
---|---|
Status |
Resolved |
Submitted |
27 Oct 2000 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
DTFR doesn't find correct material and temperature on PENDF file |
Details |
The PENDF cross sections used for the DTFR plots are not coming from the correct material and temperature. |
Comments |
A do loop was translated incorrectly in converting to block constructs without statement numbers. |
Impact |
Damages the DTFR plots, but not the multigroup tables. |
Resolution |
Fixed in up27 (27 Oct 2000) |
Number |
|
---|---|
Status |
Resolved |
Submitted |
6 Oct 2000 |
By |
Piet de Leege (p.f.a.deleege@IRA.TUDelft.NL) |
Title |
Constant nlmax used before being initialized |
Details |
In up15.13 (heatr) nlmax used before init, perhaps data statement missing. |
Comments |
The statement setting nlmax was inserted inside the "sw" conditional block instead of after both branches of the block. Therefore, it was only available to the short-word version of the code, but Piet was running at 64 bits on his alpha. |
Impact |
This only affects people using native 64-bit modes on their platform. |
Resolution |
Fixed in up30 (5 Nov 2000) |
Number |
|
---|---|
Status |
Resolved |
Submitted |
17 Oct 2000 |
By |
Cornelis Broeders (Cornelis.Broeders@iket.fzk.de) |
Title |
Problem in DLAYXS with 69 groups |
Details |
Found another bug in the routine ccccr.f. For larger group numbers (27 is OK, 69 is wrong) the output file dlayxs is wrong. After some debugging, the correction seems to be ... With this correction, my results are OK now. |
Comments |
His diagnosis was correct. The pointer in the e array for the moreio call was wrong, which caused the later parts of the data block to overwrite the first part, leading to incorrect results. |
Impact |
This only affects the generation of the CCCC delayed-neutron data file DLAYXS. |
Resolution |
Fixed in up29 (5 Nov 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
25 Oct 2000 |
By |
Andrej Trkov (A.Trkov@iaea.org) |
Ttitle |
Unresolved with no resolved for JENDL materials |
Details |
In our work on the WIMS-D library, we identified processing problems in JENDL-3.2 for the following nuclides: 44-Ru-103, 44-Ru-106, 51-Sb-125, 52-Te-127m, 61-Pm-148g, and 61-Pm-149. There may be more. These nuclides have: - Resonance range define entirely by unresolved resonance parameters, - No resolved resonance parameter ranges is defined, - For some nuclides the resonance range starts at relatively high energies (>300 eV for Ru-103). For these materials, the cross sections are zero below the resonance range. Below is a patch.... |
Comments |
We had noticed that here also, but just skipped the materials. The correction proposed by Andrej seems reasonable. |
Impact |
Only affects these JENDL fission-product evaluations (and maybe a few more not listed). There is no affect of this change for ENDF or JEF processing. |
Resolution |
Fixed in up32 (5 Nov 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
2 Nov 2000 |
By |
Chikara Konno (konno@fnshp.tokai.jaeri.go.jp) |
Title |
PDF same as CDF for law-7 charged-particle emission |
Details |
I think that acer.8863 and acer.8864 should be cut since PDF data become the same as CDF data in the angular distribution of the particle production data for law 7.... |
Comments |
The problem was an incorrect pointer for the statement that readjusted the normalization and adjusted the number of significant figures for the PDF in the routine for charged-particle emission in acer. |
Impact |
Only affects ACE files for evaluations using Law 7 for describing the energy-angle emission of secondary charged particles. The one case of this in ENDF/B-VI is alpha emission from the (n,2n)3alpha reaction in Be-9. |
Resolution |
Fixed in up28 (5 Nov 2000). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
15 Nov 2000 |
By |
Ulrich Fischer (ulrich.fischer@irs.fzk.de) |
Title |
Problems processing recent EFF-3.0 Be-9 evaluation |
Details |
We encountered some problems when tring to process the most recent Be-9 EFF-3.0 file (NMOD=3) with ACER/NJOY.14. This file is peculiar in the following sense: (a) In addition to MT=16, it has individual (n,2n) reaction channels stored in MT=875-890 (MF=3, 6, 33). (b) The file is very big (ASCII file size of 17 MB). With the following bug fixes, we succeeded in processing the file.... |
Comments |
Of the changes that Ulrich proposed, some had already been included in up20 in August before he made his report, but he also discovered that the pointer to the moreio call was not being incremented as new pages were brought in, which had the affect of losing some of the words further down in the record and clobbering the first part of the record that was read in earlier. The problem only occurs for very long records like those in the very detailed Be-9 evaluation from EFF-3.0. |
Impact |
This patch will only affect evaluations using Law7 for angle-energy distributions, and only those using very detailed tabulations. The only known example is this Be-9 evaluation from EFF-3.0. |
Resolution |
Fixed in up28 (5 Nov 2000). |
Number |
|
---|---|
Status |
Submitted |
Submitted |
15 Nov 2000 |
By |
Ulrich Fischer (ulrich.fischer@irs.fzk.de) |
Title |
Make MT=875-890 available to MCNP |
Details |
The new Be-9 evaluation for EFF-3.0 uses MT=875-890 to represent different levels of the (n,2n) reaction. Even when considering only MF=3 data using the dosimetry xs option, we cannot persuade ACER to provide them on the ACE file (of course, they are redundant) to the MF=16 data. Furthermore, we would eventually like to use the MF=6 data for MT=875-890 instead of MT=16. |
Number |
|
---|---|
Status |
Resolved |
Submitted |
2 Oct 2000 |
By |
Harm Wienke (hwienke@sckcen.be) |
Title |
Problems with multigroup runs for 150 MeV data |
Details |
Harm first reported error messages from findex when running with 150 MeV data in GROUPR. He suggested a patch to fix that, and then he had some floating-point overflows and suggested a further change. |
Comments |
It turned out that there were two problems. An index was set to 1 in an array used by tab2io, which caused the data being read in to clobber the indexes used by the dynamic storage allocation scheme. This is what caused the findex errors. The other problem was related to the "over" constant, which allows the routine to extrapolate slightly above the highest energy point in the table. It is necessary to make sure that the code doesn't use "over" at the break between the 0-20 MeV data and the 20-150 MeV data that appears when the option LTT=3 is used. Otherwise, GROUPR finds an integration panel of zero width, which leads to the overflow. |
Impact |
This will cause trouble with multigroup processing for all the Los Alamos 150-MeV evaluations that were prepared in support of the APT project. |
Resolution |
Fixed in up31 (5 Nov 2000). |
Number |
|
---|---|
Status |
Submitted |
Submitted |
23 Jan 01 |
By |
Jean-Christophe Sublet (jean-christophe.sublet@CEA.FR) |
Title |
ACER has problem opening some scratch files under Forte6 |
Details |
Using a Forte6 f95 version of NJOY99-24 or 32 for in07 and in14, ACER wants to open unit -13 or -10 through gamsum calling openz "open(nun,form=for,status=age)" but thinks that the unit (nun 13 or 10) is already connected !! |
Comments |
The Sun Forte6 f95 compiler seems to be very strict about opening a scratch file that is already open, although all other compilers used for NJOY so far were more forgiving. We will just have to be careful to close a unit used as a scratch file before reusing the unit for another purpose. |
Impact |
All ACER runs with codes compiled using Sun Forte6 f95. |
Resolution |
Fixed in up45 (8 Apr 01) and up46 (9 Apr 01). |
Number |
|
---|---|
Status |
Submitted |
Submitted |
23 Jan 01 |
By |
Jean-Christophe Sublet (jean-christophe.sublet@CEA.FR) |
Title |
Problem with 181-point tabulated angulur distributions in MF=4 |
Details |
I am trying to process a new pu-239 evaluation for deterministic code (GROUPR road) that does contain some tabulated angular distribution of secondary particles in file 4, mt 2, 51-55 with 260 to 288 energy points, each having a probability distribution f(mu,energy) given for 181 mu values included between -1 and +1. This seems to give no problem to the ACER road processing with an array size increase in ACER: <data nwmaxn/20000/ >data nwmaxn/200000/. However it fails or hangs when GROUPR is used to produces group constants. NJOY97-114 (f90 -g) fails with a core dump and a SIGFPE message. NJOY99-32 (f77 -g), goes into an infinite loop. |
Comments |
|
Impact |
|
Resolution |
|
Number |
|
---|---|
Status |
Pending |
Submitted |
5 Feb 01 |
By |
Brad Sleaford (sleaford1@llnl.gov) |
Title |
GROUPR error when running with 300 gamma groups |
Details |
When running a problem with 300 gamma groups, GROUPR stops with the message "too many groups". |
Comments |
When we increased the number of gamma groups allowed in up25, we forgot to reset the maximum value that is used by the error check. Just change the value 150 to 400 in groupr.2000. |
Impact |
Only affects GROUPR jobs with more than 150 photon groups. |
Resolution |
Fixed in up35 (8 Feb 01) |
Number |
|
---|---|
Status |
Resolved |
Submitted |
8 Feb 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Failure for Pb-208 from ENDF/B-VI Release 6 |
Details |
This failure occurred while making an MCNP4C library for 150-MeV neutron data. |
Comments |
In up17, we checked for isotropic distributions in order to use a more compact presentation. The logic misses this one special case. |
Impact |
Only this one material is currently affected. |
Resolution |
Fixed in up36 (8 Feb 01) |
Number |
|
---|---|
Status |
Resolved |
Submitted |
8 Feb 01 |
By |
Red Cullen (cullen1@llnl.gov) |
Title |
RECONR fails with 100 GeV photoatomic data |
Details |
When trying to run NJOY using the new proposed ENDF version of epdl97, the code goes into an infinite loop in RECONR. |
Comments |
NJOY currently uses 1e10 eV as its idea of an infinite energy. Thus, when the energy in RECONR hit 10 GeV, the code started to fail. We are increasing "infinity" to 1e12 eV in the places that are envolved in photoatomic processing, but leaving 1e10 elsewhere for now. With this change, a sample element from the 100 GeV library runs OK using the pattern of test/in03, and all the other standard test problems run as before. |
Impact |
None with current libraries. This patch will be needed when ENDF/B-VI Release 8 comes out. |
Resolution |
Fixed in up37 (reconr), up38 (njoy), and up39(gaminr), all 8 Feb 01. |
Number |
|
---|---|
Status |
Resolved |
Submitted |
22 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Charged-particle heating zero for some parts of Be-9 ACE file. |
Details |
The charged-particle heating for protons, deuterons, etc., is zero even above the reaction threshold, but sometimes gets to be nonzero at higher energies. Seen in the ACE file for Be-9 from ENDF/B-VI. |
Comments |
Due to a bad if clause, the contribution to heating from charged particles is not being included for MF=6, law 3 or 4. This shows up for reactions like (n,p) or (n,alpha) in Be-9. This if clause was removed, which allows the indentation level to be reduced by one notch. |
Impact |
This only affects heating in coupled neutron-proton calculations in MCNPX (the code doesn't currently handle the other charged particles). In addition, the error is extremely small because of the very small cross sections for emitting charged particles below 20 MeV. |
Resolution |
Fixed in up40 (23 Mar 01). |
Number |
114 |
---|---|
Status |
Resolved |
Submitted |
22 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
ACE interpolation-law printout |
Details |
Some sections on the ACE listing are missing the string "nr = 0" that designates them as using linear-linear interpolation for all points. |
Comments |
Easy to fix. |
Impact |
Only affects the ACE listing, not any of the numerical results. |
Resolution |
Fixed in up41 (27 Mar 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
22 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Excess timers in PURR listing. |
Details |
There are timers given for every ladder in the PURR output, and this make it hard to compare the results of successive runs using diff. There are so many hits that the real differences are hard to see and the diff listing is too long. |
Comments |
This is easy to fix, and we had already made this change in NJOY 2001. Note that when running these kinds of diffs between listings made on different computers, there will always be a large number of differences because the random number sequences will be different. This also will change with NJOY 2001. |
Impact |
Only affects the listing, not any of the numerical results. |
Resolution |
Fixed in up42 (27 Mar 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
22 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Heating zero at first energy point in ACE file. |
Details |
The heating value is zero at first energy point in ACE file, and then it has reasonable values for higher points. |
Comments |
This actually comes from HEATR where the round-up value used to make sure that the energy is inside a table when extracting values was chosen back in the old days when we had to cope with short-word machines. Now that we use double precision for energies and can routinely handle 7 or more digits, the value of "rup" needs to be made closer to unity. People who look closely at PENDF files might notice that the first energy in MF=3/MT=301 is given as 1.000010-5 instead of the correct 1.000000-5. This is why ACER thinks the heating (and damage) is zero at the first energy point. |
Impact |
Very minor for realistic problems. |
Resolution |
Fixed in up43 (27 Mar 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
22 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Allow for larger ACE files |
Details |
LANL/X-5 would like to be able to make ACE files with lengths up to 1.5 Mwords for some special files at lower temperatures. |
Comment |
It is necessary to increase the dimension of the main container array for ACE data, the number giving that size, and the integer field on the XSDIR cards that holds the file length. We hope that this doesn't adversly impact people using smaller machines. |
Impact |
No impact for runs that do not need the extra size. |
Resolution |
Fixed in up44 (29 Mar 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
13 Mar 2001 |
By |
Roberto Orsi (orsi@bologna.enea.it) |
Title |
Length of some ACE thermal files too large by one word |
Details |
The length of the data block for ACE thermal files is too large by one word in both the header of the file and the XSDIR line. Orsi was able to identify the line in NJOY 94, 97, and 99 that caused this problem. |
Comment |
This problem only occurs for materials containing incoherent elastic scattering. For endf, this is poly, h(zrh), and solid methane. In addition, it will only occur if the true length of the ACE data block is divisible by 4. this explains how we were able to get by for so many years without seeing the problem! |
Impact |
Impacts some thermal ACE files for hydrogenous solids only. |
Resolution |
Fixed in up47 (9 Apr 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
30 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
ACE parameter LANDH wrong for Be-9 |
Details |
The LANDH parameter should be zero (not -1) for isotropic subsections of MF=6 described using ACE law=3 in Be-9. |
Comments |
This occurs for the reactions (n,p0) through (n,a0) in Be-9 from ENDF/B-VI. |
Impact |
This error only affects MCNPX. The symptom is a cosine>1.0 message. |
Resolution |
Fixed in up48 (9 Apr 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
30 Mar 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Problem with ACE law3/33 for isomeric targets |
Details |
In response to some problems that Waclaw Gudowski reported with isomeric fission products, Bob Little determined that the original specifications for the ACE-format law3/33 needed to be changed to use -Q instead of abs(Q). |
Comments |
This change allows some upscatter to take place when levels are de-excited in isomeric targets. |
Impact |
This change only impacts those few evaluations that describe isomeric targets. |
Resolution |
Fixed in up49 (9 Apr 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
9 Apr 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Error message for JENDL photonuclear evaluations |
Details |
ACER crashes with an I/O error for most of the JENDL photonuclear evaluations now available from the http://iaeand.iaea.or.at/photonuclear web site. |
Comments |
These evaluations use a non-conforming format where MF=6, MT=201-207 are used to represent particle production. We are adding a clearer error message for the users' convenience. These evaluations cannot be used in NJOY or MCNPX in their current form. |
Impact |
Affects the processing of all the JENDL photonuclear evaluations except for H-2 (which has different problems with NJOY). |
Resolution |
Fixed in up50 (9 Apr 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Apr 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Two-body recoils for photonuclear evaluations |
Details |
ACER doesn't handle photonuclear evaluations that contain sections of File 6 describing two-body recoils (law=4). |
Comments |
We have added coding to handle this option. It has been tested using a modified form of the JENDL evaluation for gamma+D and a preliminary LANL evaluation for gamma+D. Full distributions are provided for both secondary neutron and secondary proton. |
Impact |
Currently only affects evaluations for gamma+D. |
Resolution |
Fixed in up51 (12 Apr 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
13 Apr 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Tabulated fission spectra for photonuclear evaluations |
Details |
ACER doesn't handle photonuclear evaluations that contain sections of File 5 using law=1 (tabulated secondary energy distribution). |
Comments |
Up until now, the code handled analytic fission spectra for photonuclear evaluations, and we have now added coding to handle tabulated fission spectra, as needed for some of the Russian photonuclear evaluations for the actinides, such as the plutonium isotopes. Note that due to other formatting problems, revised versions of the Russian evaluations are needed for processing by NJOY. These revised versions have been submitted to the IAEA photonuclear compendium. |
Impact |
Affects photonuclear processing for many of the Russian evaluations from the BOFOD library. |
Resolution |
Fixed in up52 (13 Apr 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
11 Jun 01 |
By |
Dieter Leichtle (leichtle@irs.fzk.de) |
Title |
Isotropic discrete recoils |
Details |
GROUPR gives the error message “desired energy above highest given” when computing a transfer matrix for the discreate recoil in MT=701 for Be-9 in JEFF. |
Comments |
This occurs if the File 6 distribution is fully isotropic (law=3). The getfle routine doesn't realize that when doing a discrete recoil (law=4). We add a special flag of law=-4 to pass the fact of isotropy into getfle. |
Impact |
Only affects runs that compute a transfer matrix for the recoil particle when the first particle is given as totally isotropic. One example is the recoil from the (n,t) reaction MT=701 for Be-9 from ENDF/B-VI. |
Resolution |
Fixed in up53 (11 Jun 01). |
Number |
125 |
---|---|
Status |
Resolved |
Submitted |
12 Jun 01 |
By |
Dieter Leichtle (leichtle@irs.fzk.de) |
Title |
Levels of the (n,2n) reaction in RECONR |
Details |
RECONR doesn't allow for MT=875-891, which is used to represent different levels in the (n,2n) reaction for the new JEFF Be-9. |
Comments |
The code expects that MT=16 contains the sum of MT=875-891, when present, just as MT=103 contains the sum of MT=600-649. RECONR uses MT=16 in computing the total, and it passes MT=875-891 on to subsequent modules. |
Impact |
Affects only the new JEFF Be-9 evaluation. |
Resolution |
Fixed in up54 (12 Jun 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Jun 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Levels of the (n,2n) reaction in HEATR |
Details |
HEATR doesn't allow for MT=875-891, which is used to represent different levels in the (n,2n) reaction for the new JEFF Be-9. |
Comments |
This patch matches the one described in 125 above. |
Impact |
Affects only the new JEFF Be-9 evaluation. |
Resolution |
Fixed in up55 (12 Jun 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Jun 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Histogram interpolation for law=7 |
Details |
Especially bad particle energies are observed for the discrete neutron in MT=876 for Be-9 from JEFF. |
Comments |
The integration over secondary energy for law=7 in getsix must allow for histogram interpolation. |
Impact |
Affects only the new JEFF Be-9 evaluation. |
Resolution |
Fixed in up55 (12 Jun 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Jun 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Conversion of law=7 to law=1 |
Details |
Dieter Leichtle (FZK) has reported problems processing MT=876 in Be-9 from JEFF, which contains two subsections for emitted neutrons. He proposed a possible patch, and this is another possiblity. |
Comments |
We convert subsections of File 6 using law=7 (angle-energy form) into law=1 (energy-angle form), which can then be converted into law=61 format for the ACE file. This handles multiple subsections for particle emission cleanly, but it is a little more bulky. This is only down when the new formats are requested. |
Impact |
Affects only Be-9 (ENDF/B-VI or JEFF) for MCNP4c new formats. |
Resolution |
Fixed in up56 (12 Jun 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Jun 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
ACE processing fails for charged-particle set with tabulated angular distribution. |
Details |
The code wasn't allowing for tabulated angular distributions with incident charged particles, although Legendre representations were OK. |
Comments |
Fixed bu adding a new block of coding in ACER. |
Impact |
Affects one charged-partcle evaluation so far. |
Resolution |
Fixed in up56 (12 Jun 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12 Jun 01 |
By |
Dieter Leichtle (leichtle@irs.fzk.de) |
Title |
Levels of the (n,2n) reaction in ACER |
Details |
ACER doesn't allow for MT=875-891, which is used to represent different levels in the (n,2n) reaction for the new JEFF Be-9. |
Comments |
This patch corresponds to 125 and 126 above, but for the ACER module. It causes MT=16 to be skipped when the partials are present. In addition, new reaction names are added for the partials (n,2n) cross sections, and the storage array is increased in size to handle the very large Be-9 evaluation from JEFF. |
Impact |
Affects only the new JEFF Be-9 evaluation. |
Resolution |
Fixed in up57 (12 Jun 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10 Jul 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Diffs in energy grids from BROADR |
Details |
Using diff to compare results from NJOY99 and NJOY2001 shows small differences in multigroup results and differences in the ACER energy grids. |
Comments |
The energy points where paging takes place are different in the two versions of NJOY. This causes differences less than tolerances in multigroup results for the groups that contain paging boundaries, but the shift in the grids makes diffs of ACE files pretty useless. We can alleviate this somewhat by increasing the amount of storage available for paging in BROADR. This increases the number of cases where no paging is done, and they should compare well. Some really big resonance evaluations will still show this affect. |
Impact |
This will affect comparisons of jobs run with up58 in place to jobs run previously. The actual effect of the change for real results will be negligible. But the change would be of interest to people doing things like QA comparisons between different systems. Therefore, some people may prefer to leave this patch out. |
Resolution |
Fixed in up58 (10 Jul 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
20 Jul 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
Bad vertical scale for 3-D plots. |
Details |
For charged particle emission, only a tall triange is shown for the first energy and the rest of the distribution is not shown. |
Comments |
Evaluators ofter use narrow tall triangles at the threshold for reactions. If the base of the triangle is shall, the height can be much larger than all the other data in the plot, thus causing them to be below the floor of the plot, even for log scales. We just don't use the first point in computing the z axis limits. The height of the first energy plot is cut off, but that is not important. We can then see the rest of the 3-D plot at a reasonable scale. |
Impact |
Affect the appearance of some 3-D plots. |
Resolution |
Fixed in up59 (20 Jul 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
24 Sep 01 |
By |
Stephanie Frankle (frankles@lanl.gov) |
Title |
Reactions with nonzero cross section at threshold |
Details |
ACE file threshold cross section nonzero for Gd-158 from ENDF/B-VI. |
Comments |
Occasionally, reactions are given with a nonzero cross section at threshold, even though this violates ENDF procedures. RECONR had some logic for handling this, but it was being overwritten by another change. We fix the problem by inserting an extra energy point just above the threshold and zeroing the cross section at the threshold. A diagnostic message aimed at the evaluator is provided. |
Impact |
Affects only a very few materials, one being Gd-158 from ENDF/B-VI. Even there, most results for applications won't be affected. |
Resolution |
Fixed in up60 (24 Sep 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
24 Sep 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Missing threshold for charged-particle reaction |
Details |
ACE file threshold missing for Cd-110 from ENDF/B-VI. |
Comments |
RECONR tries to define a “pseudo threshold” for very small cross sections for incident charged particles, however, this scheme seems to fail, resulting in the omission of the reaction threshold energy. This problem shows up for the (n,n't) reaction for Cd-110 in ENDF/B-VI. We aleviated this affect by reducing the constant “small” that controls this option from ee-15 to 1e-30 barns. A fuller solution would be desireable. |
Impact |
Affects only a very few materials, one being Cd-110 from ENDF/B-VI. Even there, most results for applications won't be affected. |
Resolution |
Fixed in up60 (24 Sep 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
25 Sep 01 |
By |
Kosako Kazuaki (kosako@sae.co.jp) |
Title |
Common variable coeff not set |
Details |
Compilation discovered that the common variable coeff was not being set in subroutine eval. |
Comments |
Actually routine eval is not really used in tabize anymore; it is just a leftover initialization for a(iy). The only other place it is used is in islin2, which is never called anymore. We just clean up the code by removingt these remnants. |
Impact |
No effect on calculations, but it does remove a possibly worrisome compiler message. |
Resolution |
Fixed in up61 (25 Sep 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
25 Sep 01 |
By |
Kosako Kazuaki (kosako@sae.co.jp) |
Title |
Photoatomic heating value given in eV-barns instead of MeV-barns |
Details |
The ACE files for photoatomic data have the wrong units for heating and the value for the atomic weight is not being set. |
Comments |
The current ACER treatment for photoatomic data is not complete. We were hoping that there would be an upgrade in MCNP to provide more detailed methods for handling fluorescence photons, but that has not happened as yet. We did fix the two problems noted by Kosako. People who feel heroic could try patching in some fluorescence data, but they should take care to check that the photon energies produced match the photo edges in the cross section. The key problem is to make sure that the photons appear below the edges where the cross sections are fairly low. This is not true for the old MCNP fluorescence data when used with the more recent cross sections. |
Impact |
People trying to produce photoatomic data for MCNP. |
Resolution |
Fixed in up61 (25 Sep 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
28 Sep 01 |
By |
Bob Little (rcl@lanl.gov) |
Title |
Problem with heating in probability tables |
Details |
The heating values in the ACE probability tables are not what X-5 expected. |
Comments |
We changed the calculation to get results that are consistently given as eV/reaction for lssf=0 and fluctuation factors for lssf=1. |
Impact |
Calculations of self-shielding heating in the unresolved range. |
Resolution |
Fixed in up62 (28 Sep 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
15 Oct 01 |
By |
Bob MacFarlane (ryxm@lanl.gov) |
Title |
ACER doesn't include delayed neutron data as allowed for MCNP4c |
Details |
MCNP4c added the capability to use delayed neutron fractions, time constants, and time-group spectra in both static and time dependent calculations. A library that was not produced using NJOY was provided with the 4c release (ENDF61). |
Comments |
We provided coding to take delayed data from ENDF-format evaluations and add them to the ACE class “c” files. The results were tested for some fast critical assemblies. The changes were minimal for the very hard spectrum devices like GODIVA, but we did see about an 0.2 % decrease in k-eff for the softer BIGTEN. |
Impact |
Affects people who need to go beyond the delayed data provided in ENDF61 with the MCNP4c release. |
Resolution |
Fixed in up63 (15 Oct 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
07 Dec 01 |
By |
Piet deLeege (p.f.a.deleege@iri.tudelft.nl) |
Title |
Unitialized variables introduced in up57 |
Details |
Two variables introduced in up57 are not being properly initialized in the 64-bit version. |
Comments |
The initializations for “zero” and “one” were inserted in the wrong place so that they were only available if “*set sw” was used. We moved them below the upd conditional block to make them available to both 32-bit and 64-bit versions. |
Impact |
Affects ACER for 64-bit versions from 99.57 through 99.63. |
Resolution |
Fixed in up64 (07 Dec 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
07 Dec 01 |
By |
Piet deLeege (p.f.adeleege@iri.tudelft.nl) and Roberto Orsi (orsi@bologna.enea.it) |
Title |
Undefined variables introduced in up63 |
Details |
The compilers call out some unitialized variables in ACER. |
Comments |
In constructing the new delayed neutron capability for ACER, we made some changes to common blocks. Unfortunately, two of these blocks had different names for one of the variables in different routines, and we didn't catch this at the time. In addition, the upper limit for one of the upd delete commands was too large by one, thereby deleting another common block. These three things have now been fixed. |
Impact |
Affects ACER in version 99.63 only, and probably only the block of data describing charged particle emissions. |
Resolution |
Fixed in up64 (07 Dec 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/10/01 |
By |
Jean-Christophe Sublet (CEA-Cadarache) |
Title |
Fix radionuclide production for EAF2001 |
Details |
GROUPR needs some updates to properly handle MF8 |
Comments |
For radionuclide production, NJOY uses the contents of File 8 as a directory of all the productions to be processed. |
Impact |
Affects runs using EAF2001 for nuclide production calculations. |
Resolution |
Fixed in up65 (10 Dec 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
02/12/02 |
By |
MacFarlane |
Title |
Bad MT name string in GROUPR |
Details |
The MT name string printed out by GROUPR is sometimes incorrect. |
Comments |
This occurs when automatic reactions are mixed with manual reactions where the mtname is not given. The code prints out whatever string was left from the previous case. |
Impact |
Only affects the printed output from GROUPR and not any other results. |
Resolution |
Fixed in up66 (12 Feb 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
02/12/02 |
By |
MacFarlane |
Title |
Bad Kalbach a factor for photonuclear cases |
Details |
An error in the calculation of the Kalbach a factor for the photonuclear case results in a values that are so large the sinh(a) overflows with a floating point error. |
Comments |
The energy should have been converted to MeV for this formula. |
Impact |
Affects multigroup calculations of photonuclear energy-angle distributions. |
Resolution |
Fixed in up66 (12 Feb 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
02/12/02 |
By |
MacFarlane |
Title |
Add photonuclear capability for MATXS |
Details |
Photonuclear results from GROUPR are converted to MATXS format, thus making it possible to complete neutron-photon coupled calculations in TRANSX. |
Comments |
The MATXSR input is modified to allow a unit to be specified that contains the multigroup photonuclear cross sections from GROUPR. |
Impact |
This is of interest to people who need to handle photoneutrons in multigroup calculations. |
Resolution |
Fixed in up67 (12 Feb 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
04/18/02 |
By |
MacFarlane |
Title |
Confusion of pi versus pi |
Details |
The name “pi” was being used for the constant pi and for the imaginary part of p in different places, resulting in a floating point error when calculating sqrt(pi) in PURR. |
Comments |
The name for the imaginary part of p was changed to pim. |
Impact |
Important for all users of PURR. |
Resolution |
Fixed in up68 (18 Apr 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
04/18/02 |
By |
Jean-Christophe Sublet (CEA-Cadarache) |
Title |
Temperature field in RECONR |
Details |
Allows for more digits in the temperature printout from RECONR. |
Comments |
This was needed for EAF-2001. The teperature field is not being used for resonance reconstruction, but it is passed to the output file for later use. |
Impact |
Only on people working with EAF-2001. |
Resolution |
Fixed in up68 (18 Apr 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/07/01 |
By |
LANL MCNP Effort |
Title |
MCNP Fluorescence Data |
Details |
Adds the capability to generate fluorescence data for MCNP using the existing Cashwell-Everett format with new numbers obtained from the ENDF versions of Livermore's EADL and EPDL files of atomic data. |
Comments |
The data produced with this method should give reasonable results for transport and heating for energies above 1 keV. The new evluations allow for lower incident photon energies and for more detail in photon and electron distributions from photoabsorption, and future version of NJOY and MCNP can take advantage of this. |
Impact |
Impacts MCP users of photoatomic data. This coding was generated to help with the new libraries for MCNP5. |
Resolution |
Fixed in up69 (07 Dec 01). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
02/20/02 |
By |
MacFarlane |
Title |
Delayed neutron consistency checks and plots |
Details |
Adds consistency check for ACE delayed neutron data, and adds plots for nubar and the the delayed neutron spectra. |
Comments |
This can help to avoid errors in the preparation of delayed neutron files for MCNP. |
Impact |
People who prepare MCNP data libraries. |
Resolution |
Fixed in up70 (20 Feb 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
08/25/02 |
By |
MacFarlane |
Title |
Cross sections or yields that are zero at all energies |
Details |
If a cross section or yield is zero at all energies, the union grid is spoiled. |
Comments |
We fixed lunion so the sol-called pseudo-threshold search cannot go past one less than the number of points in the section. This problems was seen in N-15 of JENDL-3.3 for MF=12/MT=104 and MF=12/MT=105. |
Impact |
Only affects a very few evaluations that should probably be fixed anyway. |
Resolution |
Fixed in up71 (25 Aug 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
09/03/02 |
By |
MacFarlane |
Title |
8-group delayed neutron data in ACER |
Details |
There is delayed neutron data in some materials of JEFF-3 that use 8 time groups rather than the more conventional 6 groups. |
Comments |
The coding was extended to allow more delayed neutron time groups, and some problems were fixed with the use of ENDF law=1 rather than the more common law=5. We also renormalize the distributions because some of the cases in JEFF-3 are not properly normalized. |
Impact |
JEFF-3 processing for MCNP. |
Resolution |
Fixed in up72 (03 Sep 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
09/03/02 |
By |
MacFarlane |
Title |
8-group delayed neutron data in GROUPR |
Details |
There is delayed neutron data in some materials of JEFF-3 that use 8 time groups rather than the more conventional 6 groups. |
Comments |
The coding was extended to allow more delayed neutron time groups. |
Impact |
JEFF-3 processing for MCNP. |
Resolution |
Fixed in up73 (03 Sep 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
09/03/02 |
By |
Skip Kahler (Bechtel Bettis) |
Title |
Editing error in LWR EPRI weight function |
Details |
Some numbers were incorrect in this weight function |
Comments |
Fixed with single and double precision versions. |
Impact |
People using the LER EPRI weight function in GROUPR. |
Resolution |
Fixed in up74 (03 Sep 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
09/03/02 |
By |
Kazuaki (SAE, Japan) |
Title |
Upscatter in getsed |
Details |
Incorrect boolean statement |
Comments |
This effected the removal of upscatter terms in subroutine getsed. |
Impact |
GROUPR users. |
Resolution |
Fixed in up74 (03 Sep 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/09/02 |
By |
MacFarlane |
Title |
Recent typographical errors |
Details |
Some typographical errors were made in recent updates. |
Comments |
Fixed. |
Impact |
ACER users. |
Resolution |
Fixed in up75 (09 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/09/02 |
By |
MacFarlane |
Title |
Increased energy resolution in PLOTR output |
Details |
More energy resolution is needed for looking at the fine details of resonance reconstruction. |
Comments |
Increased to e16.8. |
Impact |
People involved in code comparisons of resonance reconstruction. |
Resolution |
Fixed in up76 (09 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/11/02 |
By |
MacFarlane |
Title |
Plot differences or ratios |
Details |
Plots of cross section percent differences or ratios can be very useful for comparing different evaluations or testing the performance of sections of NJOY. |
Comments |
This capability previous existed for the online plotter, but now it is also implemented for direct NJOY runs. |
Impact |
People interested in the differences between evaluations or the performance of codes for such things as resonance reconstruction or Doppler broadening. |
Resolution |
Fixed in up77 (11 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/10/02 |
By |
MacFarlane |
Title |
Tighten up tolerances for integral thinning in BROADR |
Details |
In making detailed comparisons with other codes (such a PREPRO), differences are seen to to NJOY's integral thinning logic. They show up as sharp spikes up to errmax. These excursions may be a little too large for properly tracking the self-shielding in cross section minima. |
Comments |
We changed errmax to 10*errthn and also tightened up the errint value. |
Impact |
May have some small effect for capture at high energies in U-238, but this has not been quantified. It does help in reducing worries with code intercomparisons a bit. |
Resolution |
Fixed in up78 (10 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/10/02 |
By |
MacFarlane |
Title |
Cross section grids in BROADR |
Details |
In making comparisions with PREPRO, some sharp spikes of difference were observed that exceeded the expected limits. |
Comments |
We try to put the cross sections on either 7 or 9-digit printable energies, but there were some bugs in this procedure that were visible when running at 0.01% accuracty above 1 keV. |
Impact |
Basically no integrated impact, but it does cause alarm when doing comparisons at very hight precision. |
Resolution |
Fixed in up78 (10 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/10/02 |
By |
MacFarlane |
Title |
Tighten up integral thinning tolerances in RECONR |
Details |
Some deviations were higher than comfortable levels. |
Comments |
We changed errmax to 10*errthn and also tightened up the errint value. |
Impact |
Could have some small impact on things like U-238 capture at high energies, but this has not been quantified. |
Resolution |
Fixed in up79 (10 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/10/02 |
By |
MacFarlane |
Title |
Cross section grids in RECONR |
Details |
In making comparisions with PREPRO, some sharp spikes of difference were observed that exceeded the expected limits. |
Comments |
We try to put the cross sections on either 7 or 9-digit printable energies, but there were some bugs in this procedure that were visible when running at 0.01% accuracty above 1 keV. |
Impact |
Basically no integrated impact, but it does cause alarm when doing comparisons at very hight precision. |
Resolution |
Fixed in up79 (10 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
10/10/02 |
By |
MacFarlane |
Title |
RECONR printouts |
Details |
The printout for errors caused by significant figure reduction is obsolete now that we can go to 9 digits when necessary. |
Comments |
We removed the significant-figure reduction print and added a column for the error introduced into the fission integral by integral thinning. |
Impact |
Only effects the print output. |
Resolution |
Fixed in up79 (10 Oct 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/10/02 |
By |
Morgan White (LANL X-5) |
Title |
Photoatomic edge discontinuities |
Details |
The sharp discontinuities at the photoedges should have double energy points. |
Comments |
Fixed as requested. |
Impact |
Photoatomic data in MCNP. |
Resolution |
Fixed in up80 (10 Dec 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/10/02 |
By |
MacFarlane |
Title |
Particle production thresholds |
Details |
The production cross section and heating at the particle production thresholds should be zero. |
Comments |
Fixed. |
Impact |
Productions data in MCNP. |
Resolution |
Fixed in up80 (10 Dec 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/10/02 |
By |
MacFarlane |
Title |
More energy points in ACE files |
Details |
The format overflows with more than 99999 energy points. |
Comments |
Allow 6-digit energy indexes with the first column always blank. |
Impact |
ACE listings from NJOY. |
Resolution |
Fixed in up80 (10 Dec 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/10/02 |
By |
MacFarlane |
Title |
Plots with proton production present |
Details |
A high-energy U-238 case with proton production had bad plots of capture and fission resonance contours. |
Comments |
The nrt value was being overwritten when charged particle production is present. |
Impact |
ACE plots only. |
Resolution |
Fixed in up80 (10 Dec 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/02/02 |
By |
MacFarlane |
Title |
ENDF-6 covariance processing in MODER |
Details |
MODER isn't handling the ENDF-6 covariance formats properly. |
Comments |
Fixed to pass them on properly. |
Impact |
Runs that use MODER with some covariance evaluations. |
Resolution |
Fixed in up81 (02 Dec 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
01/28/03 |
By |
Morgan White (LANL X-5) |
Title |
Photoatomic heating |
Details |
Not what is needed by MCNP5. |
Comments |
Fixed using contributions from White. |
Impact |
MCNP users of photoatomic data. |
Resolution |
Fixed in up82 (28 Jan 03). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
02/12/03 |
By |
MacFarlane |
Title |
Negative percent differences |
Details |
Plotr is not printing out negative percent differences. |
Comments |
Zero percent also confuses the thinning logic, so we add a small amount ot all the percent differences. |
Impact |
People plotting percent differences. |
Resolution |
Fixed in up83 (12 Feb 03). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
09/30/02 |
By |
Bob Little (LANL X-5) |
Title |
Unresolved probability tables strange for ENDF66 |
Details |
Tables are different when starting from 77K or 300K. |
Comments |
The original logic designed to handle termperature correlations as neutrons move from one region ot another at a difference temperature is faulty. We changed things to provide a different set of total cross sections bins for each temperature. We also changed the binning logic and increased the number of samples in each pass to get more events in the sharp interference minima for materials like U-238. This takes more time, but it is possible to compensate by reducing the number of ladders calculated. |
Impact |
Users of PURR for either Monte Carlo probability tables or multigroup self-shielded cross sections. |
Resolution |
Fixed in up84 (30 Sep 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
02/25/02 |
By |
Bob Little (LANL X-5) |
Title |
Photon production eV-barns |
Details |
Would like to have the total photon energy production to adjust MCNP runs to deposit photons locally. |
Comments |
We modified HEATR to allow for MT=442, a local definition for total photon production eV-barns. This can be used by future versions of MCNP to correct for local deposition of photons for people who do not want to do complete neutron-photon-heating coupled problems. |
Impact |
None until MCNP is upgraded. |
Resolution |
Fixed in up85 (25 Feb 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
12/10/02 |
By |
Gary Russell (Manual Lujan Center, LANL) |
Title |
Coherent thermal scattering for aluminum, iron, lead |
Details |
They would like to have evaluations for cold aluminum, iron, and lead to use in simulations of their cold-neutron sources. |
Comments |
This patch provides the lattice information, structure factor calculations, and Bragg edge generation for these materials. Evaluations for cold aluminum and iron were generated and seemed to perform well. |
Impact |
Materials science people who use cold neutrons for their studies. |
Resolution |
Fixed in up86 (10 Dec 02). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
07/02/03 |
By |
MacFarlane |
Title |
ACE files for incident charged particles |
Details |
Problems observed trying to process the preENDF/B-VII light ion files from Gerry Hale (LANL T-16). |
Comments |
Some MT numbers for determining charged particle production were incorrect, and there were also failures for some cases where two identical particles are produced by File 6, such at p+t->alpha+alpha, and there was also a problem with the phase-spce option for the primary particle that showed up for the proton distribution for 3He(p,2p). |
Impact |
Only affects light isotope cases with incident charged particles, and these are not currently handled by either MCNP or MCNPX. |
Resolution |
Fixed in up87 (02 Jul 03). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
07/28/03 |
By |
MacFarlane |
Title |
MATXS self shielding and group bounds |
Details |
Self shielding data from MATXSR is incorrent, and print out of group bounds is wrong. |
Comments |
The update that introduced a photonuclear capability for MATXSR on 12Feb02 accidently removed the assignment of the scratch file used for computing the self-shielding delta cross section values. |
Impact |
MATXS files generated since 12Feb02 at LANL. Outside users didn't have the bad update and should be OK. |
Resolution |
Fixed in up88 (28 Jul 03). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
07/29/03 |
By |
MacFarlane |
Title |
Increase arrays for some evaluations |
Details |
U-233 failed at P5 in GROUPR, and W-182 and 186 had to many gammas. |
Comments |
Evaluations keep getting bigger! |
Impact |
GROUPR runs for these evaluations. |
Resolution |
Fixed in up89 (29 Jul 03). |
Number |
|
---|---|
Status |
Resolved |
Submitted |
08/06/03 |
By |
MacFarlane |
Title |
PURR zeroes two flags |
Details |
PURR providing incorrect self shielding information to GROUPR, thus giving different values for the infinitely dilute cross sections when nsigz is one or greater than one. |
Comments |
The lssf and intunr flags in the section MT=152 produced by PURR were set to zero. |
Impact |
Multigroup values for cases where PURR was run after RECONR. |
Resolution |
Fixed in up90 (06 Aug 03). |