NJOY2012.82 19 January 2017 ----------------------------------------------------------------------- This file describes updates to NJOY2012. See Readme0 for more general information on NJOY2012, installation, test problems, and so on. *ident up1 */ 21mar2013 */ - acefc */ - coding in subroutine first may read a TAB1 function without */ checking if the allocated space is sufficient. revise to */ make sure the minimum necessary space is allocated. *ident up2 */ 06may2013 */ - groupr */ - coding in njoy99.379 to replace use of an expansion for the */ complementary error function with the more accurate SLATEC */ "erfc" routine was implemented incorrectly. */ - we fix that here, plus delete some no longer needed code. */ - thanks to NRG/Petten (Hogenbirk) and LLNL (Mattoon) for */ uncovering this error. *ident up3 */ 24may2013 */ - mixr */ - need to make a number of variables global. */ - original deficiency reported by Studsvik/Idaho (Wemple). */ - while we're at it, change from hardwired dimension and do loop */ limits to parameter statements. */ - allows for easy future increase on input limits. */ - fix overloading of hollerith and real variables so User file 1 */ comment appears properly in the output file. *ident up4 */ 7jun2013 */ - groupr */ - need to make some local variables global for the flux calculator; */ - need to make some fixed arrays allocatable (fixed size was ok */ when nsigz was restricted to .le. 10). */ - these deficiencies reported by Studsvik/Idaho (Wemple). */ - need global initialization test for panel; *ident up5 */ 7aug2013 */ - groupr */ - make unr deallocation conditional (AWE, Cornock) */ - resolve flux calculator issues (Studsvik, Wemple) */ - include additional test in getfwt when searching for e on ninwt */ after 200 continue to overcome possible near duplicate energies */ at the top of the energy mesh; */ - if test revision in up4.68 needs to occur for line groupr.3883; */ - move coding to update el, wl and xl within the 210 continue loop */ to assure correct values are passed to terp1 after 230 continue. *ident up6 */ 8aug2013 */ - purr */ - correct unit number variable used in tab1io and moreio calls */ for urr energy-dependent scattering radius data (Trkov, IJS). *ident up7 */ 8aug2013 */ - samm */ - need to make a test for possible energy-dependent scattering */ radius more specific (Trkov, IJS). *ident up8 */ 8aug2013 */ - groupr */ - make sure we write the last multigroup record to ngout2 even */ when its data are zero, as has always been done in the past */ (issue reported by Martin, Sandia). */ - an improvement in moving from njoy99 to njoy2012 was a */ little too zealous is skipping over groups with zero data. *ident up9 */ 22aug2013 */ - groupr */ - provide more robust testing for initial and final multigroup number */ in getdis. previous coding could exceed the maximum multigroup and */ produce spurious results (issue reported by Sublet, CCFE). *ident up10 */ 3sept2013 */ - gaspr */ - need better checking for allocatable arrays when processing */ files with multiple temperatures (Vanhanen, Aalto U) *ident up11 */ 4sept2013 */ - groupr */ - need more storage space for tmp in getsed (Vanhanen, Aalto U). */ - new algorithm starts with a large allocation request and */ tests the allocation status, reducing the array size until */ the request is successful. */ - also tweak identical error message text to make them unique. *ident up12 */ 5sept2013 */ - covr */ - need to make sure the array we deallocate was really allocated. */ (Vanhanen, Aalto U). *ident up13 */ 19sept2013 */ - errorr */ - increase nwds in covout to match njoy99. */ - increase scr space in colaps to always be at least npage+50. */ - fixes unexplained stop reported by Vanhanen, Aalto U. *ident up14 */ 19sept2013 */ - purr */ - increase nsamp to match njoy99. */ - make sure 1 .ge. ibin .ge. nsamp in unrest *ident up15 */ 26sept2013 */ - samm */ - there is some confusion in the coding between the maximum */ Legendre order and the number of Legendre orders. This */ change also fixes the message about the maximum Legendre */ order from reconr (rem). *ident up16 */ 26sept2013 */ - reconr */ - delete coding that closes unit nscrl in rexsx. this unit holds */ elastic scattering (mf4/mt2) data that needs to be accessed if */ the logical variable, "Want_Angular_Dist" is changed to .true. */ - still need to close this unit though ... add coding in reconr */ to do so (ack). *ident up17 */ 01oct2013 */ - endf */ - allow more flexibility in findf when searching a gendf tape. *ident up18 */ 02oct2013 */ - errorr */ - if the user groupr structure extends beyond the range of the */ weight function, then the grpav routine may return with "nan" */ for some cross sections. add coding to (i) warn the user of */ the group structure/weight function mis-match, and (ii) set */ those multigroup cross sections to zero. */ - fixes "nan" issue reported by Vanhanen, Aalto U. */ - indicies used in grpav's scr aren't quite right when the number */ of union groups exceeds npage, we fix them now (ack). *ident up19 */ 03oct2013 */ - util */ - delete a couple of redundant lines in closz *ident up20 */ 08oct2013 */ - acer */ - the apparent number of photons can differ for jobs using binary */ tapes versus ascii tapes (Trumbull, KAPL). */ - loss of significant digits when working with ascii tapes makes */ this tricky, particularly when trying to extract the original */ primary photon energy from the sum of primary photon plus */ incident particle energy. This patch works for 1h which has */ a single primary photon over the 20 MeV incident neutron energy */ range while also seeing the two distinct secondary photons of */ 2.112370 MeV and 2.11238 MeV in 50V mt16. */ - NOTE: we strongly recommend use of binary tapes in all njoy */ jobs to maximize the number of significant digits that are */ saved when passing pendf data among the various modules. *ident up21 */ 09oct2013 */ - samm */ - make sure lllmax is defined. needed for array allocation if the */ "Want_Angular_Dist" logical variable is changed from its default */ false value to true. we start with a value believed to always be */ large enough, but after the real, problem dependent, value is */ determined also check that the initial allocation was sufficient. *ident up22 */ 21oct2013 */ - purr */ - need to make sure "init" in unrest is reset to zero when */ executing purr multiple times in the same job, or when */ processing multiple materials (issue reported by Koscheev, */ IPPE). */ - need additional tests to make sure we deallocate all previously */ defined allocatable arrays prior to exiting or prior to reallo- */ cating them with the current nbin, ntemp and nsig0 values. */ - rewind the scratch tape before processing each matd. */ - prior to this last patch, if multiple matd's are processed in */ the same purr job and the ntemp, nsigz and nbin values are the */ same for the first and later matd's it was the first matd's */ data that were retrieved from the scratch tape. *ident up23 */ 24oct2013 */ - samm */ - revise per updates in the ORNL SAMRML code. */ - impacts the RML resonance calculations for */ endf/b-vii.1 35Cl (rem). *ident up24 */ 18mar2014 */ - acedo */ - fix error message if desired temperature is not found on */ the input tape (Trkov, IAEA); */ - need additional logic to update the xss(lsig...) array if */ multiple sub-sections exist for a specific mf=10/mt (ack). */ - previously the code could crash or, depending upon compiler */ options, run silently to completion. */ - NOTE: If code had run to completion, the resulting ACE */ file will correctly include all sub-sections but */ the standard output will only include the first */ sub-section. */ - increase nmax (the maximum number of mt values that can be */ processed) from 100 to 350. needed when processing selected */ eaf-2010 files. */ - also check if nmax is exceeded; current code runs silently */ to completion but the .##y file contains invalid pointers */ causing an immediate mcnp crash. *ident up25 */ 07apr2014 */ - errorr */ - we failed to allocate egn for all possible multigroup options */ (reported by Becker, KIT) *ident up26 */ 07apr2014 */ - groupr */ - cm2lab is running out of space (Orsi, ENEA) */ - an if test in f6cm can be true due to round-off, */ causing epnext to be set to zero and leading to */ an infinite loop in cm2lab. Make the test less */ sensitive. *ident up27 */ 08apr2014 */ - purr */ - a simple typo when reading urr data from endf/b-vi 233u */ (the only Adler-Adler format evaluation left, we believe). *ident up28 */ 15apr2014 */ - heatr */ - need a patch in h6cm similar to what we did in groupr's */ f6cm (ack). *ident up29 */ 29apr2014 */ - groupr */ - we increased tmp array space in up11 but still need more */ (for jeff-3.2 238u, for example). */ - a difference between njoy2012 and njoy99 is that njoy2012, */ by default (ismooth=1), will expand the low energy portion */ of emission spectra defined in histogram format. we insert */ additional energy points to allow the spectrum to more closely */ follow a sqrt(E) profile. when processing delayed neutron */ spectra (mt=455 only, spectra from other mt values are ok) */ we failed to update a pointer needed to properly define the */ starting location in the tmp array for the now larger TAB1 */ function. */ - both issues reported by wemple/studsvik. */ - we also include additional array bounds tests to assure the */ new TAB1 function fits within the allocated tmp space (ack). *ident up30 */ 1may2014 */ - groupr */ - when dealing with astrophysics (1.e9 degK) temperatures */ the target nuclide za value was set to zero on the output */ gendf tape. also at these highly elevated temperatures, */ unresolved resonance data, if present, was ignored. */ (Conlin, LANL) *ident up31 */ 1may2014 */ - acefc */ - need more scratch space in acelcp (ack). *ident up32 */ 2may2014 */ - groupr */ - the "extending lin-lin as sqrt(E) ..." algorithm in getmf6 */ can produce a zero only emission spectrum under the rare */ condition of an original spectrum defined at two energy */ points, where there is a non-zero probability at zero */ emission energy and zero probability at the maximum allowed */ non-zero energy. subsequently an infinite loop condition */ may occur. */ - we expand the existing "patching low-energy distribution ..." */ algorithm to avoid this issue. (Orsi, ENEA) *ident up33 */ 27jun2014 */ - thermr */ - reduce the delta-temperature tolerance when comparing the user */ requested temperature versus the temperature being read from */ the input tape; allows for creating S(a,b) files on a finer */ temperature grid (ack). */ - a similar change follows in up34 for acefc and aceth. */ ***************************************************************** */ *** NOTE: This tighter tolerance means that test01 with its *** */ *** 300 degK input temperature will fail. We change it *** */ *** to 296 degK to better match the temperature found *** */ *** on the S(a,b) input tape. This will cause changes *** */ *** in the test01 output files. *** */ ***************************************************************** *ident up34 */ 02jul2014 */ - acefc & aceth */ - reduce the delta-temperature tolerance when comparing the user */ requested temperature versus the temperature being read from */ the input tape; allows for creating S(a,b) files on a finer */ temperature grid. */ - these changes made in coordination with up33 (ack). */ - aceth */ - current coding to get nxs(1) in continuous formatted files */ is off by one. a comparison with 99.113, when continuous kernel */ file capability was first introduced, shows a couple of code */ differences which we eliminate now. *ident up35 */ 05sept2014 */ - unresr */ - eliminate a fixed size array allocation with a calculation */ of the true needed space (issue reported by Orsi, ENEA) *ident up36 */ 05sept2014 */ - groupr */ - add missing allocation statements in gengpn (issue reported */ by Wemple, Studsvik). */ - need more space in getgyl for some JEFF-3.2 evaluations */ (issue reported by Orsi, ENEA). *ident up37 */ 05sept2014 */ - resxsr */ - fix a couple of njoy99-to-njoy2012 conversion issues */ (reported by Wemple, Studsvik). *ident up38 */ 18sept2014 */ - groupr */ - revise argument list passed in displa to be consistent with */ the declared array dimensions (resolves halt in g95 if compiled */ with array bounds checking turned on); */ - NOTE: problems run with an executable produced by the Intel */ Fortran compiler yield the same matricies with or */ without this patch, as the original array allocation */ is correct for the user's problem size. */ - expand the logic in getgyl that checks for constant production */ over low energies to account for histogram interpolation. */ - include logic in getgyl that sets econst to zero when dealing */ with file 13 photon data. *ident up39 */ 18sept2014 */ - acefc */ - coding introduced in up1 to determine the necessary allocation */ for scr when reading file 6 tab1 records is incomplete. this */ modification assures that the largest tab1 size is identified. */ - this issue identified while processing a candidate 235u cielo */ file from ornl (Leal). *ident up40 */ 14oct2014 */ - reconr */ - correct a message text typo (Koning, Petten). */ - successful execution of the cross section calculation algorithm */ in csunr2 can depend upon failing a "do while" test where we */ expect two real numbers to be equal. due to finite computer */ precision this test can remain true when the comparison is */ false. we revise the test to eliminate this very rare */ false positive condition (Wemple, Studsvik). *ident up41 */ 04dec2014 */ - acefc */ - previous coding assumes if a mf3/mt5 section exists, then light */ particle emission spectra will be found in mf6/mt5. */ - if this is not true, pointers in the ace file will be */ incorrect and a subsequent viewr job to plot the ace data */ can enter an infinite write loop (reported by Koning, Petten). */ - we provide coding revisions that define light particle flags */ while processing mt=5 to overcome this deficiency. */ - include an extra contio to up1 coding to read the end-of-section */ record when checking for the laregest tab1 function in mf6. *ident up42 */ 05jan2015 */ - acedo */ - mcnp's ace file specification for dosimetry files with multiple */ interpolation ranges says to write the NR NBT values followed */ by the NR INT values but we have always written the NR NBT,INT */ values as pairs. issue identified by Milocco (JET,UK), */ Trkov (IAEA) & Brown (LANL). */ **************************************************************** */ * NOTE: This error has existed in ALL past versions of NJOY * */ **************************************************************** */ * NOTE: All dosimetry files found in the production mcnp_data * */ * 531dos, 532dos and llldos files only contain a single * */ * interpolation interval and therefore are UNAFFECTED * */ * by this error. Also user generated dosimetry files * */ * whose endf source data only contain a single interpo- * */ * lation interval are UNAFFECTED. * */ **************************************************************** *ident up43 */ 08jan2015 - 10feb2015 */ - heatr */ - the variables used when calculating the photon energy "kick" */ are computed in multiple locations. we redefine these as */ global variables and only calculate them once (in nheat). */ - hgam102, introduced to calculated heating when the 1h capture */ data were moved from mf12 to mf6, only included the recoil */ energy. we should have included the photon "kick" to get */ more accurate heating and damage values, particularly at low */ energy, and do so now (issue reported by Konno, JAEA). */ ******************************************************** */ *** NOTE: This error has existed since endf/b-vii.0 *** */ *** files were released, but ONLY for 1h. *** */ ******************************************************** *ident up44 */ 12jan2015 */ - samm */ - reposition an integer declaration statement (mostly a cosmetic */ change as we typically declare integers prior to reals but had */ not done so in this instance; and especially when some of those */ integers appear in subsequent real declaration statements). */ - keeps the Visual Digital Fortran 6.6A compiler happy (per */ Cabellos, NEA). *ident up45 */ 12jan2015 */ - acefc */ - even though we now allocate many arrays there are still fixed */ array dimensions that become too small over time. we increase */ a couple of them here to support efforts related to the CIELO */ Project and the IAEA pfns crp (Trkov, IAEA). *ident up46 */ 27jan2015 */ - groupr */ - tweak array definition coding in getmf6 so usage is consistent */ with similar coding elsewhere in this routine. */ - revise array specification for p(7) in f6cm to be consistent */ with njoy99. depending upon compiler options, users specifying */ lord>6 in their input may have suffered a silent array overflow */ (both issues reported by Bunde, DOE/Idaho). *ident up47 */ 29jan2015 */ - heatr */ - we recently noticed that some evaluated nuclear data files */ contain photon production data for a given mt in both mf6 */ and one or more of mf12 through mf15. this causes njoy to */ incorrectly combine the "direct" and "energy-balance" methods */ when calculating heating. we set a flag when reading mf6 */ photon production data so that the corresponding mt in one or */ more of mf12 through mf15 is not processed. */ - an information message identifying the mf/mt section is */ included in njoy's standard output. *ident up48 */ 5feb2015 */ - acefc */ - in unionx we open scratch units iold and inew, and leave them */ open for subsequent use in gamsum (where they are closed). if */ there are no photon data, gamsum is not called and these units */ remain open. if a later module tries to open these unit numbers */ we crash. */ - add conditional coding in unionx to make sure these units */ are closed (issue reported by Cabellos, NEA). *ident up49 */ 10feb2015 */ - errorr */ - an array overflow error was observed when attempting to process */ the tendl-2013 83Kr file (reported by Cabellos, NEA). */ - the issue appears to be a evaluation error as the mf2/mt151 */ urr region data are given for 2 l-states while the mf32/mt151 */ urr region data are given for 3 l-states. njoy utilizes the */ mf2/mt151 l-state value when reading mf32/mt151 and so can't */ properly read this section. */ - we add a test and error stop when encountering this condition. */ - we also allow the mf32 nls be less than mf2, but advise the */ user that these covariance data are incomplete. *ident up50 */ 12feb2015 */ - acefc */ - the suite of standard ace plots is sometimes incomplete with some */ emitted neutron angular distribution and secondary neutron spectra */ plots missing (reported by Trkov, IAEA). */ - a global variable used to keep track of the number of neutron */ emitting reactions was erroneously used as a local variable */ for delayed neutron related data processing. this does not */ impact the underlying ace file since we have already written */ out the correct value and associated arrays, but did impact */ the do loop range used to control production of selected ace */ file plots. *ident up51 */ 26feb2015 * - groupr */ - correct a typo in the Kalbach-Mann formalism constants (reported */ by Huria, Westinghouse). */ - current coding uses the "k86" parameter variable to skip over */ this no longer used, and obsolete, coding. rather than */ correct the typo we'll delete these code segments. *ident up52 */ 26feb2015 */ - thermr */ - while recasting a do while loop in njoy 99 thermr/iel into a */ more traditional do loop in njoy2012 we omitted a test that */ forces a zero cross section for thermal elastic scattering at */ an energy point just beyond the user emax. */ - when later processed by groupr this yields non-zero values */ in groups beyond the thermal cutpoint. *ident up53 */ 10mar2015 */ - heatr */ - coding introduced in up47 to test for redundant photon data in */ mf6 and mf12 failed to initialize an array index. jobs that */ process multiple temperatures may suffer an array overflow */ error, depending upon the number of temperatures being processed */ and the number of mf6 sections present in the input file (issue */ reported by Orsi, ENEA). *ident up54 */ 06apr2015 */ - no code changes; but a correction to the manual. *ident up54m */ 06apr2015 */ - matxsx */ - thermal material names given in Table 21 are incomplete */ and/or don't reflect endf/b-vii nomenclature (discovered */ while debugging issues encountered by Kornreich (LANL) */ when attempting to run transx/partisn with matxsr generated */ files). *ident up55 */ - 03jun2015 */ - purr */ - the rdf2un routine doesn't properly skip over lrf=7 formatted */ resolved resonance data and so doesn't detect the subsequent */ urr data (seen in preliminary isoW files from ORNL). */ - a similar issue exists in unresr and is fixed in the */ next update. *ident up56 */ - 03jun2015 */ - unresr */ - fix rdunf2 as was done in up55 for purr's rdf2un *ident up57 */ - 04jun2015 */ - aceth */ - up34.29 exchanged one typo for another, sigh. this caused */ the thermal ace itce pointer to be wrong and coherent elastic */ data were not properly written to the thermal ace file. a */ subsequent acer check/plot (i.e., iopt=7) job fails with an */ unexpected eof error. issue reported by Zheng (Texas A&M). *ident up59 */ - 22jun2015 */ - groupr */ - level number and isomer number are used interchangably in mf6, */ mf8, mf9 and mf10. specifically, when mf8 points to mf6 we */ are dealing with "isomer number" but when mf8 points to mf9 */ or mf10 we are dealing with "level number". current groupr */ coding does not recognize this difference and if mf8 points */ to mf6 and the level number and isomer number differ we either */ process the wrong data or die with a data not found message */ (issue reported by Cornock, AWE). */ - revised coding attempts to extract the appropriate isomer */ number when auto processing (i.e., user input = 10/) by */ checking if the current 1000*Z+A read in mf8 is the same */ as the previous value. */ - for explicit user input, when requesting Xzzzaaam, treat */ "m" as isomer number for X=2 and treat "m" as a level */ number when X=3 or 4. */ - also, previous coding dealing with mfd>10000000 and lfs was */ inefficient and sometimes redundant. we clean that up now. *ident up60 */ - 14jul2015 */ - errorr */ - maxe too small error encountered when processing endf/b-vii.1 */ 58Ni (reported by Bruss, TAMU). */ - a multiplicative factor that should always be greater */ than unity can, in fact, be less than one for very small */ low energy groups. revise the calculation so that this */ factor is greater than one while guaranteeing at least */ one energy point will be sampled within this small group. */ - add a couple of timing statements so the user knows */ we're stilling working. */ - while debugging this issue we noticed the subroutine */ identified by the "maxe too small" error message could be */ incorrect and so fix those messages also. *ident up61 */ - 21aug2015 */ - plotr */ - need more space in plotr. */ - use parameter statements for arrays for easy future */ upgrade. */ - make the several "storage exceeded" error messages unique */ to better understand where this failure originates. *ident up62 */ - 21aug2015 */ - viewr */ - keep array sizes consistent with those in plotr (up61) *ident up63 */ 06sep2015 */ - mixr */ - add lrp=-1 flag in mf1/mt451 to show no mf2. */ - add awi, emax and nsub in mf1/mt451 if a version 6 file */ (missing data reported by Simakov, IAEA). */ - awi comes from the first input file; */ - emax is the largest value from the multiple input files; */ - nsub depends upon awi. */ - default is nsub=10 (incident nuetron), but recognize the */ possibility of incident photon, proton, deuteron, triton, */ 3he or alpha. *ident up64 */ 09nov2015 */ - heatr */ - coding introduced in up53 to make sure the i6p variable, introduced in */ up47, is initialized sits within an if statement that may not always be */ executed if processing multiple materials. move it to a more general */ location (and for good measure make sure we initialize the mt6yp array). */ - issue noted (and a solution offered (thank you!) by Cabellos (NEA)). *ident up65 */ 15aug2016 */ - heatr */ - array index over- or under-flow issues in h6ddx. */ - insert an additional array index test to prevent this. can occur */ when the "ep" value passed to this function is outside the energy */ range of the underlying mf6/mt secondary emission data. */ - note, coding in groupr that performs similar calculations is */ laid out differently and does not seem to suffer from the array */ index issue seen here; hence no update is planned for that module. */ - make some arrays larger. */ - delete the local "izat" variable from h6ddx since its already defined */ as a global variable, but reposition where we define it in nheat to */ make sure its defined when we need it. *ident up66 */ 17aug2016 */ - acecm */ - tighten the criteria for angular distribution reconstruction; */ - make some fixed arrays allocatable to eliminate the chance of */ a silent array overflow error. */ - acephn */ - correction to a tab1 nbt pointer value. */ - these changes are in response to concerns raised by Trkov (IAEA). *ident up67 */ 17aug2016 */ - broadr */ - broadr can seg fault if the user specifies a thnmx value */ greater than the maximum energy in the input file. insert */ a check and make a graceful exit. */ - when there are no RR parameters (lru=0 or 2), restrict */ broadening to no more than 1 MeV regardless of the file's */ EH value. cross sections are already smooth at these energies */ and there is no technical benefit to Doppler broadening. Also, */ Trkov (IAEA) reports possible issues in select files when */ attempting to broaden to higher energies. *ident up68 */ 17aug2016 */ - rpxlc2 */ - the mpid array values are dependent upon the RR formalism */ (actually a long-standing error initially reported by */ Trkov when still at the Jozef Stefan institute, but */ never corrected). */ - this update covers slbw, mlbw and rm plus an error */ message saying the other formalisms aren't treated (yet!). *ident up69 */ 18aug2016 */ - groupr */ - the dimension limit of 3 in getmf6 for iyss, izss and jjss means */ we're assuming no more than a ground state plus two isomers for */ a given ZA. Recent TENDL files may exceed this limit (issue */ reported by Sublet, CCFE). *ident up70 */ 18aug2016 */ - purr */ - legacy coding assumes the upper urr limit will not exceed 1 MeV; */ this assumption is no longer valid so let's increase that limit */ to 5 MeV (issue noted by Haeck (IRSN) and Trkov (IAEA)). *ident up71 */ 18aug2016 */ - plotr */ - a format change, the addition of a tab1 record, in tsl mf6 */ energy distribution data between njoy99/thermr and */ njoy2012/thermr wasn't accounted for in njoy2012/plotr */ (issue reported by Cabellos, NEA). */ - we add the missing call in plotr's ed3d routine now. */ - but note this is only needed when processing files produced */ by njoy2012/thermr. if a legacy njoy99 or earlier file is */ being processed, this patch will fail. */ - fix some long-standing (but harmless) typos in if-tests that */ should only cover the mt=221 to mt=250 range associated with */ tsl files (ack). *ident up72 */ 20aug2016 */ - acefc */ - we add "za" and "suff" to calculate the mcnp zaid, but if "za" */ has accidently picked up a small decimal component during */ creation of the original endf-formatted file, say for 235u */ we have 92235.01 rather than 92235.00, we won't get what we */ want. switch to using the integer za in acelod to avoid */ this error (ack). *ident up73 */ 22sept2016 */ - errorr */ - sdev and cov array allocation in rpxsamm is insufficient for */ lrf=7, lcomp=2. issue noted while debugging a Leal 16O CIELO */ candidate file. *ident up74 */ 22sept2016 */ - errorr */ - we failed to notice mt limit tests embedded in errorj when */ importing that code back into njoy several years ago were */ not up to date. */ - mt values in the 875 to 890 range are legal, and used by */ the JEFF community. */ - issue noted by Kodeli (Jozef Stefan Institute, IJS) and */ reported by Cabellos at the Spring, 2016 JEFF meeting). *ident up75 */ 15dec2016 */ - broadr */ - Re-think the logic for determining the upper energy limit for */ Doppler broadening. Too many new evaluations now have threshold */ reactions within the resolved resonance energy range and too */ often legacy default choice is not what the end expects or */ wants. */ - change logic so that if user input is thnmax=0 this means */ Doppler broadening will occur to the top of the RRR. If this */ or any other user choice causes threshold reactions to be */ broadened we know the Q-value and will zero out any spurious */ non-zero cross sections below the reaction threshold. */ - remaining options related to thnmax are unchanged. */ - caution ... up67 in this file differs from that provided to */ selected users in pre-release versions of this */ update file. *ident up76 */ 15dec2016 */ - samm */ - an ORNL patch to SAMRML that we include here to remain consistent */ with their coding. this will reduce the chances of calling the */ coulx routine which can be unstable under certain combinations of */ rho and eta, hence possibly impacting calculation of charged */ particle reaction cross sections. *ident up77 */ 16dec2016 */ - this update was actually created in nov2015 while at the IAEA but */ was overlooked during the transition from my LANL "LOFT" computer. */ - we insert it into the mainstream update sequence now ... */ - 13nov2015 */ reconr */ - reconr currently writes the user "tempr" input to the output */ pendf, but this variable is almost always zero. if the input */ tape has already been broadened to a non-zero value this */ information is lost. */ - revise reconr to copy the non-zero input endf temperature */ to the output pendf if tempr is zero. */ - speciality files, such as dosimetry, may only contain redundant */ file 3 reaction data. reconr gets confused when there is no */ file 3 data to be found after the redundant sections are removed. */ - insert upnea119 patch from Andrej Trkov (IAEA) to overcome */ this. */ - issue was discovered when attempting to process IRDFF */ v1.xx 109Ag. *ident up78 */ 16dec2016 */ - thermr */ - up78 (thermr), up79 (aceth) and up80 (leapr) are a suite of */ changes recommended by Damian (CAB) and Roubtsov (CNL) */ primarily for processing revised light and heavy water */ thermal kernels, but are ok to apply to other data sets. */ most thermr changes are related to carrying more precision */ (changes in sigfig calls) in calcem, plus a slight revision */ to the egrid array. we also use the ngrid parameter to set */ the egrid (and other) array size. we also revise the write */ format when iprint=2 to provide more digits and to appear */ similar in form as written by acer. *ident up79 */ 16dec2016 */ - aceth */ - up78 (thermr), up79 (aceth) and up80 (leapr) are a suite of */ changes recommended by Damian (CAB) and Roubtsov (CNL) */ primarily for processing revised light and heavy water */ thermal kernels, but are ok to apply to other data sets. */ most thermr changes are related to carrying more precision */ (changes in sigfig calls) in calcem, plus a slight revision */ to the egrid array. we also use the ngrid parameter to set */ the egrid (and other) array size. we also revise the write */ format when iprint=2 to provide more digits and to appear */ similar in form as written by acer. *ident up80 */ 16dec2016 */ - leapr */ - up78 (thermr), up79 (aceth) and up80 (leapr) are a suite of */ changes recommended by Damian (CAB) and Roubtsov (CNL) */ primarily for processing revised light and heavy water */ thermal kernels, but are ok to apply to other data sets. */ most thermr changes are related to carrying more precision */ (changes in sigfig calls) in calcem, plus a slight revision */ to the egrid array. we also use the ngrid parameter to set */ the egrid (and other) array size. we also revise the write */ format when iprint=2 to provide more digits and to appear */ similar in form as written by acer. *ident up81 */ 22jan2017 */ - acefc, aceth, acedo, acepa, acepn */ - when writing the xsdir stub make sure we leave a space prior */ to the len2 variable. the "i8" field used up until now has */ always included a space if the ace file length is less than */ 100M but a recent pre-endf/b-viii.0 beta file was larger. */ 22jan2017 */ - reconr */ - we fail to properly set the lssf flag for the case of energy- */ dependent unresolved resonance fission widths. this results */ in double counting of cross sections in the urr. */ - the lssf flag has been defined properly for the case of */ energy-independent urr widths or when all urr widths are */ energy-independent. */ - thanks to Andrej Trkov (IAEA) for discovering this long- */ standing error. */ - we note that the coding in the unresr and purr modules does not */ suffer from this error.