diff options
author | schnetter <schnetter@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2005-03-19 15:29:53 +0000 |
---|---|---|
committer | schnetter <schnetter@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2005-03-19 15:29:53 +0000 |
commit | b68b294b343f555efeaa085fe08ba03aa6d9cd11 (patch) | |
tree | a4e1fe73ea0410349144a375f97bf647004d819a | |
parent | df3790a6c152c338172490134561ed9e4b191747 (diff) |
Use sqrt instead of **0.5
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@88 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7
-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) |