908cd281 |
1 | SUBROUTINE ELOSS_READ |
638e45d7 |
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 |
638e45d7 |
7 | CHROOT=' ' |
8 | CALL GETENVF('ALICE_ROOT',CHROOT) |
9 | LNROOT = LNBLNK(CHROOT) |
10 | |
11 | IF(LNROOT.LE.0) THEN |
12 | FILNAM='DATAc.dat' |
13 | ELSE |
14 | FILNAM=CHROOT(1:LNROOT)//'/EVGEN/PartEloss/DATAc.dat' |
15 | ENDIF |
16 | |
17 | OPEN(UNIT=20,FILE=FILNAM,STATUS='OLD',ERR=90) |
18 | |
19 | nn = 1 |
20 | 100 read (20,*,end=110) xx(nn), ca(1,nn), ca(2,nn), ca(3,nn), |
21 | + ca(4,nn), ca(5,nn), ca(6,nn), ca(7,nn), ca(8,nn), |
22 | + ca(9,nn), ca(10,nn), ca(11,nn), ca(12,nn), ca(13,nn), |
23 | + ca(14,nn), ca(15,nn), ca(16,nn), ca(17,nn), ca(18,nn), |
24 | + ca(19,nn), ca(20,nn), ca(21,nn), ca(22,nn), ca(23,nn), |
25 | + ca(24,nn), ca(25,nn), ca(26,nn), ca(27,nn), ca(28,nn), |
26 | + ca(29,nn), ca(30,nn) |
27 | * print*, 0.005*(nn-1), ca(1,nn), ca(2,nn), ca(3,nn), |
28 | * + ca(4,nn), ca(5,nn), ca(6,nn), ca(7,nn), ca(8,nn), |
29 | * + ca(9,nn), ca(10,nn), ca(11,nn), ca(12,nn), ca(13,nn), |
30 | * + ca(14,nn), ca(15,nn), ca(16,nn), ca(17,nn), ca(18,nn), |
31 | * + ca(19,nn), ca(20,nn), ca(21,nn), ca(22,nn), ca(23,nn), |
32 | * + ca(24,nn), ca(25,nn), ca(26,nn), ca(27,nn), ca(28,nn), |
33 | * + ca(29,nn), ca(30,nn) |
34 | nn = nn + 1 |
35 | goto 100 |
36 | 110 continue |
37 | close(20) |
38 | * |
39 | IF(LNROOT.LE.0) THEN |
40 | FILNAM='DATAd.dat' |
41 | ELSE |
42 | FILNAM=CHROOT(1:LNROOT)//'/EVGEN/PartEloss/DATAd.dat' |
43 | ENDIF |
44 | |
45 | OPEN(UNIT=21,FILE=FILNAM,STATUS='OLD',ERR=90) |
46 | nn = 1 |
47 | 101 read (21,*,end=111) rrr(nn), da(nn) |
48 | * print*, rrr(nn), da(nn) |
49 | nn = nn + 1 |
50 | goto 101 |
51 | 111 continue |
52 | close(21) |
53 | * |
54 | goto 888 |
55 | 90 PRINT*, 'input - output error' |
56 | 888 continue |
57 | RETURN |
58 | END |
59 | |