35
36
37
38
39
40
41#include "implicit_f.inc"
42
43
44
45 INTEGER NUMNOD
46 INTEGER ITAB(*), ITABM1(*)
47
48
49
50 INTEGER I
51 INTEGER WORK(70000)
52 INTEGER, DIMENSION(:), ALLOCATABLE :: INDX
53
54
55
56 ALLOCATE(indx(2*numnod))
57 DO i = 1, numnod
58 indx(i) = i
59 END DO
60
62
63 IF(numnod>=1)THEN
64 itabm1(1) = itab(indx(1))
65 itabm1(numnod+1) = indx(1)
66 ENDIF
67
68 DO i = 2, numnod
69 itabm1(i) = itab(indx(i))
70 IF(itabm1(i)==itabm1(i-1))THEN
71
72
73
74
75
76 indx(i)=indx(i-1)
77 END IF
78 itabm1(numnod+i) = indx(i)
79 END DO
80 DEALLOCATE(indx)
81
82 RETURN
void my_orders(int *mode, int *iwork, int *data, int *index, int *n, int *irecl)