*/
*/ OECD/NEA compilation of NJOY updates
*/ compatible with the official patch distribution up245
*/
*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
*/ groupr A. Trkov, 24 July 2007
*/        For backward compatibility change error stop to a warning message
*d up167.22,23
            if (sl.lt.zero.or.sn.lt.zero) call mess('getunr',    
     &      ' negative cross sections found - check unresr',' ')
*/
*/ ident upnea020 is included in up235
*/
*ident upnea021
*/ groupr A. Trkov 30 Oct 2007
*/        Lahey compiler does not like initialisation of variables in
*/        common with a data statement.
*d up257.10
      ebeg=1.d-5
*d up257.12
      ebeg=1.e-5
*/
*ident upnea022
*/ errorr A. Trkov 30 Oct 2007
*/        Lahey compiler identified several errors and warnings:
*/        - Subroutine "error" must be declared external in several routines.
*/        - Extend common "ewght" similar to common "weight" in groupr,
*/          introduce dummy common "temper", initialize with temperature "tempin",
*/          both needen in subroutine "egtwtf" (a more sophisticated patch
*/          might be needed).
*/        - Replace Hollerith constants "hmt" and "uline" with character.
*/        - Define constant "small" in two subroutines.
*/        - Define constant "a33" (PLEASE CHECK)!
*i errorj.258
      external error
*i errorj.731
      external error
*i errorj.1362
      external error
*i errorj.2758
      external error
*i errorj.2996
      external error
*i errorj.3138
      external error
*i errorj.3487
      external error
*i errorj.3757
      external error
*i errorj.4099
      external error
*i errorj.5380
      external error
*i errorj.5825
      external error
*i errorj.5990
      external error
*i errorj.6360
      external error
*i errorj.6733
      external error
*i errorj.6853
      external error
*i errorj.7077
      external error
*i errorj.7142
      external error
*i errorj.7379
      external error
*i errorj.7433
      external error
*i errorj.7603
      external error
*i errorj.7655
      external error
*i errorj.7781
      external error
*/ extend common "ewght" similar to common "weight" in groupr
*/ introduce dummy common "temper", initialize with temperature "tempin"
*/ both needen in subroutine "egtwtf"
*d errorj.242
      common/ewght/iwt,jsigz,jtemp
      common/temper/temp(10),ntemp
*i errorj.351
      ntemp=1
      jtemp=1
      temp(1)=tempin
*d errorj.1356
      common/ewght/iwt,jsigz,jtemp
*d errorj.3484
      common/ewght/iwt,jsigz,jtemp
*d errorj.4090
      common/ewght/iwt,jsigz,jtemp
*d errorj.5375
      common/ewght/iwt,jsigz,jtemp
*d errorj.5568 
      common/ewght/iwt,jsigz,jtemp
*d errorj.7260
      common/ewght/iwt,jsigz,jtemp
*/ Replace Hollerith constants with character
*i errorj.7769
      character*2 hmt
      character*5 uline
*d errorj.7782 
      data hmt/'mt'/, uline/'-----'/
*i errorj.8108
      character*2 hmt
      character*5 uline
*d errorj.8119
      data hmt/'mt'/, uline/'-----'/
*/ define constant small                 (please, check !!!!!)
*i errorj.8275
      data small/1.d-10/
*i errorj.8497
      data small/1.d-10/
*/ define constant a33                   (please, check !!!!!)
*i errorj.8622
      a33 =a3*a3
*ident upnea023
*/ errorr A. Trkov 04 Dec 2007
*/      - Insufficientlength for the array at index iscr was reserved,
*/        causing the next block of data to be corrupted.
*/        All cases are affected where union grid>654 energy points.
*/      - Increase array size from 8M to 12M to process covariance
*/        matrices in 640 groups (external reference spectra for
*/        dosimetry).
*d errorj.237  
      common/estore/a(12500000)
*d errorj.272  
      namax=12500000
*d errorj.3131 
      common/estore/a(12500000)
*d errorj.2772 
      nwds=max(npage+50,nunion+1+8)
*d errorj.3499 
      nwscr=max(2*npage+50,nun1+8)
*i errorj.3520 
      if (nw.gt.nwscr)
     &   call error('colaps','storage exceeded.',' ')
*d errorj.3758 
      data nxmax/800/, irmax/60/
*ident upnea024
*/ covr A. Trkov 04 Dec 2007
*/      Increase array size from 300k to 900k to process covariance
*/      matrices in 640 groups (external reference spectra for
*/      dosimetry).
*d up111.12   
      common/storec/a(900000)
*d up111.8    
      common/storec/a(900000)
*d up111.10   
      data iamax/900000/, niad/15/, ipr/1/, ntics3/600/