EVGEN/PartEloss/readca3.f

   1       SUBROUTINE ELOSS_READCA3
   2       REAL*8           xx(400), da(30), ca(30,260), rrr(30)
   3       COMMON /data/    xx, da, ca, rrr
   4 *
   5       CHARACTER*100 CHROOT
   6       CHARACTER*100 FILNAM
   7       CHROOT=' '
   8       CALL GETENVF('ALICE_ROOT',CHROOT)
   9       LNROOT = LNBLNK(CHROOT)
  10
  11       IF(LNROOT.LE.0) THEN
  12          FILNAM='DATAcca3.dat'
  13       ELSE
  14          FILNAM=CHROOT(1:LNROOT)//'/FAST/eloss/DATAcca3.dat'
  15       ENDIF
  16
  17       OPEN(UNIT=20,FILE=FILNAM,STATUS='OLD',ERR=90)
  18       nn = 1
  19  100  read (20,*,end=110) xx(nn), ca(1,nn), ca(2,nn), ca(3,nn),
  20      +     ca(4,nn), ca(5,nn), ca(6,nn), ca(7,nn), ca(8,nn),
  21      +     ca(9,nn), ca(10,nn), ca(11,nn), ca(12,nn), ca(13,nn),
  22      +     ca(14,nn), ca(15,nn), ca(16,nn), ca(17,nn), ca(18,nn),
  23      +     ca(19,nn), ca(20,nn), ca(21,nn), ca(22,nn), ca(23,nn),
  24      +     ca(24,nn), ca(25,nn), ca(26,nn), ca(27,nn), ca(28,nn),
  25      +     ca(29,nn), ca(30,nn)
  26 *         print*, 0.005*(nn-1), ca(1,nn), ca(2,nn), ca(3,nn),
  27 *     +     ca(4,nn), ca(5,nn), ca(6,nn), ca(7,nn), ca(8,nn),
  28 *     +     ca(9,nn), ca(10,nn), ca(11,nn), ca(12,nn), ca(13,nn),
  29 *     +     ca(14,nn), ca(15,nn), ca(16,nn), ca(17,nn), ca(18,nn),
  30 *     +     ca(19,nn), ca(20,nn), ca(21,nn), ca(22,nn), ca(23,nn),
  31 *     +     ca(24,nn), ca(25,nn), ca(26,nn), ca(27,nn), ca(28,nn),
  32 *     +     ca(29,nn), ca(30,nn)
  33          nn = nn + 1
  34          goto 100
  35  110     continue
  36       close(20)
  37 *
  38       IF(LNROOT.LE.0) THEN
  39          FILNAM='DATAdca3.dat'
  40       ELSE
  41          FILNAM=CHROOT(1:LNROOT)//'/FAST/eloss/DATAdca3.dat'
  42       ENDIF
  43
  44       OPEN(UNIT=21,FILE=FILNAM,STATUS='OLD',ERR=90)
  45       nn = 1
  46  101  read (21,*,end=111) rrr(nn), da(nn)
  47 *         print*, rrr(nn), da(nn)
  48          nn = nn + 1
  49          goto 101
  50  111     continue
  51       close(21)
  52 *
  53       goto 888
  54  90   PRINT*, 'input - output error'
  55  888  continue
  56       RETURN
  57       END
  58
  59
  60
  61
  62
  63
  64
  65
  66
  67
  68