[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 148)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/08/09 01:26 AM

Modified:
 /trunk/
  param.ccl
 /trunk/src/
  GRHydro_Bondi_new.F90

Log:
 Added option to smooth metric inside horizon
 
 Signed-off-by: Bruno Coutinho Mundim <bcmsma at astro.rit.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Bondi_new.F90
Delta lines: +11 -2
===================================================================
--- trunk/src/GRHydro_Bondi_new.F90	2012-08-09 06:26:54 UTC (rev 147)
+++ trunk/src/GRHydro_Bondi_new.F90	2012-08-09 06:26:55 UTC (rev 148)
@@ -55,7 +55,7 @@
   CCTK_REAL :: M, Msq, Mdot, rs, gam, rmin_bondi, rmax_bondi, cs_sq,cs,vs_sq,vs,rhos,gtemp,hs, Kval, Qdot
   CCTK_REAL :: logrmin,dlogr,rhotmp,utmp,vtmp,rspher
   CCTK_REAL :: r_bondi(N_points), logr_bondi(N_points), rho_bondi(N_points), u_bondi(N_points), v_bondi(N_points)
-  CCTK_REAL :: drhodr, det, rhocheck, rhocheck2, riso, rnew, rsch, ucheck
+  CCTK_REAL :: drhodr, det, rhocheck, rhocheck2, riso, rnew, rsch, ucheck, psinew
   CCTK_REAL :: uiso, uisocheck, vcheck, ucheck2, vcheck2, xhat,yhat, zhat, xp, yp, zp
 
   !!$set_bondi_parameters
@@ -169,7 +169,16 @@
            yhat = yp/riso
            zhat = zp/riso
            
-           if(riso < 1.0e-7) then
+           if(use_smooth_puncture_data.ne.0 .and. riso.lt.M) then
+              psinew=2.875 - 5*(riso/M)**2 + 6*(riso/M)**4
+              gxx(i,j,k) = psinew**4
+              gyy(i,j,k) = psinew**4
+              gzz(i,j,k) = psinew**4
+              gxy(i,j,k) = 0.0
+              gxz(i,j,k) = 0.0
+              gyz(i,j,k) = 0.0
+
+           else if(riso < 1.0e-7) then
               gxx(i,j,k) = 1.0e4
               gyy(i,j,k) = 1.0e4
               gzz(i,j,k) = 1.0e4

Directory: /trunk/
==================

File [modified]: param.ccl
Delta lines: +4 -0
===================================================================
--- trunk/param.ccl	2012-08-09 06:26:54 UTC (rev 147)
+++ trunk/param.ccl	2012-08-09 06:26:55 UTC (rev 148)
@@ -97,6 +97,10 @@
   0.0:*   :: "Anything positive"
 } 1.0
 
+BOOLEAN use_smooth_puncture_data "Smooth puncture solution inside R=M using 4th order matching polynomial"
+{
+} "no"
+
 BOOLEAN change_shock_direction "Change the shock direction"
 { 
 } "no"