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
|
$NetBSD: patch-ab,v 1.1 2000/05/13 14:51:12 dmcmahill Exp $
"D" is not recognized by f2c as a comment.
--- f_source/harmd.f.orig Tue Oct 20 12:42:44 1998
+++ f_source/harmd.f Fri Apr 21 15:00:54 2000
@@ -98,5 +98,5 @@
DIMENSION INV(*),S(*)
EQUIVALENCE (N1,N(1)),(N2,N(2)),(N3,N(3))
-D DOUBLE PRECISION TIME1,TIME2,TIME3,TIME4,TIME5,TIME6,TIME7
+CD DOUBLE PRECISION TIME1,TIME2,TIME3,TIME4,TIME5,TIME6,TIME7
869 FORMAT(1H ,E13.4)
879 FORMAT(1H ,6I5)
@@ -148,5 +148,5 @@
40 KBIT=KBIT/2
KL=KBIT-2
-D CALL CLOCK(TIME1)
+CD CALL CLOCK(TIME1)
DO 50 I=1,IL1,IDIF
KLAST=KL+I
@@ -193,5 +193,5 @@
C
C DO FOR J=0
-D CALL CLOCK(TIME1)
+CD CALL CLOCK(TIME1)
DO 80 I=1,IL1,IDIF
KLAST=I+KL
@@ -256,7 +256,7 @@
A(K3+1)=AK2_1-AK3_0
80 continue
-D CALL CLOCK(TIME2)
-D WRITE(6,2000)TIME2-TIME1
-D2000 FORMAT(' AT 2 ',E12.5)
+CD CALL CLOCK(TIME2)
+CD WRITE(6,2000)TIME2-TIME1
+CD2000 FORMAT(' AT 2 ',E12.5)
IF (JLAST) 235,235,82
82 JJ=JJDIF +1
@@ -264,5 +264,5 @@
C DO FOR J=1
ILAST= IL +JJ
-D CALL CLOCK(TIME1)
+CD CALL CLOCK(TIME1)
DO 85 I = JJ,ILAST,IDIF
KLAST = KL+I
@@ -323,12 +323,12 @@
A(K3+1)=AK2_1-AK3_0
85 continue
-D CALL CLOCK(TIME2)
-D WRITE(6,3000)TIME2-TIME1
-D3000 FORMAT(' AT 3 ',E12.5)
+CD CALL CLOCK(TIME2)
+CD WRITE(6,3000)TIME2-TIME1
+CD3000 FORMAT(' AT 3 ',E12.5)
IF(JLAST-1) 235,235,90
90 JJ= JJ + JJDIF
C
C NOW DO THE REMAINING J@S
-D CALL CLOCK(TIME1)
+CD CALL CLOCK(TIME1)
DO 230 J=2,JLAST
C THIS EQUIVALENCE FORCES SOME VARIABLES INTO PERMANENT REGISTER
@@ -447,7 +447,7 @@
C
230 JJ=JJDIF+JJ
-D CALL CLOCK(TIME2)
-D WRITE(6,4000)TIME2-TIME1
-D4000 FORMAT(' AT 4 ',E12.5)
+CD CALL CLOCK(TIME2)
+CD WRITE(6,4000)TIME2-TIME1
+CD4000 FORMAT(' AT 4 ',E12.5)
C END OF J-LOOP
C
|
