NJOY 99.17 31 July 00 ----------------------------------------------------------------------- This file describes updates to NJOY99. See Readme0 for more general information on NJOY99, installation, test problems, and so on. *ident up1 */ leapr -- 04apr00 */ this fixes errors that destroy the calculation of s(alpha,beta) */ for cold para-hydrogen and for both cold deuteriums. the first */ was introduced in 97.53 while making the constants symbolic, */ and the second was introduced in the cleanup process for njoy99. We found these problems this March and April. LEAPR runs for the cold hydrogen and deuterium cases made with versions earlier than 97.53 should be OK. *ident up2 */ groupr -- 04apr00 */ there is an error in the indexing for the xmas 172-group structure */ that throws all the group bounds off by one. this structure is */ used in europe for thermal reactor calculations. This one is only important for people using this group structure. */ fix a problem introduced with the activation patch of 97.102. the */ nk parameter is only used when doing mf values for activation */ products, and it is not appropriate for fission nubar values. *d groupr.3824 */ don't strip off the upscatter groups for the neutron spectra */ coming from delayed neutron emission (mt=455). */ add a missing save statement and fix an unset variable in anased. */ these problems affect delayed neutron spectra. These changes are important for anybody doing multigroup calculations for fissionable materials. *ident up3 */ acer -- 04apr00 */ the acer consistency checks include an option to readjust */ eprime values that are greater than e, when appropriate. */ there are some problems with the logic, especially for */ cases that use histogram interpolation for the distribution. This change improves the appearance of the outputs, but it doesn't actually have any effect on calculations. MCNP only uses the differential probabilities for linear-linear interpolation, and it doesn't make any difference if they are incorrect at the upper limit for histogram cases. */ when using the old format (mcnp4b and earlier), some angle-energy */ distributions from file 6 are converted into the law 67 format, */ because these earlier versions of mcnp couldn't use all the */ file 6 representations. when converting from the cm to the lab, */ the methods used in subroutine fix6 are a little crude. they get */ confused when cm energies are so small that lab cosines of -1 */ are not reached. this patch tries to fix that in a rough way, */ but evaluations that use the cm frame in file 6 will work best */ if most of the cm energies are greater than e/(awr+1)**2 for */ each incident energy e. for mcnp4c and later, the code can */ sample directly from tabulated cm representations, and the */ approximations of the fix6 routine are avoided. this patch */ is needed for one evaluation from JEFF-3. */ increase the available storage to handle the very large */ mf6/mt16 tabulation in JEFF-3 Be-9. */ increase the space available for discontinuities in convr */ to allow for JENDL-3.2 si-nat These three sets of changes are only needed for very specific cases. We already know that there will be a few more like this coming up! *ident up4 */ reconr -- 05apr00 */ be sure to count subsections of file 12 before allocating */ storage for the elements of the new directory. otherwise, */ some materials with many sections of file 12 will overflow. */ this is a longstanding problem that we never noticed before. This is another problem that only showed up for one evaluation. It has been this way in NJOY for many years. *ident up5 */ purr -- 7may00 */ fix a problem introduced while installing the heating part */ of the probability tables. it shows up when doing elements */ that have unresolved data. also, increase the number of */ resonance sections allowed to handle the very large cd-nat */ evaluation from JENDL. PURR marks regions where the unresolved range for one isotope of an element overlaps with the resolved or smooth range from other isotopes by putting a negative sign on the energy. We have to be sure to recover the energy correctly by taking the absolute value! Note that this problem only shows up for elemental evaluations. There are only two of these in ENDF, but many more in JENDL. We also had to allow for more sections in the resonance range to handle JENDL's Cd-nat evaluation. *ident up6 */ acer -- 30may00 */ fix a typo in up3 (reported by bunde, anl) */ acer -- 30may00 */ fix problems with converting cm distributions to law=7 */ and problems reading law=7 into the ace file. these problems */ show up when running newfor=0 with njoy2000, especially on */ some materials from jef-2.2. This one is pretty important. NJOY99 was pretty thoroughly tested for producing MCNP4C libraries, but we neglected to do as much checking for MCNP4B cases (newfor=0). *ident up7 */ viewr -- 30may00 */ increase the allowed size of 3d plots. */ pushed by pb-nat from jef-2.2. In most cases, NJOY just truncates 3-D plots that are too big for viewr, but in the pb-nat case, it actually overwrote some memory. The storage available for 3-D plots has now been increased, and some plots that were truncated before without causing the code to fail will plot better now. *ident up8 */ acer -- 3jun00 */ subroutine ptleg2 does not need the dynamic array xat. */ this problem was first noted by Waclaw Gudowski for ENDF/B-VI */ si-nat. it shows up as "id xat not defined". This is a fairly isolated problem. It is a case of evolution leaving behind an obsolete organ. *ident up9 */ acer -- 07jun00 */ add the capability for processing anisotropic charged particle */ emission using tabulated legendre coefficients into the */ mcnp4c law61 format. this is needed for jeff-3 cr-52. */ allow for multiple interpolation ranges in file 6. this */ also occurs for jeff-3 cr-52. currently, the neutron */ energy-angle distribution only allows for combinations of */ histogram and linear linear interpolation, but the */ charged-particle sections allow for general combinations of */ all allowed interpolation laws. There are only one or two such cases in current evaluations, but more will probably appear in the future. The limitation to interpolation laws 1 and 2 for neutrons comes from the occasional need to handle multiple subsections in File 6 for neutron emission. The increased generality for charged-particle subsections comes about because multiple subsections are not handled. Someone will probably do that someday.... Although we did work fairly hard to handle the JEFF Cr-52 evaluation, the actual representation is fairly strange, and the treatment of the particle yields near the effective emission threshold should really should be changed. *ident up10 */ leapr -- 13jun00 */ fix two incorrect constants in leapr. one affects cases with */ diffusive effects, and it has been incorrect since njoy97.0 */ (oct 97). the other affects cold hydrogen calculations, and it */ has been incorrect since njoy97.53 (dec98). NJOY99 was exhaustively checked using ENDF/B evaluations and the newfor=1 option (to support MCNP4C). Unfortunately, we didn't test it with newfor=0 (for MCNP4B and earlier versions)! Several users have reported problems (especially Waclow Gudowski in Sweden. Some changes that affect running with newfor=0 were made above. The following changes (up11) are especially noticable with ENDF/B-VI Release 6. *ident up11 */ acer -- 26jun00 */ fix an error in determinining which reactions have to be */ converted into law=7 format when using newfor=0. because of */ overzealous code cleanup, acer is trying to convert sections */ with the kalbach representation in addition to sections with */ tabulated angular distributions. */ acer -- 26jun00 */ during the cleanup of the topfil routine, the logic to process */ sections of file 6 using law=2 (two-body distributions) into */ equally probable bins for newfor=0 was omitted. this shows up */ for evaluations that use mf6/mt51, etc., to represent inelastic */ levels. */ acer -- 26jun00 */ fix an error in the reaction naming. it affects mt=44 (n,n2p) and */ mt=45 (n,npa). this problem was introduced in may of 1995. */ examples of cases that use these reactions are later releases of */ endf/b-vi al-27. */ acer -- 26jun00 */ add missing external statement. reported by bokyun seo (kaeri) */ acer -- 26jun00 */ add missing line for the sequential (n,2n) reactions for be-9. */ this line was accidentally removed in njoy 94.19 (jan96). the */ error was continued through njoy 97 and 99. discovered by */ bob little (lanl). */ acer -- 26jun00 */ fix anisotropic photon production (e.g., endf c,n,o) *ident up12 */ njoy -- 27jul00 */ fix two typographical errors in the 64-bit version of the */ slatec math library. reported by piet de leege (delft). *ident up14 */ acer -- 20jul00 */ acer fails if you run it on a pendf tape that only has the */ single reaction mt=2 (elastic). this can happen for he-4 if you */ don't run heatr, thermr, or gaspr first. found by gudowski. */ acer -- 20jul00 */ acer fails for mf=6, law=2, lang>0 (angular distribution with */ tabulated cosines). the only known example is mt=51 for pb-208 from */ release 6 of endf/b-vi. found by waclaw gudowski. */ acer -- 21jul00 */ an error was included in up9, which added a capability to handle */ anisotropic charged-particle emission represented using legendre */ polynomials. the update disabled the case of isotropic */ charged-particle emission, which occurs in a number of important */ materials from release 6 of endf/b-vi. the symptom is a serious */ clobbering of the ace file, such that it cannot even be read into */ mcnp or even printed using acer. also reported by gudowski. *ident up15 */ heatr -- 31jul00 */ incorrect initial value found by m.mattes (u.stuttgart). */ increase the allowed number of legendre terms in h6ddx */ to handle the new jeff-3t fe-56 evaluation. */ watch for ill-defined vertical segments in distributions. these */ have been seen in zr90 from cendl3 and fe56 from jeff3. actually, */ the evaluations should be fixed to avoid such features, because */ we don't really know what y value to select in the vertical */ segment. we choose to just move the second energy of the double */ point up a little. we only print the diagnostic once to keep the */ output cleaner, but there could be more than one vertical segment. *ident up16 */ groupr -- 31jul00 */ watch for ill-defined vertical segments in distributions. these */ have been seen in zr90 from cendl3 and fe56 from jeff3. actually, */ the evaluations should be fixed to avoid such features, because */ we don't really know what y value to select in the vertical */ segment. we choose to just move the second energy of the double */ point up a little. we only print the diagnostic once to keep the */ output cleaner, but there could be more than one vertical segment. *ident up17 */ acer -- 31jul00 */ fix an incorrect index in the law=61 section for the primary */ particle. the effect of this error is to give an incorrect */ angular distribution for the energy points with scattering */ probability zero (which should be isotropic). this change */ is strictly cosmetic and shouldn't affect any results. */ acer -- 31jul00 */ we want to use the compact law=4 for isotropic charged particle */ distributions, and the more general law=61 for anisotropic cp */ distributions. unfortunately, we can't tell which is which */ without reading past the first few energies for some evaluations. */ the need for this patch was first noticed by jeff3 fe56.