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     ==================================================================
+
+      
+
+
+
+
+
+
+
+