1 c---------------------------------------------------------------------
2 c---------------------------------------------------------------------
4 subroutine set_constants
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
89 dnxm1 = 1.0d0 / dble(grid_points(1)-1)
90 dnym1 = 1.0d0 / dble(grid_points(2)-1)
91 dnzm1 = 1.0d0 / dble(grid_points(3)-1)
100 tx1 = 1.0d0 / (dnxm1 * dnxm1)
101 tx2 = 1.0d0 / (2.0d0 * dnxm1)
104 ty1 = 1.0d0 / (dnym1 * dnym1)
105 ty2 = 1.0d0 / (2.0d0 * dnym1)
108 tz1 = 1.0d0 / (dnzm1 * dnzm1)
109 tz2 = 1.0d0 / (2.0d0 * dnzm1)
130 dxmax = dmax1(dx3, dx4)
131 dymax = dmax1(dy2, dy4)
132 dzmax = dmax1(dz2, dz3)
134 dssp = 0.25d0 * dmax1(dx1, dmax1(dy1, dz1) )
136 c4dssp = 4.0d0 * dssp
137 c5dssp = 5.0d0 * dssp
146 c2dttx1 = 2.0d0*dttx1
147 c2dtty1 = 2.0d0*dtty1
148 c2dttz1 = 2.0d0*dttz1
183 xxcon1 = c3c4tx3*con43*tx3
185 xxcon3 = c3c4tx3*conz1*tx3
186 xxcon4 = c3c4tx3*con16*tx3
187 xxcon5 = c3c4tx3*c1c5*tx3
189 yycon1 = c3c4ty3*con43*ty3
191 yycon3 = c3c4ty3*conz1*ty3
192 yycon4 = c3c4ty3*con16*ty3
193 yycon5 = c3c4ty3*c1c5*ty3
195 zzcon1 = c3c4tz3*con43*tz3
197 zzcon3 = c3c4tz3*conz1*tz3
198 zzcon4 = c3c4tz3*con16*tz3
199 zzcon5 = c3c4tz3*c1c5*tz3