*/ */ OECD/NEA compilation of NJOY updates */ compatible with the official patch distribution up364 */ */ The patches include those carried-over from the set compatible */ with the official up259 but not adopted in up296. */ *ident upnea001 */ njoy - C. Broeders, 6-Oct-2006 */ Open scratch files in local disc area */ Ref.: Comments by D.E. Cullen in PrePro-2004 (Scratcha.f) on p.19 of: */ "http://www-nds.iaea.org/ndspub/endf/prepro/DOCUMENT/PDF/Overview.pdf" */ Status "scratch" does not work with Lahey compiler on Linux, */ therefore 'age" is set "unknown" *d njoy.470 age='unknown' *d njoy.471 write(fn,'(a,i2.2)') 'temp',nun open(nun,file=fn,form=for,status=age) */ *ident upnea004 */ groupr - C. Broeders, 6-Oct-2006 */ extend IWT=5 spectrum definition up to 200MeV */ ------------------------------------------- *d groupr.2145 dimension w1(92),w2(92),w3(10) *d groupr.2150 data w1/0d0,0d0,0d0,0d0,1.d0,93.d0,93.d0,5.d0, *i groupr.2175 data w3/3.d7,1.0318d-10,5.d7,6.1908d-11,1.d8,3.0954d-11, &1.5d8,2.0636d-11,2.d8,1.5477d-11/ *d groupr.2358 iw=194 *i groupr.2365 do i=1,10 a(i+183+iwght)=w3(i) enddo */ */ *ident upnea019 is included in up265 */ */ *ident upnea020 is included in up235 */ */ *ident upnea021 is included in up261 */ */ *ident upnea022 is included in up260 */ */ *ident upnea023 is included in up272 */ */ *ident upnea024 is included in up273 */ */ *ident upnea025 is included in up315 */ */ *ident upnea026 is included in up316 */ */ *ident upnea027 is included mostly in up317, except a few statements */ which are placed in upnea068 */ */ *ident upnea028 is included in up318 */ */ *ident upnea029 is included in up316 */ */ *ident upnea030 is included in up272 */ */ *ident upnea031 is included in up289 */ *ident upnea032 */ acer A. Trkov, 30 May 2008 */ Plots of discrete inelastic make sense below 10 MeV *d acer.19586 c limit discrete inelastic plots to 10 MeV if (x.gt.xmax) xmax=min(x,10.) */ */ *ident upnea033 included in up278 */ */ *ident upnea034 included in up277 */ */ *ident upnea035 included in up279 */ */ *ident upnea036 included in up279 */ */ *ident upnea037 included in up277 */ *ident upnea038 */ covr A. Trkov, 20 October 2008 * Fix diagnostic message *d covr.1793,1805 zero=0 if (abs(crit).gt.zero) write(strng, & '(''i'',i3,'' j'',i3,'' xa(i,j)='',1p,e12.4,'' xa(j,i)='', & 1p,e12.4)') i,im1,xa(im1ind,i),xa(iind,im1) if (abs(crit).gt.test .and. abs(xa(im1ind,i)).gt.1e-20) then call error('press','matrix not symmetric',strng) endif test=1 test=test/1000000 if (abs(crit).lt.test) go to 250 call mess('press','matrix not symmetric',strng) */ */ *ident upnea039 included in up272 */ *ident upnea040 */ */ GROUPR - R. Perry, 15 October 2008 */ (Original update name: upwe6) */ Problem with thermal scattering matricies in GROUPR */ There were small discrepancies between the sums thermal scatter */ matrices and the group cross-sections. This update renormalises */ the interpolated spectra (produced during integration) */ to sum to unity. */ (guarding against array overflow added by A. Trkov) *i groupr.6915 integer nsi parameter (mxnsi=15000) *i groupr.6922 dimension sint(2,mxnsi) *i groupr.7030 nsi = 0 *i groupr.7146 c c store interpolated function. c nsi = nsi+1 if(nsi.gt.mxnsi) & call error('getaed','mxnsi array limit exceeded',' ') sint(1,nsi) = ei sint(2,nsi) = fi(1) *i groupr.7155 c c sum interpolated spectrum and adjust aed. c sisum = 0.0 do il = 1, nsi - 1 sisum = sisum + (sint(2,il)+sint(2,il+1))* & (sint(1,il+1)-sint(1,il))/2 enddo c c don't re-normalise if scatter source above top energy group. c if(i .lt. ngn .or. (k1 .eq. nlo .and. k2 .eq. nhi) ) then do il = 1, i aed(1,il) = aed(1,il)/sisum enddo endif */ */ *ident upnea041 included in up282 */ */ *ident upnea042 is not included; it is redefined in upnea069 */ */ *ident upnea043 included in up282 */ */ *ident upnea044 included in up273 */ */ *ident upnea045 included in up282 */ */ *ident upnea046 withdrawn */ */ *ident upnea047 included in up287 */ */ *ident upnea048 included in up329 and extended */ */ *ident upnea049 is included in up290 */ */ *ident upnea050 is included in up289 */ */ *ident upnea051 is included in up363 */ */ *ident upnea052 is included in up361 */ */ *ident upnea053 is included in up307 */ *ident upnea054 */ errorr D.L. Aldama, July 2009 */ Fix weighting function option *d errorj.660 if (iwt.eq.1.or.iwt.eq.4.or.iwt.eq.5.or. & (iwt.ge.7.and.iwt.lt.11)) *d errorj.1492 if (iwt.eq.1.or.iwt.eq.4.or.iwt.eq.5.or. & (iwt.ge.7.and.iwt.lt.11)) */ */ *ident upnea055 is included in up314 */ */ *ident upnea056 is included in up312 */ */ *ident upnea057 is included in up305 */ */ *ident upnea058 is included in up318 */ */ *ident upnea059 is included in up205 */ */ *ident upnea060 is not needed due to alternative coding in up329 */ */ *ident upnea061 is superseeded by upnea068 */ */ *ident upnea062 is included in up324 */ */ *ident upnea063 is dealt with differently in upnea315 */ *ident upnea064 */ acer A. Trkov, February 2010 */ Particle emission spectra in the acer plots are corrupted */ for MF6 Law 1 LANG 1 when more than one interpolation range */ is specified for the distributions on incident particle */ energies. In some cases the MCNP calculations are also wrong. */ A fix is done when writing the data to a temporary file. */ The multiple ranges are suppressed. The first interpolation */ law is prescribed over the entire incident energy range */ and a message is printed. */ WARNING: */ This is a temporary patch before a proper solution is found. */ The true error probably occurs somewhere near acer.6651 */ or later, and/or possibly in MCNP. */ Implications: */ JENDL-HE files truncated to 150 MeV (as proposed for FENDL-3) */ led to strongly discrepant results in a benchmark model */ calculation representative of the ITER device. The discrepancy */ is greatly reduced with the use of this patch. The assumption */ in the patch has no influence on the evaluated data because */ the interpolation law changes above 20 MeV where the yield */ drops to zero due to the reaction representation in MT5. *d acer.2373,2374 call tab2io(nin,0,0,b,nb,nw) nr=n1h ne=nint(b(6)) if(nr.gt.1) then write(string,'(a,i3)') & 'multiple interp. ranges for mf6, mt',mt call mess('topfil',string & ,'first law applied everywhere') nr=1 b(5)=nr b(7)=ne end if call tab2io(0,nout,0,b,nb,nw) */ */ *ident upnea065 is included in up341 */ */ *ident upnea066 is included in up342 */ */ *ident upnea067 is included in up363 */ *ident upnea068 */ acer A. Trkov, January 2010 */ (Revision of upnea061, noticed by Skip Kahler) *d acer.119,120 c * 7 read fast ace files to print or edit * c * 8 read thermal ace files to print or edit * c * 9 read dosimetry ace files to print or edit * *d acer.407 else if (iopt.eq.6.or.iopt.gt.9) then *d acer.435 else if (iopt.ge.7.and.iopt.le.9) then *i acer.17639 c c ***read type 3 ace format file else if (itype.eq.3) then if (mcnpx.eq.0) then read(nin) hz(1:10),aw0,tz,hd,hko,hm,(izo(i),awo(i),i=1,16), & (nxs(i),i=1,16),(jxs(i),i=1,32) else read(nin) hz(1:13),aw0,tz,hd,hko,hm,(izo(i),awo(i),i=1,16), & (nxs(i),i=1,16),(jxs(i),i=1,32) endif len2=nxs(1) n=(len2+ner-1)/ner l=0 do i=1,n max=len2-l if (max.gt.ner) max=ner read (nin) (xss(l+j),j=1,max) l=l+max enddo call closz(-nin) c flag incident particle "undefined" izai=-1 awi =-1 */ */ *ident upnea069 is included in up362 */ */ *ident upnea070 is included in up363 */ */ *ident upnea071 is included in up363 */ *ident upnea073 */ acer A. Trkov, January 2011 */ Error message lenth exceeds string size *d up360.91,93 write(strng,'(''mtd='',i3,'' mt='',i6, & '' ie='',i4,i5,'' nd='',3i4,'' ed='', & 1p,e9.3e1)') */ *ident upnea074 */ acer A. Trkov, June 2011 */ zaid is re-defined incorrectly when editing fast ace files if original */ suffix is greater than 0.50 *d acer.17652 iza=int(zaid+0.001) */ */ *ident upnea075 is included in up363 */ *ident upnea076 */ */ HEATR - K. Vignitchouk, 06/09/11 */ Calculate mean neutron energy for law 12-energy */ dependent fission neutron spectrum (madland-nix) */ (originaly submitted as upwd20) */ *i heatr.2169 c c ***law 12--energy dependent fission neutron spectrum (madland-nix) else if (lf.eq.12) then s=((a(lnext)+a(lnext+1))/2)+(4*theta/3) */ *ident upnea077 */ broadr A. Trkov, October 2011 (see also upnea100) */ Increase the maximum number of non-threshold cross sections */ (case: NNDC evaluation of Tm-168) *d up175.13 parameter (nttmax=48) *d up175.21 parameter (nttmax=48) *d up175.25 parameter (nttmax=48) *d up175.30 parameter (nttmax=48) *d up175.34 parameter (nttmax=48) */ *ident upnea078 */ ACER - D. L. Aldama (IAEA/NDS Consultant) */ Nov./2011 */ Correction to NJOY update up360 for processing */ MF 6 with different number of discrete */ photon energies by sections */ (needed to process the FENDL-3 library) *d up360.17 nwordg=2000 *d up360.94 a(idise)=ep *i up360.95 else if (abs(ep).lt.abs(a(idise+nd0-1))) then a(idise+nd0)=ep nd0=nd0+1 *d up360.123 & nwordg,nd0 */ *ident upnea079 */ HEATR - D. L. Aldama (IAEA/NDS Consultant) */ Nov./2011 */ Correction for large MF4/MT=2 sections */ (needed to process the FENDL-3 library) *i heatr.3832 if (nb.ne.0) then l=iraw do while (nb.ne.0) l=l+nwc call moreio(nin,0,0,c(l),nb,nwc) enddo endif *i heatr.3837 if (nb.ne.0) then l=iraw do while (nb.ne.0) l=l+nwc call moreio(nin,0,0,c(l),nb,nwc) enddo endif *i heatr.3877 if (nb.ne.0) then l=iraw do while (nb.ne.0) l=l+nwc call moreio(nin,0,0,c(l),nb,nwc) enddo endif */ *ident upnea080 */ errorr - A. Trkov, December 2011 */ On multiple reaction mf34 processing errorr tries to read */ spectrum specifications for IWT=4 for each reaction. */ Flag eb negative to suppress reading for multiple reactions. */ (A similar patch might be needed in the genwtf routine */ in groupr). *i errorj.253 common/iwt4/eb,tb,ab,ec,tc,ac c *i errorj.278 eb=-1 *d errorj.5481 if(eb.lt.0) read(nsysi,*) eb,tb,ec,tc */ *ident upnea081 */ errorr - O. Cabellos, January 2012 */ In order to compile with VDF6.6A in WindowsXP, some dimensions */ in ERRORR module are changed without effect in this processing */ work. In NJOY99.364 the dimension of this vectors is 50 000 000, */ in this work, it is reduced to 9 200 000. */ (The patch is compiler-dependent, so it was made inactive, */ but can be activated by deleting the "*/" string (i.e the */ first two columns of the patch. (Implemented by A. Trkov)) */*d up344.6 */ common/estore/a(9200000) */*d up344.8 */ namax=9200000 */*d up344.10 */ common/estore/a(9200000) */ *ident upnea082 */ gaspr - O. Cabellos, January 2012 */ Typos in up353 for the GASPR module are identified and corrected *d up353.17 if (mth.eq.152) izr=izr-5 *d up353.27 if (mth.eq.162) izr=izr-1006 *d up353.37 if (mth.eq.172) izr=izr-1006 *d up353.47 if (mth.eq.182) izr=izr-2005 *d up353.57 if (mth.eq.192) izr=izr-3005 */ *ident upnea083 */ matxsr - O. Cabellos, nov2011 */ processing JEFF-3.1.2 with VDF6.6A on WindowsXP *i matxsr.1965 integer*4 jg1lo,jg1hi *i matxsr.2066 integer*4 jg1lo,jg1hi */ *ident upnea084 */ acer - O. Cabellos, nov2011 */ processing JEFF-3.1.2 (H3, He3, Y90, etc.) with VDF6.6A on WindowsXP *i up45.9 if (ngmt.eq.0) call closz(-iold) if (ngmt.eq.0) call closz(-inew) */ *ident upnea085 */ errorr - O. Cabellos, January 2012 */ In order to compile with VDF6.6A in WindowsXP, some dimensions */ in ERRORR module are changed without effect in this processing */ work. In NJOY99.364 the length of nwds is 10 000 000, */ in this work, it is reduced to 1 000 000. */ (The patch is compiler-dependent, so it was made inactive, */ but can be activated by deleting the "*/" string (i.e the */ first two columns of the patch. (Implemented by A. Trkov)) */*d up363.106 */ nwds=1000000 */ */ *ident upnea086 */ */ groupr - O. Cabellos, nov2011 */ */ Fix the input request for specific lfs */ */ processing JEFF-3.1.2 with VDF6.6A in WindowsXP */ */ (This ident interferes with automatic processing */ */ of data in MF 10; A. Trkov) */ *i up347.64 */ if (mfd.gt.10000000) lfs=mod(mfd,10) */ *ident upnea087 */ gaminr - O. Cabellos, nov2011 */ processing JEFF-3.1.2 with IVF in Windows7 */ To compile NJOY with Intel Visual Fortran, Windows7 (64b), */ it is necessary to deactivate the compiler option */ "Check Routine Interfaces". This compiler shows a problem */ in GAMINR soubroutine. The variable z(10) should be */ dimensioned to z(100) for consistency with the dimension */ declaration in gtflx. *d gaminr.92 dimension z(100) */ *ident upnea088 */ matxsr - Go Chiba, January 2012 */ increase maxw in matxsr *d matxsr.512 maxw=50000 */ *ident upnea089 */ errorj - Go Chiba, January 2012 */ Fix problem when the width of a user energy grid is quite thin. *i up272.271 y1=sig(i,j) y2=sig(i-1,j) *d up272.275 gsig(j,ig)=(yl*zl+yr*zr)/sumde */ *ident upnea090 */ broadr - O. Cabellos, June 2009 (renamed ident janis2) */ Problems in some calculations at 293.6 K */ processing JEFF-3.0/A and EAF-2007 *d broadr.520 if(inutot.eq.0) then fnu=0.0 else call terpa(fnu,therm,xnext,idis,a(inutot),ip,ir) endif *d broadr.541 if(inutot.eq.0) then fnu=0.0 else call terpa(fnu,tev,xnext,idis,a(inutot),ip,ir) endif *d broadr.557 if(inutot.eq.0) then fnu=0.0 else call terpa(fnu,enow,xnext,idis,a(inutot),ip,ir) endif */ *ident upnea091 */ reconr - O. Cabellos, June 2009 (renamed ident janis4) */ Increase of necessary array size */ processing Cr-54 from BROND-2.2 *d reconr.2580 dimension sigp(4), sigj(16,2), gj(16), a(*) */ *ident upnea092 */ broadr - O. Cabellos, June 2009 (renamed ident janis5) */ Increase of necessary array size */ (no details given) *d broadr.140 nwscr=1500 */ *ident upnea093 */ njoy - O. Cabellos, June 2009 (renamed ident janis7) */ increase page size (see also upnea109) *d njoy.149 npage=612 */ *ident upnea094 */ heatr - O. Cabellos, June 2009 (renamed ident janis7) */ increase array size of elist (JENDL-3.3) */ increase dimension of array a (janis8) */ (patches for specific materials in janis8 not included) *d heatr.91 common/heat2/qa(3),tempr,efirst,elast,za,awr,elist(200) *d heatr.410 common/heat2/qa(3),tempr,efirst,elast,za,awr,elist(200) *d heatr.781 common/heat2/qa(3),tempr,efirst,elast,za,awr,elist(200) *d heatr.4405 common/heat2/qa(3),tempr,efirst,elast,za,awr,elist(200) *d heatr.4997 common/heat2/qa(3),tempr,efirst,elast,za,awr,elist(200) *d heatr.114 ilmax=200 *d heatr.94 dimension a(150000) *d heatr.110 namax=150000 */ *ident upnea095 */ heatr - Kazuaki Kosako, Shimizu Corporation (renamed ident janis9) */ kinematic method */ (patch received from O. Cabellos) *i heatr.55 c * kkerma 0/1=total (mt301) is energy balance / kinematic * c * (default=0) * *i heatr.90 common/heat6/kkerma *d heatr.139 kkerma=0 read(nsysi,*) matd,npk,nqa,ntemp,local,iprint,break,kkerma *i heatr.151 if (kkerma.eq.1) then if (npk.gt.0) then mtkk=0 do i=1,npk if (mtk(i).eq.443) mtkk=i enddo if (mtkk.eq.0) then npk=npk+1 if (kchk.eq.1) then npkk=3*npk+7 else npkk=npk+3 endif if (npkk.gt.npkmax) call error('heatr', & 'requested too many kerma mt-s (6+mt301 allowed).',' ') mtk(npk)=443 endif else npk=1 mtk(npk)=443 endif endif *i heatr.4996 common/heat6/kkerma *i heatr.5045 if (kkerma.eq.1) then do i=2,npk if (mtp(i).eq.443) mtkk=i enddo endif *i heatr.5121 if (mt.eq.301.and.kkerma.eq.1) a(ibase+i)=c(mtkk) *i heatr.5146 if (kkerma.eq.1) write(nsyso,'(/'' total kerma (mt=301) was '', & ''replaced to kinematic kerma (mt=443).'')') */ *ident upnea096 */ gaspr - O. Cabellos, June 2009 (renamed ident janis10) */ Fix residual ZA *d gaspr.522 izr=izr-3006 */ *ident upnea097 */ gaspr - O. Cabellos, June 2009 (renamed ident janis11) */ Additional updates to take into accoutn MF=10 branching reactions */ (WARNING: If mf3 and mf10 are both present, check if there is */ double counting) *d up244.47,56 if (mfi.eq.3.or.mfi.eq.10.and.mti.eq.103) mt103=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.eq.104) mt104=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.eq.105) mt105=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.eq.106) mt106=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.eq.107) mt107=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.ge.mpmin.and.mti.le.mpmax) & mt600=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.ge.mdmin.and.mti.le.mdmax) & mt650=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.ge.mtmin.and.mti.le.mtmax) & mt700=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.ge.m3min.and.mti.le.m3max) & mt750=1 if (mfi.eq.3.or.mfi.eq.10.and.mti.ge.m4min.and.mti.le.m4max) & mt800=1 *d gaspr.102 if (mfi.gt.10) then *d up244.82 if (mth.gt.m4max.or.math.eq.-1.or.math.eq.0) go to 250 if (mth.eq.0) go to 240 if (.not.(mfh.eq.3).and. .not.(mfh.eq.10)) goto 240 *i gaspr.346 if10=0 nf10=n1h 2700 continue *d up244.108 if (mth.gt.m4max.or.math.eq.-1.or.math.eq.0) go to 330 if (mth.eq.0) go to 270 if (.not.(mfh.eq.3).and. .not.(mfh.eq.10)) goto 310 *i gaspr.574 if(mfh.eq.10) then if10=if10+1 if (if10.eq.nf10) goto 310 goto 2700 endif */ *ident upnea098 */ gaspr - O. Cabellos, June 2009 (renamed ident janis12) *d gaspr.205,208 *d up244.61,62 *d gaspr.211,217 98 call contio(npend,0,0,b(1),nb,nw) mfb=mfh mtb=mth if (math.gt.0) then if (mfh.eq.3.or.mfh.eq.10) then call contio(0,nscr1,0,b(1),nb,nw) call tosend(npend,nscr1,0,a(1)) elseif (mfh.eq.0.and.mth.eq.0) then call afend(0,nscr1) elseif (mfh.ne.0) then call tosend(npend,0,0,a(1)) endif goto 98 end if call afend(0,nscr1) */ *ident upnea099 */ gaspr - O. Cabellos, June 2009 (renamed ident janis3) */ Fix problem processing SLTs from ENDF/B-VII.0 *d plotr.330 dimension a(11200),z(15),aa(100000) *d plotr.349 data max/2500/, nwamax/11200/ */ *ident upnea100 */ broadr -- O. Cabellos, November 2011 */ Further increase of array size introduced in upnea077 */ (processing JEFF-3.1.2 Ti46, Ti47, Ge70, Fe57, Fe56, etc.) *d upnea077.6 parameter (nttmax=160) *d upnea077.8 parameter (nttmax=160) *d upnea077.10 parameter (nttmax=160) *d upnea077.12 parameter (nttmax=160) *d upnea077.14 parameter (nttmax=160) */ *ident upnea101 */ errorr -- A. Trkov, March 2012 */ The patch up363 (former upnea051) did not provide for B-W. */ I believe the coding would not work for lrf=4,7 either. */ An extension is implemented (e.g. Na-23 from ENDF/B-VII.1). *( An error message is added. *d up279.19 dimension mpid(6),mpidbw(6),mpidrm(6) *d up363.28 data mpidrm /1,3,4,5,6,0/ data mpidbw /1,4,5,6,2,0/ *i errorj.6735 do i=1,6 if(lrf.le.2) then mpid(i)=mpidbw(i) else mpid(i)=mpidrm(i) end if end do *i up279.24 if(lrf.gt.3) then write(strng1,'(''No coding for lrf'',i4)') lrf strng2=' processing terminated' call error('Resprx_RRR_Lcomp2',strng1,strng2) end if */ *ident upnea102 */ errorr -- R. Perel, June 2012 */ Correct an inconsistent update in upnea040, which does */ not work for 64-bit compilers *d upnea040.12,13 *b groupr.6917 integer nsi parameter (mxnsi=15000) */ *ident upnea103 */ groupr -- A. Trkov, September 2012 */ Correct an inconsistent update in up339, upnea004, which do */ not work for 64-bit compilers *d up339.66 *b groupr.1026 common/argcom1/lfs *d groupr.2213 data w1/0.e0,0.e0,0.e0,0.e0,1.e0,93.e0,93.e0,5.e0, *i groupr.2238 data w3/3.e7,1.0318e-10,5.e7,6.1908e-11,1.e8,3.0954e-11, &1.5e8,2.0636e-11,2.e8,1.5477e-11/ */ *ident upnea104 */ matxsr -- A. Trkov, September 2012 */ Correct an inconsistent update in upnea083, which does */ not work for 64-bit compilers *d upnea083.5 *b matxsr.1967 integer*4 jg1lo,jg1hi *d upnea083.7 *b matxsr.2068 integer*4 jg1lo,jg1hi */ *ident upnea105 */ acer -- A. Trkov, September 2012 */ Correct an inconsistent update in up327, which does */ not work for 64-bit compilers *d up327.17 *b acer.8048 character*60 strng */ *ident upnea106 */ purr -- A. Trkov, September 2012 */ Correct an inconsistent update in up332, which does */ not work for 64-bit compilers *d up332.7 *b purr.1101 parameter (mxns0=100) *d up332.17 *b purr.1498 parameter (mxns0=100) *d up332.24 *b purr.1619 parameter (mxns0=100) *d up332.31 *b up84.30 parameter (mxns0=100) */ *ident upnea107 */ leapr -- A. Trkov, September 2012 */ Correct an inconsistent update in up304, which does */ not work for 64-bit compilers *d up304.14,15 *b leapr.1106 character*60 strng external error */ *ident upnea108 */ errorr -- W.F.G. van Rooijen, September 2012 */ MAXB is too small when processing covariances of U-238 */ with JENDL-4.0. Increase from 30000 to 40000. *d up363.110 parameter (maxe=600000,mxnpar=7000,maxb=40000) */ *ident upnea109 */ njoy -- Revert upnea093 to original form */ (Problem with this patch identified by C. Broeders) *d upnea093.5 npage=306 */