


















































































































 


































































































































































 






























































































































































































 



















































 










 








































































































































































































































 










                                
                                

                                
                                

                                
                                

                                





































c

      SUBROUTINE unphysical_to_physical(grr,grt,grp,gtt,gtp,gpp,dgtt,
     &     dgtp,dgpp,Psis,dPsis,Psi_power)

c
c
c
c
c
c
c

      IMPLICIT NONE

c
      INTEGER :: 
     &    Psi_power
      REAL*8,INTENT(INOUT),DIMENSION (:,:) :: 
     &    grr,grt,grp,gtt,gtp,gpp,dgtt,dgtp,dgpp
      REAL*8,INTENT(IN),DIMENSION (:,:) :: 
     &    Psis,dPsis

c
c

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

      








