1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
|
subroutine addevs(tevts,evtspt,nevts,pointi,t,evtnb,ierr)
C
c Copyright INRIA
double precision tevts(nevts)
integer evtspt(nevts)
integer nevts
integer pointi
double precision t
integer evtnb
integer ierr
C
C.. Local Scalars ..
integer i,j
C
C
ierr = 0
if (evtspt(evtnb).ne.-1) then
ierr=1
return
else
evtspt(evtnb)=0
tevts(evtnb)=t
endif
if(pointi.eq.0) then
pointi=evtnb
return
endif
if (t.le.tevts(pointi)) then
evtspt(evtnb)=pointi
pointi=evtnb
return
endif
i = pointi
itest=0
C
100 itest=itest+1
if(itest.gt.nevts) then
ierr=1
return
endif
if(evtspt(i).eq.0) then
evtspt(i)=evtnb
return
endif
if (t.gt.tevts(evtspt(i))) then
j = evtspt(i)
if(evtspt(j).eq.0) then
evtspt(j)=evtnb
return
endif
i=j
goto 100
else
evtspt(evtnb)=evtspt(i)
evtspt(i)=evtnb
endif
end
subroutine putevs(tevts,evtspt,nevts,pointi,t,evtnb,ierr)
C
double precision tevts(nevts)
integer evtspt(nevts)
integer nevts
integer pointi
double precision t
integer evtnb
integer ierr
C
C.. Local Scalars ..
integer i,j
C
C
ierr = 0
if (evtspt(evtnb).ne.-1) then
ierr=1
return
else
evtspt(evtnb)=0
tevts(evtnb)=t
endif
if(pointi.eq.0) then
pointi=evtnb
return
endif
evtspt(evtnb)=pointi
pointi=evtnb
return
end
|
