9.5.4} RESNUCLE output

 ----------------------

 Option RESNUCLE produces either formatted or unformatted output (for the
 latter, see RESNUCLE description for instructions on how to read it).
 The formatted output begins with the run title and run time, followed by
 a short information about:
   1) Total number of primary particles followed
   2) Total weight of the primaries
   3) Number and name of the residual nuclei detector, type of reactions
      considered (high energy or low energy only, or all), region number
   4) Detector volume in cm3
   5) Range of Z and N-Z printed
   6) Tabulation format
 The information reported in 3), 4) and 5) is printed also in the expanded input
 summary on main output. For instance:

     Res. nuclei n.   1  "Pb-Region " , "high" energy products, region n.     3
        detector volume:  1.0000E+00 cm**3
        Max. Z:  86, Max. N-Z:  49 Min. N-Z: -4
        data will be printed on unit   21 (unformatted if < 0)

 On the formatted RESNUCLE output, the above text is followed by one additional
 line explaining how to read the result matrix which follows:

      Data follow in a matrix A(z,n-z-k), k: -5 format (1(5x,1p,10(1x,e11.4)))

 Here is an example of a simple program which can be used to display the same
 results in a more plain way:

      PROGRAM READRN
      CHARACTER*125 LINE, FILINP, FILOUT
      PARAMETER (MAXZ = 86, MINNMZ = -4, MAXNMZ = 49, K = -5)
      DIMENSION RESULT(MAXZ, MINNMZ-K:MAXNMZ-K)

      WRITE(*,*) "Filename?"
      READ(*,'(A)') FILINP
      OPEN(UNIT=1, FILE=FILINP, STATUS='OLD')
      LQ = INDEX(FILINP,' ') - 1
      FILOUT = FILINP(1:LQ)//'.rn'
      OPEN(UNIT=2, FILE=FILOUT, STATUS='UNKNOWN')

      DO 1 I = 1, 14
         READ(1,'(A)') LINE    !  skip header lines
  1   CONTINUE

      READ(1,100,END=4) RESULT
  4   CONTINUE

      WRITE(2,'(A)')   '   Z   A    Residual nuclei'
      WRITE(2,'(A,/)') '         per cm**3 per primary'
      DO 2 I = 1, MAXZ
         DO 3 J = MINNMZ-K, MAXNMZ-K
            IF(RESULT(I,J) .GT. 0.D0)
     &       WRITE(2,'(2I4,1P, G15.6)') I, J+K+2*I, RESULT(I,J)
  3      CONTINUE
  2   CONTINUE
 100  FORMAT(1(5X,1P,10(1X,E11.4)))
      END

 ------------------------------------------------------------------------

Previous Index Next