diff options
Diffstat (limited to 'src/constants.F90')
| -rw-r--r-- | src/constants.F90 | 24 |
1 files changed, 23 insertions, 1 deletions
diff --git a/src/constants.F90 b/src/constants.F90 index 35555d3..5f162dc 100644 --- a/src/constants.F90 +++ b/src/constants.F90 @@ -8,7 +8,7 @@ module constants public delta2, delta3, delta4 public eta4 - public epsilon2, epsilon3 + public epsilon2, epsilon3, epsilon4 public pi @@ -43,6 +43,28 @@ module constants & 0,0,-1, 0,0,0, +1,0,0, & & 0,+1,0, -1,0,0, 0,0,0 /), (/3,3,3/)) + CCTK_REAL, parameter :: epsilon4(0:3,0:3,0:3,0:3) & + = reshape ((/ 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, & + & 0,0,0,0, 0,0,0,0, 0,0,0,+1, 0,0,-1,0, & + & 0,0,0,0, 0,0,0,-1, 0,0,0,0, 0,+1,0,0, & + & 0,0,0,0, 0,0,+1,0, 0,-1,0,0, 0,0,0,0, & + + & 0,0,0,0, 0,0,0,0, 0,0,0,-1, 0,0,+1,0, & + & 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, & + & 0,0,0,+1, 0,0,0,0, 0,0,0,0, -1,0,0,0, & + & 0,0,-1,0, 0,0,0,0, +1,0,0,0, 0,0,0,0, & + + & 0,0,0,0, 0,0,0,+1, 0,0,0,0, 0,-1,0,0, & + & 0,0,0,-1, 0,0,0,0, 0,0,0,0, +1,0,0,0, & + & 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, & + & 0,+1,0,0, -1,0,0,0, 0,0,0,0, 0,0,0,0, & + + & 0,0,0,0, 0,0,-1,0, 0,+1,0,0, 0,0,0,0, & + & 0,0,+1,0, 0,0,0,0, -1,0,0,0, 0,0,0,0, & + & 0,-1,0,0, +1,0,0,0, 0,0,0,0, 0,0,0,0, & + & 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0 /), & + & (/4,4,4,4/)) + CCTK_REAL, parameter :: pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068_rk |