1      SUBROUTINE eomccsdtq_density1(d_d1,d_i0,d_t1,d_t2,d_t3,d_t4,d_x0,d
2     &_x1,d_x2,d_x3,d_x4,d_y0,d_y1,d_y2,d_y3,d_y4,k_d1_offset,k_i0_offse
3     &t,k_t1_offset,k_t2_offset,k_t3_offset,k_t4_offset,k_x0_offset,k_x1
4     &_offset,k_x2_offset,k_x3_offset,k_x4_offset,k_y0_offset,k_y1_offse
5     &t,k_y2_offset,k_y3_offset,k_y4_offset)
6C     $Id$
7C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
8C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
9C     i0 ( )_yxd + = 1 * Sum ( h2 p1 ) * d ( p1 h2 )_d * i1 ( h2 p1 )_yx
10C         i1 ( h2 p1 )_yx + = 1 * x ( )_x * y ( h2 p1 )_y
11C         i1 ( h2 p1 )_yx + = 1 * Sum ( h4 p3 ) * x ( p3 h4 )_x * y ( h2 h4 p1 p3 )_y
12C         i1 ( h2 p1 )_yx + = 1/4 * Sum ( h6 h5 p4 p3 ) * x ( p3 p4 h5 h6 )_x * y ( h2 h5 h6 p1 p3 p4 )_y
13C         i1 ( h2 p1 )_yx + = 1/36 * Sum ( h8 h7 h6 p5 p4 p3 ) * x ( p3 p4 p5 h6 h7 h8 )_x * y ( h2 h6 h7 h8 p1 p3 p4 p5 )_y
14C     i0 ( )_dxy + = 1 * y ( )_y * i1 ( )_dx
15C         i1 ( )_dx + = 1 * Sum ( h1 p2 ) * d ( h1 p2 )_d * x ( p2 h1 )_x
16C         i1 ( )_dtx + = 1 * x ( )_x * i2 ( )_dt
17C             i2 ( )_dt + = 1 * Sum ( h1 p2 ) * d ( h1 p2 )_d * t ( p2 h1 )_t
18C     i0 ( )_yxd + = -1 * Sum ( h2 h1 ) * d ( h1 h2 )_d * i1 ( h2 h1 )_yx
19C         i1 ( h2 h1 )_yx + = 1 * Sum ( p3 ) * x ( p3 h1 )_x * y ( h2 p3 )_y
20C         i1 ( h2 h1 )_yx + = 1/2 * Sum ( h5 p4 p3 ) * x ( p3 p4 h1 h5 )_x * y ( h2 h5 p3 p4 )_y
21C         i1 ( h2 h1 )_yx + = 1/12 * Sum ( h7 h6 p5 p4 p3 ) * x ( p3 p4 p5 h1 h6 h7 )_x * y ( h2 h6 h7 p3 p4 p5 )_y
22C         i1 ( h2 h1 )_yx + = 1/144 * Sum ( h9 h8 h7 p6 p5 p4 p3 ) * x ( p3 p4 p5 p6 h1 h7 h8 h9 )_x * y ( h2 h7 h8 h9 p3 p4 p5 p6 )_y
23C         i1 ( h2 h1 )_ytx + = 1 * x ( )_x * i2 ( h2 h1 )_yt
24C             i2 ( h2 h1 )_yt + = 1 * Sum ( p3 ) * t ( p3 h1 )_t * y ( h2 p3 )_y
25C             i2 ( h2 h1 )_yt + = 1/2 * Sum ( h5 p4 p3 ) * t ( p3 p4 h1 h5 )_t * y ( h2 h5 p3 p4 )_y
26C             i2 ( h2 h1 )_yt + = 1/12 * Sum ( h7 h6 p5 p4 p3 ) * t ( p3 p4 p5 h1 h6 h7 )_t * y ( h2 h6 h7 p3 p4 p5 )_y
27C             i2 ( h2 h1 )_yt + = 1/144 * Sum ( h9 h8 h7 p6 p5 p4 p3 ) * t ( p3 p4 p5 p6 h1 h7 h8 h9 )_t * y ( h2 h7 h8 h9 p3 p4 p5 p6 )_y
28C         i1 ( h2 h1 )_yxt + = 1 * Sum ( p3 ) * t ( p3 h1 )_t * i2 ( h2 p3 )_yx
29C             i2 ( h2 p3 )_yx + = 1 * Sum ( h5 p4 ) * x ( p4 h5 )_x * y ( h2 h5 p3 p4 )_y
30C             i2 ( h2 p3 )_yx + = 1/4 * Sum ( h7 h6 p5 p4 ) * x ( p4 p5 h6 h7 )_x * y ( h2 h6 h7 p3 p4 p5 )_y
31C             i2 ( h2 p3 )_yx + = 1/36 * Sum ( h9 h8 h7 p6 p5 p4 ) * x ( p4 p5 p6 h7 h8 h9 )_x * y ( h2 h7 h8 h9 p3 p4 p5 p6 )_y
32C         i1 ( h2 h1 )_yxt + = 1/2 * Sum ( h5 p3 p4 ) * t ( p3 p4 h1 h5 )_t * i2 ( h2 h5 p3 p4 )_yx
33C             i2 ( h2 h5 p3 p4 )_yx + = 1 * Sum ( h7 p6 ) * x ( p6 h7 )_x * y ( h2 h5 h7 p3 p4 p6 )_y
34C             i2 ( h2 h5 p3 p4 )_yx + = 1/4 * Sum ( h9 h8 p7 p6 ) * x ( p6 p7 h8 h9 )_x * y ( h2 h5 h8 h9 p3 p4 p6 p7 )_y
35C         i1 ( h2 h1 )_yxt + = 1/12 * Sum ( h6 h7 p3 p4 p5 ) * t ( p3 p4 p5 h1 h6 h7 )_t * i2 ( h2 h6 h7 p3 p4 p5 )_yx
36C             i2 ( h2 h6 h7 p3 p4 p5 )_yx + = 1 * Sum ( h9 p8 ) * x ( p8 h9 )_x * y ( h2 h6 h7 h9 p3 p4 p5 p8 )_y
37C     i0 ( )_dxy + = 1 * Sum ( p1 h3 ) * y ( h3 p1 )_y * i1 ( p1 h3 )_dx
38C         i1 ( p1 h3 )_dx + = 1 * Sum ( p2 ) * d ( p1 p2 )_d * x ( p2 h3 )_x
39C         i1 ( p1 h3 )_dtx + = 1 * x ( )_x * i2 ( p1 h3 )_dt
40C             i2 ( p1 h3 )_dt + = 1 * Sum ( p2 ) * d ( p1 p2 )_d * t ( p2 h3 )_t
41C     i0 ( )_yxd + = 1 * Sum ( p12 h11 ) * d ( h11 p12 )_d * i1 ( p12 h11 )_yx
42C         i1 ( p12 h11 )_yx + = 1 * Sum ( h4 p3 ) * x ( p3 p12 h4 h11 )_x * y ( h4 p3 )_y
43C         i1 ( p12 h11 )_yx + = 1/4 * Sum ( h6 h5 p4 p3 ) * x ( p3 p4 p12 h5 h6 h11 )_x * y ( h5 h6 p3 p4 )_y
44C         i1 ( p12 h11 )_yx + = 1/36 * Sum ( h8 h7 h6 p5 p4 p3 ) * x ( p3 p4 p5 p12 h6 h7 h8 h11 )_x * y ( h6 h7 h8 p3 p4 p5 )_y
45C         i1 ( p12 h11 )_yxt + = -1 * Sum ( h2 ) * t ( p12 h2 )_t * i2 ( h2 h11 )_yx
46C             i2 ( h2 h11 )_yx + = 1 * Sum ( p4 ) * x ( p4 h11 )_x * y ( h2 p4 )_y
47C             i2 ( h2 h11 )_yx + = -1/2 * Sum ( h6 p5 p4 ) * x ( p4 p5 h6 h11 )_x * y ( h2 h6 p4 p5 )_y
48C             i2 ( h2 h11 )_yx + = 1/12 * Sum ( h8 h7 p6 p5 p4 ) * x ( p4 p5 p6 h7 h8 h11 )_x * y ( h2 h7 h8 p4 p5 p6 )_y
49C             i2 ( h2 h11 )_yx + = -1/144 * Sum ( h10 h9 h8 p7 p6 p5 p4 ) * x ( p4 p5 p6 p7 h8 h9 h10 h11 )_x * y ( h2 h8 h9 h10 p4 p5 p6 p7 )_y
50C             i2 ( h2 h11 )_ytx + = 1 * x ( )_x * i3 ( h2 h11 )_yt
51C                 i3 ( h2 h11 )_yt + = 1 * Sum ( p3 ) * t ( p3 h11 )_t * y ( h2 p3 )_y
52C                 i3 ( h2 h11 )_yt + = -1/2 * Sum ( h5 p4 p3 ) * t ( p3 p4 h5 h11 )_t * y ( h2 h5 p3 p4 )_y
53C                 i3 ( h2 h11 )_yt + = 1/12 * Sum ( h7 h6 p5 p4 p3 ) * t ( p3 p4 p5 h6 h7 h11 )_t * y ( h2 h6 h7 p3 p4 p5 )_y
54C                 i3 ( h2 h11 )_yt + = -1/144 * Sum ( h9 h8 h7 p6 p5 p4 p3 ) * t ( p3 p4 p5 p6 h7 h8 h9 h11 )_t * y ( h2 h7 h8 h9 p3 p4 p5 p6 )_y
55C             i2 ( h2 h11 )_yxt + = 1 * Sum ( p3 ) * t ( p3 h11 )_t * i3 ( h2 p3 )_yx
56C                 i3 ( h2 p3 )_yx + = 1 * Sum ( h6 p5 ) * x ( p5 h6 )_x * y ( h2 h6 p3 p5 )_y
57C                 i3 ( h2 p3 )_yx + = 1/4 * Sum ( h8 h7 p6 p5 ) * x ( p5 p6 h7 h8 )_x * y ( h2 h7 h8 p3 p5 p6 )_y
58C                 i3 ( h2 p3 )_yx + = 1/36 * Sum ( h10 h9 h8 p7 p6 p5 ) * x ( p5 p6 p7 h8 h9 h10 )_x * y ( h2 h8 h9 h10 p3 p5 p6 p7 )_y
59C             i2 ( h2 h11 )_yxt + = 1/2 * Sum ( h5 p3 p4 ) * t ( p3 p4 h5 h11 )_t * i3 ( h2 h5 p3 p4 )_yx
60C                 i3 ( h2 h5 p3 p4 )_yx + = -1 * Sum ( h8 p7 ) * x ( p7 h8 )_x * y ( h2 h5 h8 p3 p4 p7 )_y
61C                 i3 ( h2 h5 p3 p4 )_yx + = -1/4 * Sum ( h10 h9 p8 p7 ) * x ( p7 p8 h9 h10 )_x * y ( h2 h5 h9 h10 p3 p4 p7 p8 )_y
62C             i2 ( h2 h11 )_yxt + = 1/12 * Sum ( h6 h7 p3 p4 p5 ) * t ( p3 p4 p5 h6 h7 h11 )_t * i3 ( h2 h6 h7 p3 p4 p5 )_yx
63C                 i3 ( h2 h6 h7 p3 p4 p5 )_yx + = 1 * Sum ( h10 p9 ) * x ( p9 h10 )_x * y ( h2 h6 h7 h10 p3 p4 p5 p9 )_y
64C         i1 ( p12 h11 )_ytx + = -1 * Sum ( h1 ) * x ( p12 h1 )_x * i2 ( h1 h11 )_yt
65C             i2 ( h1 h11 )_yt + = 1 * Sum ( p3 ) * t ( p3 h11 )_t * y ( h1 p3 )_y
66C             i2 ( h1 h11 )_yt + = -1/2 * Sum ( h5 p4 p3 ) * t ( p3 p4 h5 h11 )_t * y ( h1 h5 p3 p4 )_y
67C             i2 ( h1 h11 )_yt + = 1/12 * Sum ( h7 h6 p5 p4 p3 ) * t ( p3 p4 p5 h6 h7 h11 )_t * y ( h1 h6 h7 p3 p4 p5 )_y
68C             i2 ( h1 h11 )_yt + = -1/144 * Sum ( h9 h8 h7 p6 p5 p4 p3 ) * t ( p3 p4 p5 p6 h7 h8 h9 h11 )_t * y ( h1 h7 h8 h9 p3 p4 p5 p6 )_y
69C         i1 ( p12 h11 )_yxt + = 1 * t ( p12 h11 )_t * i2 ( )_yx
70C             i2 ( )_yx + = 1 * Sum ( h4 p3 ) * x ( p3 h4 )_x * y ( h4 p3 )_y
71C             i2 ( )_yx + = 1/4 * Sum ( h6 h5 p4 p3 ) * x ( p3 p4 h5 h6 )_x * y ( h5 h6 p3 p4 )_y
72C             i2 ( )_yx + = 1/36 * Sum ( h8 h7 h6 p5 p4 p3 ) * x ( p3 p4 p5 h6 h7 h8 )_x * y ( h6 h7 h8 p3 p4 p5 )_y
73C             i2 ( )_yx + = 1/576 * Sum ( h10 h9 h8 h7 p6 p5 p4 p3 ) * x ( p3 p4 p5 p6 h7 h8 h9 h10 )_x * y ( h7 h8 h9 h10 p3 p4 p5 p6 )_y
74C         i1 ( p12 h11 )_ytx + = -1/2 * Sum ( h1 h2 p13 ) * x ( p12 p13 h1 h2 )_x * i2 ( h1 h2 h11 p13 )_yt
75C             i2 ( h1 h2 h11 p13 )_yt + = 1 * Sum ( p3 ) * t ( p3 h11 )_t * y ( h1 h2 p3 p13 )_y
76C             i2 ( h1 h2 h11 p13 )_yt + = 1/2 * Sum ( h5 p4 p3 ) * t ( p3 p4 h5 h11 )_t * y ( h1 h2 h5 p3 p4 p13 )_y
77C             i2 ( h1 h2 h11 p13 )_yt + = 1/12 * Sum ( h7 h6 p5 p4 p3 ) * t ( p3 p4 p5 h6 h7 h11 )_t * y ( h1 h2 h6 h7 p3 p4 p5 p13 )_y
78C         i1 ( p12 h11 )_ytx + = -1/12 * Sum ( h13 h6 h7 p1 p2 ) * x ( p1 p2 p12 h6 h7 h13 )_x * i2 ( h6 h7 h13 h11 p1 p2 )_yt
79C             i2 ( h6 h7 h13 h11 p1 p2 )_yt + = 1 * Sum ( p3 ) * t ( p3 h11 )_t * y ( h6 h7 h13 p1 p2 p3 )_y
80C             i2 ( h6 h7 h13 h11 p1 p2 )_yt + = 1/2 * Sum ( h5 p4 p3 ) * t ( p3 p4 h5 h11 )_t * y ( h5 h6 h7 h13 p1 p2 p3 p4 )_y
81C         i1 ( p12 h11 )_ytx + = 1/144 * Sum ( h7 h8 h9 h10 p4 p5 p6 ) * x ( p4 p5 p6 p12 h7 h8 h9 h10 )_x * i2 ( h7 h8 h9 h10 h11 p4 p5 p6 )_yt
82C             i2 ( h7 h8 h9 h10 h11 p4 p5 p6 )_yt + = 1 * Sum ( p3 ) * t ( p3 h11 )_t * y ( h7 h8 h9 h10 p3 p4 p5 p6 )_y
83C         i1 ( p12 h11 )_yxt + = 1 * Sum ( h4 p3 ) * t ( p3 p12 h4 h11 )_t * i2 ( h4 p3 )_yx
84C             i2 ( h4 p3 )_yx + = 1 * x ( )_x * y ( h4 p3 )_y
85C             i2 ( h4 p3 )_yx + = 1 * Sum ( h6 p5 ) * x ( p5 h6 )_x * y ( h4 h6 p3 p5 )_y
86C             i2 ( h4 p3 )_yx + = 1/4 * Sum ( h8 h7 p6 p5 ) * x ( p5 p6 h7 h8 )_x * y ( h4 h7 h8 p3 p5 p6 )_y
87C             i2 ( h4 p3 )_yx + = 1/36 * Sum ( h10 h9 h8 p7 p6 p5 ) * x ( p5 p6 p7 h8 h9 h10 )_x * y ( h4 h8 h9 h10 p3 p5 p6 p7 )_y
88C         i1 ( p12 h11 )_yxt + = -1/2 * Sum ( h2 h1 p13 ) * t ( p12 p13 h1 h2 )_t * i2 ( h1 h2 h11 p13 )_yx
89C             i2 ( h1 h2 h11 p13 )_yx + = 1 * Sum ( p6 ) * x ( p6 h11 )_x * y ( h1 h2 p6 p13 )_y
90C             i2 ( h1 h2 h11 p13 )_yx + = 1/2 * Sum ( h8 p7 p6 ) * x ( p6 p7 h8 h11 )_x * y ( h1 h2 h8 p6 p7 p13 )_y
91C             i2 ( h1 h2 h11 p13 )_yx + = 1/12 * Sum ( h10 h9 p8 p7 p6 ) * x ( p6 p7 p8 h9 h10 h11 )_x * y ( h1 h2 h9 h10 p6 p7 p8 p13 )_y
92C             i2 ( h1 h2 h11 p13 )_yxt + = -5/18 * Sum ( h7 p5 p6 ) * t ( p5 p6 h7 h11 )_t * i3 ( h1 h2 h7 p5 p6 p13 )_yx
93C                 i3 ( h1 h2 h7 p5 p6 p13 )_yx + = 1 * Sum ( h10 p9 ) * x ( p9 h10 )_x * i4 ( h1 h2 h7 h10 p5 p6 p9 p13 )_y
94C                     i4 ( h1 h2 h7 h10 p5 p6 p9 p13 )_y + = 1 * y ( h1 h2 h7 h10 p5 p6 p9 p13 )_y
95C                     i4 ( h1 h2 h7 h10 p5 p6 p9 p13 )_y + = 4/5 * y ( h1 h2 h7 h10 p5 p6 p9 p13 )_y
96C         i1 ( p12 h11 )_ytx + = 1/4 * x ( )_x * i2 ( p12 h11 )_yt
97C             i2 ( p12 h11 )_yt + = 1 * Sum ( h6 h5 p4 p3 ) * t ( p3 p4 p12 h5 h6 h11 )_t * y ( h5 h6 p3 p4 )_y
98C             i2 ( p12 h11 )_yt + = 1/9 * Sum ( h8 h7 h6 p5 p4 p3 ) * t ( p3 p4 p5 p12 h6 h7 h8 h11 )_t * y ( h6 h7 h8 p3 p4 p5 )_y
99C             i2 ( p12 h11 )_ytt + = 2 * Sum ( h13 h2 p8 ) * t ( p8 p12 h2 h13 )_t * i3 ( h2 h13 h11 p8 )_yt
100C                 i3 ( h2 h13 h11 p8 )_yt + = 1 * Sum ( p6 ) * t ( p6 h11 )_t * y ( h2 h13 p6 p8 )_y
101C                 i3 ( h2 h13 h11 p8 )_yt + = -5/18 * Sum ( h7 p6 p5 ) * t ( p5 p6 h7 h11 )_t * i4 ( h2 h7 h13 p5 p6 p8 )_y
102C                     i4 ( h2 h7 h13 p5 p6 p8 )_y + = 1 * y ( h2 h7 h13 p5 p6 p8 )_y
103C                     i4 ( h2 h7 h13 p5 p6 p8 )_y + = 4/5 * y ( h2 h7 h13 p5 p6 p8 )_y
104C                 i3 ( h2 h13 h11 p8 )_yt + = 1/12 * Sum ( h7 h6 p5 p4 p3 ) * t ( p3 p4 p5 h6 h7 h11 )_t * y ( h2 h6 h7 h13 p3 p4 p5 p8 )_y
105C             i2 ( p12 h11 )_ytt + = -1/3 * Sum ( h5 h6 h7 p3 p4 ) * t ( p3 p4 p12 h5 h6 h7 )_t * i3 ( h5 h6 h7 h11 p3 p4 )_yt
106C                 i3 ( h5 h6 h7 h11 p3 p4 )_yt + = 1 * Sum ( p8 ) * t ( p8 h11 )_t * y ( h5 h6 h7 p3 p4 p8 )_y
107C                 i3 ( h5 h6 h7 h11 p3 p4 )_yt + = -1/2 * Sum ( h10 p9 p8 ) * t ( p8 p9 h10 h11 )_t * y ( h5 h6 h7 h10 p3 p4 p8 p9 )_y
108C             i2 ( p12 h11 )_ytt + = -1/36 * Sum ( h6 h7 h8 h9 p3 p4 p5 ) * t ( p3 p4 p5 p12 h6 h7 h8 h9 )_t * i3 ( h6 h7 h8 h9 h11 p3 p4 p5 )_yt
109C                 i3 ( h6 h7 h8 h9 h11 p3 p4 p5 )_yt + = 1 * Sum ( p10 ) * t ( p10 h11 )_t * y ( h6 h7 h8 h9 p3 p4 p5 p10 )_y
110C         i1 ( p12 h11 )_yxt + = -1/12 * Sum ( h5 h6 h7 p3 p4 ) * t ( p3 p4 p12 h5 h6 h7 )_t * i2 ( h5 h6 h7 h11 p3 p4 )_yx
111C             i2 ( h5 h6 h7 h11 p3 p4 )_yx + = 1 * Sum ( p8 ) * x ( p8 h11 )_x * y ( h5 h6 h7 p3 p4 p8 )_y
112C             i2 ( h5 h6 h7 h11 p3 p4 )_yx + = -1/2 * Sum ( h10 p9 p8 ) * x ( p8 p9 h10 h11 )_x * y ( h5 h6 h7 h10 p3 p4 p8 p9 )_y
113C         i1 ( p12 h11 )_yxt + = 1/4 * Sum ( h5 h6 p3 p4 ) * t ( p3 p4 p12 h5 h6 h11 )_t * i2 ( h5 h6 p3 p4 )_yx
114C             i2 ( h5 h6 p3 p4 )_yx + = 1 * Sum ( h8 p7 ) * x ( p7 h8 )_x * y ( h5 h6 h8 p3 p4 p7 )_y
115C             i2 ( h5 h6 p3 p4 )_yx + = 1/4 * Sum ( h10 h9 p8 p7 ) * x ( p7 p8 h9 h10 )_x * y ( h5 h6 h9 h10 p3 p4 p7 p8 )_y
116C         i1 ( p12 h11 )_yxt + = -1/144 * Sum ( h6 h7 h8 h9 p3 p4 p5 ) * t ( p3 p4 p5 p12 h6 h7 h8 h9 )_t * i2 ( h6 h7 h8 h9 h11 p3 p4 p5 )_yx
117C             i2 ( h6 h7 h8 h9 h11 p3 p4 p5 )_yx + = 1 * Sum ( p10 ) * x ( p10 h11 )_x * y ( h6 h7 h8 h9 p3 p4 p5 p10 )_y
118C         i1 ( p12 h11 )_yxt + = 1/36 * Sum ( h6 h7 h8 p3 p4 p5 ) * t ( p3 p4 p5 p12 h6 h7 h8 h11 )_t * i2 ( h6 h7 h8 p3 p4 p5 )_yx
119C             i2 ( h6 h7 h8 p3 p4 p5 )_yx + = 1 * Sum ( h10 p9 ) * x ( p9 h10 )_x * y ( h6 h7 h8 h10 p3 p4 p5 p9 )_y
120C         i1 ( p12 h11 )_yxtt + = -1/2 * Sum ( p1 ) * t ( p1 h11 )_t * i2 ( p12 p1 )_yxt
121C             i2 ( p12 p1 )_yxt + = 1 * Sum ( h4 h5 p3 ) * t ( p3 p12 h4 h5 )_t * i3 ( h4 h5 p1 p3 )_yx
122C                 i3 ( h4 h5 p1 p3 )_yx + = -1 * Sum ( h8 p7 ) * x ( p7 h8 )_x * y ( h4 h5 h8 p1 p3 p7 )_y
123C                 i3 ( h4 h5 p1 p3 )_yx + = -1/4 * Sum ( h10 h9 p8 p7 ) * x ( p7 p8 h9 h10 )_x * y ( h4 h5 h9 h10 p1 p3 p7 p8 )_y
124C             i2 ( p12 p1 )_yxt + = 1/6 * Sum ( h5 h6 h7 p3 p4 ) * t ( p3 p4 p12 h5 h6 h7 )_t * i3 ( h5 h6 h7 p1 p3 p4 )_yx
125C                 i3 ( h5 h6 h7 p1 p3 p4 )_yx + = 1 * Sum ( h10 p9 ) * x ( p9 h10 )_x * y ( h5 h6 h7 h10 p1 p3 p4 p9 )_y
126C     i0 ( )_yxd + = -1/2 * Sum ( p2 p1 ) * d ( p1 p2 )_d * i1 ( p2 p1 )_yx
127C         i1 ( p2 p1 )_yx + = -1 * Sum ( h5 h4 p3 ) * x ( p2 p3 h4 h5 )_x * y ( h4 h5 p1 p3 )_y
128C         i1 ( p2 p1 )_yx + = -1/6 * Sum ( h7 h6 h5 p4 p3 ) * x ( p2 p3 p4 h5 h6 h7 )_x * y ( h5 h6 h7 p1 p3 p4 )_y
129C         i1 ( p2 p1 )_yx + = -1/72 * Sum ( h9 h8 h7 h6 p5 p4 p3 ) * x ( p2 p3 p4 p5 h6 h7 h8 h9 )_x * y ( h6 h7 h8 h9 p1 p3 p4 p5 )_y
130C         i1 ( p2 p1 )_yxt + = -2 * Sum ( h3 ) * t ( p2 h3 )_t * i2 ( h3 p1 )_yx
131C             i2 ( h3 p1 )_yx + = 1 * Sum ( h5 p4 ) * x ( p4 h5 )_x * y ( h3 h5 p1 p4 )_y
132C             i2 ( h3 p1 )_yx + = 1/4 * Sum ( h7 h6 p5 p4 ) * x ( p4 p5 h6 h7 )_x * y ( h3 h6 h7 p1 p4 p5 )_y
133C             i2 ( h3 p1 )_yx + = 1/36 * Sum ( h9 h8 h7 p6 p5 p4 ) * x ( p4 p5 p6 h7 h8 h9 )_x * y ( h3 h7 h8 h9 p1 p4 p5 p6 )_y
134C         i1 ( p2 p1 )_ytx + = 1 * x ( )_x * i2 ( p2 p1 )_yt
135C             i2 ( p2 p1 )_yt + = -1 * Sum ( h5 h4 p3 ) * t ( p2 p3 h4 h5 )_t * y ( h4 h5 p1 p3 )_y
136C             i2 ( p2 p1 )_yt + = -1/6 * Sum ( h7 h6 h5 p4 p3 ) * t ( p2 p3 p4 h5 h6 h7 )_t * y ( h5 h6 h7 p1 p3 p4 )_y
137C             i2 ( p2 p1 )_yt + = -1/72 * Sum ( h9 h8 h7 h6 p5 p4 p3 ) * t ( p2 p3 p4 p5 h6 h7 h8 h9 )_t * y ( h6 h7 h8 h9 p1 p3 p4 p5 )_y
138C         i1 ( p2 p1 )_yxt + = -1 * Sum ( h4 h5 p3 ) * t ( p2 p3 h4 h5 )_t * i2 ( h4 h5 p1 p3 )_yx
139C             i2 ( h4 h5 p1 p3 )_yx + = 1 * Sum ( h7 p6 ) * x ( p6 h7 )_x * y ( h4 h5 h7 p1 p3 p6 )_y
140C             i2 ( h4 h5 p1 p3 )_yx + = 1/4 * Sum ( h9 h8 p7 p6 ) * x ( p6 p7 h8 h9 )_x * y ( h4 h5 h8 h9 p1 p3 p6 p7 )_y
141C         i1 ( p2 p1 )_yxt + = -1/6 * Sum ( h5 h6 h7 p3 p4 ) * t ( p2 p3 p4 h5 h6 h7 )_t * i2 ( h5 h6 h7 p1 p3 p4 )_yx
142C             i2 ( h5 h6 h7 p1 p3 p4 )_yx + = 1 * Sum ( h9 p8 ) * x ( p8 h9 )_x * y ( h5 h6 h7 h9 p1 p3 p4 p8 )_y
143      IMPLICIT NONE
144#include "global.fh"
145#include "mafdecls.fh"
146#include "util.fh"
147#include "errquit.fh"
148#include "tce.fh"
149      INTEGER d_i0
150      INTEGER k_i0_offset
151      INTEGER d_d1
152      INTEGER k_d1_offset
153      INTEGER d_i1
154      INTEGER k_i1_offset
155      INTEGER d_y0
156      INTEGER k_y0_offset
157      INTEGER d_y1
158      INTEGER k_y1_offset
159      INTEGER l_i1_offset
160      INTEGER d_x0
161      INTEGER k_x0_offset
162      INTEGER size_i1
163      INTEGER d_x1
164      INTEGER k_x1_offset
165      INTEGER d_y2
166      INTEGER k_y2_offset
167      INTEGER d_x2
168      INTEGER k_x2_offset
169      INTEGER d_y3
170      INTEGER k_y3_offset
171      INTEGER d_x3
172      INTEGER k_x3_offset
173      INTEGER d_y4
174      INTEGER k_y4_offset
175      INTEGER d_i2
176      INTEGER k_i2_offset
177      INTEGER l_i2_offset
178      INTEGER d_t1
179      INTEGER k_t1_offset
180      INTEGER size_i2
181      INTEGER d_x4
182      INTEGER k_x4_offset
183      INTEGER d_t2
184      INTEGER k_t2_offset
185      INTEGER d_t3
186      INTEGER k_t3_offset
187      INTEGER d_t4
188      INTEGER k_t4_offset
189      INTEGER d_i3
190      INTEGER k_i3_offset
191      INTEGER l_i3_offset
192      INTEGER size_i3
193      INTEGER d_i4
194      INTEGER k_i4_offset
195      INTEGER l_i4_offset
196      INTEGER size_i4
197      CHARACTER*255 filename
198c#if !defined(CCSDTQ)
199c      CALL ERRQUIT('CCSDTQ not compiled',0,CAPMIS_ERR)
200c      RETURN
201c      END
202c#else
203      CALL OFFSET_eomccsdtq_density1_1_1(l_i1_offset,k_i1_offset,size_i1
204     &)
205      CALL TCE_FILENAME('eomccsdtq_density1_1_1_i1',filename)
206      CALL CREATEFILE(filename,d_i1,size_i1)
207      CALL eomccsdtq_density1_1_1(d_x0,k_x0_offset,d_y1,k_y1_offset,d_i1
208     &,k_i1_offset)
209      CALL eomccsdtq_density1_1_2(d_x1,k_x1_offset,d_y2,k_y2_offset,d_i1
210     &,k_i1_offset)
211      CALL eomccsdtq_density1_1_3(d_x2,k_x2_offset,d_y3,k_y3_offset,d_i1
212     &,k_i1_offset)
213      CALL eomccsdtq_density1_1_4(d_x3,k_x3_offset,d_y4,k_y4_offset,d_i1
214     &,k_i1_offset)
215      CALL RECONCILEFILE(d_i1,size_i1)
216      CALL eomccsdtq_density1_1(d_d1,k_d1_offset,d_i1,k_i1_offset,d_i0,k
217     &_i0_offset)
218      CALL DELETEFILE(d_i1)
219      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('eomccsdtq_densit
220     &y1',-1,MA_ERR)
221      CALL OFFSET_eomccsdtq_density1_2_1(l_i1_offset,k_i1_offset,size_i1
222     &)
223      CALL TCE_FILENAME('eomccsdtq_density1_2_1_i1',filename)
224      CALL CREATEFILE(filename,d_i1,size_i1)
225      CALL eomccsdtq_density1_2_1(d_d1,k_d1_offset,d_x1,k_x1_offset,d_i1
226     &,k_i1_offset)
227      CALL OFFSET_eomccsdtq_density1_2_2_1(l_i2_offset,k_i2_offset,size_
228     &i2)
229      CALL TCE_FILENAME('eomccsdtq_density1_2_2_1_i2',filename)
230      CALL CREATEFILE(filename,d_i2,size_i2)
231      CALL eomccsdtq_density1_2_2_1(d_d1,k_d1_offset,d_t1,k_t1_offset,d_
232     &i2,k_i2_offset)
233      CALL RECONCILEFILE(d_i2,size_i2)
234      CALL eomccsdtq_density1_2_2(d_x0,k_x0_offset,d_i2,k_i2_offset,d_i1
235     &,k_i1_offset)
236      CALL DELETEFILE(d_i2)
237      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
238     &y1',-1,MA_ERR)
239      CALL RECONCILEFILE(d_i1,size_i1)
240      CALL eomccsdtq_density1_2(d_y0,k_y0_offset,d_i1,k_i1_offset,d_i0,k
241     &_i0_offset)
242      CALL DELETEFILE(d_i1)
243      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('eomccsdtq_densit
244     &y1',-1,MA_ERR)
245      CALL OFFSET_eomccsdtq_density1_3_1(l_i1_offset,k_i1_offset,size_i1
246     &)
247      CALL TCE_FILENAME('eomccsdtq_density1_3_1_i1',filename)
248      CALL CREATEFILE(filename,d_i1,size_i1)
249      CALL eomccsdtq_density1_3_1(d_x1,k_x1_offset,d_y1,k_y1_offset,d_i1
250     &,k_i1_offset)
251      CALL eomccsdtq_density1_3_2(d_x2,k_x2_offset,d_y2,k_y2_offset,d_i1
252     &,k_i1_offset)
253      CALL eomccsdtq_density1_3_3(d_x3,k_x3_offset,d_y3,k_y3_offset,d_i1
254     &,k_i1_offset)
255      CALL eomccsdtq_density1_3_4(d_x4,k_x4_offset,d_y4,k_y4_offset,d_i1
256     &,k_i1_offset)
257      CALL OFFSET_eomccsdtq_density1_3_5_1(l_i2_offset,k_i2_offset,size_
258     &i2)
259      CALL TCE_FILENAME('eomccsdtq_density1_3_5_1_i2',filename)
260      CALL CREATEFILE(filename,d_i2,size_i2)
261      CALL eomccsdtq_density1_3_5_1(d_t1,k_t1_offset,d_y1,k_y1_offset,d_
262     &i2,k_i2_offset)
263      CALL eomccsdtq_density1_3_5_2(d_t2,k_t2_offset,d_y2,k_y2_offset,d_
264     &i2,k_i2_offset)
265      CALL eomccsdtq_density1_3_5_3(d_t3,k_t3_offset,d_y3,k_y3_offset,d_
266     &i2,k_i2_offset)
267      CALL eomccsdtq_density1_3_5_4(d_t4,k_t4_offset,d_y4,k_y4_offset,d_
268     &i2,k_i2_offset)
269      CALL RECONCILEFILE(d_i2,size_i2)
270      CALL eomccsdtq_density1_3_5(d_x0,k_x0_offset,d_i2,k_i2_offset,d_i1
271     &,k_i1_offset)
272      CALL DELETEFILE(d_i2)
273      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
274     &y1',-1,MA_ERR)
275      CALL OFFSET_eomccsdtq_density1_3_6_1(l_i2_offset,k_i2_offset,size_
276     &i2)
277      CALL TCE_FILENAME('eomccsdtq_density1_3_6_1_i2',filename)
278      CALL CREATEFILE(filename,d_i2,size_i2)
279      CALL eomccsdtq_density1_3_6_1(d_x1,k_x1_offset,d_y2,k_y2_offset,d_
280     &i2,k_i2_offset)
281      CALL eomccsdtq_density1_3_6_2(d_x2,k_x2_offset,d_y3,k_y3_offset,d_
282     &i2,k_i2_offset)
283      CALL eomccsdtq_density1_3_6_3(d_x3,k_x3_offset,d_y4,k_y4_offset,d_
284     &i2,k_i2_offset)
285      CALL RECONCILEFILE(d_i2,size_i2)
286      CALL eomccsdtq_density1_3_6(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i1
287     &,k_i1_offset)
288      CALL DELETEFILE(d_i2)
289      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
290     &y1',-1,MA_ERR)
291      CALL OFFSET_eomccsdtq_density1_3_7_1(l_i2_offset,k_i2_offset,size_
292     &i2)
293      CALL TCE_FILENAME('eomccsdtq_density1_3_7_1_i2',filename)
294      CALL CREATEFILE(filename,d_i2,size_i2)
295      CALL eomccsdtq_density1_3_7_1(d_x1,k_x1_offset,d_y3,k_y3_offset,d_
296     &i2,k_i2_offset)
297      CALL eomccsdtq_density1_3_7_2(d_x2,k_x2_offset,d_y4,k_y4_offset,d_
298     &i2,k_i2_offset)
299      CALL RECONCILEFILE(d_i2,size_i2)
300      CALL eomccsdtq_density1_3_7(d_t2,k_t2_offset,d_i2,k_i2_offset,d_i1
301     &,k_i1_offset)
302      CALL DELETEFILE(d_i2)
303      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
304     &y1',-1,MA_ERR)
305      CALL OFFSET_eomccsdtq_density1_3_8_1(l_i2_offset,k_i2_offset,size_
306     &i2)
307      CALL TCE_FILENAME('eomccsdtq_density1_3_8_1_i2',filename)
308      CALL CREATEFILE(filename,d_i2,size_i2)
309      CALL eomccsdtq_density1_3_8_1(d_x1,k_x1_offset,d_y4,k_y4_offset,d_
310     &i2,k_i2_offset)
311      CALL RECONCILEFILE(d_i2,size_i2)
312      CALL eomccsdtq_density1_3_8(d_t3,k_t3_offset,d_i2,k_i2_offset,d_i1
313     &,k_i1_offset)
314      CALL DELETEFILE(d_i2)
315      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
316     &y1',-1,MA_ERR)
317      CALL RECONCILEFILE(d_i1,size_i1)
318      CALL eomccsdtq_density1_3(d_d1,k_d1_offset,d_i1,k_i1_offset,d_i0,k
319     &_i0_offset)
320      CALL DELETEFILE(d_i1)
321      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('eomccsdtq_densit
322     &y1',-1,MA_ERR)
323      CALL OFFSET_eomccsdtq_density1_4_1(l_i1_offset,k_i1_offset,size_i1
324     &)
325      CALL TCE_FILENAME('eomccsdtq_density1_4_1_i1',filename)
326      CALL CREATEFILE(filename,d_i1,size_i1)
327      CALL eomccsdtq_density1_4_1(d_d1,k_d1_offset,d_x1,k_x1_offset,d_i1
328     &,k_i1_offset)
329      CALL OFFSET_eomccsdtq_density1_4_2_1(l_i2_offset,k_i2_offset,size_
330     &i2)
331      CALL TCE_FILENAME('eomccsdtq_density1_4_2_1_i2',filename)
332      CALL CREATEFILE(filename,d_i2,size_i2)
333      CALL eomccsdtq_density1_4_2_1(d_d1,k_d1_offset,d_t1,k_t1_offset,d_
334     &i2,k_i2_offset)
335      CALL RECONCILEFILE(d_i2,size_i2)
336      CALL eomccsdtq_density1_4_2(d_x0,k_x0_offset,d_i2,k_i2_offset,d_i1
337     &,k_i1_offset)
338      CALL DELETEFILE(d_i2)
339      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
340     &y1',-1,MA_ERR)
341      CALL RECONCILEFILE(d_i1,size_i1)
342      CALL eomccsdtq_density1_4(d_y1,k_y1_offset,d_i1,k_i1_offset,d_i0,k
343     &_i0_offset)
344      CALL DELETEFILE(d_i1)
345      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('eomccsdtq_densit
346     &y1',-1,MA_ERR)
347      CALL OFFSET_eomccsdtq_density1_5_1(l_i1_offset,k_i1_offset,size_i1
348     &)
349      CALL TCE_FILENAME('eomccsdtq_density1_5_1_i1',filename)
350      CALL CREATEFILE(filename,d_i1,size_i1)
351      CALL eomccsdtq_density1_5_1(d_x2,k_x2_offset,d_y1,k_y1_offset,d_i1
352     &,k_i1_offset)
353      CALL eomccsdtq_density1_5_2(d_x3,k_x3_offset,d_y2,k_y2_offset,d_i1
354     &,k_i1_offset)
355      CALL eomccsdtq_density1_5_3(d_x4,k_x4_offset,d_y3,k_y3_offset,d_i1
356     &,k_i1_offset)
357      CALL OFFSET_eomccsdtq_density1_5_4_1(l_i2_offset,k_i2_offset,size_
358     &i2)
359      CALL TCE_FILENAME('eomccsdtq_density1_5_4_1_i2',filename)
360      CALL CREATEFILE(filename,d_i2,size_i2)
361      CALL eomccsdtq_density1_5_4_1(d_x1,k_x1_offset,d_y1,k_y1_offset,d_
362     &i2,k_i2_offset)
363      CALL eomccsdtq_density1_5_4_2(d_x2,k_x2_offset,d_y2,k_y2_offset,d_
364     &i2,k_i2_offset)
365      CALL eomccsdtq_density1_5_4_3(d_x3,k_x3_offset,d_y3,k_y3_offset,d_
366     &i2,k_i2_offset)
367      CALL eomccsdtq_density1_5_4_4(d_x4,k_x4_offset,d_y4,k_y4_offset,d_
368     &i2,k_i2_offset)
369      CALL OFFSET_eomccsdtq_density1_5_4_5_1(l_i3_offset,k_i3_offset,siz
370     &e_i3)
371      CALL TCE_FILENAME('eomccsdtq_density1_5_4_5_1_i3',filename)
372      CALL CREATEFILE(filename,d_i3,size_i3)
373      CALL eomccsdtq_density1_5_4_5_1(d_t1,k_t1_offset,d_y1,k_y1_offset,
374     &d_i3,k_i3_offset)
375      CALL eomccsdtq_density1_5_4_5_2(d_t2,k_t2_offset,d_y2,k_y2_offset,
376     &d_i3,k_i3_offset)
377      CALL eomccsdtq_density1_5_4_5_3(d_t3,k_t3_offset,d_y3,k_y3_offset,
378     &d_i3,k_i3_offset)
379      CALL eomccsdtq_density1_5_4_5_4(d_t4,k_t4_offset,d_y4,k_y4_offset,
380     &d_i3,k_i3_offset)
381      CALL RECONCILEFILE(d_i3,size_i3)
382      CALL eomccsdtq_density1_5_4_5(d_x0,k_x0_offset,d_i3,k_i3_offset,d_
383     &i2,k_i2_offset)
384      CALL DELETEFILE(d_i3)
385      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
386     &y1',-1,MA_ERR)
387      CALL OFFSET_eomccsdtq_density1_5_4_6_1(l_i3_offset,k_i3_offset,siz
388     &e_i3)
389      CALL TCE_FILENAME('eomccsdtq_density1_5_4_6_1_i3',filename)
390      CALL CREATEFILE(filename,d_i3,size_i3)
391      CALL eomccsdtq_density1_5_4_6_1(d_x1,k_x1_offset,d_y2,k_y2_offset,
392     &d_i3,k_i3_offset)
393      CALL eomccsdtq_density1_5_4_6_2(d_x2,k_x2_offset,d_y3,k_y3_offset,
394     &d_i3,k_i3_offset)
395      CALL eomccsdtq_density1_5_4_6_3(d_x3,k_x3_offset,d_y4,k_y4_offset,
396     &d_i3,k_i3_offset)
397      CALL RECONCILEFILE(d_i3,size_i3)
398      CALL eomccsdtq_density1_5_4_6(d_t1,k_t1_offset,d_i3,k_i3_offset,d_
399     &i2,k_i2_offset)
400      CALL DELETEFILE(d_i3)
401      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
402     &y1',-1,MA_ERR)
403      CALL OFFSET_eomccsdtq_density1_5_4_7_1(l_i3_offset,k_i3_offset,siz
404     &e_i3)
405      CALL TCE_FILENAME('eomccsdtq_density1_5_4_7_1_i3',filename)
406      CALL CREATEFILE(filename,d_i3,size_i3)
407      CALL eomccsdtq_density1_5_4_7_1(d_x1,k_x1_offset,d_y3,k_y3_offset,
408     &d_i3,k_i3_offset)
409      CALL eomccsdtq_density1_5_4_7_2(d_x2,k_x2_offset,d_y4,k_y4_offset,
410     &d_i3,k_i3_offset)
411      CALL RECONCILEFILE(d_i3,size_i3)
412      CALL eomccsdtq_density1_5_4_7(d_t2,k_t2_offset,d_i3,k_i3_offset,d_
413     &i2,k_i2_offset)
414      CALL DELETEFILE(d_i3)
415      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
416     &y1',-1,MA_ERR)
417      CALL OFFSET_eomccsdtq_density1_5_4_8_1(l_i3_offset,k_i3_offset,siz
418     &e_i3)
419      CALL TCE_FILENAME('eomccsdtq_density1_5_4_8_1_i3',filename)
420      CALL CREATEFILE(filename,d_i3,size_i3)
421      CALL eomccsdtq_density1_5_4_8_1(d_x1,k_x1_offset,d_y4,k_y4_offset,
422     &d_i3,k_i3_offset)
423      CALL RECONCILEFILE(d_i3,size_i3)
424      CALL eomccsdtq_density1_5_4_8(d_t3,k_t3_offset,d_i3,k_i3_offset,d_
425     &i2,k_i2_offset)
426      CALL DELETEFILE(d_i3)
427      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
428     &y1',-1,MA_ERR)
429      CALL RECONCILEFILE(d_i2,size_i2)
430      CALL eomccsdtq_density1_5_4(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i1
431     &,k_i1_offset)
432      CALL DELETEFILE(d_i2)
433      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
434     &y1',-1,MA_ERR)
435      CALL OFFSET_eomccsdtq_density1_5_5_1(l_i2_offset,k_i2_offset,size_
436     &i2)
437      CALL TCE_FILENAME('eomccsdtq_density1_5_5_1_i2',filename)
438      CALL CREATEFILE(filename,d_i2,size_i2)
439      CALL eomccsdtq_density1_5_5_1(d_t1,k_t1_offset,d_y1,k_y1_offset,d_
440     &i2,k_i2_offset)
441      CALL eomccsdtq_density1_5_5_2(d_t2,k_t2_offset,d_y2,k_y2_offset,d_
442     &i2,k_i2_offset)
443      CALL eomccsdtq_density1_5_5_3(d_t3,k_t3_offset,d_y3,k_y3_offset,d_
444     &i2,k_i2_offset)
445      CALL eomccsdtq_density1_5_5_4(d_t4,k_t4_offset,d_y4,k_y4_offset,d_
446     &i2,k_i2_offset)
447      CALL RECONCILEFILE(d_i2,size_i2)
448      CALL eomccsdtq_density1_5_5(d_x1,k_x1_offset,d_i2,k_i2_offset,d_i1
449     &,k_i1_offset)
450      CALL DELETEFILE(d_i2)
451      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
452     &y1',-1,MA_ERR)
453      CALL OFFSET_eomccsdtq_density1_5_6_1(l_i2_offset,k_i2_offset,size_
454     &i2)
455      CALL TCE_FILENAME('eomccsdtq_density1_5_6_1_i2',filename)
456      CALL CREATEFILE(filename,d_i2,size_i2)
457      CALL eomccsdtq_density1_5_6_1(d_x1,k_x1_offset,d_y1,k_y1_offset,d_
458     &i2,k_i2_offset)
459      CALL eomccsdtq_density1_5_6_2(d_x2,k_x2_offset,d_y2,k_y2_offset,d_
460     &i2,k_i2_offset)
461      CALL eomccsdtq_density1_5_6_3(d_x3,k_x3_offset,d_y3,k_y3_offset,d_
462     &i2,k_i2_offset)
463      CALL eomccsdtq_density1_5_6_4(d_x4,k_x4_offset,d_y4,k_y4_offset,d_
464     &i2,k_i2_offset)
465      CALL RECONCILEFILE(d_i2,size_i2)
466      CALL eomccsdtq_density1_5_6(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i1
467     &,k_i1_offset)
468      CALL DELETEFILE(d_i2)
469      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
470     &y1',-1,MA_ERR)
471      CALL OFFSET_eomccsdtq_density1_5_7_1(l_i2_offset,k_i2_offset,size_
472     &i2)
473      CALL TCE_FILENAME('eomccsdtq_density1_5_7_1_i2',filename)
474      CALL CREATEFILE(filename,d_i2,size_i2)
475      CALL eomccsdtq_density1_5_7_1(d_t1,k_t1_offset,d_y2,k_y2_offset,d_
476     &i2,k_i2_offset)
477      CALL eomccsdtq_density1_5_7_2(d_t2,k_t2_offset,d_y3,k_y3_offset,d_
478     &i2,k_i2_offset)
479      CALL eomccsdtq_density1_5_7_3(d_t3,k_t3_offset,d_y4,k_y4_offset,d_
480     &i2,k_i2_offset)
481      CALL RECONCILEFILE(d_i2,size_i2)
482      CALL eomccsdtq_density1_5_7(d_x2,k_x2_offset,d_i2,k_i2_offset,d_i1
483     &,k_i1_offset)
484      CALL DELETEFILE(d_i2)
485      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
486     &y1',-1,MA_ERR)
487      CALL OFFSET_eomccsdtq_density1_5_8_1(l_i2_offset,k_i2_offset,size_
488     &i2)
489      CALL TCE_FILENAME('eomccsdtq_density1_5_8_1_i2',filename)
490      CALL CREATEFILE(filename,d_i2,size_i2)
491      CALL eomccsdtq_density1_5_8_1(d_t1,k_t1_offset,d_y3,k_y3_offset,d_
492     &i2,k_i2_offset)
493      CALL eomccsdtq_density1_5_8_2(d_t2,k_t2_offset,d_y4,k_y4_offset,d_
494     &i2,k_i2_offset)
495      CALL RECONCILEFILE(d_i2,size_i2)
496      CALL eomccsdtq_density1_5_8(d_x3,k_x3_offset,d_i2,k_i2_offset,d_i1
497     &,k_i1_offset)
498      CALL DELETEFILE(d_i2)
499      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
500     &y1',-1,MA_ERR)
501      CALL OFFSET_eomccsdtq_density1_5_9_1(l_i2_offset,k_i2_offset,size_
502     &i2)
503      CALL TCE_FILENAME('eomccsdtq_density1_5_9_1_i2',filename)
504      CALL CREATEFILE(filename,d_i2,size_i2)
505      CALL eomccsdtq_density1_5_9_1(d_t1,k_t1_offset,d_y4,k_y4_offset,d_
506     &i2,k_i2_offset)
507      CALL RECONCILEFILE(d_i2,size_i2)
508      CALL eomccsdtq_density1_5_9(d_x4,k_x4_offset,d_i2,k_i2_offset,d_i1
509     &,k_i1_offset)
510      CALL DELETEFILE(d_i2)
511      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
512     &y1',-1,MA_ERR)
513      CALL OFFSET_eomccsdtq_density1_5_10_1(l_i2_offset,k_i2_offset,size
514     &_i2)
515      CALL TCE_FILENAME('eomccsdtq_density1_5_10_1_i2',filename)
516      CALL CREATEFILE(filename,d_i2,size_i2)
517      CALL eomccsdtq_density1_5_10_1(d_x0,k_x0_offset,d_y1,k_y1_offset,d
518     &_i2,k_i2_offset)
519      CALL eomccsdtq_density1_5_10_2(d_x1,k_x1_offset,d_y2,k_y2_offset,d
520     &_i2,k_i2_offset)
521      CALL eomccsdtq_density1_5_10_3(d_x2,k_x2_offset,d_y3,k_y3_offset,d
522     &_i2,k_i2_offset)
523      CALL eomccsdtq_density1_5_10_4(d_x3,k_x3_offset,d_y4,k_y4_offset,d
524     &_i2,k_i2_offset)
525      CALL RECONCILEFILE(d_i2,size_i2)
526      CALL eomccsdtq_density1_5_10(d_t2,k_t2_offset,d_i2,k_i2_offset,d_i
527     &1,k_i1_offset)
528      CALL DELETEFILE(d_i2)
529      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
530     &y1',-1,MA_ERR)
531      CALL OFFSET_eomccsdtq_density1_5_11_1(l_i2_offset,k_i2_offset,size
532     &_i2)
533      CALL TCE_FILENAME('eomccsdtq_density1_5_11_1_i2',filename)
534      CALL CREATEFILE(filename,d_i2,size_i2)
535      CALL eomccsdtq_density1_5_11_1(d_x1,k_x1_offset,d_y2,k_y2_offset,d
536     &_i2,k_i2_offset)
537      CALL eomccsdtq_density1_5_11_2(d_x2,k_x2_offset,d_y3,k_y3_offset,d
538     &_i2,k_i2_offset)
539      CALL eomccsdtq_density1_5_11_3(d_x3,k_x3_offset,d_y4,k_y4_offset,d
540     &_i2,k_i2_offset)
541      CALL OFFSET_eomccsdtq_density1_5_11_4_1(l_i3_offset,k_i3_offset,si
542     &ze_i3)
543      CALL TCE_FILENAME('eomccsdtq_density1_5_11_4_1_i3',filename)
544      CALL CREATEFILE(filename,d_i3,size_i3)
545      CALL OFFSET_eomccsdtq_density1_5_11_4_1_1(l_i4_offset,k_i4_offset,
546     &size_i4)
547      CALL TCE_FILENAME('eomccsdtq_density1_5_11_4_1_1_i4',filename)
548      CALL CREATEFILE(filename,d_i4,size_i4)
549      CALL eomccsdtq_density1_5_11_4_1_1(d_y4,k_y4_offset,d_i4,k_i4_offs
550     &et)
551      CALL eomccsdtq_density1_5_11_4_1_2(d_y4,k_y4_offset,d_i4,k_i4_offs
552     &et)
553      CALL RECONCILEFILE(d_i4,size_i4)
554      CALL eomccsdtq_density1_5_11_4_1(d_x1,k_x1_offset,d_i4,k_i4_offset
555     &,d_i3,k_i3_offset)
556      CALL DELETEFILE(d_i4)
557      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('eomccsdtq_densit
558     &y1',-1,MA_ERR)
559      CALL RECONCILEFILE(d_i3,size_i3)
560      CALL eomccsdtq_density1_5_11_4(d_t2,k_t2_offset,d_i3,k_i3_offset,d
561     &_i2,k_i2_offset)
562      CALL DELETEFILE(d_i3)
563      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
564     &y1',-1,MA_ERR)
565      CALL RECONCILEFILE(d_i2,size_i2)
566      CALL eomccsdtq_density1_5_11(d_t2,k_t2_offset,d_i2,k_i2_offset,d_i
567     &1,k_i1_offset)
568      CALL DELETEFILE(d_i2)
569      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
570     &y1',-1,MA_ERR)
571      CALL OFFSET_eomccsdtq_density1_5_12_1(l_i2_offset,k_i2_offset,size
572     &_i2)
573      CALL TCE_FILENAME('eomccsdtq_density1_5_12_1_i2',filename)
574      CALL CREATEFILE(filename,d_i2,size_i2)
575      CALL eomccsdtq_density1_5_12_1(d_t3,k_t3_offset,d_y2,k_y2_offset,d
576     &_i2,k_i2_offset)
577      CALL eomccsdtq_density1_5_12_2(d_t4,k_t4_offset,d_y3,k_y3_offset,d
578     &_i2,k_i2_offset)
579      CALL OFFSET_eomccsdtq_density1_5_12_3_1(l_i3_offset,k_i3_offset,si
580     &ze_i3)
581      CALL TCE_FILENAME('eomccsdtq_density1_5_12_3_1_i3',filename)
582      CALL CREATEFILE(filename,d_i3,size_i3)
583      CALL eomccsdtq_density1_5_12_3_1(d_t1,k_t1_offset,d_y2,k_y2_offset
584     &,d_i3,k_i3_offset)
585      CALL OFFSET_eomccsdtq_density1_5_12_3_2_1(l_i4_offset,k_i4_offset,
586     &size_i4)
587      CALL TCE_FILENAME('eomccsdtq_density1_5_12_3_2_1_i4',filename)
588      CALL CREATEFILE(filename,d_i4,size_i4)
589      CALL eomccsdtq_density1_5_12_3_2_1(d_y3,k_y3_offset,d_i4,k_i4_offs
590     &et)
591      CALL eomccsdtq_density1_5_12_3_2_2(d_y3,k_y3_offset,d_i4,k_i4_offs
592     &et)
593      CALL RECONCILEFILE(d_i4,size_i4)
594      CALL eomccsdtq_density1_5_12_3_2(d_t2,k_t2_offset,d_i4,k_i4_offset
595     &,d_i3,k_i3_offset)
596      CALL DELETEFILE(d_i4)
597      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('eomccsdtq_densit
598     &y1',-1,MA_ERR)
599      CALL eomccsdtq_density1_5_12_3_3(d_t3,k_t3_offset,d_y4,k_y4_offset
600     &,d_i3,k_i3_offset)
601      CALL RECONCILEFILE(d_i3,size_i3)
602      CALL eomccsdtq_density1_5_12_3(d_t2,k_t2_offset,d_i3,k_i3_offset,d
603     &_i2,k_i2_offset)
604      CALL DELETEFILE(d_i3)
605      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
606     &y1',-1,MA_ERR)
607      CALL OFFSET_eomccsdtq_density1_5_12_4_1(l_i3_offset,k_i3_offset,si
608     &ze_i3)
609      CALL TCE_FILENAME('eomccsdtq_density1_5_12_4_1_i3',filename)
610      CALL CREATEFILE(filename,d_i3,size_i3)
611      CALL eomccsdtq_density1_5_12_4_1(d_t1,k_t1_offset,d_y3,k_y3_offset
612     &,d_i3,k_i3_offset)
613      CALL eomccsdtq_density1_5_12_4_2(d_t2,k_t2_offset,d_y4,k_y4_offset
614     &,d_i3,k_i3_offset)
615      CALL RECONCILEFILE(d_i3,size_i3)
616      CALL eomccsdtq_density1_5_12_4(d_t3,k_t3_offset,d_i3,k_i3_offset,d
617     &_i2,k_i2_offset)
618      CALL DELETEFILE(d_i3)
619      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
620     &y1',-1,MA_ERR)
621      CALL OFFSET_eomccsdtq_density1_5_12_5_1(l_i3_offset,k_i3_offset,si
622     &ze_i3)
623      CALL TCE_FILENAME('eomccsdtq_density1_5_12_5_1_i3',filename)
624      CALL CREATEFILE(filename,d_i3,size_i3)
625      CALL eomccsdtq_density1_5_12_5_1(d_t1,k_t1_offset,d_y4,k_y4_offset
626     &,d_i3,k_i3_offset)
627      CALL RECONCILEFILE(d_i3,size_i3)
628      CALL eomccsdtq_density1_5_12_5(d_t4,k_t4_offset,d_i3,k_i3_offset,d
629     &_i2,k_i2_offset)
630      CALL DELETEFILE(d_i3)
631      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
632     &y1',-1,MA_ERR)
633      CALL RECONCILEFILE(d_i2,size_i2)
634      CALL eomccsdtq_density1_5_12(d_x0,k_x0_offset,d_i2,k_i2_offset,d_i
635     &1,k_i1_offset)
636      CALL DELETEFILE(d_i2)
637      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
638     &y1',-1,MA_ERR)
639      CALL OFFSET_eomccsdtq_density1_5_13_1(l_i2_offset,k_i2_offset,size
640     &_i2)
641      CALL TCE_FILENAME('eomccsdtq_density1_5_13_1_i2',filename)
642      CALL CREATEFILE(filename,d_i2,size_i2)
643      CALL eomccsdtq_density1_5_13_1(d_x1,k_x1_offset,d_y3,k_y3_offset,d
644     &_i2,k_i2_offset)
645      CALL eomccsdtq_density1_5_13_2(d_x2,k_x2_offset,d_y4,k_y4_offset,d
646     &_i2,k_i2_offset)
647      CALL RECONCILEFILE(d_i2,size_i2)
648      CALL eomccsdtq_density1_5_13(d_t3,k_t3_offset,d_i2,k_i2_offset,d_i
649     &1,k_i1_offset)
650      CALL DELETEFILE(d_i2)
651      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
652     &y1',-1,MA_ERR)
653      CALL OFFSET_eomccsdtq_density1_5_14_1(l_i2_offset,k_i2_offset,size
654     &_i2)
655      CALL TCE_FILENAME('eomccsdtq_density1_5_14_1_i2',filename)
656      CALL CREATEFILE(filename,d_i2,size_i2)
657      CALL eomccsdtq_density1_5_14_1(d_x1,k_x1_offset,d_y3,k_y3_offset,d
658     &_i2,k_i2_offset)
659      CALL eomccsdtq_density1_5_14_2(d_x2,k_x2_offset,d_y4,k_y4_offset,d
660     &_i2,k_i2_offset)
661      CALL RECONCILEFILE(d_i2,size_i2)
662      CALL eomccsdtq_density1_5_14(d_t3,k_t3_offset,d_i2,k_i2_offset,d_i
663     &1,k_i1_offset)
664      CALL DELETEFILE(d_i2)
665      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
666     &y1',-1,MA_ERR)
667      CALL OFFSET_eomccsdtq_density1_5_15_1(l_i2_offset,k_i2_offset,size
668     &_i2)
669      CALL TCE_FILENAME('eomccsdtq_density1_5_15_1_i2',filename)
670      CALL CREATEFILE(filename,d_i2,size_i2)
671      CALL eomccsdtq_density1_5_15_1(d_x1,k_x1_offset,d_y4,k_y4_offset,d
672     &_i2,k_i2_offset)
673      CALL RECONCILEFILE(d_i2,size_i2)
674      CALL eomccsdtq_density1_5_15(d_t4,k_t4_offset,d_i2,k_i2_offset,d_i
675     &1,k_i1_offset)
676      CALL DELETEFILE(d_i2)
677      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
678     &y1',-1,MA_ERR)
679      CALL OFFSET_eomccsdtq_density1_5_16_1(l_i2_offset,k_i2_offset,size
680     &_i2)
681      CALL TCE_FILENAME('eomccsdtq_density1_5_16_1_i2',filename)
682      CALL CREATEFILE(filename,d_i2,size_i2)
683      CALL eomccsdtq_density1_5_16_1(d_x1,k_x1_offset,d_y4,k_y4_offset,d
684     &_i2,k_i2_offset)
685      CALL RECONCILEFILE(d_i2,size_i2)
686      CALL eomccsdtq_density1_5_16(d_t4,k_t4_offset,d_i2,k_i2_offset,d_i
687     &1,k_i1_offset)
688      CALL DELETEFILE(d_i2)
689      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
690     &y1',-1,MA_ERR)
691      CALL OFFSET_eomccsdtq_density1_5_17_1(l_i2_offset,k_i2_offset,size
692     &_i2)
693      CALL TCE_FILENAME('eomccsdtq_density1_5_17_1_i2',filename)
694      CALL CREATEFILE(filename,d_i2,size_i2)
695      CALL OFFSET_eomccsdtq_density1_5_17_1_1(l_i3_offset,k_i3_offset,si
696     &ze_i3)
697      CALL TCE_FILENAME('eomccsdtq_density1_5_17_1_1_i3',filename)
698      CALL CREATEFILE(filename,d_i3,size_i3)
699      CALL eomccsdtq_density1_5_17_1_1(d_x1,k_x1_offset,d_y3,k_y3_offset
700     &,d_i3,k_i3_offset)
701      CALL eomccsdtq_density1_5_17_1_2(d_x2,k_x2_offset,d_y4,k_y4_offset
702     &,d_i3,k_i3_offset)
703      CALL RECONCILEFILE(d_i3,size_i3)
704      CALL eomccsdtq_density1_5_17_1(d_t2,k_t2_offset,d_i3,k_i3_offset,d
705     &_i2,k_i2_offset)
706      CALL DELETEFILE(d_i3)
707      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
708     &y1',-1,MA_ERR)
709      CALL OFFSET_eomccsdtq_density1_5_17_2_1(l_i3_offset,k_i3_offset,si
710     &ze_i3)
711      CALL TCE_FILENAME('eomccsdtq_density1_5_17_2_1_i3',filename)
712      CALL CREATEFILE(filename,d_i3,size_i3)
713      CALL eomccsdtq_density1_5_17_2_1(d_x1,k_x1_offset,d_y4,k_y4_offset
714     &,d_i3,k_i3_offset)
715      CALL RECONCILEFILE(d_i3,size_i3)
716      CALL eomccsdtq_density1_5_17_2(d_t3,k_t3_offset,d_i3,k_i3_offset,d
717     &_i2,k_i2_offset)
718      CALL DELETEFILE(d_i3)
719      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('eomccsdtq_densit
720     &y1',-1,MA_ERR)
721      CALL RECONCILEFILE(d_i2,size_i2)
722      CALL eomccsdtq_density1_5_17(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i
723     &1,k_i1_offset)
724      CALL DELETEFILE(d_i2)
725      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
726     &y1',-1,MA_ERR)
727      CALL RECONCILEFILE(d_i1,size_i1)
728      CALL eomccsdtq_density1_5(d_d1,k_d1_offset,d_i1,k_i1_offset,d_i0,k
729     &_i0_offset)
730      CALL DELETEFILE(d_i1)
731      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('eomccsdtq_densit
732     &y1',-1,MA_ERR)
733      CALL OFFSET_eomccsdtq_density1_6_1(l_i1_offset,k_i1_offset,size_i1
734     &)
735      CALL TCE_FILENAME('eomccsdtq_density1_6_1_i1',filename)
736      CALL CREATEFILE(filename,d_i1,size_i1)
737      CALL eomccsdtq_density1_6_1(d_x2,k_x2_offset,d_y2,k_y2_offset,d_i1
738     &,k_i1_offset)
739      CALL eomccsdtq_density1_6_2(d_x3,k_x3_offset,d_y3,k_y3_offset,d_i1
740     &,k_i1_offset)
741      CALL eomccsdtq_density1_6_3(d_x4,k_x4_offset,d_y4,k_y4_offset,d_i1
742     &,k_i1_offset)
743      CALL OFFSET_eomccsdtq_density1_6_4_1(l_i2_offset,k_i2_offset,size_
744     &i2)
745      CALL TCE_FILENAME('eomccsdtq_density1_6_4_1_i2',filename)
746      CALL CREATEFILE(filename,d_i2,size_i2)
747      CALL eomccsdtq_density1_6_4_1(d_x1,k_x1_offset,d_y2,k_y2_offset,d_
748     &i2,k_i2_offset)
749      CALL eomccsdtq_density1_6_4_2(d_x2,k_x2_offset,d_y3,k_y3_offset,d_
750     &i2,k_i2_offset)
751      CALL eomccsdtq_density1_6_4_3(d_x3,k_x3_offset,d_y4,k_y4_offset,d_
752     &i2,k_i2_offset)
753      CALL RECONCILEFILE(d_i2,size_i2)
754      CALL eomccsdtq_density1_6_4(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i1
755     &,k_i1_offset)
756      CALL DELETEFILE(d_i2)
757      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
758     &y1',-1,MA_ERR)
759      CALL OFFSET_eomccsdtq_density1_6_5_1(l_i2_offset,k_i2_offset,size_
760     &i2)
761      CALL TCE_FILENAME('eomccsdtq_density1_6_5_1_i2',filename)
762      CALL CREATEFILE(filename,d_i2,size_i2)
763      CALL eomccsdtq_density1_6_5_1(d_t2,k_t2_offset,d_y2,k_y2_offset,d_
764     &i2,k_i2_offset)
765      CALL eomccsdtq_density1_6_5_2(d_t3,k_t3_offset,d_y3,k_y3_offset,d_
766     &i2,k_i2_offset)
767      CALL eomccsdtq_density1_6_5_3(d_t4,k_t4_offset,d_y4,k_y4_offset,d_
768     &i2,k_i2_offset)
769      CALL RECONCILEFILE(d_i2,size_i2)
770      CALL eomccsdtq_density1_6_5(d_x0,k_x0_offset,d_i2,k_i2_offset,d_i1
771     &,k_i1_offset)
772      CALL DELETEFILE(d_i2)
773      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
774     &y1',-1,MA_ERR)
775      CALL OFFSET_eomccsdtq_density1_6_6_1(l_i2_offset,k_i2_offset,size_
776     &i2)
777      CALL TCE_FILENAME('eomccsdtq_density1_6_6_1_i2',filename)
778      CALL CREATEFILE(filename,d_i2,size_i2)
779      CALL eomccsdtq_density1_6_6_1(d_x1,k_x1_offset,d_y3,k_y3_offset,d_
780     &i2,k_i2_offset)
781      CALL eomccsdtq_density1_6_6_2(d_x2,k_x2_offset,d_y4,k_y4_offset,d_
782     &i2,k_i2_offset)
783      CALL RECONCILEFILE(d_i2,size_i2)
784      CALL eomccsdtq_density1_6_6(d_t2,k_t2_offset,d_i2,k_i2_offset,d_i1
785     &,k_i1_offset)
786      CALL DELETEFILE(d_i2)
787      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
788     &y1',-1,MA_ERR)
789      CALL OFFSET_eomccsdtq_density1_6_7_1(l_i2_offset,k_i2_offset,size_
790     &i2)
791      CALL TCE_FILENAME('eomccsdtq_density1_6_7_1_i2',filename)
792      CALL CREATEFILE(filename,d_i2,size_i2)
793      CALL eomccsdtq_density1_6_7_1(d_x1,k_x1_offset,d_y4,k_y4_offset,d_
794     &i2,k_i2_offset)
795      CALL RECONCILEFILE(d_i2,size_i2)
796      CALL eomccsdtq_density1_6_7(d_t3,k_t3_offset,d_i2,k_i2_offset,d_i1
797     &,k_i1_offset)
798      CALL DELETEFILE(d_i2)
799      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('eomccsdtq_densit
800     &y1',-1,MA_ERR)
801      CALL RECONCILEFILE(d_i1,size_i1)
802      CALL eomccsdtq_density1_6(d_d1,k_d1_offset,d_i1,k_i1_offset,d_i0,k
803     &_i0_offset)
804      CALL DELETEFILE(d_i1)
805      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('eomccsdtq_densit
806     &y1',-1,MA_ERR)
807      RETURN
808      END
809c#endif
810