diff options
author | allen <allen@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 1999-11-01 14:56:51 +0000 |
---|---|---|
committer | allen <allen@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 1999-11-01 14:56:51 +0000 |
commit | 661c56c847f4334feb53fd8d94551b2ab1efdb9a (patch) | |
tree | 3dbdeebf2c0b8f0f0fe30218e9a58fc0cd1638a1 /src/unphysical_to_physical.F | |
parent | d8cf32e113a96b24744d231013e3d0b81bf55c30 (diff) |
This commit was generated by cvs2svn to compensate for changes in r2, which
included commits to RCS files with non-trunk default branches.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Extract/trunk@3 5301f0c2-dbc4-4cee-b2f5-8d7afba4d129
Diffstat (limited to 'src/unphysical_to_physical.F')
-rw-r--r-- | src/unphysical_to_physical.F | 56 |
1 files changed, 56 insertions, 0 deletions
diff --git a/src/unphysical_to_physical.F b/src/unphysical_to_physical.F new file mode 100644 index 0000000..74b4ad8 --- /dev/null +++ b/src/unphysical_to_physical.F @@ -0,0 +1,56 @@ + +#include "cctk.h" + +c ================================================================== + + SUBROUTINE unphysical_to_physical(grr,grt,grp,gtt,gtp,gpp,dgtt, + & dgtp,dgpp,Psis,dPsis,Psi_power) + +c ------------------------------------------------------------------ +c +c Convert unphysical metric components and radial (eta) derivatives +c on the 2-sphere into physical quantities, using the conformal factor +c on the sphere +c +c ------------------------------------------------------------------ + + IMPLICIT NONE + +c Input variables + INTEGER :: + & Psi_power + CCTK_REAL,INTENT(INOUT),DIMENSION (:,:) :: + & grr,grt,grp,gtt,gtp,gpp,dgtt,dgtp,dgpp + CCTK_REAL,INTENT(IN),DIMENSION (:,:) :: + & Psis,dPsis + +c Output variables +c WE ARE CHANGING THE INPUT VARIABLES !!!! + +c ------------------------------------------------------------------ + + dgtt = Psis**4*dgtt + 4.0d0*Psis**3*dPsis*gtt + dgtp = Psis**4*dgtp + 4.0d0*Psis**3*dPsis*gtp + dgpp = Psis**4*dgpp + 4.0d0*Psis**3*dPsis*gpp + + grr = Psis**4*grr + grt = Psis**4*grt + grp = Psis**4*grp + gtt = Psis**4*gtt + gtp = Psis**4*gtp + gpp = Psis**4*gpp + + + END SUBROUTINE unphysical_to_physical + +c ================================================================== + + + + + + + + + + |