2 c---------------------------------------------------------------------
3 c---------------------------------------------------------------------
5 subroutine set_constants
7 c---------------------------------------------------------------------
8 c---------------------------------------------------------------------
90 dnxm1 = 1.0d0 / dble(grid_points(1)-1)
91 dnym1 = 1.0d0 / dble(grid_points(2)-1)
92 dnzm1 = 1.0d0 / dble(grid_points(3)-1)
101 tx1 = 1.0d0 / (dnxm1 * dnxm1)
102 tx2 = 1.0d0 / (2.0d0 * dnxm1)
105 ty1 = 1.0d0 / (dnym1 * dnym1)
106 ty2 = 1.0d0 / (2.0d0 * dnym1)
109 tz1 = 1.0d0 / (dnzm1 * dnzm1)
110 tz2 = 1.0d0 / (2.0d0 * dnzm1)
131 dxmax = dmax1(dx3, dx4)
132 dymax = dmax1(dy2, dy4)
133 dzmax = dmax1(dz2, dz3)
135 dssp = 0.25d0 * dmax1(dx1, dmax1(dy1, dz1) )
137 c4dssp = 4.0d0 * dssp
138 c5dssp = 5.0d0 * dssp
147 c2dttx1 = 2.0d0*dttx1
148 c2dtty1 = 2.0d0*dtty1
149 c2dttz1 = 2.0d0*dttz1
184 xxcon1 = c3c4tx3*con43*tx3
186 xxcon3 = c3c4tx3*conz1*tx3
187 xxcon4 = c3c4tx3*con16*tx3
188 xxcon5 = c3c4tx3*c1c5*tx3
190 yycon1 = c3c4ty3*con43*ty3
192 yycon3 = c3c4ty3*conz1*ty3
193 yycon4 = c3c4ty3*con16*ty3
194 yycon5 = c3c4ty3*c1c5*ty3
196 zzcon1 = c3c4tz3*con43*tz3
198 zzcon3 = c3c4tz3*conz1*tz3
199 zzcon4 = c3c4tz3*con16*tz3
200 zzcon5 = c3c4tz3*c1c5*tz3