diff options
| -rw-r--r-- | src/planewaves.F77 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/planewaves.F77 b/src/planewaves.F77 index 4a40688..e33ece2 100644 --- a/src/planewaves.F77 +++ b/src/planewaves.F77 @@ -108,12 +108,12 @@ c precalc kx = 2*pi*sin(the)*cos(phi)/wave ky = 2*pi*sin(the)*sin(phi)/wave kz = 2*pi*cos(the)/wave - w = (kx*kx+ky*ky+kz*kz)**0.5 + w = sqrt(kx*kx+ky*ky+kz*kz) kxp = 2*pi*sin(the)*cos(phi)/wavep kyp = 2*pi*sin(the)*sin(phi)/wavep kzp = 2*pi*cos(the)/wavep - wp = (kxp*kxp+kyp*kyp+kzp*kzp)**0.5 + wp = sqrt(kxp*kxp+kyp*kyp+kzp*kzp) c *************** plane waves ******************** do k=1,cctk_lsh(3) |