1      SUBROUTINE beta_1(d_i0,d_t1,d_tra1,d_tra2,d_trb1,d_trb2,d_trc1,d_t
2     &rc2,d_v2,d_y1,d_y2,k_i0_offset,k_t1_offset,k_tra1_offset,k_tra2_of
3     &fset,k_trb1_offset,k_trb2_offset,k_trc1_offset,k_trc2_offset,k_v2_
4     &offset,k_y1_offset,k_y2_offset)
5C     $Id$
6C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
7C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
8C     i0 ( )_vtratrctrby + = -1 * Sum ( p9 h8 ) * y ( h8 p9 )_y * i1 ( p9 h8 )_vtratrctrb
9C         i1 ( p9 h8 )_vtratrctrb + = 1 * Sum ( h4 ) * trb ( p9 h4 )_trb * i2 ( h4 h8 )_vtratrc
10C             i2 ( h4 h8 )_vtratrc + = 1 * Sum ( p5 ) * trc ( p5 h8 )_trc * i3 ( h4 p5 )_vtra
11C                 i3 ( h4 p5 )_vtra + = 1 * Sum ( h2 p1 ) * tra ( p1 h2 )_tra * v ( h2 h4 p1 p5 )_v
12C             i2 ( h4 h8 )_vtrctra + = 1 * Sum ( p1 ) * tra ( p1 h8 )_tra * i3 ( h4 p1 )_vtrc
13C                 i3 ( h4 p1 )_vtrc + = 1 * Sum ( h6 p5 ) * trc ( p5 h6 )_trc * v ( h4 h6 p1 p5 )_v
14C         i1 ( p9 h8 )_vtrbtrctra + = -1 * Sum ( h2 ) * tra ( p9 h2 )_tra * i2 ( h2 h8 )_vtrbtrc
15C             i2 ( h2 h8 )_vtrbtrc + = 1 * Sum ( p5 ) * trc ( p5 h8 )_trc * i3 ( h2 p5 )_vtrb
16C                 i3 ( h2 p5 )_vtrb + = 1 * Sum ( h4 p3 ) * trb ( p3 h4 )_trb * v ( h2 h4 p3 p5 )_v
17C             i2 ( h2 h8 )_vtrctrb + = -1 * Sum ( p3 ) * trb ( p3 h8 )_trb * i3 ( h2 p3 )_vtrc
18C                 i3 ( h2 p3 )_vtrc + = 1 * Sum ( h6 p5 ) * trc ( p5 h6 )_trc * v ( h2 h6 p3 p5 )_v
19C         i1 ( p9 h8 )_vtratrbtrc + = 1 * Sum ( h6 ) * trc ( p9 h6 )_trc * i2 ( h6 h8 )_vtratrb
20C             i2 ( h6 h8 )_vtratrb + = 1 * Sum ( p3 ) * trb ( p3 h8 )_trb * i3 ( h6 p3 )_vtra
21C                 i3 ( h6 p3 )_vtra + = 1 * Sum ( h2 p1 ) * tra ( p1 h2 )_tra * v ( h2 h6 p1 p3 )_v
22C             i2 ( h6 h8 )_vtrbtra + = -1 * Sum ( p1 ) * tra ( p1 h8 )_tra * i3 ( h6 p1 )_vtrb
23C                 i3 ( h6 p1 )_vtrb + = 1 * Sum ( h4 p3 ) * trb ( p3 h4 )_trb * v ( h4 h6 p1 p3 )_v
24C     i0 ( )_ytrbtratrcv + = -1 * Sum ( h7 p8 h10 h11 ) * i1 ( h7 p8 h10 h11 )_ytrbtratrc * v ( h10 h11 h7 p8 )_v
25C         i1 ( h7 p8 h10 h11 )_ytrbtratrc + = 1 * Sum ( h6 ) * trc ( p8 h6 )_trc * i2 ( h6 h7 h10 h11 )_ytrbtra
26C             i2 ( h6 h7 h10 h11 )_ytrbtra + = 1 * Sum ( p1 ) * tra ( p1 h10 )_tra * i3 ( h6 h7 h11 p1 )_ytrb
27C                 i3 ( h6 h7 h11 p1 )_ytrb + = 1 * Sum ( p3 ) * trb ( p3 h11 )_trb * y ( h6 h7 p1 p3 )_y
28C         i1 ( h7 p8 h10 h11 )_ytrctratrb + = 1 * Sum ( h4 ) * trb ( p8 h4 )_trb * i2 ( h4 h7 h10 h11 )_ytrctra
29C             i2 ( h4 h7 h10 h11 )_ytrctra + = 1 * Sum ( p1 ) * tra ( p1 h10 )_tra * i3 ( h4 h7 h11 p1 )_ytrc
30C                 i3 ( h4 h7 h11 p1 )_ytrc + = 1 * Sum ( p5 ) * trc ( p5 h11 )_trc * y ( h4 h7 p1 p5 )_y
31C         i1 ( h7 p8 h10 h11 )_ytrctrbtra + = 1 * Sum ( h2 ) * tra ( p8 h2 )_tra * i2 ( h2 h7 h10 h11 )_ytrctrb
32C             i2 ( h2 h7 h10 h11 )_ytrctrb + = 1 * Sum ( p3 ) * trb ( p3 h10 )_trb * i3 ( h2 h7 h11 p3 )_ytrc
33C                 i3 ( h2 h7 h11 p3 )_ytrc + = 1 * Sum ( p5 ) * trc ( p5 h11 )_trc * y ( h2 h7 p3 p5 )_y
34C     i0 ( )_vtrctrbtray + = -1 * Sum ( p12 p9 h10 h11 ) * y ( h10 h11 p9 p12 )_y * i1 ( p9 p12 h10 h11 )_vtrctrbtra
35C         i1 ( p9 p12 h10 h11 )_vtrctrbtra + = 1 * Sum ( h2 ) * tra ( p12 h2 )_tra * i2 ( h2 p9 h10 h11 )_vtrctrb
36C             i2 ( h2 p9 h10 h11 )_vtrctrb + = 1 * Sum ( p3 ) * trb ( p3 h10 )_trb * i3 ( h2 p9 h11 p3 )_vtrc
37C                 i3 ( h2 p9 h11 p3 )_vtrc + = -1 * Sum ( p5 ) * trc ( p5 h11 )_trc * v ( h2 p9 p3 p5 )_v
38C             i2 ( h2 p9 h10 h11 )_vtrctrb + = -1 * Sum ( h6 p4 ) * trb ( p4 p9 h6 h10 )_trb * i3 ( h2 h6 h11 p4 )_vtrc
39C                 i3 ( h2 h6 h11 p4 )_vtrc + = 1 * Sum ( p7 ) * trc ( p7 h11 )_trc * v ( h2 h6 p4 p7 )_v
40C             i2 ( h2 p9 h10 h11 )_vtrbtrc + = -1/4 * Sum ( h8 ) * trc ( p9 h8 )_trc * i3 ( h2 h8 h10 h11 )_vtrb
41C                 i3 ( h2 h8 h10 h11 )_vtrb + = -1 * Sum ( p3 p4 ) * trb ( p3 p4 h10 h11 )_trb * v ( h2 h8 p3 p4 )_v
42C             i2 ( h2 p9 h10 h11 )_vtrctrb + = -1/2 * Sum ( p4 ) * trb ( p4 p9 h10 h11 )_trb * i3 ( h2 p4 )_vtrc
43C                 i3 ( h2 p4 )_vtrc + = 1 * Sum ( h8 p7 ) * trc ( p7 h8 )_trc * v ( h2 h8 p4 p7 )_v
44C             i2 ( h2 p9 h10 h11 )_vtrbtrc + = -1 * Sum ( h8 p6 ) * trc ( p6 p9 h8 h11 )_trc * i3 ( h2 h8 h10 p6 )_vtrb
45C                 i3 ( h2 h8 h10 p6 )_vtrb + = 1 * Sum ( p3 ) * trb ( p3 h10 )_trb * v ( h2 h8 p3 p6 )_v
46C             i2 ( h2 p9 h10 h11 )_vtrctrb + = -1/4 * Sum ( h4 ) * trb ( p9 h4 )_trb * i3 ( h2 h4 h10 h11 )_vtrc
47C                 i3 ( h2 h4 h10 h11 )_vtrc + = -1 * Sum ( p5 p6 ) * trc ( p5 p6 h10 h11 )_trc * v ( h2 h4 p5 p6 )_v
48C             i2 ( h2 p9 h10 h11 )_vtrbtrc + = 1/2 * Sum ( p6 ) * trc ( p6 p9 h10 h11 )_trc * i3 ( h2 p6 )_vtrb
49C                 i3 ( h2 p6 )_vtrb + = 1 * Sum ( h4 p3 ) * trb ( p3 h4 )_trb * v ( h2 h4 p3 p6 )_v
50C             i2 ( h2 p9 h10 h11 )_vtrctrbt + = 1/2 * Sum ( h4 ) * t ( p9 h4 )_t * i3 ( h2 h4 h10 h11 )_vtrctrb
51C                 i3 ( h2 h4 h10 h11 )_vtrctrb + = 1 * Sum ( p5 ) * trb ( p5 h10 )_trb * i4 ( h2 h4 h11 p5 )_vtrc
52C                     i4 ( h2 h4 h11 p5 )_vtrc + = 1 * Sum ( p7 ) * trc ( p7 h11 )_trc * v ( h2 h4 p5 p7 )_v
53C         i1 ( p9 p12 h10 h11 )_vtrctratrb + = 1 * Sum ( h4 ) * trb ( p12 h4 )_trb * i2 ( h4 p9 h10 h11 )_vtrctra
54C             i2 ( h4 p9 h10 h11 )_vtrctra + = 1 * Sum ( p1 ) * tra ( p1 h10 )_tra * i3 ( h4 p9 h11 p1 )_vtrc
55C                 i3 ( h4 p9 h11 p1 )_vtrc + = -1 * Sum ( p5 ) * trc ( p5 h11 )_trc * v ( h4 p9 p1 p5 )_v
56C         i1 ( p9 p12 h10 h11 )_vtrbtratrc + = 1 * Sum ( h6 ) * trc ( p12 h6 )_trc * i2 ( h6 p9 h10 h11 )_vtrbtra
57C             i2 ( h6 p9 h10 h11 )_vtrbtra + = 1 * Sum ( p1 ) * tra ( p1 h10 )_tra * i3 ( h6 p9 h11 p1 )_vtrb
58C                 i3 ( h6 p9 h11 p1 )_vtrb + = -1 * Sum ( p3 ) * trb ( p3 h11 )_trb * v ( h6 p9 p1 p3 )_v
59C         i1 ( p9 p12 h10 h11 )_vtrbtratrc + = -1/2 * Sum ( h8 ) * trc ( p9 p12 h8 h11 )_trc * i2 ( h8 h10 )_vtrbtra
60C             i2 ( h8 h10 )_vtrbtra + = 1 * Sum ( p1 ) * tra ( p1 h10 )_tra * i3 ( h8 p1 )_vtrb
61C                 i3 ( h8 p1 )_vtrb + = 1 * Sum ( h4 p3 ) * trb ( p3 h4 )_trb * v ( h4 h8 p1 p3 )_v
62C         i1 ( p9 p12 h10 h11 )_vtrctrbtra + = 1/4 * Sum ( h3 h4 ) * tra ( p9 p12 h3 h4 )_tra * i2 ( h3 h4 h10 h11 )_vtrctrb
63C             i2 ( h3 h4 h10 h11 )_vtrctrb + = 1 * Sum ( p5 ) * trb ( p5 h10 )_trb * i3 ( h3 h4 h11 p5 )_vtrc
64C                 i3 ( h3 h4 h11 p5 )_vtrc + = -1 * Sum ( p7 ) * trc ( p7 h11 )_trc * v ( h3 h4 p5 p7 )_v
65C         i1 ( p9 p12 h10 h11 )_vtrbtratrc + = -1 * Sum ( h8 ) * trc ( p9 h8 )_trc * i2 ( h8 p12 h10 h11 )_vtrbtra
66C             i2 ( h8 p12 h10 h11 )_vtrbtra + = 1 * Sum ( h4 p2 ) * tra ( p2 p12 h4 h10 )_tra * i3 ( h4 h8 h11 p2 )_vtrb
67C                 i3 ( h4 h8 h11 p2 )_vtrb + = 1 * Sum ( p5 ) * trb ( p5 h11 )_trb * v ( h4 h8 p2 p5 )_v
68C             i2 ( h8 p12 h10 h11 )_vtrbtra + = 1/2 * Sum ( p2 ) * tra ( p2 p12 h10 h11 )_tra * i3 ( h8 p2 )_vtrb
69C                 i3 ( h8 p2 )_vtrb + = 1 * Sum ( h6 p5 ) * trb ( p5 h6 )_trb * v ( h6 h8 p2 p5 )_v
70C         i1 ( p9 p12 h10 h11 )_vtrctrbtra + = -1/2 * Sum ( h4 ) * tra ( p9 p12 h4 h10 )_tra * i2 ( h4 h11 )_vtrctrb
71C             i2 ( h4 h11 )_vtrctrb + = 1 * Sum ( p5 ) * trb ( p5 h11 )_trb * i3 ( h4 p5 )_vtrc
72C                 i3 ( h4 p5 )_vtrc + = 1 * Sum ( h8 p7 ) * trc ( p7 h8 )_trc * v ( h4 h8 p5 p7 )_v
73C             i2 ( h4 h11 )_vtrbtrc + = -1 * Sum ( p7 ) * trc ( p7 h11 )_trc * i3 ( h4 p7 )_vtrb
74C                 i3 ( h4 p7 )_vtrb + = 1 * Sum ( h6 p5 ) * trb ( p5 h6 )_trb * v ( h4 h6 p5 p7 )_v
75C         i1 ( p9 p12 h10 h11 )_vtrctrbtrat + = 1/2 * Sum ( h2 ) * t ( p12 h2 )_t * i2 ( h2 p9 h10 h11 )_vtrctrbtra
76C             i2 ( h2 p9 h10 h11 )_vtrctrbtra + = 1 * Sum ( h4 ) * tra ( p9 h4 )_tra * i3 ( h2 h4 h10 h11 )_vtrctrb
77C     i0 ( )_ytrbtrctrav + = -1/4 * Sum ( p11 p12 h10 h13 ) * i1 ( p11 p12 h10 h13 )_ytrbtrctra * v ( h10 h13 p11 p12 )_v
78C         i1 ( p11 p12 h10 h13 )_ytrbtrctra + = 1 * Sum ( h2 ) * tra ( p11 h2 )_tra * i2 ( h2 p12 h10 h13 )_ytrbtrc
79C             i2 ( h2 p12 h10 h13 )_ytrbtrc + = 1 * Sum ( h8 ) * trc ( p12 h8 )_trc * i3 ( h2 h8 h10 h13 )_ytrb
80C                 i3 ( h2 h8 h10 h13 )_ytrb + = -1 * Sum ( p3 p4 ) * trb ( p3 p4 h10 h13 )_trb * y ( h2 h8 p3 p4 )_y
81C                 i3 ( h2 h8 h10 h13 )_ytrbt + = -4 * Sum ( p3 ) * t ( p3 h10 )_t * i4 ( h2 h8 h13 p3 )_ytrb
82C                     i4 ( h2 h8 h13 p3 )_ytrb + = 1 * Sum ( p5 ) * trb ( p5 h13 )_trb * y ( h2 h8 p3 p5 )_y
83C             i2 ( h2 p12 h10 h13 )_ytrctrb + = 4 * Sum ( h5 p3 ) * trb ( p3 p12 h5 h10 )_trb * i3 ( h2 h5 h13 p3 )_ytrc
84C                 i3 ( h2 h5 h13 p3 )_ytrc + = 1 * Sum ( p7 ) * trc ( p7 h13 )_trc * y ( h2 h5 p3 p7 )_y
85C             i2 ( h2 p12 h10 h13 )_ytrctrb + = 1 * Sum ( h1 ) * trb ( p12 h1 )_trb * i3 ( h1 h2 h10 h13 )_ytrc
86C                 i3 ( h1 h2 h10 h13 )_ytrc + = 1 * Sum ( p5 p6 ) * trc ( p5 p6 h10 h13 )_trc * y ( h1 h2 p5 p6 )_y
87C                 i3 ( h1 h2 h10 h13 )_ytrct + = -4 * Sum ( p3 ) * t ( p3 h10 )_t * i4 ( h1 h2 h13 p3 )_ytrc
88C                     i4 ( h1 h2 h13 p3 )_ytrc + = -1 * Sum ( p7 ) * trc ( p7 h13 )_trc * y ( h1 h2 p3 p7 )_y
89C             i2 ( h2 p12 h10 h13 )_ytrbtrc + = 4 * Sum ( h7 p5 ) * trc ( p5 p12 h7 h13 )_trc * i3 ( h2 h7 h10 p5 )_ytrb
90C                 i3 ( h2 h7 h10 p5 )_ytrb + = 1 * Sum ( p3 ) * trb ( p3 h10 )_trb * y ( h2 h7 p3 p5 )_y
91C             i2 ( h2 p12 h10 h13 )_ytrctrbt + = -2 * Sum ( h4 ) * t ( p12 h4 )_t * i3 ( h2 h4 h10 h13 )_ytrctrb
92C                 i3 ( h2 h4 h10 h13 )_ytrctrb + = 1 * Sum ( p5 ) * trb ( p5 h10 )_trb * i4 ( h2 h4 h13 p5 )_ytrc
93C                     i4 ( h2 h4 h13 p5 )_ytrc + = 1 * Sum ( p7 ) * trc ( p7 h13 )_trc * y ( h2 h4 p5 p7 )_y
94C         i1 ( p11 p12 h10 h13 )_ytrbtratrc + = 2 * trc ( p12 h13 )_trc * i2 ( p11 h10 )_ytrbtra
95C             i2 ( p11 h10 )_ytrbtra + = 1 * Sum ( h2 ) * tra ( p11 h2 )_tra * i3 ( h2 h10 )_ytrb
96C                 i3 ( h2 h10 )_ytrb + = -1 * Sum ( h5 p3 p4 ) * trb ( p3 p4 h5 h10 )_trb * y ( h2 h5 p3 p4 )_y
97C             i2 ( p11 h10 )_ytrbtra + = -1 * Sum ( h3 h4 p1 ) * tra ( p1 p11 h3 h4 )_tra * i3 ( h3 h4 h10 p1 )_ytrb
98C                 i3 ( h3 h4 h10 p1 )_ytrb + = -1 * Sum ( p5 ) * trb ( p5 h10 )_trb * y ( h3 h4 p1 p5 )_y
99C         i1 ( p11 p12 h10 h13 )_ytrbtrctra + = -2 * tra ( p11 h10 )_tra * i2 ( p12 h13 )_ytrbtrc
100C             i2 ( p12 h13 )_ytrbtrc + = 1 * Sum ( h8 ) * trc ( p12 h8 )_trc * i3 ( h8 h13 )_ytrb
101C                 i3 ( h8 h13 )_ytrb + = -1 * Sum ( h5 p3 p4 ) * trb ( p3 p4 h5 h13 )_trb * y ( h5 h8 p3 p4 )_y
102C             i2 ( p12 h13 )_ytrctrb + = 1 * Sum ( h5 h6 p3 ) * trb ( p3 p12 h5 h6 )_trb * i3 ( h5 h6 h13 p3 )_ytrc
103C                 i3 ( h5 h6 h13 p3 )_ytrc + = -1 * Sum ( p7 ) * trc ( p7 h13 )_trc * y ( h5 h6 p3 p7 )_y
104C             i2 ( p12 h13 )_ytrctrb + = -1 * Sum ( h4 ) * trb ( p12 h4 )_trb * i3 ( h4 h13 )_ytrc
105C                 i3 ( h4 h13 )_ytrc + = -1 * Sum ( h7 p5 p6 ) * trc ( p5 p6 h7 h13 )_trc * y ( h4 h7 p5 p6 )_y
106C             i2 ( p12 h13 )_ytrbtrc + = -1 * Sum ( h7 h8 p5 ) * trc ( p5 p12 h7 h8 )_trc * i3 ( h7 h8 h13 p5 )_ytrb
107C                 i3 ( h7 h8 h13 p5 )_ytrb + = -1 * Sum ( p3 ) * trb ( p3 h13 )_trb * y ( h7 h8 p3 p5 )_y
108C         i1 ( p11 p12 h10 h13 )_ytrbtratrc + = 4 * Sum ( h8 ) * trc ( p12 h8 )_trc * i2 ( h8 p11 h10 h13 )_ytrbtra
109C             i2 ( h8 p11 h10 h13 )_ytrbtra + = 1 * Sum ( h3 p1 ) * tra ( p1 p11 h3 h10 )_tra * i3 ( h3 h8 h13 p1 )_ytrb
110C                 i3 ( h3 h8 h13 p1 )_ytrb + = 1 * Sum ( p5 ) * trb ( p5 h13 )_trb * y ( h3 h8 p1 p5 )_y
111C         i1 ( p11 p12 h10 h13 )_ytrctrbtra + = 1 * Sum ( h3 h4 ) * tra ( p11 p12 h3 h4 )_tra * i2 ( h3 h4 h10 h13 )_ytrctrb
112C             i2 ( h3 h4 h10 h13 )_ytrctrb + = 1 * Sum ( p5 ) * trb ( p5 h10 )_trb * i3 ( h3 h4 h13 p5 )_ytrc
113C                 i3 ( h3 h4 h13 p5 )_ytrc + = -1 * Sum ( p7 ) * trc ( p7 h13 )_trc * y ( h3 h4 p5 p7 )_y
114C         i1 ( p11 p12 h10 h13 )_ytrbtratrct + = -4 * Sum ( h2 ) * t ( p11 h2 )_t * i2 ( h2 p12 h10 h13 )_ytrbtratrc
115C             i2 ( h2 p12 h10 h13 )_ytrbtratrc + = 1 * Sum ( h8 ) * trc ( p12 h8 )_trc * i3 ( h2 h8 h10 h13 )_ytrbtra
116C                 i3 ( h2 h8 h10 h13 )_ytrbtra + = 1 * Sum ( p3 ) * tra ( p3 h10 )_tra * i4 ( h2 h8 h13 p3 )_ytrb
117C             i2 ( h2 p12 h10 h13 )_ytrctratrb + = 1 * Sum ( h6 ) * trb ( p12 h6 )_trb * i3 ( h2 h6 h10 h13 )_ytrctra
118C                 i3 ( h2 h6 h10 h13 )_ytrctra + = 1 * Sum ( p3 ) * tra ( p3 h10 )_tra * i4 ( h2 h6 h13 p3 )_ytrc
119C                     i4 ( h2 h6 h13 p3 )_ytrc + = 1 * Sum ( p7 ) * trc ( p7 h13 )_trc * y ( h2 h6 p3 p7 )_y
120C             i2 ( h2 p12 h10 h13 )_ytrctrbtra + = 1/2 * Sum ( h4 ) * tra ( p12 h4 )_tra * i3 ( h2 h4 h10 h13 )_ytrctrb
121      IMPLICIT NONE
122#include "global.fh"
123#include "mafdecls.fh"
124#include "util.fh"
125#include "errquit.fh"
126#include "tce.fh"
127      INTEGER d_i0
128      INTEGER k_i0_offset
129      INTEGER d_y1
130      INTEGER k_y1_offset
131      INTEGER d_i1
132      INTEGER k_i1_offset
133      INTEGER d_v2
134      INTEGER k_v2_offset
135      INTEGER d_y2
136      INTEGER k_y2_offset
137      INTEGER l_i1_offset
138      INTEGER d_trb1
139      INTEGER k_trb1_offset
140      INTEGER d_i2
141      INTEGER k_i2_offset
142      INTEGER size_i1
143      INTEGER d_tra1
144      INTEGER k_tra1_offset
145      INTEGER d_trc1
146      INTEGER k_trc1_offset
147      INTEGER l_i2_offset
148      INTEGER d_i3
149      INTEGER k_i3_offset
150      INTEGER size_i2
151      INTEGER l_i3_offset
152      INTEGER size_i3
153      INTEGER d_trc2
154      INTEGER k_trc2_offset
155      INTEGER d_tra2
156      INTEGER k_tra2_offset
157      INTEGER d_t1
158      INTEGER k_t1_offset
159      INTEGER d_trb2
160      INTEGER k_trb2_offset
161      INTEGER d_i4
162      INTEGER k_i4_offset
163      INTEGER l_i4_offset
164      INTEGER size_i4
165      CHARACTER*255 filename
166      CALL OFFSET_beta_1_1_1(l_i1_offset,k_i1_offset,size_i1)
167      CALL TCE_FILENAME('beta_1_1_1_i1',filename)
168      CALL CREATEFILE(filename,d_i1,size_i1)
169      CALL OFFSET_beta_1_1_1_1(l_i2_offset,k_i2_offset,size_i2)
170      CALL TCE_FILENAME('beta_1_1_1_1_i2',filename)
171      CALL CREATEFILE(filename,d_i2,size_i2)
172      CALL OFFSET_beta_1_1_1_1_1(l_i3_offset,k_i3_offset,size_i3)
173      CALL TCE_FILENAME('beta_1_1_1_1_1_i3',filename)
174      CALL CREATEFILE(filename,d_i3,size_i3)
175      CALL beta_1_1_1_1_1(d_tra1,k_tra1_offset,d_v2,k_v2_offset,d_i3,k_i
176     &3_offset)
177      CALL RECONCILEFILE(d_i3,size_i3)
178      CALL beta_1_1_1_1(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
179     &offset)
180      CALL DELETEFILE(d_i3)
181      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
182     &R)
183      CALL OFFSET_beta_1_1_1_2_1(l_i3_offset,k_i3_offset,size_i3)
184      CALL TCE_FILENAME('beta_1_1_1_2_1_i3',filename)
185      CALL CREATEFILE(filename,d_i3,size_i3)
186      CALL beta_1_1_1_2_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
187     &3_offset)
188      CALL RECONCILEFILE(d_i3,size_i3)
189      CALL beta_1_1_1_2(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
190     &offset)
191      CALL DELETEFILE(d_i3)
192      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
193     &R)
194      CALL RECONCILEFILE(d_i2,size_i2)
195      CALL beta_1_1_1(d_trb1,k_trb1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
196     &fset)
197      CALL DELETEFILE(d_i2)
198      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
199     &R)
200      CALL OFFSET_beta_1_1_2_1(l_i2_offset,k_i2_offset,size_i2)
201      CALL TCE_FILENAME('beta_1_1_2_1_i2',filename)
202      CALL CREATEFILE(filename,d_i2,size_i2)
203      CALL OFFSET_beta_1_1_2_1_1(l_i3_offset,k_i3_offset,size_i3)
204      CALL TCE_FILENAME('beta_1_1_2_1_1_i3',filename)
205      CALL CREATEFILE(filename,d_i3,size_i3)
206      CALL beta_1_1_2_1_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
207     &3_offset)
208      CALL RECONCILEFILE(d_i3,size_i3)
209      CALL beta_1_1_2_1(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
210     &offset)
211      CALL DELETEFILE(d_i3)
212      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
213     &R)
214      CALL OFFSET_beta_1_1_2_2_1(l_i3_offset,k_i3_offset,size_i3)
215      CALL TCE_FILENAME('beta_1_1_2_2_1_i3',filename)
216      CALL CREATEFILE(filename,d_i3,size_i3)
217      CALL beta_1_1_2_2_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
218     &3_offset)
219      CALL RECONCILEFILE(d_i3,size_i3)
220      CALL beta_1_1_2_2(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
221     &offset)
222      CALL DELETEFILE(d_i3)
223      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
224     &R)
225      CALL RECONCILEFILE(d_i2,size_i2)
226      CALL beta_1_1_2(d_tra1,k_tra1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
227     &fset)
228      CALL DELETEFILE(d_i2)
229      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
230     &R)
231      CALL OFFSET_beta_1_1_3_1(l_i2_offset,k_i2_offset,size_i2)
232      CALL TCE_FILENAME('beta_1_1_3_1_i2',filename)
233      CALL CREATEFILE(filename,d_i2,size_i2)
234      CALL OFFSET_beta_1_1_3_1_1(l_i3_offset,k_i3_offset,size_i3)
235      CALL TCE_FILENAME('beta_1_1_3_1_1_i3',filename)
236      CALL CREATEFILE(filename,d_i3,size_i3)
237      CALL beta_1_1_3_1_1(d_tra1,k_tra1_offset,d_v2,k_v2_offset,d_i3,k_i
238     &3_offset)
239      CALL RECONCILEFILE(d_i3,size_i3)
240      CALL beta_1_1_3_1(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
241     &offset)
242      CALL DELETEFILE(d_i3)
243      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
244     &R)
245      CALL OFFSET_beta_1_1_3_2_1(l_i3_offset,k_i3_offset,size_i3)
246      CALL TCE_FILENAME('beta_1_1_3_2_1_i3',filename)
247      CALL CREATEFILE(filename,d_i3,size_i3)
248      CALL beta_1_1_3_2_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
249     &3_offset)
250      CALL RECONCILEFILE(d_i3,size_i3)
251      CALL beta_1_1_3_2(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
252     &offset)
253      CALL DELETEFILE(d_i3)
254      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
255     &R)
256      CALL RECONCILEFILE(d_i2,size_i2)
257      CALL beta_1_1_3(d_trc1,k_trc1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
258     &fset)
259      CALL DELETEFILE(d_i2)
260      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
261     &R)
262      CALL RECONCILEFILE(d_i1,size_i1)
263      CALL beta_1_1(d_y1,k_y1_offset,d_i1,k_i1_offset,d_i0,k_i0_offset)
264      CALL DELETEFILE(d_i1)
265      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
266     &R)
267      CALL OFFSET_beta_1_2_1(l_i1_offset,k_i1_offset,size_i1)
268      CALL TCE_FILENAME('beta_1_2_1_i1',filename)
269      CALL CREATEFILE(filename,d_i1,size_i1)
270      CALL OFFSET_beta_1_2_1_1(l_i2_offset,k_i2_offset,size_i2)
271      CALL TCE_FILENAME('beta_1_2_1_1_i2',filename)
272      CALL CREATEFILE(filename,d_i2,size_i2)
273      CALL OFFSET_beta_1_2_1_1_1(l_i3_offset,k_i3_offset,size_i3)
274      CALL TCE_FILENAME('beta_1_2_1_1_1_i3',filename)
275      CALL CREATEFILE(filename,d_i3,size_i3)
276      CALL beta_1_2_1_1_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i3,k_i
277     &3_offset)
278      CALL RECONCILEFILE(d_i3,size_i3)
279      CALL beta_1_2_1_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
280     &offset)
281      CALL DELETEFILE(d_i3)
282      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
283     &R)
284      CALL RECONCILEFILE(d_i2,size_i2)
285      CALL beta_1_2_1(d_trc1,k_trc1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
286     &fset)
287      CALL DELETEFILE(d_i2)
288      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
289     &R)
290      CALL OFFSET_beta_1_2_2_1(l_i2_offset,k_i2_offset,size_i2)
291      CALL TCE_FILENAME('beta_1_2_2_1_i2',filename)
292      CALL CREATEFILE(filename,d_i2,size_i2)
293      CALL OFFSET_beta_1_2_2_1_1(l_i3_offset,k_i3_offset,size_i3)
294      CALL TCE_FILENAME('beta_1_2_2_1_1_i3',filename)
295      CALL CREATEFILE(filename,d_i3,size_i3)
296      CALL beta_1_2_2_1_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i3,k_i
297     &3_offset)
298      CALL RECONCILEFILE(d_i3,size_i3)
299      CALL beta_1_2_2_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
300     &offset)
301      CALL DELETEFILE(d_i3)
302      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
303     &R)
304      CALL RECONCILEFILE(d_i2,size_i2)
305      CALL beta_1_2_2(d_trb1,k_trb1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
306     &fset)
307      CALL DELETEFILE(d_i2)
308      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
309     &R)
310      CALL OFFSET_beta_1_2_3_1(l_i2_offset,k_i2_offset,size_i2)
311      CALL TCE_FILENAME('beta_1_2_3_1_i2',filename)
312      CALL CREATEFILE(filename,d_i2,size_i2)
313      CALL OFFSET_beta_1_2_3_1_1(l_i3_offset,k_i3_offset,size_i3)
314      CALL TCE_FILENAME('beta_1_2_3_1_1_i3',filename)
315      CALL CREATEFILE(filename,d_i3,size_i3)
316      CALL beta_1_2_3_1_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i3,k_i
317     &3_offset)
318      CALL RECONCILEFILE(d_i3,size_i3)
319      CALL beta_1_2_3_1(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
320     &offset)
321      CALL DELETEFILE(d_i3)
322      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
323     &R)
324      CALL RECONCILEFILE(d_i2,size_i2)
325      CALL beta_1_2_3(d_tra1,k_tra1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
326     &fset)
327      CALL DELETEFILE(d_i2)
328      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
329     &R)
330      CALL RECONCILEFILE(d_i1,size_i1)
331      CALL beta_1_2(d_i1,k_i1_offset,d_v2,k_v2_offset,d_i0,k_i0_offset)
332      CALL DELETEFILE(d_i1)
333      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
334     &R)
335      CALL OFFSET_beta_1_3_1(l_i1_offset,k_i1_offset,size_i1)
336      CALL TCE_FILENAME('beta_1_3_1_i1',filename)
337      CALL CREATEFILE(filename,d_i1,size_i1)
338      CALL OFFSET_beta_1_3_1_1(l_i2_offset,k_i2_offset,size_i2)
339      CALL TCE_FILENAME('beta_1_3_1_1_i2',filename)
340      CALL CREATEFILE(filename,d_i2,size_i2)
341      CALL OFFSET_beta_1_3_1_1_1(l_i3_offset,k_i3_offset,size_i3)
342      CALL TCE_FILENAME('beta_1_3_1_1_1_i3',filename)
343      CALL CREATEFILE(filename,d_i3,size_i3)
344      CALL beta_1_3_1_1_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
345     &3_offset)
346      CALL RECONCILEFILE(d_i3,size_i3)
347      CALL beta_1_3_1_1(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
348     &offset)
349      CALL DELETEFILE(d_i3)
350      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
351     &R)
352      CALL OFFSET_beta_1_3_1_2_1(l_i3_offset,k_i3_offset,size_i3)
353      CALL TCE_FILENAME('beta_1_3_1_2_1_i3',filename)
354      CALL CREATEFILE(filename,d_i3,size_i3)
355      CALL beta_1_3_1_2_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
356     &3_offset)
357      CALL RECONCILEFILE(d_i3,size_i3)
358      CALL beta_1_3_1_2(d_trb2,k_trb2_offset,d_i3,k_i3_offset,d_i2,k_i2_
359     &offset)
360      CALL DELETEFILE(d_i3)
361      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
362     &R)
363      CALL OFFSET_beta_1_3_1_3_1(l_i3_offset,k_i3_offset,size_i3)
364      CALL TCE_FILENAME('beta_1_3_1_3_1_i3',filename)
365      CALL CREATEFILE(filename,d_i3,size_i3)
366      CALL beta_1_3_1_3_1(d_trb2,k_trb2_offset,d_v2,k_v2_offset,d_i3,k_i
367     &3_offset)
368      CALL RECONCILEFILE(d_i3,size_i3)
369      CALL beta_1_3_1_3(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
370     &offset)
371      CALL DELETEFILE(d_i3)
372      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
373     &R)
374      CALL OFFSET_beta_1_3_1_4_1(l_i3_offset,k_i3_offset,size_i3)
375      CALL TCE_FILENAME('beta_1_3_1_4_1_i3',filename)
376      CALL CREATEFILE(filename,d_i3,size_i3)
377      CALL beta_1_3_1_4_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
378     &3_offset)
379      CALL RECONCILEFILE(d_i3,size_i3)
380      CALL beta_1_3_1_4(d_trb2,k_trb2_offset,d_i3,k_i3_offset,d_i2,k_i2_
381     &offset)
382      CALL DELETEFILE(d_i3)
383      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
384     &R)
385      CALL OFFSET_beta_1_3_1_5_1(l_i3_offset,k_i3_offset,size_i3)
386      CALL TCE_FILENAME('beta_1_3_1_5_1_i3',filename)
387      CALL CREATEFILE(filename,d_i3,size_i3)
388      CALL beta_1_3_1_5_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
389     &3_offset)
390      CALL RECONCILEFILE(d_i3,size_i3)
391      CALL beta_1_3_1_5(d_trc2,k_trc2_offset,d_i3,k_i3_offset,d_i2,k_i2_
392     &offset)
393      CALL DELETEFILE(d_i3)
394      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
395     &R)
396      CALL OFFSET_beta_1_3_1_6_1(l_i3_offset,k_i3_offset,size_i3)
397      CALL TCE_FILENAME('beta_1_3_1_6_1_i3',filename)
398      CALL CREATEFILE(filename,d_i3,size_i3)
399      CALL beta_1_3_1_6_1(d_trc2,k_trc2_offset,d_v2,k_v2_offset,d_i3,k_i
400     &3_offset)
401      CALL RECONCILEFILE(d_i3,size_i3)
402      CALL beta_1_3_1_6(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
403     &offset)
404      CALL DELETEFILE(d_i3)
405      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
406     &R)
407      CALL OFFSET_beta_1_3_1_7_1(l_i3_offset,k_i3_offset,size_i3)
408      CALL TCE_FILENAME('beta_1_3_1_7_1_i3',filename)
409      CALL CREATEFILE(filename,d_i3,size_i3)
410      CALL beta_1_3_1_7_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
411     &3_offset)
412      CALL RECONCILEFILE(d_i3,size_i3)
413      CALL beta_1_3_1_7(d_trc2,k_trc2_offset,d_i3,k_i3_offset,d_i2,k_i2_
414     &offset)
415      CALL DELETEFILE(d_i3)
416      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
417     &R)
418      CALL OFFSET_beta_1_3_1_8_1(l_i3_offset,k_i3_offset,size_i3)
419      CALL TCE_FILENAME('beta_1_3_1_8_1_i3',filename)
420      CALL CREATEFILE(filename,d_i3,size_i3)
421      CALL OFFSET_beta_1_3_1_8_1_1(l_i4_offset,k_i4_offset,size_i4)
422      CALL TCE_FILENAME('beta_1_3_1_8_1_1_i4',filename)
423      CALL CREATEFILE(filename,d_i4,size_i4)
424      CALL beta_1_3_1_8_1_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i4,k
425     &_i4_offset)
426      CALL RECONCILEFILE(d_i4,size_i4)
427      CALL beta_1_3_1_8_1(d_trb1,k_trb1_offset,d_i4,k_i4_offset,d_i3,k_i
428     &3_offset)
429      CALL DELETEFILE(d_i4)
430      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
431     &R)
432      CALL RECONCILEFILE(d_i3,size_i3)
433      CALL beta_1_3_1_8(d_t1,k_t1_offset,d_i3,k_i3_offset,d_i2,k_i2_offs
434     &et)
435      CALL DELETEFILE(d_i3)
436      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
437     &R)
438      CALL RECONCILEFILE(d_i2,size_i2)
439      CALL beta_1_3_1(d_tra1,k_tra1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
440     &fset)
441      CALL DELETEFILE(d_i2)
442      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
443     &R)
444      CALL OFFSET_beta_1_3_2_1(l_i2_offset,k_i2_offset,size_i2)
445      CALL TCE_FILENAME('beta_1_3_2_1_i2',filename)
446      CALL CREATEFILE(filename,d_i2,size_i2)
447      CALL OFFSET_beta_1_3_2_1_1(l_i3_offset,k_i3_offset,size_i3)
448      CALL TCE_FILENAME('beta_1_3_2_1_1_i3',filename)
449      CALL CREATEFILE(filename,d_i3,size_i3)
450      CALL beta_1_3_2_1_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
451     &3_offset)
452      CALL RECONCILEFILE(d_i3,size_i3)
453      CALL beta_1_3_2_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
454     &offset)
455      CALL DELETEFILE(d_i3)
456      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
457     &R)
458      CALL RECONCILEFILE(d_i2,size_i2)
459      CALL beta_1_3_2(d_trb1,k_trb1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
460     &fset)
461      CALL DELETEFILE(d_i2)
462      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
463     &R)
464      CALL OFFSET_beta_1_3_3_1(l_i2_offset,k_i2_offset,size_i2)
465      CALL TCE_FILENAME('beta_1_3_3_1_i2',filename)
466      CALL CREATEFILE(filename,d_i2,size_i2)
467      CALL OFFSET_beta_1_3_3_1_1(l_i3_offset,k_i3_offset,size_i3)
468      CALL TCE_FILENAME('beta_1_3_3_1_1_i3',filename)
469      CALL CREATEFILE(filename,d_i3,size_i3)
470      CALL beta_1_3_3_1_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
471     &3_offset)
472      CALL RECONCILEFILE(d_i3,size_i3)
473      CALL beta_1_3_3_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
474     &offset)
475      CALL DELETEFILE(d_i3)
476      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
477     &R)
478      CALL RECONCILEFILE(d_i2,size_i2)
479      CALL beta_1_3_3(d_trc1,k_trc1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
480     &fset)
481      CALL DELETEFILE(d_i2)
482      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
483     &R)
484      CALL OFFSET_beta_1_3_4_1(l_i2_offset,k_i2_offset,size_i2)
485      CALL TCE_FILENAME('beta_1_3_4_1_i2',filename)
486      CALL CREATEFILE(filename,d_i2,size_i2)
487      CALL OFFSET_beta_1_3_4_1_1(l_i3_offset,k_i3_offset,size_i3)
488      CALL TCE_FILENAME('beta_1_3_4_1_1_i3',filename)
489      CALL CREATEFILE(filename,d_i3,size_i3)
490      CALL beta_1_3_4_1_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
491     &3_offset)
492      CALL RECONCILEFILE(d_i3,size_i3)
493      CALL beta_1_3_4_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
494     &offset)
495      CALL DELETEFILE(d_i3)
496      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
497     &R)
498      CALL RECONCILEFILE(d_i2,size_i2)
499      CALL beta_1_3_4(d_trc2,k_trc2_offset,d_i2,k_i2_offset,d_i1,k_i1_of
500     &fset)
501      CALL DELETEFILE(d_i2)
502      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
503     &R)
504      CALL OFFSET_beta_1_3_5_1(l_i2_offset,k_i2_offset,size_i2)
505      CALL TCE_FILENAME('beta_1_3_5_1_i2',filename)
506      CALL CREATEFILE(filename,d_i2,size_i2)
507      CALL OFFSET_beta_1_3_5_1_1(l_i3_offset,k_i3_offset,size_i3)
508      CALL TCE_FILENAME('beta_1_3_5_1_1_i3',filename)
509      CALL CREATEFILE(filename,d_i3,size_i3)
510      CALL beta_1_3_5_1_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
511     &3_offset)
512      CALL RECONCILEFILE(d_i3,size_i3)
513      CALL beta_1_3_5_1(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
514     &offset)
515      CALL DELETEFILE(d_i3)
516      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
517     &R)
518      CALL RECONCILEFILE(d_i2,size_i2)
519      CALL beta_1_3_5(d_tra2,k_tra2_offset,d_i2,k_i2_offset,d_i1,k_i1_of
520     &fset)
521      CALL DELETEFILE(d_i2)
522      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
523     &R)
524      CALL OFFSET_beta_1_3_6_1(l_i2_offset,k_i2_offset,size_i2)
525      CALL TCE_FILENAME('beta_1_3_6_1_i2',filename)
526      CALL CREATEFILE(filename,d_i2,size_i2)
527      CALL OFFSET_beta_1_3_6_1_1(l_i3_offset,k_i3_offset,size_i3)
528      CALL TCE_FILENAME('beta_1_3_6_1_1_i3',filename)
529      CALL CREATEFILE(filename,d_i3,size_i3)
530      CALL beta_1_3_6_1_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
531     &3_offset)
532      CALL RECONCILEFILE(d_i3,size_i3)
533      CALL beta_1_3_6_1(d_tra2,k_tra2_offset,d_i3,k_i3_offset,d_i2,k_i2_
534     &offset)
535      CALL DELETEFILE(d_i3)
536      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
537     &R)
538      CALL OFFSET_beta_1_3_6_2_1(l_i3_offset,k_i3_offset,size_i3)
539      CALL TCE_FILENAME('beta_1_3_6_2_1_i3',filename)
540      CALL CREATEFILE(filename,d_i3,size_i3)
541      CALL beta_1_3_6_2_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
542     &3_offset)
543      CALL RECONCILEFILE(d_i3,size_i3)
544      CALL beta_1_3_6_2(d_tra2,k_tra2_offset,d_i3,k_i3_offset,d_i2,k_i2_
545     &offset)
546      CALL DELETEFILE(d_i3)
547      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
548     &R)
549      CALL RECONCILEFILE(d_i2,size_i2)
550      CALL beta_1_3_6(d_trc1,k_trc1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
551     &fset)
552      CALL DELETEFILE(d_i2)
553      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
554     &R)
555      CALL OFFSET_beta_1_3_7_1(l_i2_offset,k_i2_offset,size_i2)
556      CALL TCE_FILENAME('beta_1_3_7_1_i2',filename)
557      CALL CREATEFILE(filename,d_i2,size_i2)
558      CALL OFFSET_beta_1_3_7_1_1(l_i3_offset,k_i3_offset,size_i3)
559      CALL TCE_FILENAME('beta_1_3_7_1_1_i3',filename)
560      CALL CREATEFILE(filename,d_i3,size_i3)
561      CALL beta_1_3_7_1_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i3,k_i
562     &3_offset)
563      CALL RECONCILEFILE(d_i3,size_i3)
564      CALL beta_1_3_7_1(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
565     &offset)
566      CALL DELETEFILE(d_i3)
567      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
568     &R)
569      CALL OFFSET_beta_1_3_7_2_1(l_i3_offset,k_i3_offset,size_i3)
570      CALL TCE_FILENAME('beta_1_3_7_2_1_i3',filename)
571      CALL CREATEFILE(filename,d_i3,size_i3)
572      CALL beta_1_3_7_2_1(d_trb1,k_trb1_offset,d_v2,k_v2_offset,d_i3,k_i
573     &3_offset)
574      CALL RECONCILEFILE(d_i3,size_i3)
575      CALL beta_1_3_7_2(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
576     &offset)
577      CALL DELETEFILE(d_i3)
578      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
579     &R)
580      CALL RECONCILEFILE(d_i2,size_i2)
581      CALL beta_1_3_7(d_tra2,k_tra2_offset,d_i2,k_i2_offset,d_i1,k_i1_of
582     &fset)
583      CALL DELETEFILE(d_i2)
584      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
585     &R)
586      CALL OFFSET_beta_1_3_8_1(l_i2_offset,k_i2_offset,size_i2)
587      CALL TCE_FILENAME('beta_1_3_8_1_i2',filename)
588      CALL CREATEFILE(filename,d_i2,size_i2)
589      CALL OFFSET_beta_1_3_8_1_1(l_i3_offset,k_i3_offset,size_i3)
590      CALL TCE_FILENAME('beta_1_3_8_1_1_i3',filename)
591      CALL CREATEFILE(filename,d_i3,size_i3)
592      CALL OFFSET_beta_1_3_8_1_1_1(l_i4_offset,k_i4_offset,size_i4)
593      CALL TCE_FILENAME('beta_1_3_8_1_1_1_i4',filename)
594      CALL CREATEFILE(filename,d_i4,size_i4)
595      CALL beta_1_3_8_1_1_1(d_trc1,k_trc1_offset,d_v2,k_v2_offset,d_i4,k
596     &_i4_offset)
597      CALL RECONCILEFILE(d_i4,size_i4)
598      CALL beta_1_3_8_1_1(d_trb1,k_trb1_offset,d_i4,k_i4_offset,d_i3,k_i
599     &3_offset)
600      CALL DELETEFILE(d_i4)
601      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
602     &R)
603      CALL RECONCILEFILE(d_i3,size_i3)
604      CALL beta_1_3_8_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
605     &offset)
606      CALL DELETEFILE(d_i3)
607      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
608     &R)
609      CALL RECONCILEFILE(d_i2,size_i2)
610      CALL beta_1_3_8(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i1,k_i1_offset
611     &)
612      CALL DELETEFILE(d_i2)
613      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
614     &R)
615      CALL RECONCILEFILE(d_i1,size_i1)
616      CALL beta_1_3(d_y2,k_y2_offset,d_i1,k_i1_offset,d_i0,k_i0_offset)
617      CALL DELETEFILE(d_i1)
618      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
619     &R)
620      CALL OFFSET_beta_1_4_1(l_i1_offset,k_i1_offset,size_i1)
621      CALL TCE_FILENAME('beta_1_4_1_i1',filename)
622      CALL CREATEFILE(filename,d_i1,size_i1)
623      CALL OFFSET_beta_1_4_1_1(l_i2_offset,k_i2_offset,size_i2)
624      CALL TCE_FILENAME('beta_1_4_1_1_i2',filename)
625      CALL CREATEFILE(filename,d_i2,size_i2)
626      CALL OFFSET_beta_1_4_1_1_1(l_i3_offset,k_i3_offset,size_i3)
627      CALL TCE_FILENAME('beta_1_4_1_1_1_i3',filename)
628      CALL CREATEFILE(filename,d_i3,size_i3)
629      CALL beta_1_4_1_1_1(d_trb2,k_trb2_offset,d_y2,k_y2_offset,d_i3,k_i
630     &3_offset)
631      CALL OFFSET_beta_1_4_1_1_2_1(l_i4_offset,k_i4_offset,size_i4)
632      CALL TCE_FILENAME('beta_1_4_1_1_2_1_i4',filename)
633      CALL CREATEFILE(filename,d_i4,size_i4)
634      CALL beta_1_4_1_1_2_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i4,k
635     &_i4_offset)
636      CALL RECONCILEFILE(d_i4,size_i4)
637      CALL beta_1_4_1_1_2(d_t1,k_t1_offset,d_i4,k_i4_offset,d_i3,k_i3_of
638     &fset)
639      CALL DELETEFILE(d_i4)
640      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
641     &R)
642      CALL RECONCILEFILE(d_i3,size_i3)
643      CALL beta_1_4_1_1(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
644     &offset)
645      CALL DELETEFILE(d_i3)
646      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
647     &R)
648      CALL OFFSET_beta_1_4_1_2_1(l_i3_offset,k_i3_offset,size_i3)
649      CALL TCE_FILENAME('beta_1_4_1_2_1_i3',filename)
650      CALL CREATEFILE(filename,d_i3,size_i3)
651      CALL beta_1_4_1_2_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i3,k_i
652     &3_offset)
653      CALL RECONCILEFILE(d_i3,size_i3)
654      CALL beta_1_4_1_2(d_trb2,k_trb2_offset,d_i3,k_i3_offset,d_i2,k_i2_
655     &offset)
656      CALL DELETEFILE(d_i3)
657      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
658     &R)
659      CALL OFFSET_beta_1_4_1_3_1(l_i3_offset,k_i3_offset,size_i3)
660      CALL TCE_FILENAME('beta_1_4_1_3_1_i3',filename)
661      CALL CREATEFILE(filename,d_i3,size_i3)
662      CALL beta_1_4_1_3_1(d_trc2,k_trc2_offset,d_y2,k_y2_offset,d_i3,k_i
663     &3_offset)
664      CALL OFFSET_beta_1_4_1_3_2_1(l_i4_offset,k_i4_offset,size_i4)
665      CALL TCE_FILENAME('beta_1_4_1_3_2_1_i4',filename)
666      CALL CREATEFILE(filename,d_i4,size_i4)
667      CALL beta_1_4_1_3_2_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i4,k
668     &_i4_offset)
669      CALL RECONCILEFILE(d_i4,size_i4)
670      CALL beta_1_4_1_3_2(d_t1,k_t1_offset,d_i4,k_i4_offset,d_i3,k_i3_of
671     &fset)
672      CALL DELETEFILE(d_i4)
673      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
674     &R)
675      CALL RECONCILEFILE(d_i3,size_i3)
676      CALL beta_1_4_1_3(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
677     &offset)
678      CALL DELETEFILE(d_i3)
679      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
680     &R)
681      CALL OFFSET_beta_1_4_1_4_1(l_i3_offset,k_i3_offset,size_i3)
682      CALL TCE_FILENAME('beta_1_4_1_4_1_i3',filename)
683      CALL CREATEFILE(filename,d_i3,size_i3)
684      CALL beta_1_4_1_4_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i3,k_i
685     &3_offset)
686      CALL RECONCILEFILE(d_i3,size_i3)
687      CALL beta_1_4_1_4(d_trc2,k_trc2_offset,d_i3,k_i3_offset,d_i2,k_i2_
688     &offset)
689      CALL DELETEFILE(d_i3)
690      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
691     &R)
692      CALL OFFSET_beta_1_4_1_5_1(l_i3_offset,k_i3_offset,size_i3)
693      CALL TCE_FILENAME('beta_1_4_1_5_1_i3',filename)
694      CALL CREATEFILE(filename,d_i3,size_i3)
695      CALL OFFSET_beta_1_4_1_5_1_1(l_i4_offset,k_i4_offset,size_i4)
696      CALL TCE_FILENAME('beta_1_4_1_5_1_1_i4',filename)
697      CALL CREATEFILE(filename,d_i4,size_i4)
698      CALL beta_1_4_1_5_1_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i4,k
699     &_i4_offset)
700      CALL RECONCILEFILE(d_i4,size_i4)
701      CALL beta_1_4_1_5_1(d_trb1,k_trb1_offset,d_i4,k_i4_offset,d_i3,k_i
702     &3_offset)
703      CALL DELETEFILE(d_i4)
704      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
705     &R)
706      CALL RECONCILEFILE(d_i3,size_i3)
707      CALL beta_1_4_1_5(d_t1,k_t1_offset,d_i3,k_i3_offset,d_i2,k_i2_offs
708     &et)
709      CALL DELETEFILE(d_i3)
710      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
711     &R)
712      CALL RECONCILEFILE(d_i2,size_i2)
713      CALL beta_1_4_1(d_tra1,k_tra1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
714     &fset)
715      CALL DELETEFILE(d_i2)
716      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
717     &R)
718      CALL OFFSET_beta_1_4_2_1(l_i2_offset,k_i2_offset,size_i2)
719      CALL TCE_FILENAME('beta_1_4_2_1_i2',filename)
720      CALL CREATEFILE(filename,d_i2,size_i2)
721      CALL OFFSET_beta_1_4_2_1_1(l_i3_offset,k_i3_offset,size_i3)
722      CALL TCE_FILENAME('beta_1_4_2_1_1_i3',filename)
723      CALL CREATEFILE(filename,d_i3,size_i3)
724      CALL beta_1_4_2_1_1(d_trb2,k_trb2_offset,d_y2,k_y2_offset,d_i3,k_i
725     &3_offset)
726      CALL RECONCILEFILE(d_i3,size_i3)
727      CALL beta_1_4_2_1(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
728     &offset)
729      CALL DELETEFILE(d_i3)
730      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
731     &R)
732      CALL OFFSET_beta_1_4_2_2_1(l_i3_offset,k_i3_offset,size_i3)
733      CALL TCE_FILENAME('beta_1_4_2_2_1_i3',filename)
734      CALL CREATEFILE(filename,d_i3,size_i3)
735      CALL beta_1_4_2_2_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i3,k_i
736     &3_offset)
737      CALL RECONCILEFILE(d_i3,size_i3)
738      CALL beta_1_4_2_2(d_tra2,k_tra2_offset,d_i3,k_i3_offset,d_i2,k_i2_
739     &offset)
740      CALL DELETEFILE(d_i3)
741      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
742     &R)
743      CALL RECONCILEFILE(d_i2,size_i2)
744      CALL beta_1_4_2(d_trc1,k_trc1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
745     &fset)
746      CALL DELETEFILE(d_i2)
747      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
748     &R)
749      CALL OFFSET_beta_1_4_3_1(l_i2_offset,k_i2_offset,size_i2)
750      CALL TCE_FILENAME('beta_1_4_3_1_i2',filename)
751      CALL CREATEFILE(filename,d_i2,size_i2)
752      CALL OFFSET_beta_1_4_3_1_1(l_i3_offset,k_i3_offset,size_i3)
753      CALL TCE_FILENAME('beta_1_4_3_1_1_i3',filename)
754      CALL CREATEFILE(filename,d_i3,size_i3)
755      CALL beta_1_4_3_1_1(d_trb2,k_trb2_offset,d_y2,k_y2_offset,d_i3,k_i
756     &3_offset)
757      CALL RECONCILEFILE(d_i3,size_i3)
758      CALL beta_1_4_3_1(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
759     &offset)
760      CALL DELETEFILE(d_i3)
761      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
762     &R)
763      CALL OFFSET_beta_1_4_3_2_1(l_i3_offset,k_i3_offset,size_i3)
764      CALL TCE_FILENAME('beta_1_4_3_2_1_i3',filename)
765      CALL CREATEFILE(filename,d_i3,size_i3)
766      CALL beta_1_4_3_2_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i3,k_i
767     &3_offset)
768      CALL RECONCILEFILE(d_i3,size_i3)
769      CALL beta_1_4_3_2(d_trb2,k_trb2_offset,d_i3,k_i3_offset,d_i2,k_i2_
770     &offset)
771      CALL DELETEFILE(d_i3)
772      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
773     &R)
774      CALL OFFSET_beta_1_4_3_3_1(l_i3_offset,k_i3_offset,size_i3)
775      CALL TCE_FILENAME('beta_1_4_3_3_1_i3',filename)
776      CALL CREATEFILE(filename,d_i3,size_i3)
777      CALL beta_1_4_3_3_1(d_trc2,k_trc2_offset,d_y2,k_y2_offset,d_i3,k_i
778     &3_offset)
779      CALL RECONCILEFILE(d_i3,size_i3)
780      CALL beta_1_4_3_3(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
781     &offset)
782      CALL DELETEFILE(d_i3)
783      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
784     &R)
785      CALL OFFSET_beta_1_4_3_4_1(l_i3_offset,k_i3_offset,size_i3)
786      CALL TCE_FILENAME('beta_1_4_3_4_1_i3',filename)
787      CALL CREATEFILE(filename,d_i3,size_i3)
788      CALL beta_1_4_3_4_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i3,k_i
789     &3_offset)
790      CALL RECONCILEFILE(d_i3,size_i3)
791      CALL beta_1_4_3_4(d_trc2,k_trc2_offset,d_i3,k_i3_offset,d_i2,k_i2_
792     &offset)
793      CALL DELETEFILE(d_i3)
794      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
795     &R)
796      CALL RECONCILEFILE(d_i2,size_i2)
797      CALL beta_1_4_3(d_tra1,k_tra1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
798     &fset)
799      CALL DELETEFILE(d_i2)
800      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
801     &R)
802      CALL OFFSET_beta_1_4_4_1(l_i2_offset,k_i2_offset,size_i2)
803      CALL TCE_FILENAME('beta_1_4_4_1_i2',filename)
804      CALL CREATEFILE(filename,d_i2,size_i2)
805      CALL OFFSET_beta_1_4_4_1_1(l_i3_offset,k_i3_offset,size_i3)
806      CALL TCE_FILENAME('beta_1_4_4_1_1_i3',filename)
807      CALL CREATEFILE(filename,d_i3,size_i3)
808      CALL beta_1_4_4_1_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i3,k_i
809     &3_offset)
810      CALL RECONCILEFILE(d_i3,size_i3)
811      CALL beta_1_4_4_1(d_tra2,k_tra2_offset,d_i3,k_i3_offset,d_i2,k_i2_
812     &offset)
813      CALL DELETEFILE(d_i3)
814      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
815     &R)
816      CALL RECONCILEFILE(d_i2,size_i2)
817      CALL beta_1_4_4(d_trc1,k_trc1_offset,d_i2,k_i2_offset,d_i1,k_i1_of
818     &fset)
819      CALL DELETEFILE(d_i2)
820      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
821     &R)
822      CALL OFFSET_beta_1_4_5_1(l_i2_offset,k_i2_offset,size_i2)
823      CALL TCE_FILENAME('beta_1_4_5_1_i2',filename)
824      CALL CREATEFILE(filename,d_i2,size_i2)
825      CALL OFFSET_beta_1_4_5_1_1(l_i3_offset,k_i3_offset,size_i3)
826      CALL TCE_FILENAME('beta_1_4_5_1_1_i3',filename)
827      CALL CREATEFILE(filename,d_i3,size_i3)
828      CALL beta_1_4_5_1_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i3,k_i
829     &3_offset)
830      CALL RECONCILEFILE(d_i3,size_i3)
831      CALL beta_1_4_5_1(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
832     &offset)
833      CALL DELETEFILE(d_i3)
834      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
835     &R)
836      CALL RECONCILEFILE(d_i2,size_i2)
837      CALL beta_1_4_5(d_tra2,k_tra2_offset,d_i2,k_i2_offset,d_i1,k_i1_of
838     &fset)
839      CALL DELETEFILE(d_i2)
840      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
841     &R)
842      CALL OFFSET_beta_1_4_6_1(l_i2_offset,k_i2_offset,size_i2)
843      CALL TCE_FILENAME('beta_1_4_6_1_i2',filename)
844      CALL CREATEFILE(filename,d_i2,size_i2)
845      CALL OFFSET_beta_1_4_6_1_1(l_i3_offset,k_i3_offset,size_i3)
846      CALL TCE_FILENAME('beta_1_4_6_1_1_i3',filename)
847      CALL CREATEFILE(filename,d_i3,size_i3)
848      CALL OFFSET_beta_1_4_6_1_1_1(l_i4_offset,k_i4_offset,size_i4)
849      CALL TCE_FILENAME('beta_1_4_6_1_1_1_i4',filename)
850      CALL CREATEFILE(filename,d_i4,size_i4)
851      CALL beta_1_4_6_1_1_1(d_trb1,k_trb1_offset,d_y2,k_y2_offset,d_i4,k
852     &_i4_offset)
853      CALL RECONCILEFILE(d_i4,size_i4)
854      CALL beta_1_4_6_1_1(d_tra1,k_tra1_offset,d_i4,k_i4_offset,d_i3,k_i
855     &3_offset)
856      CALL DELETEFILE(d_i4)
857      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
858     &R)
859      CALL RECONCILEFILE(d_i3,size_i3)
860      CALL beta_1_4_6_1(d_trc1,k_trc1_offset,d_i3,k_i3_offset,d_i2,k_i2_
861     &offset)
862      CALL DELETEFILE(d_i3)
863      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
864     &R)
865      CALL OFFSET_beta_1_4_6_2_1(l_i3_offset,k_i3_offset,size_i3)
866      CALL TCE_FILENAME('beta_1_4_6_2_1_i3',filename)
867      CALL CREATEFILE(filename,d_i3,size_i3)
868      CALL OFFSET_beta_1_4_6_2_1_1(l_i4_offset,k_i4_offset,size_i4)
869      CALL TCE_FILENAME('beta_1_4_6_2_1_1_i4',filename)
870      CALL CREATEFILE(filename,d_i4,size_i4)
871      CALL beta_1_4_6_2_1_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i4,k
872     &_i4_offset)
873      CALL RECONCILEFILE(d_i4,size_i4)
874      CALL beta_1_4_6_2_1(d_tra1,k_tra1_offset,d_i4,k_i4_offset,d_i3,k_i
875     &3_offset)
876      CALL DELETEFILE(d_i4)
877      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
878     &R)
879      CALL RECONCILEFILE(d_i3,size_i3)
880      CALL beta_1_4_6_2(d_trb1,k_trb1_offset,d_i3,k_i3_offset,d_i2,k_i2_
881     &offset)
882      CALL DELETEFILE(d_i3)
883      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
884     &R)
885      CALL OFFSET_beta_1_4_6_3_1(l_i3_offset,k_i3_offset,size_i3)
886      CALL TCE_FILENAME('beta_1_4_6_3_1_i3',filename)
887      CALL CREATEFILE(filename,d_i3,size_i3)
888      CALL OFFSET_beta_1_4_6_3_1_1(l_i4_offset,k_i4_offset,size_i4)
889      CALL TCE_FILENAME('beta_1_4_6_3_1_1_i4',filename)
890      CALL CREATEFILE(filename,d_i4,size_i4)
891      CALL beta_1_4_6_3_1_1(d_trc1,k_trc1_offset,d_y2,k_y2_offset,d_i4,k
892     &_i4_offset)
893      CALL RECONCILEFILE(d_i4,size_i4)
894      CALL beta_1_4_6_3_1(d_trb1,k_trb1_offset,d_i4,k_i4_offset,d_i3,k_i
895     &3_offset)
896      CALL DELETEFILE(d_i4)
897      IF (.not.MA_POP_STACK(l_i4_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
898     &R)
899      CALL RECONCILEFILE(d_i3,size_i3)
900      CALL beta_1_4_6_3(d_tra1,k_tra1_offset,d_i3,k_i3_offset,d_i2,k_i2_
901     &offset)
902      CALL DELETEFILE(d_i3)
903      IF (.not.MA_POP_STACK(l_i3_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
904     &R)
905      CALL RECONCILEFILE(d_i2,size_i2)
906      CALL beta_1_4_6(d_t1,k_t1_offset,d_i2,k_i2_offset,d_i1,k_i1_offset
907     &)
908      CALL DELETEFILE(d_i2)
909      IF (.not.MA_POP_STACK(l_i2_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
910     &R)
911      CALL RECONCILEFILE(d_i1,size_i1)
912      CALL beta_1_4(d_i1,k_i1_offset,d_v2,k_v2_offset,d_i0,k_i0_offset)
913      CALL DELETEFILE(d_i1)
914      IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('beta_1',-1,MA_ER
915     &R)
916      RETURN
917      END
918